- %建模数据集,123×8
- mg_data=[];
- %多元线性回归建模
- n=length(mg_data);
- X=[ones(123,1),mg_data(:,1:8)];
- %
- %定义期望输出Y(主导变量:出池时的淀粉量)
- Y=mg_data(:,8);
- %求取回归系数距阵A
- A=XY;
- %估计值yy
- yy=X*A;
- %绝对误差(残差)
- e=Y-yy;
- %定义误差性能函数,评价预测模型的拟合能力
- YSSE=e'*e %误差平方和函数 YSSE =
- YMSE=YSSE/n %均方误差函数 YMSE =
- %绘制出实际输出与预测输出的拟合曲线
- figure;
- plot(g_data0(:,8),'r');
- hold on;
- plot(yy,'--');
- legend('化验值','估计值');
- title('多元线性回归模型输出拟合曲线');
- xlabel('输入样本点');
- ylabel('淀粉利用率');