【2017年整理】单纯形法.doc

PAGE \* MERGEFORMAT19 运筹学实验报告 题目：用单纯形法求解线性规划问题 姓 名 顾文远 学 号 2012436034 年级专业 12数电类2班 指导教师 苏珂 2013年11月14日 TOC \o 1-2 \h \z \u HYPERLINK \l _Toc372323318 一．实验目的 PAGEREF _Toc372323318 \h 2 HYPERLINK \l _Toc372323319 二． 实验环境 PAGEREF _Toc372323319 \h 3 HYPERLINK \l _Toc372323320 三． 实验内容 PAGEREF _Toc372323320 \h 4 HYPERLINK \l _Toc372323321 1.线性问题 PAGEREF _Toc372323321 \h 4 HYPERLINK \l _Toc372323322 2.函数调用形式： PAGEREF _Toc372323322 \h 4 HYPERLINK \l _Toc372323323 3.参数介绍： PAGEREF _Toc372323323 \h 4 HYPERLINK \l _Toc372323324 5.实验原理及步骤： PAGEREF _Toc372323324 \h 5 HYPERLINK \l _Toc372323325 四．函数代码及注释： PAGEREF _Toc372323325 \h 12 HYPERLINK \l _Toc372323326 五．数据测试： PAGEREF _Toc372323326 \h 16 HYPERLINK \l _Toc372323327 六、 实际问题求解 PAGEREF _Toc372323327 \h 19 一．实验目的 掌握两阶段法求线性规划问题的最优解，并完成 matlab求解单纯形算法求线性规划问题的最优可行解的程序。加深对单纯形算法的理解，掌握matlab的使用技巧. 实验环境 Matlab 7.1 系统：microsoft windows XP Professional CPU:Inter（R）Core（TM）2 Duo CPU E7500 @ 2.93GHz 3.24GB的内存 实验内容 1.线性问题 用matlab程序两阶段法解下面形式的方程组的最优解： min f=cT*x s.t. A*x=b x0 2.函数调用形式： [x,minf,flag,cpt]=dcxsf(A,b,c) 3.参数介绍： A为约束条件的系数矩阵。 b为约束条件的常数列向量。 c为目标函数的系数行向量。 4.函数输出介绍： X：为线性规划的解向量，当函数有可行解并且有最优值时输出最优解；当有可行解但无最优值或没有可行解时输出x[]。 minf：为目标函数的最优值，当函数有可行解并且有最优值时输出最优值；当有可行解但没有最优值时minf=-inf，输出-1/0；当没有可行解时输出minf=[]。 flag：为记录函数解的情况，当函数有可行解并且有最优值时，flag=1；当函数有可行解但没有最优值时，flag=0；当函数没有可行解时，flag=-1. cpt：为单纯形表，当函数有可行解并且有最优值时输出最优解对应的单纯形表；当有可行解但没有最优值时cpt=[]。 5.实验原理及步骤： 程序由用户输入约束方程组的系数矩阵、常数列向量和目标函数的系数行向量。通过调用函数的形式，采用两阶段法求解线性规划问题，并输出相关结果。 程序的第一阶段判断问题是否有可行解，若没有，程序停止，输出结果；如果有可行解，按第一阶段的方法可以求得一个初始的基本可行解，使运算进入第二阶段。第二阶段是从这个初始的基本可行解开始，使用单纯形法或者判定线性规划问题无界，或者求得一个最优解。 具体运算如下：対原问题增加m个人工变量Xα=(Xn+1,…,Xn+m)T后，用 单纯形法解如下问题： min g= s.t. Ax+xa=b x≥0, xa≥0 如上形式的辅助问题，它所对应的目标函数值为g0=。则该辅助线性规划问题一定有可行解及最优解，由于人工变量对应的m列构成单位矩阵且有b0，则辅助问题的第一个基本可行解为x=0,xa=b，该解对应的目标函数值，使用单纯型算法进行迭代，直到求得最优值。 构造辅助问题的单纯形表时，将检验向量与最优值合并作为单纯形表的第一行g。最后得到的基变量中如果有人工变量，则从该