环境计算题缩印版.docVIP

例1：空气污染治理最优化问题 如图有区域A，在不同功能区设m个控制点（j，m），有n个点源排污（i，n），造成部分控制点超标j的污染物监测值Cj超过允许的标准值Sj，即Cj≥ Sj。 设：Mi为i源产污量（kg/h），xi为i源治理去除量（kg/h），Mi—xi则为i源排放量（源强，kg/h）。由于采用不同的治理方案，各源单位治理成本Ci不同（元/kg），Ci xi为各源治理费用（元），区域治理总费用为（元）。 根据高斯模式，i源对j点的浓度传递系数aij可写出（对确定区域、确定源和确定控制点为确定值，在各控制点达标的前提下，考虑xi的分配方案，使区域治理总费用最小。 按题意可写出如下线性规划模型： 目标函数：minZ i 1,2,3,…,约束条件： i 1,2,3,…, j 1,2,3,…,m 非负约束：xi ≥ 0约束条件也可写为：，其中bj 。 设C [c1,c2,…c]，为单位成本矩阵（常数） X [x1,x2,…x]T，为决策变量阵（方案组合） B [b1, b 2,…b]T，约束条件阵（常数），[aij]n×m，为传递系数（常数）上述线性规划问题可写为：目标函数：minZ C X 约束条件：AX≥B 非负约束：X≥ 0 例2：河流污染治理最优化问题设3个污染源的污染物去处率分别为x1、x2、x3，则相应的处理费用为100 000*0.4 x1、50 000*0.4 x2、200 000*0.4 x3（元/ d），总费用为40 000 x1+20 000 x2+80 000x3（元/ d）。 据此，可确定目标函数：minZ 40 000 x1+20 000 x2+80 000x3 非负约束：x1、x2、x3≥ 0 根据整个河段（3个断面）不超标写出3个约束条件： 0.2*5*105+10*105 1- x1 / 5*105+105 1 即x1≥ 0.5 0.2*5+10 1- x1 exp -0.03*10 +20*0.5 1- x2 / 6+0.5 1 即exp -0.3 x1+x2≥ 0.35+1.1 exp -0.3 0.2*5+10 1- x1 exp -0.03*10 +20*0.5 1- x2 exp -0.03*15 +8*2 1- x3 / 6.5+2 1 即exp -0.75 x1+ exp -0.45 x2+1.6 x3≥ 0.75+1.1 exp -0.45 + exp -0.75 河流污染治理最优化目标函数：minZ 40 000 x1+20 000 x2+80 000x3 约束条件：x1≥ 0.5 exp -0.3 x1+x2≥ 0.35+1.1 exp -0.3 exp -0.75 x1+ exp -0.45 x2+1.6 x3≥ 0.75+1.1 exp -0.45 + exp -0.75 非负约束：x1、x2、x3 ≥ 0 也可表示为矩阵形式： 目标函数：minZ C X 约束条件：AX ≥ B 非负约束：X ≥ 0更一般数学规划问题为： 目标函数：min max Z f X 约束条件：G X ≥ B 非负约束：X ≥ 0 如何将线性规划化为标准型？ 将非标准型化为标准型的方法： minZ max（-Z） 约束条件为≥ ≤ 号的，可在不等式左边加（减）一个非负的新变量（松弛变量或剩余变量），同时修改目标函数； bj为负值的，两边同乘以（-1）。 例：LP1：maxZ 6x1+4x2 2x1+3x2≤100 4x1+2x2≤120 x1 14 x2≥22 x1，x2≥0 可化为LP2 ： maxZ 6x1+4x2+0x3+0x4+0x5 2x1+3x2+x3 100 4x1+2x2 + x4 120 x1 14 x2 -x5 22 x1，x2，x3，x4，x5≥0 上述约束方程也可写成矩阵形式： 能满足LP2的X* （x1，x2，x3，x4，x5）T，其对应的X* （x1，x2）T也能满足LP1，因此，可以认为，LP2在特定条件下是与LP1等价的线性规划。 4.怎样把握动态规划的建模与求解方法？（p92） fk,n Ck Sk,xk +fk+1,n Sk+1,…,Sn+1 1 常见的指标函数是取各段指标和的形式： 2 例7 求解图15-6的最短路径问题。 解：根据（1）式和（2）式，首先得到f4 S4 有两个允许决策: 路径分别为D1→E和D2→Ek 3时, 最短路线为C1→D2→E。 其最短路线为C2→D2→E。f3 C3 7，最短路线为C3→D2→E 当k 2时， 最短路线为B1→C1→D2→E。，f2 B2