[求助]正交实验的参数估计怎么做？
 

我用正交试验设计了一个3个参数的实验，获得了9组数据，看能不能用最小二乘法进行参数估计？

具体的实验数据如下：

1 1 1 1 0.9
2 1 2 2 1.2
3 1 3 3 1.0
4 2 1 2 1.1
5 2 2 3 1.0
6 2 3 1 0.4
7 3 1 3 1.4
8 3 2 1 0.9
9 3 3 2 1.1


第一列是实验序号，前面是水平，最后一列是数据！

看各位大侠有没有这方面的经验，具体水平下的指标我就不列入了！

希望大侠们能够给点最高指示！

谢谢！

回复: [求助]正交实验的参数估计怎么做？
 

请运行下列小程序：
clear,clc
x1=[1 2 3 1 2 3 1 2 3]';%原第二列
x2=[1 2 3 2 3 1 3 1 2]';%原第三列
x3=x1.^2;
y=[.9, 1.2, 1, 1.1, 1, .4, 1.4, .9, 1.1]';
X=[ones(9,1),x2,x3];
b=X\y
SSy=var(y)*8
RSS=y'*y-b'*X'*y
MSe=RSS/(9-3)
R2=(SSy-RSS)/SSy
scatter3(x1,x2,y)
stem3(x1,x2,y,'filled')
str=num2str([1:9]');
text(x1,x2,y+.05,str,'fontsize',12)
hold on
[x11,x22]=meshgrid(0.5:.05:3.5);
y1=b(1)+b(2)*x22+b(3)*x11.^2;
surf(x11,x22,y1)
axis tight

%总体上数据点较少，自变数的作用不太明显，这种试验可能不大能说明问题。