灰色模型matlab程序(1).doc

《灰色模型matlab程序(1).doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、灰色模型matlab程序%alpha是包含α、μ值的参数矩阵；ago是预测后累加值矩阵；var是还原预测值矩阵；error是残差矩阵；c是后验差比值clc%清屏，以使结果独立显示formatlong;%设置计算精度x=[10.90,11.57,11.36,11.49,11.82];iflength(x(:,1))==1%对输入矩阵进行判断，如不是一维列矩阵，进行转置变换    x=x';endn=length(x);%取输入数据的样本量z=0;fori=1:n%计算累加值，并将值赋予矩阵be    z=z+x(i,:

2、);    be(i,:)=z;endfori=2:n%对原始数列平行移位    y(i-1,:)=x(i,:);endfori=1:n-1%计算数据矩阵B的第一列数据    c(i,:)=-0.5*(be(i,:)+be(i+1,:));endforj=1:n-1%计算数据矩阵B的第二列数据    e(j,:)=1;endfori=1:n-1%构造数据矩阵B    B(i,1)=c(i,:);    B(i,2)=e(i,:);endalpha=inv(B'*B)*B'*y;%计算参数矩阵fori=1:n+1%计算

3、数据估计值的累加数列，如改为n+1为n+m可预测后m-1个值    ago(i,:)=(x(1,:)-alpha(2,:)/alpha(1,:))*exp(-alpha(1,:)*(i-1))+alpha(2,:)/alpha(1,:);endvar(1,:)=ago(1,:)fori=1:n%如改n为n+m-1，可预测后m-1个值    var(i+1,:)=ago(i+1,:)-ago(i,:);%估计值的累加数列的还原，并计算出下一预测值endfori=1:n    error(i,:)=var(i,:)-x(

4、i,:);%计算残差endc=std(error)/std(x)%调用统计工具箱的标准差函数计算后验差的比值c