实验报告(单纯形法的matlab程序)_第1页
实验报告(单纯形法的matlab程序)_第2页
实验报告(单纯形法的matlab程序)_第3页
实验报告(单纯形法的matlab程序)_第4页
实验报告(单纯形法的matlab程序)_第5页
全文预览已结束

下载本文档

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

文档简介

南昌航空大学数学与信息科学学院实验报告 实验一：线性规划单纯形算法一、实验目的通过实验熟悉单纯形法的原理，掌握Matlab循环语句的应用,提高编程的能力和技巧。二、算法 对于一般的标准形式线性规划问题(求极小问题),首先给定一个初始基本可行解。设初始基为B,然后执行如下步骤：(1).解,求得，(2).计算单纯形乘子, ,得到,对于非基变量，计算判别数,令 ,R为非基变量集合若判别数 ,则得到一个最优基本可行解，运算结束；否则，转到下一步(3).解,得到；若,即的每个分量均非正数，则停止计算，问题不存在有限最优解，否则，进行步骤(4).(4).确定下标r,使为离基变量。为进基变量，用替换，得到新的基矩阵，返回步骤（1）。对于极大化问题，可以给出完全类似的步骤，只是确定进基变量的准则不同。对于极大化问题，应令 四、计算框图 开始 初始可行解令计算单纯形乘子，计算判别数（非基变量）令是得到最优解解方程，得到。否是不存在有限最优解确定下标，是否为进基变量，用替换，得到新的基矩阵五、计算程序function x,f=zuiyouhua(A,b,c)size(A)=m,n;i=n+1:n+m;%基变量集合,后面m个松弛变量为初始基变量;N=1:n;%初始非基变量;B=eye(m,m);xb=b;xn=zeros(m,1);f1=0;w=zeros(1,m);z=-c;%初始判别数;flag=1;while(1) a,k=max(z);%x(k)为进基变量; if a=0 flag=0; break else y=inv(B)*A(:,k) if y0); a,r1=min(b1(t)./y(t) r=t(r1); %基变量中第r个变量为退基变量; i(:,r)=k B(:,r)=A(:,k);%换基，即将原基中第个变量换成第个变量; cb=c(:,i);%新的价值系数; xb=inv(B)*b; b0=xb; x=zeros(1,n+m) x(:,i)=xb f=cb*xb z=cb*inv(B)*A-c;%可用z=cb*(BA)-c,判别数. end end六、数值实验及结果分析求解线性规划问题:在工作区输入:A=3,3,1,0;-4,-4,0,1;2,-1,0,0;b=30,16,12;c=-3,1,0,0; x,f=zuiyouhua(A,b,c)x =
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

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

最新文档

评论

0/150

提交评论