%---------------------------------------------%
%					??????%
%??????????工作室提供代做matlab仿真	??????%
%					??????%
%??詳情請(qǐng)?jiān)L問：http://cn.mikecrm.com/5k6v1DP??%
%					??????%
%---------------------------------------------%

%使用卡爾曼濾波方法對(duì)飛行航班進(jìn)行軌跡預(yù)測(cè)
%數(shù)據(jù)來源：FlightAware（https://zh.flightaware.com）
%航????班：CES9937?????寧波櫟社國際機(jī)場(chǎng)飛往成都雙流國際機(jī)場(chǎng)
%飛行時(shí)間：2018-07-17??07:07-09:55
%說明：取起飛后前20組數(shù)據(jù)作為實(shí)驗(yàn)數(shù)據(jù)。對(duì)時(shí)間點(diǎn)進(jìn)行近似取值，假設(shè)每隔30s上報(bào)一次數(shù)據(jù)
clear;
clc;
%采樣點(diǎn)的個(gè)數(shù)
N=228;
%測(cè)試數(shù)據(jù)：緯度
latitude=load（‘C:\workspace\matlab?workspace\CES9937\latitude.txt‘）;
%真實(shí)維度值
lat=latitude;
%卡爾曼濾波處理的狀態(tài)，即估計(jì)值
lat_kf=zeros（1N）;
%測(cè)報(bào)值
lat_z=zeros（1N）;
P=zeros（1N）;
%初始緯度值
lat（1）=29.8131;
%初始值的協(xié)方差
P（1）=0.09;
%初始測(cè)報(bào)值
lat_z（1）=29.8027;
%初始估計(jì)狀態(tài)。假設(shè)和初始測(cè)報(bào)值相同
lat_kf（1）=lat_z（1）;
%噪聲方差
%系統(tǒng)噪聲方差
Q=0.1;
%測(cè)量噪聲方差
R=0.001;
%方差決定噪聲大小
W=sqrt

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????153993??2018-07-24?16:28??卡爾曼濾波算法實(shí)現(xiàn)飛行物體運(yùn)動(dòng)軌跡預(yù)測(cè)\CES9937\CES9937.docx

?????文件???????1825??2018-07-17?15:07??卡爾曼濾波算法實(shí)現(xiàn)飛行物體運(yùn)動(dòng)軌跡預(yù)測(cè)\CES9937\latitude.txt

?????文件???????2051??2018-07-17?15:15??卡爾曼濾波算法實(shí)現(xiàn)飛行物體運(yùn)動(dòng)軌跡預(yù)測(cè)\CES9937\longitude.txt

?????文件????????410??2018-07-17?11:16??卡爾曼濾波算法實(shí)現(xiàn)飛行物體運(yùn)動(dòng)軌跡預(yù)測(cè)\CES9937\飛行數(shù)據(jù)-20組.txt

?????文件???????2227??2018-09-19?08:09??卡爾曼濾波算法實(shí)現(xiàn)飛行物體運(yùn)動(dòng)軌跡預(yù)測(cè)\QA_KF_lat.m

?????文件?????????55??2018-08-27?11:33??卡爾曼濾波算法實(shí)現(xiàn)飛行物體運(yùn)動(dòng)軌跡預(yù)測(cè)\【源碼使用必讀】.url

?????目錄??????????0??2018-08-14?20:18??卡爾曼濾波算法實(shí)現(xiàn)飛行物體運(yùn)動(dòng)軌跡預(yù)測(cè)\CES9937

?????目錄??????????0??2018-09-16?11:01??卡爾曼濾波算法實(shí)現(xiàn)飛行物體運(yùn)動(dòng)軌跡預(yù)測(cè)

-----------?---------??----------?-----??----

???????????????160561????????????????????8