灰色算法预测模型matlab程序.pdfVIP

灰色预测模型matlab程序(2007-05-3111:02:49)

：灰色模型gm(11)二次拟合matlab分类：滴

%byallen@红嘴海鸥

%灰色模型预测是在数据不呈现一定规律下可以采取的一种建模和预测方法，其预测数据与

原始数据存在一定的规律相似性

%下面程序是灰色模型GM(1,1)程序二次拟合和等维新陈代谢改进预测程序,matlab6.5,使

用本程序请注明，程序为gm1.m

%x=[5999,5903,5848,5700,7884];gm1(x);测试数据

%二次拟合预测GM(1,1)模型

functiongmcal=gm1(x)

sizexd2=size(x,2);

%求数组长度

k=0;

fory1=x

k=k+1;

ifk1

x1(k)=x1(k-1)+x(k);

%累加生成

z1(k-1)=-0.5*(x1(k)+x1(k-1));

%z1维数减1，用于计算B

yn1(k-1)=x(k);

else

x1(k)=x(k);

end

end

%x1,z1,k,yn1

sizez1=size(z1,2);

%size(yn1);

z2=z1;

z3=ones(1,sizez1);

YN=yn1;%转置

%YN

B=[z2z3];

au0=inv(B*B)*B*YN;

au=au0;

%B,au0,au

afor=au(1);

ufor=au(2);

ua=au(2)./au(1);

%afor,ufor,ua

%输出预测的au和u/a的值

constant1=x(1)-ua;

afor1=-afor;

x1t1=x1(t+1);

estr=exp;

tstr=t;

leftbra=(;

rightbra=);

%constant1,afor1,x1t1,estr,tstr,leftbra,rightbra

strcat(x1t1,=,num2str(constant1),estr,leftbra,num2str(afor1),tstr,rightbra,+

,leftbra,num2str(ua),rightbra)

%输出时间响应方程

%******************************************************

%二次拟合

k2=0;

fory2=x1

k2=k2+1;

ifk2k

else

ze1(k2)=exp(-(k2-1)*afor);

end

end

%ze1

sizeze1=size(ze1,2);

z4=ones(1,sizeze1);

G=[ze1z4];

X1=x1;

au20=inv(G*G)*G*X1;

au2=au20;

%z4,X1,G,au20

Aval=au2(1);

Bval=au2(2);

%Aval,Bval

%输出预测的A,B的值

strcat(x1t1,=,num2str(Aval),estr,leftbra,num2str(afor1),tstr,rightbra,+,lef

tbra,num2str(Bval),rightbra)

%输出时间响应方程

nfinal=sizexd2-1+1;

%决定预测的步骤数5这个步骤可以通过函数传入

%nfinal=sizexd2-1+1;

%预测的步骤数1

fork3=1