正规单纯形程序.m

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

n=input('please input n:')
p0=(sqrt(n+1)+n-1)/(n*sqrt(2))
q0=(sqrt(n+1)-1)/(n*sqrt(2))
x=zeros(n+1,n);
a=zeros(n,1);
t0=zeros(n,1);
te=zeros(n,1);
for i=1:n
    i,
    t0(i)=input('please input t0:')
    te(i)=input('please input te:')
end
for i=1:n
    x(1,i)=input('please input x1:')
end
p=zeros(n,1);
q=zeros(n,1);
for i=1:n
    a(i)=2*(te(i)-t0(i))/3
    p(i)=p0*a(i)
    q(i)=q0*a(i)
end
for i=2:n+1
    for j=1:n
        if j==i-1
            x(i,j)=x(1,j)+p(j)
        else
            x(i,j)=x(1,j)+q(j)
        end
    end
end
y=zeros(n+1,1);
for i=1:n+1
    y(i)=input('please input y:')
end
simp=[x y]
a=simpx(n,y,x)


    




.