%原始數據導入部分
clearclc;
fin=fopen（‘data.txt‘‘r‘）;
y（1）=str2num（fgetl（fin））;
counter=2;
while?feof（fin）==0
????y（counter）=str2num（fgetl（fin））;
????counter=counter+1;
end
n=size（y2）;
%畫出原始數據圖像
plot（n-5:ny（end-5:end）‘pg‘）;
plot（n-5:ny（end-5:end）‘-g‘）;
hold?on;
%求其一階差分
dy=y（2:n）-y（1:n-1）;
%計算dy的均值和方差
u=mean（dy）;
sigma2=mean（（dy-u）.^2）;
%核函數句柄
Kfun=@（u）1/（2*pi）^0.5*exp（-u.^2/2）;
%一階所以p=1
p=1;
%計算窗寬h
h=（4/（p+2）/n）^（1/（p+4））;
%對差分的數據進行預測（只預測未來的6個點）
for?i=7:-1:2
????k=1;
????for?j=1:（n-1）-i
????????%計算滿足|x-xi|????????if（abs（dy（j）-dy（n-i））????????????yg（k）=dy（j）;
????????????k=k+1;
????????end
????end
????%用核函數進行預測估計
????K=Kfun（（dy（n-i）-yg）/h）;
????yp（8-i）=dot（K/sum（K）yg）;
end
%還原數據
yp（1）=y（n-7）+yp（1）;
for?i=2:6
????yp（i）=yp（i-1）+yp（i）;
end
%加上隨機誤差項
yp=yp+u-sqrt（sigma2）+2*sqrt（sigma2）*rand（16）;
%繪制預測點及折線
plot（n-5:nyp‘*r‘）;
plot（n-5:nyp‘-r‘）;
%計算MSE和MAE
MSE=sqrt（sum（（yp-y（n-5:n））.^2））/6
MAE=sum（abs（yp-y（n-5:n）））/6

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

?????文件???????7889??2011-06-21?21:06??時間序列\data.txt

?????文件???????1097??2011-06-21?21:06??時間序列\TS.m

?????文件?????412485??2011-06-21?14:22??時間序列\人口增長率的非參數自回歸預測模型.pdf

?????文件?????185344??2011-06-21?22:24??時間序列\時間序列的數據分析.doc

?????文件?????194584??2011-06-21?21:06??時間序列\畢業論文?李吟.docx

?????文件????2215936??2011-06-21?12:23??時間序列\非參數回歸模型.doc

?????目錄??????????0??2011-06-30?09:25??時間序列

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

??????????????3017335????????????????????7