用matlab做灰色预测的代码.docVIP

用matlab做灰色预测的代码 1 程序代码: %普通的灰色预测GM1.m clear X=input(请输入原始数据:,s);%原始数据(可以多行，每一行为一类原始数据，即可多类原始数据) X=str2num(X); [m1 m2]=size(X);%m1和m2分别表示X的行数和列数 k0=input(请输入所要预测的阶数:); %GM(1,1)模型 for i=1:m1 n=i; x0=X(i,:);%将原始数据X中的第i行数据赋给x0，即取出一类原始数据 disp(1.原始数据:); Y=; for z=1:m2 Y=strcat(Y,(,num2str(x0(z)),)); end disp(Y); % 1. 利用一次累加(1-AGO)生成新数列 E=triu(ones(m2));%E表示元素为1右上三角阵 x1=x0*E;%对原始数据进行一次累加(1-AGO)生成新数列x1 disp(2.一次累加(1-AGO)生成的数据:); Y=; for z=1:m2 Y=strcat(Y,(,num2str(x1(z)),)); end disp(Y); % 2. 计算出发展系数a，灰作用量u b1=x1;b1(1)=[]; b2=x1;b2(m2)=[]; b=-0.5*(b1+b2);%生成B的第一行 B=[b;ones