灰色模型Matlab程序.docVIP

%by allen @ 红嘴海鸥 %灰色模型预测是在数据不呈现一定规律下可以采取的一种建模和预测方法，其预测数据与原始数据存在一定的规律相似性 %下面程序是灰色模型GM(1,1)程序二次拟合和等维新陈代谢改进预测程序,matlab6.5 ,使用本程序请注明，程序存储为gm1.m %x =[5999,5903,5848,5700,7884];gm1(x); 测试数据 %二次拟合预测GM(1,1)模型function gmcal=gm1(x)sizexd2 = size(x,2);%求数组长度 k=0;for y1=xk=k+1;ifk1x1(k)=x1(k-1)+x(k);%累加生成z1(k-1)=-0.5*(x1(k)+x1(k-1)); %z1维数减1，用于计算Byn1(k-1)=x(k);elsex1(k)=x(k);endend%x1,z1,k,yn1 sizez1=size(z1,2);%size(yn1);z2 = z1;z3 = ones(1,sizez1); YN = yn1; %转置%YN B=[z2 z3];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 %输出预测的 a u 和 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;for y2 = x1k2 = k2 +1;if k2 k elseze1(k2) =exp(-(k2-1)*afor); endend%ze1 sizeze1 = size(ze1,2);z4 = ones(1,sizeze1);G=[ze1 z4];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,+,leftbra,num2str(Bval),rightbra)%输出时间响应方程 nfinal = sizexd2-1 + 1;%决定预测的步骤数5 这个步骤可以通过函数传入 %nfinal = sizexd2 - 1 + 1;%预测的步骤数 1 for k3=1:nfinalx3fcast(k3)= constant1*exp(afor1*k3)+ua;end%x3fcast%一次拟合累加值 for k31=nfinal:-1:0ifk311x31fcast(k31+1) = x3fcast(k31)-x3fcast(k31-1);elseif k310x31fcast(k31+1) = x3fcast(k31)-x(1);elsex31fcast(k31+1) = x(1);endendendx31fcast%一次拟合预测值 for k4=1:nfinalx4fcast(k4)= Aval*exp(afor1*k4)+Bval;end%x4fcast for k41=nfinal:-1:0ifk411x41fcast(k41+1) = x4fcast(k41)-x4fcast(k41-1);elseif k410x41fcast(k41+1) = x4fcast(k41)-x(1);elsex41fcast(k41+1) = x(1);endendendx41fcast,x%二次拟合预测值 %***精度检验pC************//////////////////////////////////k5 = 0;for y5 = xk5 = k5 +1;i