MATLAB 灰色模型分析代码示例

% 输入数据 x = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]; y = [10.6, 14.2, 18.2, 22.4, 27.2, 32.4, 38.2, 44.4, 51.2, 58.4];

% 灰色预测模型 n = length(x); X = cumsum(y); Y = y; for i = 2:n Y(i) = (y(i) + y(i-1))/2; end B = [-X', ones(n,1)]; a = inv(B'*B)*B'*Y'; a1 = a(1); a2 = a(2); x0 = (Y(1)-a2/a1); y_predict = (y(1)-x0)exp(-a1(x-x(1))) + x0;

% 灰色预测结果可视化 plot(x, y, 'r', x, y_predict, 'b'); xlabel('时间'); ylabel('数据'); title('灰色预测模型'); legend('原始数据', '灰色预测结果');