时间序列地matlab实例编程

1、%时间序列ex1 * %x-72个源数据%各个图像需要单独的操作完成手工进行x=90078106 8928 913710017 10826 11317 10744 9713 9938 9161 8927 77506981 8038 8422 8714 9512 10120 9823 8743 9192 8710 8680 8162 7306 8124 78709387 9556 10093 96208285 8433 8160 8034 7717 74617776 7925 8634894510078 91798037 84887874 8647 77926957 7726 8106 8890

2、 9299 1062593028314 8850 8265 8796 7836 6892 7791 8129 9115 9434 10484 9827 9110 9070 8633 9240；Tl=zeros(1,72);t=965175 8723.75 8585.8+1/30 8396.75 8576.8+1/30 8796.75；亡人数for i=1:6T1(12*k11):(12*i)=t(i);%趋势项 T1 的赋值end%据图和分段趋势% plot(x)：% hold on% plot仃仁T);% hold off%m ethodi -分段趋势法轻彩対Snze6s(h72)八Em2

3、隘血国S13関戾 aa 占 S 二江 2三+2*5)nsumx(+l 2、oe,rl (+2WO5DU®endRlMXHlwr 或 ss FM 加游®s?pot(s)% ho-d on%PO(R5"% ho_d offmeihod，ramniyh。nes(2-2)八Y(2Trf7J?All.rw(YWY茨TFA()+ASV»三72=治强鲨召T23彈曲2%pot(©% ho-d 03w% po_m% ho_d off s2Hzeros?72x蛹淞文案fori=1:12 海节项S2的赋值S2(i:12Xi+12*5)jsum(x(i+12.*IO:

4、5)-T2(i+12.*0:5)/6; endH2=xT2S2； %协財R2及绘图 % plot(S2);% hold on% plot(R2；r,);% hold off%forecastfcT=A( 1 )+A(2).*73:84I； %趋势项预测fcS=S2(1:12);%季节项预测fcx=fcT+fcS;%死亡人数预测function asb,sig2arma2_4(rsk)%功能:根据自协方差函数列求ARMA(2,2)模型%输入:卜自协方差函数列，k矩阵的阶数%输出:a-模型的数值项系数；b-噪声项系数;sig2礫声项的方差%2021-4-17,Designed by luli. a

5、temp=r(3) r(2);r(4) r(3)r(4)：r(5); %系数 a 的计算a=-1 ;atemp(1 ):atemp(2)j;ry=zeros(1,3); %变换后的MA(2)模型的自协方差函数列梢彩文档foir-l -3ryflra丄srFl) r(i+2)=3ax(3>)sr(+)m(4)3ax(3L) endRuzerospki)；RmTry; aho1Q乂 cxlok omaHzeros(2>八 oma (二TrysjomauMrrys oma(2lrry (3)一 r2 丄ry(2)Tyw)k gmaHzerosk) o二1汶?!*gmawrRfirbs (

6、工4才endendpr(oma/gmaroma Sig2ury (二，c.mci b=(r2AP<cusig? Vfunction gammak=arma2 5(L,Max)精彩文档%功能:计算ARMA(2,2)的自协方差因数列%输入心变量的计算下标;Max-i+W自协方差的最大次数(默认 10000)% 输出:gamak<1»effiK%2021 -4-17> Designed by lulla1 =0.0894;%初始参数的赋值a2=-0.6265; % bO=1;b1h0.3334;b2=0.8158;sig2=4.0119； psy=eye(1Max); %

7、wold 系数的初值psy(2)=b1+a1;psy (3)=b2+a24a1 *psy ；for j=4:Max%wo Id系数列的计算psy( j)=a1 *psy0-1)+a2*psyG-2);end自协方差函数的计算gammak=sig2*sum(psy(1:Max-L).*psy(L+1:Max); %fun ctio n v=MA2(r1 ,k)%功能:根据自协方差函数求MA(2)模型%输入:门自协方差因数序矩阵的阶数 %输出M模型的参数组合%2021-4-17 Designed by luliFzeroscl>去=O芦 口辿！ i ;F?汁 gYrliAxo 1o0=CX1

8、s;Omauze63(2>= oma?lnn (2)j oma?2Trl (3)； oma(2lTn(3)i r2 占二 25(3 汙gmauzeH3s(k-k)jorlr'-l -ksegma 131orlr-pkgma(LTr(abs (工)+才endendpro3a/gmaoma； Sig2=r1(f prp b«r2>卫6)/sigw vx片b- 2g2K £爵圧迓禅 funaonITAR2simll-i nk(mNMJOWJheia)农常>RS3翟oc J Hl獅腿文案%输出：HH的第一行是样本数向最N；H的第二行是佶计的样本均值Avem

9、u;H的第三行是估计的噪%声项均值Avee;H的第四行是估计的样本方差Stdmu ;H的第五行是估计的噪声项方差 Stde.%2021 -5-1,Designed by luli.format short gH=zeros(5je ngth(N);%各参数的初始for j=1 :le ngth(N)%不同观測数的循环模拟Avemu=zeros(1 ,le ngth(N);Avee=zeros(1,le ngth(N);Stdmu=zeros(1 Je ngth(N);Stde=zeros(1,le ngth(N)；.for k=1:M%屡次模拟计算Y=zeros(1,m+N(j);terr=z

10、eros(1 E 十 N(j)；mu=zeros(1,M);%序列值的计郭e=zeros(1tM);for l=3:m+N(j)terr(l)=ra ndn;Yl)-2Vow*cos(theta)*Y(l-1)-rowA2*Y(l-2)+terr(l);end%各参数的赋值 l f w »* 轻影文档XHY?m4pm+N0) errHs-Iru -m+Ner mu(kTmea n( x)i e(knmea n( err).if MlAvemu(TmLGAveeolle ssmuuTp srdeeHoie_serlJ % A ve mAus 匸 m3u) * M 八AveejTsum(

11、e=M 八srdrnuuHsqmsum(muAvemu§A2)、(Ma= sfdeoTsqrt(su3(?Avee( j)>2)、？) j endHETr N( j)>vemu(j5vee( DcoidmuesQlegi6?Endfun cfo'nTTsimu-i nk3zM6WJh2.a)HHZe6s(5Je ngthN)Kfor j=1 :le ngth(N)for k=1:MY=zeros(1,m-FN(j);terr=zeros(1fm+N( );X=zeros(1,N(j);err=zeros(1,N(j);for l=3:m+N(j)terr(l)=ra ndn;Y(l)=2*row*cos(theta)*Y(l-1)-rowA2*Y(l-2)+terr(l); endX=Y(1,m4-1:m+N);err=terr(1 ,m+1 :m+N);mu(k)=mea n( X);e(k)=mea n( err);endif M=1Av