二次指数平滑及其时间序列预测matlab代码

function[Y,S1,S2,a,b]=expsmooth2(Yt,alpha,t)%Yt:原时间序列；alpha:平滑系数；t:预测时长%Y:预测值；S1/S2:一次/二次指数平滑值；a/b:预测公式参数n=length(Yt);%计算一次指数平滑值S1(1)=Yt(1);fori=2:nS1(i)=alpha*Yt(i)+(1-alpha)*S1(i-1);end%计算二次指数平滑值S2(1)=S1(1);fori=2:nS2(i)=alpha*S1(i)+(1-alpha)*S2(i-1);end%计算参数a和bfori=1:na(i)=2*S1(i)-S2(i);b(i)=alpha/(1-alpha)*(S1(i)-S2(i));end%计算预测值Yfori=1:tY(i)=a(n)+b(n)*i;end%绘图plot(1:n,Yt,(n+1):(n+t),Y,'*');end