线性规划单纯形法matlab解法_第1页
线性规划单纯形法matlab解法_第2页
线性规划单纯形法matlab解法_第3页
线性规划单纯形法matlab解法_第4页
免费预览已结束，剩余1页可下载查看

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

线性规划单纯形法matlab解法线性规划单纯形法matlab解法线性规划单纯形法matlab解法资料仅供参考文件编号：2022年4月线性规划单纯形法matlab解法版本号：A修改号：1页次：1.0审核：批准:发布日期：线性规划单纯形法matlab解法%单纯形法matlab程序-ssimplex%求解标准型线性规划:maxc*x;.A*x=b;x>=0%本函数中的A是单纯初始表，包括:最后一行是初始的检验数，最后一列是资源向量b%N是初始的基变量的下标%输出变量sol是最优解,其中松弛变量（或剩余变量）可能不为0%输出变量val是最优目标值，kk是迭代次数%例：max2*x1+3*x2%.x1+2*x2<=8%

4*x1<=16%

4*x2<=12%

x1,x2>=0%加入松驰变量，化为标准型，得到%A=[121008;%

4001016;%

0400112;%

230000];%N=[345];%[sol,val,kk]=ssimplex(A,N)%然后执行[sol,val,kk]=ssimplex(A,N)就可以了。function[sol,val,kk]=ssimplex(A,N)[mA,nA]=size(A);kk=0;%迭代次数flag=1;whileflagkk=kk+1;

ifA(mA,:)<=0%已找到最优解

flag=0;

sol=zeros(1,nA-1);

fori=1:mA-1

sol(N(i))=A(i,nA);

end

val=-A(mA,nA);

else

fori=1:nA-1

ifA(mA,i)>0&A(1:mA-1,i)<=0%问题有无界解

disp('haveinfinitesolution!');

flag=0;

break;

end

end

ifflag%还不是最优表，进行转轴运算

temp=0;

fori=1:nA-1

ifA(mA,i)>temp

temp=A(mA,i);

inb=i;%进基变量的下标

end

end

sita=zeros(1,mA-1);

fori=1:mA-1

ifA(i,inb)>0

sita(i)=A(i,nA)/A(i,inb);

end

end

temp=inf;

fori=1:mA-1

ifsita(i)>0&sita(i)<temp

temp=sita(i);

outb=i;%出基变量下标

end

end

%以下更新N

fori=1:mA-1

ifi==outb

N(i)=inb;

end

end

%以下进行转轴运算

A(outb,:)=A(outb,:)/A(outb,inb);

fori=1:mA

ifi~=outb

人人文库> 全部分类> 教育资料 > 备课教案

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论