ssimplex.m

一个单纯形法的程序 可以显示出最优解 最佳值 和迭代次数

function [sol,val,kk]=ssimplex(A,N)
[mA,nA]=size(A);
kk=0;     %  迭代次数
flag=1;
while flag
    kk=kk+1;
    if A(mA,:)<=0      %  已找到最优解
        flag=0;
        sol=zeros(1,nA-1);
        for i=1:mA-1
            sol(N(i))=A(i,nA);
        end
        val=-A(mA,nA);
    else
        for i=1:nA-1
            if A(mA,i)>0&A(1:mA-1,i)<=0    %  问题有无界解
                disp('have infinite solution!');
                flag=0;
                break;
            end
        end
        if flag        %    还不是最优表，进行转轴运算
            temp=0;
            for i=1:nA-1
                if A(mA,i)>temp
                    temp=A(mA,i);
                    inb=i;   %   进基变量的下标
                end
            end
            sita=zeros(1,mA-1);
            for i=1:mA-1
                if A(i,inb)>0
                    sita(i)=A(i,nA)/A(i,inb);
                end
            end
            temp=inf;
            for i=1:mA-1
                if sita(i)>0&sita(i)<temp
                    temp=sita(i);
                    outb=i;   %  出基变量下标
                end
            end
            %  以下更新N
            for i=1:mA-1
                if i==outb
                    N(i)=inb;
                end
            end
            %  以下进行转轴运算
            A(outb,:)=A(outb,:)/A(outb,inb);
            for i=1:mA
                if i~=outb
                    A(i,:)=A(i,:)-A(outb,:)*A(i,inb);
                end
            end
        end
    end
end