%原始數(shù)據(jù)導(dǎo)入部分
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）;
%畫出原始數(shù)據(jù)圖像
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）;
%計(jì)算dy的均值和方差
u=mean（dy）;
sigma2=mean（（dy-u）.^2）;
%核函數(shù)句柄
Kfun=@（u）1/（2*pi）^0.5*exp（-u.^2/2）;
%一階所以p=1
p=1;
%計(jì)算窗寬h
h=（4/（p+2）/n）^（1/（p+4））;
%對(duì)差分的數(shù)據(jù)進(jìn)行預(yù)測(cè)（只預(yù)測(cè)未來的6個(gè)點(diǎn)）
for?i=7:-1:2
????k=1;
????for?j=1:（n-1）-i
????????%計(jì)算滿足|x-xi|????????if（abs（dy（j）-dy（n-i））????????????yg（k）=dy（j）;
????????????k=k+1;
????????end
????end
????%用核函數(shù)進(jìn)行預(yù)測(cè)估計(jì)
????K=Kfun（（dy（n-i）-yg）/h）;
????yp（8-i）=dot（K/sum（K）yg）;
end
%還原數(shù)據(jù)
yp（1）=y（n-7）+yp（1）;
for?i=2:6
????yp（i）=yp（i-1）+yp（i）;
end
%加上隨機(jī)誤差項(xiàng)
yp=yp+u-sqrt（sigma2）+2*sqrt（sigma2）*rand（16）;
%繪制預(yù)測(cè)點(diǎn)及折線
plot（n-5:nyp‘*r‘）;
plot（n-5:nyp‘-r‘）;
%計(jì)算MSE和MAE
MSE=sqrt（sum（（yp-y（n-5:n））.^2））/6
MAE=sum（abs（yp-y（n-5:n）））/6

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

?????文件???????7889??2011-06-21?21:06??時(shí)間序列\(zhòng)data.txt

?????文件???????1097??2011-06-21?21:06??時(shí)間序列\(zhòng)TS.m

?????文件?????412485??2011-06-21?14:22??時(shí)間序列\(zhòng)人口增長(zhǎng)率的非參數(shù)自回歸預(yù)測(cè)模型.pdf

?????文件?????185344??2011-06-21?22:24??時(shí)間序列\(zhòng)時(shí)間序列的數(shù)據(jù)分析.doc

?????文件?????194584??2011-06-21?21:06??時(shí)間序列\(zhòng)畢業(yè)論文?李吟.docx

?????文件????2215936??2011-06-21?12:23??時(shí)間序列\(zhòng)非參數(shù)回歸模型.doc

?????目錄??????????0??2011-06-30?09:25??時(shí)間序列

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

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