matlab求序列输出相应,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软件做这些有关时间序列的，简单的就不要再浪费时间了