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

[slgu]

请运行下列小程序：
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

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