响应面法matlab代码

以下是一个简单的响应面法的 Matlab 代码示例：

% 定义实验设计矩阵
X = [1 1 -1 -1;
1 -1 1 -1;
1 1 1 1;
1 -1 -1 1;
1 0 0 0;
1 0 0 0;
1 0 0 0;
1 0 0 0;
1 0 0 0;
1 0 0 0;
1 0 0 0;
1 0 0 0;
1 0 0 0;
1 0 0 0;
1 0 0 0;
1 0 0 0];

% 定义响应变量
y = [60 50 70 65 75 70 80 85 90 95 100 105 110 115 120 125];

% 计算回归系数
b = inv(X'*X)*X'*y;

% 定义预测变量
x1 = linspace(-1,1,101);
x2 = linspace(-1,1,101);

% 计算预测响应
[Y1,Y2] = meshgrid(x1,x2);
Y = b(1) + b(2)*Y1 + b(3)*Y2 + b(4)*Y1.*Y2;

% 绘制响应面
surf(Y1,Y2,Y);
xlabel('x_1');
ylabel('x_2');
zlabel('y');
title('响应面');