单纯形法的MATLAB实现

单纯形函数定义如下:

function [y,T]=Optimize(A,b,c)

%作者:MOK
%创作时间:2018/10/10
%A为系数矩阵
%b为常数约束矩阵
%c为目标函数系数矩阵
%y为最优函数值,T为最优化时的单纯形表
%以上条件默认优化问题是标准形式的


l=size(A,2);
r=rank(A);

Choose=nchoosek(1:l,r);
for i=size(Choose,1):-1:1
    choose=Choose(i,:);
    if (det(A(:,choose))~=0)
        B=A(:,choose);
        tem=A;
        tem(:,choose)=[];
        N=tem;
        cB=c(choose);
        tem=c;
        tem(choose)=[];
        cN=tem;
        break
     end
end

T=[(cB)*inv(B)*b,zeros(1,r),cB*(B)^(-1)*N-cN;(B^-1)*b,diag(ones(r,1)),(B^-1)*N];


flag=0;
while(any(T(1,2:end)>0))

frow=T(1,:);
[~,loc1]=max(frow);
fcolumn=T(2:end,loc1);
base=T(2:end,1);
cita=base./fcolumn;
[~,loc2]=min(cita);
loc2=loc2+1;
T(loc2,:)=T(loc2,:)/T(loc2,loc1);
%if (all(T(2:end,loc1)<=0)&&(T(1,loc1)>0))
%    flag=1;
%    break
%end

tem=T(2:end,loc1);
  • 11
    点赞
  • 109
    收藏
    觉得还不错? 一键收藏
  • 12
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 12
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值