論文 最優泄洪方案(定稿) 論文題目最優泄洪方案隊員指導教師:日期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