用MATLAB写一段灰色预测代码

由于灰色预测需要对原始数据进行处理，因此需要先用MATLAB对原始数据进行处理，然后再进行灰色预测。

以下是一段MATLAB灰色预测的代码示例：

% 原始数据
data = [10, 15, 20, 25, 30, 35, 40];

% 累加生成序列
data_cumsum = cumsum(data);

% 累加平均生成序列
data_cummean = data_cumsum ./ [1:length(data)];

% 灰色预测
alpha = 0.5; % 灰色预测中的发展系数
data_predict = zeros(size(data));
data_predict(1) = data(1);
for i = 2:length(data)
data_predict(i) = (data(1) - data_cummean(i-1)) * exp(-alpha*(i-1)) + data_cummean(i-1);
end

% 绘制原始数据和预测数据的图像
figure
plot([1:length(data)], data, '-o', 'LineWidth', 2)
hold on
plot([1:length(data)], data_predict, '-x', 'LineWidth', 2)
legend('原始数据', '预测数据')
xlabel('时间')
ylabel('数据值')
title('灰色预测')