给我一段Matlab代码，功能是实现灰色预测模型GM（1,1）

以下是一个简单的Matlab代码，用于实现灰色预测模型GM（1,1）：

% 输入矩阵x为原始数据，n为数据长度
% 输出y为预测结果，a为灰色预测模型GM（1,1）的参数

function [y,a] = grey_model(x,n)

% 累加生成数列
for i = 2:n
x1(i) = x1(i-1) + x(i);
end

% 紧邻均值生成数列
for i = 2:n
z(i) = (x(i) + x(i-1)) / 2;
end

% 灰色预测模型GM（1,1）求解
B = [-z(2:n)', ones(n-1,1)];
Y = x(2:n)';
a = B\Y;
u = x(1) - (a(2)/a(1));
y(1) = x(1);
for i = 2:n
y(i) = (x(1)-u)*exp(-a(1)*(i-1)) + u + a(2)/a(1);
end

end

注意：此代码仅供参考，实际应用中需要根据具体情况进行修改和优化。