時間序列MATLAB程序

1、.時間序列移動平均法clc,cleary=533.8 574.6 606.9 649.8 705.1 772.0 816.4 892.7 963.9 1015.1 1102.7;m=length(y);n=4,5;%n 為移動平均的項數(shù)fori=1:length(n)%由于 n 的取值不同，下面使用了細胞數(shù)組forj=1:m-n(i)+1.yhati(j)=sum(y(j:j+n(i)-1)/n(i);endy12(i)=yhati(end);%提出第 12 月份的預測值s(i)=sqrt(mean(y(n(i)+1:end)-yhati(1:end-1).2);%求預測的標準誤差endy12

2、, s%分別顯示兩種方法的預測值和預測的標準誤差指數(shù)平滑一次指數(shù)平滑程序：clc,clearyt=load(dianqi.txt);%讀取 dianqi.txt的數(shù)據(jù)，n=length(yt);%求 yt的長度alpha=0.2 0.5 0.8;%輸入 a 的值m=length(alpha)yhat(1,1:m)=(yt(1)+yt(2)/2;%求第一個預測值索引fori=2:nyhat(i,:)=alpha*yt(i-1)+(1-alpha).*yhat(i-1,:);.endyhat%求預測值err=sqrt(mean(repmat(yt,1,m)-yhat).2)%求預測的標準誤差xls

3、write(dianqi.xls,yhat) %把預測數(shù)據(jù)寫到Excel文件，準備在word表格中使用yhat1988=alpha*yt(n)+(1-alpha).*yhat(n,:)%求 1988的預測值二次指數(shù)平滑.程序clc,clearyt=load(劍魔 fadian.txt);%原始發(fā)電總量數(shù)據(jù)以列向量的方式存放在純文本文件中n=length(yt);alpha=0.3;st1(1)=yt(1);st2(1)=yt(1);fori=2:nst1(i)=alpha*yt(i)+(1-alpha)*st1(i-1);st2(i)=alpha*st1(i)+(1-alpha)*st2(i-

4、1);endxlswrite( fadian.xls ,st1,st2) %把數(shù)據(jù)寫入表單 Sheet1 中的前兩列 at=2*st1-st2;bt=alpha/(1-alpha)*(st1-st2);yhat=at+bt;%最后的一個分量為1986 年的預測值xlswrite(fadian.xls,yhat,Sheet1, C2 ) %把預測值寫入第3 列str= C,int2str(n+2);%準備寫 1987年預測值位置的字符串xlswrite(fadian.xls,at(n)+2*bt(n),Sheet1,str) %把 1987年預測值寫到相應位置yt=load(劍魔 fadian.

5、txt);%原始發(fā)電總量數(shù)據(jù)以列向量的方式存放在純文本文件中n=length(yt);alpha=0.3;st1(1)=yt(1);.st2(1)=yt(1);fori=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 年的預測值xlswrite(fadian.xls,yhat,Sheet1, C2 ) %把預測值寫入第3 列str= C,int2str(n+2);%