時間序列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ù)for i=1:length(n) %由于n的取值不同，下面使用了細胞數(shù)組 for j=1:m-n(i)+1 yhati(j)=sum(y(j:j+n(i)-1)/n(i); end y12(i)=yhati(end); %提出第12月份的預測值 s(i)=sqrt(mean(y(n(i)+1:end)-yhati(1:end-1).2); %求預測的標準誤

2、差endy12, 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; %求第一個預測值 索引for i=2:n yhat(i,:)=alpha*yt(i-1)+(1-alpha).*yhat(i-1,:);endyhat %求預測值err=sqrt(mean(repmat(yt,1

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

4、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=&

5、#39;C',int2str(n+2); %準備寫1987年預測值位置的字符串xlswrite('fadian.xls',at(n)+2*bt(n),'Sheet1',str)%把1987年預測值寫到相應位置 yt=load('c:UsersasusDesktop劍魔fadian.txt'); %原始發(fā)電總量數(shù)據(jù)以列向量的方式存放在純文本文件中n=length(yt);alpha=0.3;st1(1)=yt(1); st2(1)=yt(1);for i=2:n st1(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); %準備寫1