运筹学单纯形法matlab程序
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
function [xx,fm]=myprgmh(m,n,A,b,c)
B0=A(:,1:m);
cb=c(:,1:m);
xx=1:n;
sgm=c-cb*B0^-1*A;
h=-1;
sta=ones(m,1);
for i=m+1:n
if sgm(i)>0
h=1;
end
end
while h>0
[msg,mk]=max(sgm);
for i=1:m
sta(i)=b(i)/A(i,mk);
end
[mst,mr]=min(sta);
zy=A(mr,mk);
for i=1:m
if i==mr
for j=1:n
A(i,j)=A(i,j)/zy;
end
b(i)=b(i)/zy;
end
end
for i=1:m
if i~=mr
for j=1:n
A(i,j)=A(i,j)-A(i,mk)*A(mr,j);
end
b(i)=b(i)-A(i,mk)*b(mr);
end
end
B1=A(:,1:m);
cb(mr)=c(mk);
xx(mr)=mk;
sgm=c-cb*B1*A;
for i=m+1:n
if sgm(i)>0
h=1;
end
end end
fm=c*xx;