MATLAB - + - 灰色预测程序，数学建模[1]

B=[B' ones(n-1,1)]; %生成矩阵 Y

for i=1:(n-1);

y(i)=x_ori(i+1); ed

Y=y'; %计算系数 a=au(1) u=au(2)

au=(inv(B'*B))*(B'*Y); %--------------------------------------------------------

符号

coef1=au(2)/au(1); coef2=x_or (1)-coef1; co3=0-au(1);

costr1=nm2str(coef1);

costr2=numstr(abs(coef2)); costr3=ntr(coef3);

eq=strcat(ctr1,'+',costr2,'e^',costr3,'*(t-1))');

for t=1:(n+predict)

mcv(t)=co1+coef2*exp(coef3*(t-1)); end

x_mcv0=diff(mcv);

x_mcve=[x_orig(1) x_mcv0]

x_c_error=x_orig_n-x_mcv;

x_errr=mn(abs(x_c_error./x_orig_n));

if x_error>0.2 disp('del disqualification!'); elseif x_error>0.1

dip('model check out'); disp('model is perfect!');

end plot(1:n,x_orig,'o',1:n+predict,x_mcve);

p=x_mcve(end-predict+1:end); 始数列的点

xlabel('年份(从第一个数据年份起)'); ylabel('产水量（万吨）'); tie('灰度模型 GM(1,1)'); grid on y=eq;

e=x_error;

p=x_mcve(end-predict+1:end);

%把huise模型公式转换成 %计算每一个值 %输出图形中的各点 %相对误差的均值 %画出预测模型和初