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

文章来源: yunyaniu.blog.csdn.net，作者：一个处女座的程序猿，版权归原作者所有，如需转载，请联系作者。

原文链接：yunyaniu.blog.csdn.net/article/details/80380384