clear?all
clc
%緯度
latitude?=?load（‘C:\workspace\matlab?workspace\CES9937\latitude.txt‘）;
%經(jīng)度
longitude?=?load（‘C:\workspace\matlab?workspace\CES9937\longitude.txt‘）;
%每五對數(shù)據(jù)點做一次數(shù)據(jù)擬合
%以第六對數(shù)據(jù)點的緯度值作為基準點，計算經(jīng)度值
%計算真實經(jīng)度值和通過擬合曲線計算出的經(jīng)度值之間的差值
%=================================================================
%真實值與擬合值之差
err_fit?=?[];
%用來做數(shù)據(jù)擬合的點集
data_fit_lat?=?[];
data_fit_lon?=?[];
%數(shù)據(jù)擬合得到的結(jié)果集
data_fit?=?[];
%對應的真實數(shù)據(jù)集合
data_real?=?[];
%擬合預測過程
for?i?=?1:length（latitude）
%當前用來做數(shù)據(jù)擬合的集合中有幾組數(shù)據(jù)
????n?=?length（data_fit_lat）;
????%當集合中的數(shù)據(jù)不滿足擬合需要的最少數(shù)據(jù)點時，繼續(xù)向集合中添加數(shù)據(jù)
????if?n<=10
????????data_fit_lat（n+1）?=?latitude（i）;
????????data_fit_lon（n+1）?=?longitude（i）;
????????continue
????end
????%當集合中的數(shù)據(jù)滿足擬合需要的最少數(shù)據(jù)點時
????%曲線擬合
????if?i????????%最小二乘擬合-二階擬合
????????line_fit?=?polyfit（data_fit_latdata_fit_lon2）;
????????%通過擬合函數(shù)求下一時刻經(jīng)度值
????????lon_fit?=?polyval（line_fitlatitude（i+1））;
????????%保存結(jié)果
????????df?=?length（data_fit）;
????????data_fit（df+1）?=?lon_fit;
????????%計算預測值與真實值之間的差值
????????m?=?length（err_fit）;
????????err_fit（m+1）?=?abs（longitude（i+1）?-?lon_fit）;
????????%記錄真實值
????????dr?=?length（data_real）;
????????data_real（dr+1）?=?longitude（i+1）;
????????%更新擬合數(shù)據(jù)
????????data_fit_lat（1）?=?[];
????????data_fit_lat（10）?=?latitude（i+1）;
????????data_fit_lon（1）?=?[];
????????data_fit_lon（10）?=?longitude（i+1）;
????end
end
%=================================================================
%繪制誤差圖
k?=?1:length（err_fit）;
figure
hold?on;
grid?on;
plot（kerr_fit‘k.‘）;
xlabel（‘X‘）;
ylabel（‘誤差值‘）;
title（‘擬合‘）;
%繪制預測值與真實值曲線圖
b?=?1:length（data_real）;
figure
hold?on;
grid?on;
plot（bdata_real‘g.-‘bdata_fit‘r.-‘）;
xlabel（‘X‘）;
ylabel（‘誤差值‘）;
title（‘擬合‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????153993??2018-07-24?16:28??CES9937\CES9937.docx

?????文件???????1825??2018-07-17?15:07??CES9937\latitude.txt

?????文件???????2051??2018-07-17?15:15??CES9937\longitude.txt

?????文件????????410??2018-07-17?11:16??CES9937\飛行數(shù)據(jù)-20組.txt

?????文件???????2000??2018-08-14?13:46??re_fit_4.m

?????目錄??????????0??2018-08-14?20:18??CES9937

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

???????????????160279????????????????????6