第四章到第六章地下水系统.pptVIP

第六章 线性规划 第四节 单纯形法 线性规划问题的最优解一定是某个基本可行解（可行域的一个顶点）。单纯形法的基本思想是：根据问题的标准型，从可行域中任意找出一个基本可行解（称为初始基本可行解），从这个解出发，转换到另一个基本可行解（顶点），并保证目标函数的值逐渐增大，当目标函数达到最大值时，就得到了问题的最优解。 第六章 线性规划 单纯形法的优点：只通过几次迭代就可以得到最优解。若采用枚举法，由于可行域有可能会有多个顶点，所以需要多次比较才能确定最优解。 单纯形法的计算整个过程可以归结为：选定初始基本可行解，最优性检验，作基变换求一组新的基本可行解，重复进行直到求出最优解。 第六章 线性规划 用图表进行单纯形法计算的步骤： 1 将线性规划问题进行标准式转换； 2 按表格形式填写有关变量及系数：第一行填决策变量及松弛变量下标；第二行填系数CJ;第三行填常数项b，以及决策变量和松弛变量；对应于b项以下填入约束条件的B列向量。第一次迭代时，以松弛变量为基变量，依次填入该列中。最下的zj-cj行各元素为对应列zj值与cj值之差，它是检验迭代是否继续进行下去的指标。 3 在zj-cj＜0的各行中取绝对值最大的所在列（称主列）为换入基变量。 4 主列大于零的系数分别除以常数项，取最小值所在行（称主行）的基变量作为换出基变量。 第六章 线性规划 5 主行与主列相交的系数称为主系数（又称主元）； 6 通过消元法，把主列各系数（主系数除外）化为零； 7 若表中有任一个（zj-cj ）＜0，且该列有aij＞0，则表明该目标函数仍有改善的余地，应继续迭代下去，直至表中所有（zj-cj ）≥0，即达到目标函数的极大值，迭代客终止。此时，b列所对应的基变量及zj值，为所求变量及目标函数值。 引进松弛变量，列出初始单纯形表 计算非基本变量各列检验数δj δj=zj-cj0 选择δk绝对值最大的k列作为主列，相应该列的变量为调入变量 主列各行ai,k0 计算出最小检验比θr确定r主行和主系数，ar,k主行变量为调出变量 计算出最小检验比θr确定r主行和主系数，ar,k主行变量为调出变量 修改单纯形表 1 以调入变量替换调出变量 2 计算主行新系数 3 计算其它各行新系数 计算最优解 结 束 否 是 问题无解 结 束 是 否 单纯形法计算框图 第六章 线性规划 第五节 具有人工变量的单纯形计算 一 ﹣M法（罚款法） 当线性规划数学模型含有“≥”或“=”约束方程时，应在其左端加一非负的人工变量yi，构成单位矩阵。但加入yi的矩阵，就违背了原来的约束条件，所以必须保证在最后的解中， yi =0，为此在目标函数式中，给人工变量一个很大的系数。对极大问题，系数用-M表示，对极小问题，系数用M表示（M本身为正数），只有当yi =0时，才能使M yi =0，目标函数才达到最优化。由于yi具有很大的系数而得到严格控制。 加入人工变量yi后，即可以yi 作为初始基本解，按单纯形法计算。 第六章 线性规划 二 二阶段法 第一阶段是判断原线性规划问题是否有解，并寻求一个初始基本解。为此，要构造一个辅助规划，这个辅助规划有一个单位矩阵，其目标函数以人工变量之和的负值为极大值，或以人工变量之和的正值为极小值，按单纯形法计算。若此辅助规划的最优解使其目标函数等于零，则说明没有一个人工变量在基变量内取正值，从而可得到原问题的一个基本可行解，否则原题无解。 第二阶段是求原问题的最优解。在第一阶段最后单纯形表的基础上，剔除人工变量，以第一阶段求得的最优解作为第一个基本可行解，以原问题的目标函数，继续用单纯形法进行迭代，直到达到最优解为止。 第六章 线性规划 第六节 利用线性规划模型建立 地下水管理模型的方法 地下水数学模型一般分为预报模型和管理模型。 用数学方程来描述地下水系统内部各变量及参数间的关系，从而对地下水系统的水头、水量或水质进行预报的模型称预报模型； 而在一定约束条件下，通过对地下水系统中某些状态变量的操纵，使系统按既定目标达到最优化，则称为地下水管理模型。 第六章 线性规划 地下水资源管理模型包括二部分： 第一部分是目标函数； 第二部分是约束条件，包括二个方面： 1 限量和需求约束； 2 水均衡约束。 第六章 线性规划 嵌入法 在线性变化的地下水系统预测模型中，输出可表现为输入的组合。在应用有限差分或有限单元法解地下水流动问题的数值法中，最终也归结为线性方程组以求解地下