运筹学-大M法或两阶段法的上机实验.docx

实验报告 实验课程名称实验项目名称大M法或两阶段法的上机实验年级专业学生姓名学号00学院 实验时间：年月日 姓名 学 号 实验组 实验时间 指导教 师 成绩 实验项目名称 大M法或两阶段法的上机实验 实验目的及要求： 实验目的： 学会用Tora软件或Lindo软件求解线性规划问题， 理解每一步迭代计算中进基与出基变量等，了解大M法或两段法的上机实验。 实验要求： 完成作业P97页第6题及第7题（4）。 实验（或算法）原理： 1.大M法思路： 在单纯形法的基础上，为了使解线性规划有一个统一的解法，我们把所有求目标函数最小值的问题化为求目标函数最大值的问题。只要把目标函数乘以-1，就可以把原来求目标函数最小值的问题化为求目标函数最大值问题。为了找到一个满足条件的单位向量（非负），就需要加人工变量，注意人工变量与松弛变量和剩余变量是不同的，松弛变量和人工变量可以取零值也可以取正值，而人工变量只可以取零值，否则就会不等价。我们规定人工变量在目标函数中的系数为-M，M为任意大的数，这样只要人工变量大于零，所求的目标函数就是一个任意小的数，为了使目标函数最大，就必须将人工变量从基变量中换出。如果一直到最后，人工变量仍不能从基变量中换出，也就是说人工变量仍不为零，则该问题无可行解。像这样，为了构造初始可行基得到初始可行解，把人工变量”强行”的加到原来的约束方程中去，又为了尽力地把人工变量从基变量中替换出来就令人工变量在求最大的目标函数里的系数为-M的方法叫做大M法，M叫做罚因子。 两阶段法原理： 两阶段法是处理人工变量的另一种方法，这种方法是将加入人工变量后的线性规划问题分两阶段求解。第一阶段：要判断原线性规划问题是否有基本可行解，保持线性规划问题的约束条件原线性规划问题一样，而目标是求人工变量的相反数之和的最大值，如果此值大于零，即说明不存在使所有人工变量都为零的可行解，即原问题无可行解，因停止计算。如果此值为零，即说明存在一个可行解使得所有的人工变量都为零。第二阶段：将第一阶段的最终单纯形表中的人工变量（都是非基变量）取消，将目标函数换为原来的目标函数，把此可行解作为初始解进行计算，接下来的计算和单纯形法计算原理是一样的。 实验硬件及软件平台： PC机，Tora软件，Internet网。 实验步骤： 大M法步骤： 打开TORA命令窗口；选择Linearprogramming-Selectinputmode-Gotoinputscreen;输入待解的方程组-〉Slolvemenu-Solveproblem-Algebraic-Iterations-〉M-method-〉输入值-〉点击GoToOutputFormatScreen-〉点击GoToOutputScreen-〉点击Alllterations。 得出运行结果。 改变3步骤中的值（例100改为100000），再按之后的步骤运行，得出结果。 观察对比结果。 两阶段法步骤： 1）打开TORA命令窗口；2）选择Linearprogramming-Selectinputmode-Gotoinputscreen；3）输入待解的方程组-〉Slolvemenu-Solveproblem-Algebraic-Iterations-〉Two-phasemethod-〉点击GoToOutputFormatScreen-）点击Alllterations；4）得出运行结果。 实验内容（包括实验具体内容、算法分析、源代码等等）： 1.书上P97页第6题：用大M法和两阶段法求解下列线性规划问题。 maxz=5x+x+3x； 约束条件：x1+4x2+2x310，x-2x+x16. A：大M法TORA.LINEAR图1?10FileEditGrid TORA. LINEAR 图1?1 0FileEditGrid B：两阶段法 Phase1(Iter1 Basic xl x3 Sx4 Rx5 sx6 Solution z(min) 1(10 400 2DO -1一叩 o.oa 0.00 10.00 Hk5 1.DD 4.DD 2.DO 4.叩 I.Od 0.00 10.00 sk6 10。 NDD 1前 。一叩 o.oa 1.00 16.00 LowerBound D-DD 壮皿 (LOO UpperBound infini如 infinity infinity Unrestrd(jf/d]? n n Phase1(Iter2 Basic k1 k2 x3 Sx4 RxS sxG Solution z(min] 0.00 0.00 0.00 0.00 -1.0Q 0.00 0.00 x2 。一西 1一而 []一成 -。一应 0.25 0.00 25。 sxG 2一而 0.50 1.00 21*。 LowerBoun