《论文—— 最优泄洪方案(定稿)》.doc

论文 最优泄洪方案(定稿) 论文题目最优泄洪方案队员指导教师:日期xx年8月31日最优泄洪方案摘要关键字matlab 一、问题的提出有一条河流由于河床泥沙淤积，每当上游发生洪水时，就会破堤淹没两岸，造成人员和财产的损失,为减少总的损失，人们采取破堤泄洪方法。 图1是该河一侧河岸区域的信息示意图。 在该区域边界上有很高的山，使该区域为封闭区域。 区域内分为17个小区，每个小区内标有三个数字，分别表示该小区的海拔高度h(m)，面积S(km万元)。 2)和被完全淹没时的土地、房屋和财产等损失总数k(百河流h4.8,S4.8,k2.03.35.85.4,9.3,5.7,7.0,5.0,8.4,7.04.6,6.1,1.45.4,9.1,1.66.0,1.0,4.06.00.83.23.5,8.5,6.04.2,0.9,0.94.3,3.6,9.44.8,0.6,2.14.8,0.7,2.44.8,8.3,5.33.4,2.3,4.14.5,1.5,4.14.0,4.6,3.0我们作以下假设1.各小区间有相对高度为1.3m的小堤相互隔离，例如左上方第一块和第二块小区间事实上有海拔6.7m的小堤。 2.当洪水淹没一个小区且水位高于该小区高度p(m)时，该小区的损失为该小区的k和p的函数损失=,01,1kppkp?3.假设决堤口可选在大堤和*小堤的任何地方，决堤口数目不受限制。 但一经决口，就不能再补合。 从河流经大堤决口流入小区的洪水量按决口数比例分配。 如在小区之间小堤开一决口，则假设该两小区之间的这段小堤不复存在。 若水位高过小堤，则将自动向邻近最低的一个小区泄洪。 若这样的小区有几块时，则平均泄洪。 需要解决的问题: （1）整个区域全部受损失的最小泄洪量maxQ （2）当洪水量为max213Q,max3Q时，分别制定泄洪方案，使总损失最小（在一种方案中，决堤同时进行）。 须计算出该方案的损失数。 二、模型的假设问题假设 1、各小区间有相对高度为1.3m的小堤相互隔离； 2、假设决堤口可选在大堤和小堤的任何地方，决堤口数目不受限制；符号说明h小区的海拔高度S小区的面积k被完全淹没时的土地、房屋和财产等损失总数 三、模型的建立与求解问题一的求解在观察了17个小区的划分区域后，我们将这17个小区按矩阵排列方式分别加以命名。 下面的ij b表示第i行第j列的小区，可得到各小区的海拔高矩阵h,面积矩阵S。 4.85.45.75.04.60海拔高矩阵h=5.46.06.03.54.24.34.84.84.83.44.54.04.89.37.08.46.10面积矩阵S=9.11.00.88.50.93.60.60.78.32.31.54.6由海拔高矩阵h可得，该17个区域中海拔最高的区域为22b和23b，即maxh=6.0（m），由题设可知，22b和23b区周围小堤的最小海拔高H=6.0+1.3=7.3（m）为求得最小泄洪量，我们可以想象，当17个小区完全受损时，泄洪区的水平面必须高于海拔最高小区的最小海拔高h=7.3（m）。 因此，全部区域受损的最小泄洪量minQ=i=i6131（ij hH?）ijS-0.3*22S-0.3*23S通过matlab和数据处理，求解得到当泄洪量为整个区域全部受损失的最小泄洪量时，每个小区对应的泄洪高度矩阵h12.51.91.62.32.70h1=1.91.31.33.83.13.02.52.52.53.92.83.3以及当泄洪高度为h1时，每个小区对应的泄洪量矩阵Q如下12.0017.6711.2019.3216.470Q=h1*S=17.291.301.0432.302.7910.801.501.7520.758.974.xx.18经过计算可得最小泄洪量minQ为810*9345.1（3m）。 问题二的求解 四、模型的结果分析求解及算法的主要步骤，使用的数学软件等 五、模型检验结果表示、分析与检验，误差分析等 六、模型评价本模型的特点，优缺点，改进方法 七、参考文献1献按正文中的引用次序列出，其中书籍的表述方式为编号作者，书名，出版地出版社，出版年。 参考文献中期刊杂志论文的表述方式为编号作者，论文名，杂志名，卷期号起止页码，出版年。 参考文献中网上资源的表述方式为编号作者，资源标题，网址，访问时间（年月日）。 2陈永春MATLAB M语言高级编程清华大学出版社 八、附录附录1在现有条件下，公司安排生产使获利最大时的源程序max=(2100-1400-400-135000*1.4/(1.4*m+n)-85000/(m+2*n)-9000/m)*m+(2000-1225-425-135000/(1.4*m+n)-85000*2/(m+2*n)-75000/n)*n;(1.4*m+n)3500;计算结果附录2通过引入“外包加工”费用，使公司获利最大时源程序如下max=(2100*x+465075-139563000/(1.4*(333