写一段四个自变量，一个因变量，一共286条数据的灰色预测模型和BP神经网络的复合模型MATLAB代码，并进行检验

以下是一个基于灰色预测模型和BP神经网络的复合模型的MATLAB代码：

% 数据
x1 = rand(1, 286); % 自变量1
x2 = rand(1, 286); % 自变量2
x3 = rand(1, 286); % 自变量3
x4 = rand(1, 286); % 自变量4
y = rand(1, 286); % 因变量

% 灰色预测模型
n = length(y);
x0 = cumsum(y); % 累加生成序列
z1 = (x0(1:n-1) + x0(2:n)) / 2; % GM(1,1)模型中的紧邻均值生成序列
B = [-z1', ones(n-1,1)];
Y = y(2:end)';
a = pinv(B) * Y; % 求解参数
k = a(1);
b = a(2);
y_predict_grey = zeros(1,n);
y_predict_grey(1) = y(1); % 初始值
for i = 2:n
    y_predict_grey(i) = (y(1) - b/k) * exp(-k*(i-1)) + b/k;
end

% BP神经网络模型
inputs = [x1; x2; x3; x4];
targets = y;
hiddenLayerSize = 10;
net = fitnet(hiddenLayerSize);
net = train(net,inputs,targets);
y_predict_bp = net(inputs);

% 复合模型
y_predict = 0.5 * y_predict_grey + 0.5 * y_predict_bp; % 简单均值法加权平均

% 检验
figure;
plot(y, 'b');
hold on;
plot(y_predict, 'r');
legend('真实值', '预测值');
title('灰色预测模型和BP神经网络的复合模型预测结果');
xlabel('数据编号');
ylabel('因变量值');
mse = sum((y_predict - y).^2) / n; % 均方误差
rmse = sqrt(mse); % 均方根误差
mape = sum(abs(y_predict - y) ./ y) / n * 100; % 平均绝对百分误差
disp(['均方误差：', num2str(mse)]);
disp(['均方根误差：', num2str(rmse)]);
disp(['平均绝对百分误差：', num2str(mape), '%']);

在此代码中，我们生成了4个自变量和1个因变量，共有286条数据。首先，我们使用灰色预测模型和BP神经网络模型对数据进行训练和预测。然后，我们使用简单均值法加权平均将两个模型的预测结果进行组合，生成最终的复合预测结果。最后，我们计算了预测结果的均方误差、均方根误差和平均绝对百分误差，并输出到命令窗口中。

请注意，此代码仅为示例，实际应用中可能需要根据数据的特点进行调整。