danchun.m

来自「学会MATLAB编程； 学会用单纯型法求解线性规划； 学会用对偶单纯型法求解」· M 代码 · 共 31 行

function [y,A]=danchun(A,x,y)
[m,n]=size(A);
if min(A(1,1:n-1))<0
    flag=0;
else
    flag=1;
end
while flag==0
    [h1,j]=min(A(1,1:n-1));
    for p=2:m
        if A(p,j)<=0 | A(p,n)==0
            q(p-1)=inf;
        else
            q(p-1)=A(p,n)./A(p,j);
        end
    end
    [h2,i]=min(q);
    y(i)=x(j);
    i=i+1;
    A(i,:)=A(i,:)./A(i,j);
    for k=1:m
        if k~=i
            A(k,:)=A(k,:)+(-A(k,j)).*A(i,:);
        end
    end
    if min(A(1,1:n-1))<0
        flag=0;
    else
        flag=1;
    end
end