時間序列MATLAB程序

1、時間序列移動平均法例I拄金業(yè)I rII 口份的悄俾收入時間序列朗衣I禮試用一報怖單滑動平 均注揪測E 12月份的梢密收入.農(nóng)1上禮稱害收人12j4L3GS7t 6tkib.生6 ME7C5. 772. Dmtdti9IUHH1&. 1M2 7業(yè)a鼻MIB 110- 7解；仃刑1U V = 4+ Ar = 5的價朗瓷式禽=甘兒嚴(yán)先r冷r當(dāng)N4時.ftfflffii-f? =993.6料利的標(biāo)準(zhǔn)謀差為11Jf-iAR=5時* ffiMftj4伍泅的乩準(zhǔn)址丫力r I*iHT姑黑總明./V = 4時俺則的杠朮茂差較小，所口逸取N =4.鶴測站JJ份的 椚裁我人為怕J機(jī)clc,cleary=533、8

2、574、6 606、9 649、8 705、1 772、0 816、4 892、7 963、9 10151 1102、7;m=le ngth(y);n=4,5;%n為移動平均的項數(shù)for i=1:length(n)%由于n的取值不同，下面使用了細(xì)胞數(shù)組for j=1:m-n(i)+1yhati(j)=sum(y(j：j+n( i)-1)/n(i);endy12(i)=yhati(end);%提岀第12月份的預(yù)測值s(i)=sqrt(mean(y(n(i)+1:end)-yhati(1:end-1)、A2);%求預(yù)測的標(biāo)準(zhǔn)誤差endy12, s%分別顯示兩種方法的預(yù)測值與預(yù)測的標(biāo)準(zhǔn)誤差指數(shù)平滑一

3、次指數(shù)平滑*ih年某種電器銷竹幽山!占電所示.試廉測lK年詒電器銷杵解 劇II描數(shù)平滯叢片分別lUo = 0.2,0,50,a行計寡創(chuàng)貽值-212厶“即2亍產(chǎn)J5f九Ji和-cc)元計尊齊期協(xié)麗沱，列于衣4中”瑕4莫榊電器常昔*肚拆蝕半淆懺測曲訃算決 t華他 片尤)t側(cè) til ifa * OJ師越悅xa = 0-S禹睡QltN76SflSISI51山廠525IISSOS50.2rmJ175I.CM.2SSIM應(yīng)79AS15OJ349. B47.yl49刃腫沖06和孫WJ64M50.11W訂.2S1S149.6J:、is 11際；翱P 49.95149J50.454KJ7 M44礬衛(wèi)3W410

4、5246.47船梅輛LI514S.774畑加9bIW1：盟49J2W 12孫49程序：clc,clearyt=load(dianqi、txt) ;%讀取 dianqi、txt 的數(shù)據(jù),n=length(yt);%求 yt 的長度alpha=O 、2 0、5 0、8 ;%輸入 a 的值m=le ngth(alpha)yhat(1,1:m)=(yt(1)+yt(2)/2;%求第一個預(yù)測值索引for i=2: nyhat(i,:)=alpha*yt(i-1)+(1-alpha) end、*yhat(i-1,:);yhat%求預(yù)測值err=sqrt(mea n(repmat(yt,1,m)-yhat)

5、 xlswrite(dianqi、xls,yhat) %yhat1988=alpha*yt( n)+(1-alpha)、A2)%求預(yù)測的標(biāo)準(zhǔn)誤差把預(yù)測數(shù)據(jù)寫到Excel文件，準(zhǔn)備在word表格中使用、*yhat(n,:)%求 1988 的預(yù)測值二次指數(shù)平滑解Kia = 03.削晞值顯厠顯澤極取岸列的片頊敵fib EH5；l, = XJ，=&76-S?/ = 3S23.L= 3032.6由藝式【冷幾Fftf = 21時4 = i?y/=4ou）2H24Nr= ?i ： ritjjV/j 用力九.=4013,7+ 210 J4m預(yù)測19K6訃ftl I9S7夠的發(fā)電直量治 單也 億度） BH421

6、19S-4434.19卻J承*屈的枳擬他町將成（W 優(yōu)入口純乜協(xié)蟆星（1）.y = (25；-5(u,) + 各“)1 or即令21乂“20由企式Q0町求出各期的HUflfll.計奪端果見糞心tt = 0.3,擁齢値耳利顯澤肚取序畀的育fllttflh B|I5；1, = V=676 -汁ff5/5；n列嚴(yán)農(nóng)酩御外S?/ = 3523,1.* = 3032.6由舍式 19).可帑2$1 -1d2l =15-5 = 4013.7,垢=(卑卉 h21QQ41 -cr兀* =4013.7+2IL24m預(yù)測19K6年 1M7年的KA量為（僅 億度） n=length(yt);=Kt = 4223.95

7、-4434.19為f承4劇的模KltfL可轄式（19）代入fl純越孫模壁（IBK并專的壬1.程序clc,clearyt=load（的方式存放在純文本文件中n tz劍魔fadian 、txt ）;%原始發(fā)電總量數(shù)據(jù)以列向量alpha=0 、 3;st1(1)=yt(1);st2(1)=yt(1);for i=2:nst1(i)=alpha*yt(i)+(1-alpha)*st1(i-1);st2(i)=alpha*st1(i)+(1-alpha)*st2(i-1); endxlswrite( fadian 、 xls ,st1,st2)%把數(shù)據(jù)寫入表單 Sheet1 中的前兩列at=2*st1-

8、st2;bt=alpha/(1-alpha)*(st1-st2);yhat=at+bt;%最后的一個分量為 1986 年的預(yù)測值xlswrite(fadian、 xls,yhat, Sheet1, C2 )%把預(yù)測值寫入第 3 列str= C ,int2str(n+2);%準(zhǔn)備寫 1987 年預(yù)測值位置的字符串xlswrite(fadian、xls,at(n)+2*bt(n),Sheet1,str) %把 1987 年預(yù)測值寫到相應(yīng)位置yt=load(劍魔 fadian、 txt ); %原始發(fā)電總量數(shù)據(jù)以列向量的方式存放在純文本文件中n=length(yt); alpha=0 、 3;st1(1)=yt(1); st2(1)=yt(1);for i=2:nst1(i)=alpha*yt(i)+(1-alpha)*st1(i-1); st2(i)=alpha*st1(i)+(1-alpha)*st2(i-1);endxlswrite( fadian 、 xls ,st1,st2)%把數(shù)據(jù)寫入表單 Sheet1 中的前兩列at=2*st1-st2;bt=alpha/(1-alpha)*(st1-st2);yhat=at+bt;%最后的一個分量為 1986 年的預(yù)測值xlswrite(fadian、 x