灰色预测模型matlab程序精确版.docxVIP

灰色预测模型matlab程序精确 版 %x=[1019,1088,1324,1408,1601];gm1(x); 测 试数据 %二次拟合预测GM(1,1)模型 fun cti on gmcal=gm1(x) if nargin==0 x=[1019,1088,1324,1408,1601] end format long g sizex=le ngth(x); %求数组长度 k=0; for y1=x k=k+1; if k1 x1(k)=x1(k-1)+x(k); % 累加生成 z1(k-1)=-0.5*(x1(k)+x1(k-1)); %z1 维数减1,用于计算B yn 1(k-1)=x(k); else x1(k)=x(k); end end %x1,z1,k,y n1 sizez1=le ngth(zl); %size(y n1); z2 = z1; z3 = on es(1,sizez1); YN = yn 1; % 转置 %YN B=[z2 z3]; au0=i nv(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的值 con sta nt1 = x(1)-ua; afor1 = -afor; x1t1