正交试验设计程序.m

"试验设计"是很多工程人员尤其是化学工作者经常要面对的话题

x=[1	1	1		1
1	1	2	2
1	2	1	2
1	2	2	1
2	1	1	2
2	1	2	1
2	2	1	1
2	2	2	2]
y=[56.5
78.9
57.2
61.8
88.9
93.5
69.9
92.3]
x1=x2fx(x)
b=regress(y,x1)
yc=x1*b
plot(y,yc,'+')
hold on
linex=50:100/100:100
liney=linex
plot(liney,linex,'-')
xa=[2 1 2 2];
ya=95.4
y=[y;ya]
x1=[x;xa]
x1=x2fx(x1)
b=regress(y,x1)
yc=x1*b
plot(y,yc,'.')
hold on
linex=50:100/100:100
liney=linex
plot(liney,linex,'-')