线性规划单纯形解法实验报告

线性规划单纯形解法实验报告

ID:17455401

大小:255.23 KB

页数:31页

时间:2018-09-01

线性规划单纯形解法实验报告_第1页
线性规划单纯形解法实验报告_第2页
线性规划单纯形解法实验报告_第3页
线性规划单纯形解法实验报告_第4页
线性规划单纯形解法实验报告_第5页
资源描述:

《线性规划单纯形解法实验报告》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、运筹学实验报告题目:线性规划单纯形解法实验单纯形法求解线性规划运行环境:电脑型号方正操作系统WindowsXP32位处理器英特尔酷睿2双核内存4GBMatlab7.1一.实验目的:进一步熟练掌握单纯形法求解线性规划。二.实验内容:单纯形法求解线性规划4个三.实验原理:线性规划单纯形法(线性规划解有四种情形,唯一最优解,无穷多个最解,无界解,无可行解)两阶段法:第一阶段是判断线性规划是否有可行解,如果没有可行解,当然就没有基本可行解,计算停止;如果有可行解,按第一阶段的方法可以求得一个初始的基本可行解,使运

2、算进入第二阶段。第二阶段是从这个初始的基本可行解开始,使用单纯形方法或者判定线性规划问题无届,或者求得一个最优解。添加松弛变量、人工变量列出初始单纯形表计算非基变量各列的检验数σj唯一最优解某非基变量检验数为0基变量中有非零的人工变量所有σj≤0是否否否否无穷多最优解无可行解是是对任一σj>0有pj≤0无界解是否迭代运算1.用非基变量xk替换基变量xl2.对主元素行(第l行)令bl=bl/alk;alj=alj/alk3.对主元素列(第k列)令alk=1;其他元素=04.表中其他行列元素令alj=aij-

3、(ali/alk)*aikbi=bi-(bl/alk)*aik令σk=max(σj)对所有σik0计算θ=bi/aik令θl=min(θi)Xl为换出变量aik为主元素四.实验步骤1要求上机实验前先编写出程序代码2编辑录入程序3调试程序并记录调试过程中出现的问题及修改程序的过程4经反复调试后,运行程序并验证程序运行是否正确。5记录运行时的输入和输出。6.解决实际问题五.例题1.唯一最优解:A=[1,0,0,1,2;0,1,0,4,0;0,0,1,0,0]b=[8;16;12]c=[0,0,0,2,3]>>

4、A=[1,0,0,1,2;0,1,0,4,0;0,0,1,0,0];>>b=[8;16;12];>>c=[0,0,0,2,3];>>[x,minf,zuiyoubiao,flag]=linp(A,c,b)list=000-2-30000111520003610012100801040010160010000112list=000-2-3000001140-1002810012100801040010160010000112list=000-2-3000000100-1-101210012100801040

5、010160010000112list=000-2-3000000000-1-1-1010012100801040010160010000112已找到最优解!最优可行解为:x=8161200最优值为:minf=0最优解对应的单纯形表为:zuiyoubiao=000-2-3010012801040160010012flag=1x=8161200minf=0flag=1>>2.无穷多最优解:A=[1,0,0,1,2;0,1,0,4,0;0,0,1,0,0]b=[8;16;12]c=[0,0,0,2,4]>>A

6、=[1,0,0,1,2;0,1,0,4,0;0,0,1,0,0];>>b=[8;16;12];>>c=[0,0,0,2,4];>>[x,minf,zuiyoubiao,flag]=linp(A,c,b)list=000-2-40000111520003610012100801040010160010000112list=000-2-4000001140-1002810012100801040010160010000112list=000-2-4000000100-1-1012100121008010400

7、10160010000112list=000-2-4000000000-1-1-1010012100801040010160010000112已找到最优解!最优可行解为:x=8161200最优值为:minf=0最优解对应的单纯形表为:zuiyoubiao=000-2-4010012801040160010012flag=1x=8161200minf=0flag=1>>3.无界解:A=[1,0,-2,1;0,1,1,-1]b=[4;2]c=[0,0,1,1]>>A=[1,0,-2,1;0,1,1,-1];>

8、>b=[4;2];>>c=[0,0,1,1];>>[x,minf,zuiyoubiao,flag]=linp(A,c,b)list=00-1-100011-1000610-21104011-1012list=00-1-1000011-1-10210-21104011-1012list=00-1-10000000-1-1010-21104011-1012已找到最优解!最优可行解为:x=4200最优值为:minf=0最优解对应的单

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。