Matlab：序列分析法MATLAB代码

Matlab：序列分析法MATLAB代码

---

输出结果

更新……

设计代码

###下面所有带代码中的n值需要以自己输入的数据为准###

1、简单一次滑动平均法预测MATLAB程序代码
y=[_______];

n=length(y);

Sum=cumsum(y);% 求累积和

mt=(Sum(4:n)-[0 temp(1:n-4)])/4;

y12=mt(end)

ythat=mt(1:end-1);%预测的值

fangcha=mean((y(5:n)-ythat).^2);

sigma=sqrt(fangcha)

2、一次指数平滑法预测MATLAB程序代码
alpha=____;%平滑常数；当时间数列相对平稳时，可取较大的a；当时间数列波动较大时，应取较小的a.

y=[________];
n=length(y)l

s1(1)=y(1);

for i=2:n

s1(i)=alpha*y(i)+(1-alpha)*s1(i-1);

end

yhat(n+1)=s1(end)

sigma=sqrt(mean((s1(1:end-1)-y(2:end)).^2))

3、二次指数平滑法预测MATLAB程序代码
clc,clear
alpha= _____ ;%平滑常数
y=[_______];
n=length(y);
s1(1)=y(1);
for i=2:n
s1(i)=alpha*y(i)+(1-alpha)*s1(i-1);
end
s2=y(1);
for i=2:n
s2(i)=alpha*s1(i)+(1-alpha)*s2(i-1);
end
an=2*s1(n)-s2(n) %n值需要修改，下同。
bn=alpha/(1-alpha)*(s1(n)-s2(n))
yhatn+1=an+bn
yhat(1)=y(1)
for i=2:n
yhat(i)=s1(i-1)+1/(1-alpha)*(s1(i-1)-s2(i-1));
end
temp=sum((yhat-y).^2);
sigma=sqrt(temp/(n-2))%标准误差

插一句，还是推荐SPSS软件做这些有关时间序列的，简单的就不要再浪费时间了