线性规划单纯形法matlab解法(共4页)

1、精选优质文档-倾情为你奉上线性规划单纯形法matlab解法%单纯形法matlab程序-ssimplex% 求解标准型线性规划:max c*x; s.t. A*x=b; x>=0% 本函数中的A是单纯初始表，包括:最后一行是初始的检验数，最后一列是资源向量b% N是初始的基变量的下标% 输出变量sol是最优解, 其中松弛变量（或剩余变量）可能不为0% 输出变量val是最优目标值，kk是迭代次数% 例：max 2*x1+3*x2% s.t. x1+2*x2<=8%      4*x1<=16%   &

2、#160;  4*x2<=12%      x1,x2>=0% 加入松驰变量，化为标准型，得到% A=1 2 1 0 0 8; %    4 0 0 1 0 16; %    0 4 0 0 1 12; %    2 3 0 0 0 0;% N=3 4 5;% sol,val,kk=ssimplex(A,N)% 然后执行 sol,val,kk=ssimplex(A,N)就可以了。function sol,val,kk=ssimplex

3、(A,N) mA,nA=size(A); kk=0; % 迭代次数flag=1; while flagkk=kk+1;    if A(mA,:)<=0 % 已找到最优解        flag=0;        sol=zeros(1,nA-1);        for i=1:mA-1     &

4、#160;      sol(N(i)=A(i,nA);        end        val=-A(mA,nA);    else        for i=1:nA-1           

5、 if A(mA,i)>0&A(1:mA-1,i)<=0 % 问题有无界解                disp('have infinite solution!');                flag=0;    

6、            break;            end        end        if flag % 还不是最优表，进行转轴运算       

7、60;    temp=0;            for i=1:nA-1                if A(mA,i)>temp              

8、;      temp=A(mA,i);                    inb=i; % 进基变量的下标                end     

9、       end            sita=zeros(1,mA-1);            for i=1:mA-1                if A

10、(i,inb)>0                    sita(i)=A(i,nA)/A(i,inb);                end         

11、   end            temp=inf;            for i=1:mA-1                if sita(i)>0&sita(i)<temp

12、60;                   temp=sita(i);                    outb=i; % 出基变量下标       &#

13、160;        end            end            % 以下更新N            for i=1:mA-1    

14、60;           if i=outb                    N(i)=inb;                end 

15、60;          end            % 以下进行转轴运算            A(outb,:)=A(outb,:)/A(outb,inb);            for i=1:mA                if i=outb                    A(i,:)=A(i,:)-A(outb,:)*A(