灰色ga-bp代码MATLAB

以下是一份简单的灰色神经网络（GM(1,1)）的MATLAB代码实现：

% GM(1,1)灰色神经网络算法
function [y_pred, coef] = gm11(x)
    n = length(x);
    % 累加生成序列
    X1 = cumsum(x);
    % 紧邻均值生成序列
    for i = 2:n
        X2(i) = (X1(i-1) + X1(i)) / 2;
    end
    % 构造数据矩阵B和数据向量Y
    B = [-X2(2:n)', ones(n-1,1)];
    Y = x(2:end)';
    % 计算GM(1,1)微分方程的参数a和u
    coef = pinv(B) * Y;
    a = coef(1);
    u = coef(2);
    % 对原始数据进行预测
    y_pred = zeros(1,n);
    y_pred(1) = x(1);
    for i = 2:n
        y_pred(i) = (x(1) - u/a) * exp(-a*(i-1)) + u/a;
    end
    % 可视化预测结果
    plot(1:n, x, 'b-o', 1:n, y_pred, 'r-*');
    legend('真实值', '预测值');
end

代码中的`x`为输入的原始数据序列，输出为预测结果序列`y_pred`和GM(1,1)微分方程的参数`coef`（包括系数`a`和截距`u`）。代码中还包括可视化预测结果的部分，可以直观地观察预测效果。