还有一个多变量的灰色预测题(带原程序)。。。
近年来,许多文献从不 同的角度给出了非等间距灰色预测模型的改进,然而实 际的社会 、经济系统 中往往包含多个变量,各变量相互影响 、相互作用 。
为了解决多变量等间距原始数据的灰色预测问题 ,文献 分别讨论 了多变量等间距和非等间距预测模型的建模方法和应用 。多变量灰色模型 的参数估计 、模型预测以及模型的精度检验都需要一些 比较复杂的计算。
部分源程序
- %输入原始序列X0
- X0=[11.42
- 12.78
- 15.23
- 18.66
- 23.57
- 40
- 46];
- %对原始序列X0一次累加生成序列X1
- [n,m]=size(X0);
- for i=1:n c=0;
- for j=1:i
- c=X0(j,m)+c;
- end
- X1(i,m)=c;
- end
- disp('X1=');
- disp(X1);
- %计算数据矩阵B
- for i=2:n
- Z1(i,m)=(X1(i-1,m)+X1(i,m))/2;
- end
- disp('Z1=');
- disp(Z1);
- B=[-Z1 ones(n,1)];
- B(1,:)=[];
- disp('B=');
- disp(B);
- %计算数据矩阵Y及a,b的值
- X0(1,:)=[];
- Y=X0;
- disp('Y=');
- disp(Y);
- Q=inv(B'*B)*B'*Y;
- disp('Q=');
- disp(Q);