均匀试验设计程序(1).m

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

clear 
a=[200,26,2800,8
350,50,2600,7
500,26,3000,5
650,50,2700,4
800,18,2500,9
950,42,2900,7
1100,18,2600,6
1250,42,3000,4
1400,10,2800,9
1550,34,2500,8
1700,10,2900,6
1900,34,2700,5]
y=[0.151,0.113,0.199,0.116,0.091,0.142,0.099,0.135,0.128,0.029,0.116,0.016]'
n=12
m=4

t=2*m+1
x=zeros(n,t);
for i=1:n
    x(i,1)=1;
    for j=1:m
        x(i,j+1)=a(i,j);
        x(i,j+1+m)=a(i,j)^2;
    end
end
b=regress(y,x)
yc=x*b
y1=0;
y2=0;
for i=1:n
    y1=y1+y(i);
    y2=y2+y(i)*y(i);
end
lyy=y2-y1*y1/n;
bnew=zeros(t,1);
for j=1:t
    x1=0;
    x2=0;
    for i=1:n
        x1=x1+x(i,j);
        x2=x2+x(i,j)*x(i,j);
    end
    lxx=x2-x1*x1/n;
    bnew(j)=b(j)*sqrt(lxx/lyy);
end
plot(y,yc,'+')
hold on
linex=0:0.2/10:0.2
liney=linex
plot(liney,linex,'-')