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