clear?all
clc
%緯度
latitude?=?load（‘C:\workspace\matlab?workspace\CES9937\latitude.txt‘）;
%經度
longitude?=?load（‘C:\workspace\matlab?workspace\CES9937\longitude.txt‘）;
%每五對數據點做一次數據擬合
%以第六對數據點的緯度值作為基準點，計算經度值
%計算真實經度值和通過擬合曲線計算出的經度值之間的差值
%=================================================================
%真實值與擬合值之差
err_fit?=?[];
%用來做數據擬合的點集
data_fit_lat?=?[];
data_fit_lon?=?[];
%數據擬合得到的結果集
data_fit?=?[];
%對應的真實數據集合
data_real?=?[];
%擬合預測過程
for?i?=?1:length（latitude）
%當前用來做數據擬合的集合中有幾組數據
????n?=?length（data_fit_lat）;
????%當集合中的數據不滿足擬合需要的最少數據點時，繼續向集合中添加數據
????if?n<=10
????????data_fit_lat（n+1）?=?latitude（i）;
????????data_fit_lon（n+1）?=?longitude（i）;
????????continue
????end
????%當集合中的數據滿足擬合需要的最少數據點時
????%曲線擬合
????if?i????????%最小二乘擬合-二階擬合
????????line_fit?=?polyfit（data_fit_latdata_fit_lon2）;
????????%通過擬合函數求下一時刻經度值
????????lon_fit?=?polyval（line_fitlatitude（i+1））;
????????%保存結果
????????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）;
????????%更新擬合數據
????????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\飛行數據-20組.txt

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

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

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

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