单纯形法基本算法并求解线性规划问题解说.docVIP

数学与计算科学学院 实 验 报 告 实验项目名称 单纯形法的基本算法并求解线性规划问题 所属课程名称 最优化方法 实 验 类 型 算法编程 实 验 日 期 班 级 信息与计算科学1101班 学 号 201153100124 姓 名 陈 兵 成 绩 一、实验概述： 【实验目的】 （1）掌握Matlab的各种数据类型及基本运算; （2）有关矩阵基本运算的编程练习； （3）向量函数的梯度及Ｈｅｓｓｅ矩阵编程练习； （4）单纯形法算法得Ｍatlab编程源代码并进行机器求解； 【实验原理】 过复习Matlab数据的类型、定义、赋值和使用进行有关矩阵基本运算；编程练习,对单纯形法给出Matlab源代码进行机器求解。 【实验环境】 1. 2.Windows7旗舰版; Matlab 7 二、实验内容： 【实验方案】 一、单纯形法计算步骤 第一步：求初始基可行解，列出初始单纯形表。 对非标准型的线性规划问题首先要化为便准形式。由于总可以设法使越是方程的系数矩阵中包含一个单位矩阵，以此作为基求出问题的一个初始可行解。 为检验一个基可行解是否最优，需要将其目标函数值与相邻基可行解的目标函数值进行比较。为了书写规范和便于计算，对单纯型法的计算设计了一种专门的表格，称为单纯形表（如下表）。迭代计算中每找出一个新的基可行解时，就重画一张单纯形表。含初始基可行解的单纯形表成为初始单纯形表，含最优解的单纯形表称为最终单纯形表。 　 　 　 　 　 　 　 　 　 基 　 　 　 　 　 　 　 　 　 　 　 1 0 0　 　 0 0 0　 　 　 　 　　 0　 　 　0 　 　 　 　 单纯形表结构为：表的第2~3列列出基可行解中的基变量及其取值。接下来列出问题中所有变量，基变量下面列是单位矩阵，非基变量下面数字是该变量系数向量表为基向量线性组合时的系数。因是单位向量，故有 上表最上端的一行数是各变量在目标函数中的系数值，最左端一列数是与各基变量对应的目标函数中的系数值 对只要将它下面这一列数字与中同行的数字分别相乘，再用它上端的值减去上述乘积之和有 （1） 对j=1,…,n，将分别按上式求得的检验数，或写成（）记入表的最下面一行。 第二步：最优性检验 如表中所有检验数，且基变量中不含有人工变量时，表中的基可行解即为最优解，计算结束。当表中存在时，如有，则问题为无界解，计算结束；否则转下一步。 第三步：从一个基可行解转换到相邻的目标函数值更大的基可行解，列出新的单纯行表。 1、确定换入基的变量，只要有检验数，对应的变量???就可作为换入基的变量，当有一个以上检验数大于零时，一般从中找出最大一个 其对应的变量???作为换入基的变量(简称换入变量)。 2、确定换出基的变量。根据上一节中确定的规则，对列由式计算得到 确定是换出基的变量(简称换入变量)。元素决定了从一个基可行解到相邻基可行解的转移去向，取名主元素。 3、用换入变量替换基变量中的换出变量???，得到一个新的基。对应这个基可以找出一个新的基可行解，并相应的画出一个新的单纯行表 在这个新的表中的基仍应是单位矩阵，即应变换成单位向量。为此在表中进行行的处等变换，并将运算结果填入表相应格中。 （1）将主元素所在l行数字除以主元素，即有 (2) 将下表中刚计算得到的第l行数字乘上()加到上表的第i行数字上，记入下表的相应行。即有 　 　 　 　 　 　 　 　 　 基 　 　 　 　 　 　 　 　 　 　 　 1 0 0 　 　 0 0 1　 　 　 　 0 1 1 　　 0　 　 0 　 　 0　 （3）上表中各变量的检验数求法见（1）式。其中 第四步：重复第二、三步，一直到计算结束为止。 二、单纯形法计算步骤流程图 三、单纯形法标准化总结 线性规划模型 化标准形式 变量 取值无约束 不变 令，则 令，其中， 约束条件 右端项 不变 约束条件两端乘“-1” 形式 目标函数 极大或极小 不变 令，化为求 变量前的系数 加松弛变量时 加人工变量时 【实验过程】（实验步骤、记录、数据、分析） 一、手工计算，用单纯形法求解线性规划问题 解：将问题化为标准形式有 其约束条件系数矩阵的增广矩阵为