美赛（matlab自学）时间序列

确定性时间序列分析方法

平均移动法

例题

分别取N=4，N=5；

clc,clear
y=[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为移动平均的项数
for i =1:length(n)%由于n的取值不同,因此下面使用了细胞数组.
for j =1:m-n(i) +1
yhat{i}(j) =sum(y(j:j +n(i) -1))/n(i);
end 
y12(i) =yhat{i}(end);%提出12月的预测值
s(i) = sqrt( mean((y(n(i) +1:end) -yhat{i}(1:end-1)).^2));%求预测的标准误差
end
y12, s%分别显示两种方法的预测值和预测的标准误差

指数平滑法

一次指数平滑法

例题

clc,clear
yt = load('dianqi .txt');%实际销售额数据以列向量的方式存放在纯文本文件中
n= length(yt); alpha=[0.2 0.5 0.8]; 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,:);
end
yhat
err = sqrt(mean( ( repmat(yt,1 ,m) -yhat).^2))
xlswrite( 'dianqi .xls' ,yhat)%把预测数据写到Excel文件,准备在Word表格中使用
yhat1988 =alpha*yt(n) +(1 -alpha). *yhat(n,:)

二次指数平滑法

三次指数平滑法

差分指数平滑法

一阶差分指数平滑法

模板代码

clc,clear
yt = load( 'ranliao.txt');%实际燃料消耗量数据以列向量的方式存放在纯文本文件中.
n= length(yt); alpha=0.4;
dyt =diff(yt); %求yt的一阶向前差分
dyt =[0;dyt]; %这里使用的是一阶向后差分，加“0”补位
dyhat(2) =dyt(2); %指数平滑值的初始值
for i =2:n
dyhat(i +1) =alpha *dyt(i) +(1 - alpha) * dyhat(i);
end
for i =1:n
yhat(i +1) =dyhat(i +1) +yt(i);  
end
yhat
xlswrite( 'ranliao .xls' ,[yt ,dyt])
xlswrite( 'ranliao .x1s',[ dyhat' ,yhat'], 'Sheet1','C1')

二阶差分指数平滑法

具有季节性特点的时间序列预测