抗旱方案的制定題目及答案

1、題目A題 抗旱方案的制定位于我國西南地區(qū)的某個偏遠貧困村，年平均降水量不足20mm，是典型的缺水地區(qū)。過去村民的日常生活和農(nóng)業(yè)生產(chǎn)用水一方面靠的是每家每戶自行建造的小蓄水池，用來屯積每逢下雨時獲得的雨水，另一方面是利用村里現(xiàn)有的四口水井。由于近年來環(huán)境破壞，經(jīng)常是一連數(shù)月滴雨不下，這些小蓄水池的功能完全喪失。而現(xiàn)有的四口水井經(jīng)過多年使用后，年產(chǎn)水量也在逐漸減少，在表1中給出它們在近9年來的產(chǎn)水量粗略統(tǒng)計數(shù)字。2009年以來，由于水井的水遠遠不能滿足需要，不僅各種農(nóng)業(yè)生產(chǎn)全部停止，而且大量的村民每天要被迫翻山越嶺到相隔十幾里外去背水來維持日常生活。為此，今年政府打算著手幫助該村解決用水難的問題。

2、從兩方面考慮，一是地質(zhì)專家經(jīng)過勘察，在該村附近又找到了8個可供打井的位置，它們的地質(zhì)構造不同，因而每個位置打井的費用和預計的年產(chǎn)水量也不同，詳見表2，而且預計每口水井的年產(chǎn)水量還會以平均每年10%左右的速率減少。二是從長遠考慮，可以通過鋪設管道的辦法從相隔20公里外的地方把河水引入該村。鋪設管道的費用為 （萬元），其中 表示每年的可供水量（萬噸/年）， 表示管道長度（公里）。鋪設管道從開工到完成需要三年時間，且每年投資鋪設管道的費用為萬元的整數(shù)倍。要求完成之后，每年能夠通過管道至少提供100萬噸水。政府從2010年開始，連續(xù)三年，每年最多可提供60萬元用于該村打井和鋪設管道，為了保證該村從20

3、10至2014年這五年間每年分別能至少獲得150、160、170、180、190萬噸水，請作出一個從2010年起三年的打井和鋪設管道計劃，以使整個計劃的總開支盡量節(jié)?。ú豢紤]小蓄水池的作用和利息的因素在內(nèi)）。 表1 現(xiàn)有各水井在近幾年的產(chǎn)水量（萬噸）年份產(chǎn)水量編號2001200220032004200520062007200820091號井32.231.329.728.627.526.125.323.722.72號井21.515.911.88.76.54.83.52.62.03號井27.925.823.821.619.517.415.513.311.24號井46.232.626.723.020

4、.018.917.516.3表2 10個位置打井費用（萬元）和當年產(chǎn)水量（萬噸）編號12345678打井費用57546553當年產(chǎn)水2536321531282212A題 抗旱方案的制定摘要 我國西部地區(qū)是我國干旱災害較為嚴重的地區(qū)之一，隨著經(jīng)濟的迅速的發(fā)展、人口增長及由此引起的以氣候變暖為標志的全球氣候變化的發(fā)生,干旱災害有進一步加重的趨勢。為了使我國西南地區(qū)的人們能過上正常的生活，制定抗旱方案有其重要的意義。是的，水資源日趨緊張,如何合理開發(fā)、充分利用水資源,并能可持續(xù)供給,已成為人們研究、探求、關注的焦點。本文運用數(shù)學建模的相關知識,針對我國西南地區(qū)的某個偏遠貧困村的干旱情況,提出了一些如

5、何開發(fā)利用水資源,使人們能夠有充足的水資源，并且使整個抗旱方案的總開支盡量節(jié)省的方案。通過認真分析我國西北該貧困村的干旱情況，結合政府為解決旱災所出臺的政策，我們提出了兩個抗旱方案，并據(jù)此建立了三個數(shù)學模型。方案一：假設只打井，不鋪設管道，由此建立模型一：“只打井”模型。方案二：既打井又鋪設管道引水。對于方案二，我們根據(jù)假設的不同建立了模型二：“既打井又鋪設管道引水”模型和模型三：“既打井又鋪設管道引水”模型。假設一為：該村每年打井和鋪設管道的費用控制在60萬元之內(nèi)，即政府每年的撥款不累積且打井和鋪設管道的費用不拖欠；假設二為：政府提供給該村三年打井和鋪設管道的總費用為180萬元，打井和鋪設管

6、道的費用在三年后竣工時一次付清。通過建立模型，我們知道模型一的建立是失敗的。“只打井”的模型是無法滿足該村5年內(nèi)的需水量。從側面反映出該地區(qū)水資源相當匱乏。而模型二和模型三是可行的。模型二所體現(xiàn)的抗旱方案如下：三年的最小費用為172萬元；第一年打2、3、6、7號井，鋪設管道4.34公里，共花費52萬元；第二年打5號井，鋪設管道7.81公里，共花費60萬元；第三年打1號井，鋪設管道7.96公里，共花費60萬元。模型三所體現(xiàn)的抗旱方案如下：三年的最小費用為172萬元；第一年打1、2、3、5、7號井，鋪設管道19.82公里，共花費165萬元；第二年不打井，鋪設管道0.29公里，共花費2萬元；第三年打

7、6號井，鋪設管道0公里，共花費5萬元。模型二和模型三所呈現(xiàn)的優(yōu)化結果是一樣的，即三年的最小總費用均為172。但模型二和模型三還是有區(qū)別的，我們要根據(jù)現(xiàn)實的實際情況而采納不同的模型。關鍵詞:打井 鋪設管道 總開支 優(yōu)化問題 matlab lingo 一、問題重述和分析（一）、問題重述1、 問題的背景 位于我國西南地區(qū)的某個偏遠貧困村，年平均降水量不足20mm，是典型的缺水地區(qū)。過去村民的日常生活和農(nóng)業(yè)生產(chǎn)用水一方面靠的是每家每戶自行建造的小蓄水池，用來屯積每逢下雨時獲得的雨水，另一方面是利用村里現(xiàn)有的四口水井。由于近年來環(huán)境破壞，經(jīng)常是一連數(shù)月滴雨不下，這些小蓄水池的功能完全喪失。而現(xiàn)有的四口水

8、井經(jīng)過多年使用后，年產(chǎn)水量也在逐漸減少，在表1中給出它們在近9年來的產(chǎn)水量粗略統(tǒng)計數(shù)字。2009年以來，由于水井的水遠遠不能滿足需要，不僅各種農(nóng)業(yè)生產(chǎn)全部停止，而且大量的村民每天要被迫翻山越嶺到相隔十幾里外去背水來維持日常生活。 表1 現(xiàn)有各水井在近幾年的產(chǎn)水量（萬噸）年份產(chǎn)水量編號2001200220032004200520062007200820091號井32.231.329.728.627.526.125.323.722.72號井21.515.911.88.76.54.83.52.62.03號井27.925.823.821.619.517.415.513.311.24號井46.232.6

9、26.723.020.018.917.516.32、 相關信息我國西南地區(qū)干旱的主要原因：一是青藏高原的高原熱力因素。云南、四川省受“世界屋脊”的影響較大，青藏高原去年冬積雪少，高原熱力因素好，影響省內(nèi)出現(xiàn)干旱。二是厄爾尼諾現(xiàn)象影響。在此影響下，太平洋表層的熱流向東走，將熱帶水汽帶走，導致我國東部降水較多。與之相反，熱帶水汽減少，由南方進入到西南的水汽也偏少了。去年入秋以來，西南氣流偏弱、水汽輸送少造成降水偏少、氣溫偏高，這是云南等地產(chǎn)生持續(xù)氣象干旱的直接原因，就目前的干旱趨勢分析表明，西南等地出現(xiàn)冬春連旱的可能性較大。三是冷空氣勢力較弱。雖然今冬有數(shù)次冷空氣影響我國，但大多數(shù)冷空氣在北方開始

10、偏東移動，加之，受西南地形因素影響，冷空氣影響不到西南地區(qū)。四是云南、貴州特殊的喀斯特地貌，形成了雨水蓄不住，地下水用不上的狀況。五是抗旱基礎設施薄弱的問題依然突出。云南省的水資源是比較豐沛的，在全國排在第三位，但水資源利用率只有6%。貴州省已建成的17893處蓄水工程中，中型水庫僅34座，99.9%都是小型水庫，總蓄水量不到20億立方米。因為缺乏必要的水利設施，特別是缺少大中型骨干水利工程。3、問題的提出干旱已經(jīng)嚴重影響人們的正常生活，為此，今年政府打算著手幫助該村解決用水難的問題。從兩方面考慮，一是地質(zhì)專家經(jīng)過勘察，在該村附近又找到了8個可供打井的位置，它們的地質(zhì)構造不同，因而每個位置打井

11、的費用和預計的年產(chǎn)水量也不同，詳見表2，而且預計每口水井的年產(chǎn)水量還會以平均每年10%左右的速率減少。二是從長遠考慮，可以通過鋪設管道的辦法從相隔20公里外的地方把河水引入該村。鋪設管道從開工到完成需要三年時間，且每年投資鋪設管道的費用為萬元的整數(shù)倍。要求完成之后，每年能夠通過管道至少提供100萬噸水。政府從2010年開始，連續(xù)三年，每年最多可提供60萬元用于該村打井和鋪設管道，為了保證該村從2010至2014年這五年間每年分別能至少獲得150、160、170、180、190萬噸水，請作出一個從2010年起三年的打井和鋪設管道計劃，以使整個計劃的總開支盡量節(jié)省。表2 10個位置打井費用（萬元）

12、和當年產(chǎn)水量（萬噸）編號12345678打井費用57546553當年產(chǎn)水2536321531282212（2） 、問題分析由相關資料表明，該地區(qū)在未來的20102014年間的降水量依舊很少，光靠天然降水量是無法緩解該地區(qū)的旱情。所以制定和實施抗旱方案有其現(xiàn)實性和緊迫性。通過認真分析我國西北該貧困村的干旱情況，結合政府為解決旱災所出臺的政策，我們提出了兩個抗旱方案，并據(jù)此建立了三個數(shù)學模型。方案一：假設只打井，不鋪設管道。由于經(jīng)專家勘測得該地區(qū)有8個可供打井的位置，所以我們據(jù)此可以建立模型一：“只打井”模型。方案二：既打井又鋪設管道引水。該村20公里外豐富的河水為鋪設管道引水提供了條件，再加上8

13、個可供打井的位置，我們提出了方案二。對于方案二，我們根據(jù)假設的不同可以建立模型二和模型三。由于政府從2010年開始，連續(xù)三年，每年最多可提供60萬元用于該村打井和鋪設管道，所以我們假設一：該村每年打井和鋪設管道的費用控制在60萬元之內(nèi)，即政府每年的撥款不累積且打井和鋪設管道的費用不拖欠。并據(jù)此建立了模型二：“既打井又鋪設管道引水”模型；假設二：政府提供給該村三年打井和鋪設管道的總費用為180萬元，打井和鋪設管道的費用在三年后竣工時一次付清。并據(jù)此建立了模型三：“既打井又鋪設管道引水”模型。為了模型的建立，我們利用matlab軟件對該村原有的四口井進行分析，畫出20012009年原有四口井的年產(chǎn)

14、水量隨時間變化的散點圖（如圖1所示），通過計算機擬合，可得這四口井年產(chǎn)水量隨時間變化的函數(shù)關系式，從而計算出2010、2011、2012、2013、2014年這四口井的年產(chǎn)水量（產(chǎn)水量應大于零），計算結果如表4所示。針對模型一，我們要考慮分析3點： 1、20102012年這三年里如何計劃打井； 2、只需供應該村5年的需水量，那么打井的數(shù)量該是多少，第幾年打； 3、要怎樣才能讓總費用最少。對此，我們列出目標函數(shù)及約束條件，并利用lingo進行求解。針對模型二和模型三，我們要在模型一的基礎上進行進一步地分析。模型二和模型三均為優(yōu)化模型。我們通過討論、分析、總結、可得出該模型的目標函數(shù)和約束條件，最

15、后用lingo進行求解，即可得出該優(yōu)化問題的最優(yōu)解。模型中的決策變量是：（1）每年投入到管道鋪設的費用（必須為整數(shù)）；（2）是否在8 個位置中的某一位置打井；（3）在某一位置打井的時間；（4）管道鋪設成功之后每年提供水的量；而題目的約束條件為：（1）模型二：每年的總費用投入在60 萬元以內(nèi)；模型三：三年的總費用投入在180 萬元以內(nèi)；（2）管道鋪設完之后，每年至少提供100 萬噸水；（3）2010 年至2014 年必須保證每年提供的水量至少等于該年要求的供水量；二、模型的假設約定（1）假定打完井即可供水；（2）不考慮管道的滲漏因素；（3）打井所需的時間忽略不計；（4）不考慮人口增長對供水量的影

16、響；（5）不考慮小蓄水池的作用和利息的因素；（6）原有的四口井的年產(chǎn)水量滿足某種函數(shù)關系；（7）8 個位置的供水量會服從專家們預計的方式變化；（8） 假設管道鋪設后每年通過的水量是穩(wěn)定的，均為Q；（9） 規(guī)定每年鋪設的管道公里數(shù)保留兩位小數(shù)。三、符號說明及名詞定義（一）、符號及說明表3 符號說明符號定義單位表示三年打井及鋪設管道的最小花費；萬元表示原有四口井每年的產(chǎn)水量；萬噸表示專家勘察的8口井中第口井在第年打井；表示第年打井及修管道所花的費用；萬元表示八口井每年的供水量；萬噸表示勘測的八口井在挖掘當年的產(chǎn)水量；萬噸表示打第口井所需的費用萬元表示第年鋪設的管道長度公里表示第年鋪設管道的費用萬元

17、（二）、名詞定義散點圖：在回歸分析中，數(shù)據(jù)點在直角坐系平面上的分布圖。表示因變量隨自變量而變化的大致趨勢，據(jù)此可以選擇合適的函數(shù)對數(shù)據(jù)點進行擬合。 擬合：是指已知某函數(shù)的若干離散函數(shù)值f1,f2,fn，通過調(diào)整該函數(shù)中若干待定系數(shù)f(1, 2,n)，使得該函數(shù)與已知點集的差別(最小二乘意義)最小。源代碼（也稱源程序）：是指未編譯的按照一定的程序設計語言規(guī)范書寫的文本文件，是指一系列人類可讀的計算機語言指令。數(shù)學模型：根據(jù)對研究對象所觀察到的現(xiàn)象及實踐經(jīng)驗，歸結成的一套反映其內(nèi)部因素數(shù)量關系的數(shù)學公式、邏輯準則和具體算法。用以描述和研究客觀現(xiàn)象的運動規(guī)律。四、模型的建立和求解一、由表1所提供的原

18、有各水井在近幾年的年產(chǎn)水量的數(shù)據(jù)，利用matlab可得四口井的年產(chǎn)水量與時間的變化圖（程序見附錄一），如下圖：圖一 20012009年原有的四口井的年產(chǎn)水量隨時間變化的散點圖通過散點圖可以看出，1 號井、3 號井的變化函數(shù)接近一次函數(shù)，2 號井、4 號井的變化函數(shù)接近三次函數(shù)，再由計算機擬合，可得四口井的年產(chǎn)水量隨時間變化的函數(shù)關系：由函數(shù)關系得原有四口井在未來五年內(nèi)的年產(chǎn)水量，如下表表4 20102014年原有四口井年產(chǎn)水量年份產(chǎn)水量20102011201220132014第一口井21.4620.2619.0617.8616.66第二口井2.441.730.6400第三口井9.0974.91

19、2.820.73第四口井8.810.52000四口井總和41.829.5124.6120.6817.39二、建立模型模型一：“只打井”模型。 由于若8口井同時修，所需費用為40萬，則不考慮費用的約束條件。先求8口井每年的總供水量：2010年供水量： 2011年供水量：2012年供水量：2013年供水量：2014年供水量： 目標函數(shù)為： 約束條件為：（1）每口井最多只能打一次 ； （2）為了保證該村從2010至2014年5年間每年分別能至少獲得150、160、170、180、190萬噸水 用lingo求解（源代碼見附錄二）。lingo求解的結果是無解。由此我們可以得出不管怎樣打井都無法供應該村近

20、五年內(nèi)的需水量，因此此模型不成立。模型二：“既打井又鋪設管道引水”模型。 目標函數(shù)為：約束條件為：（1）每口井最多只能打一次 （2）每年的總費用為打井和鋪設管道的和 （3）題目要求三年鋪設的管道長度不小于20公里 （4）題目要求完成之后管道每年至少提供100萬噸水 （5）由假設一，從2010年起，連續(xù)三年，政府每年最多可以提供60萬元用于該村打井和鋪設管道 （6）為了保證該村從2010至2014年5年間每年分別能至少獲得150、160、170、180、190萬噸水 用lingo求解（源代碼見附錄三）得： 表4 三年打井及鋪設管道的優(yōu)化方案三年打井及鋪設管道的最小費用172萬第一年打井及鋪設管道

21、的花費52萬打井的號數(shù)2、3、6、7號井鋪設管道的長度4.34公里第二年打井及鋪設管道的花費60萬打井的號數(shù)5號井鋪設管道的長度7.81公里第三年打井及鋪設管道的花費60萬打井的號數(shù)1號井鋪設管道的長度7.96公里管道每年可提供的水量模型三：“既打井又鋪設管道引水”模型。目標函數(shù)為：約束條件為：（1）每口井最多只能打一次 （2）每年的總費用為打井和鋪設管道的和 （3）題目要求三年鋪設的管道長度不小于20公里 （4）題目要求完成之后管道每年至少提供100萬噸水 （5）由假設二，從2010年起，連續(xù)三年，政府每年最多可以提供60萬元用于該村打井和鋪設管道，則三年累計共180萬元，打井和鋪設管道的費

22、用在三年后竣工再結算 （6）為了保證該村從2010至2014年5年間每年分別能至少獲得150、160、170、180、190萬噸水 用lingo求解（源代碼見附錄四），得：表5 三年打井及鋪設管道的優(yōu)化方案三年打井及鋪設管道的最小費用172萬第一年打井及鋪設管道的花費165萬打井的號數(shù)1、2、3、5、7號井鋪設管道的長度19.82公里第二年打井及鋪設管道的花費2萬打井的號數(shù)無鋪設管道的長度0.29公里第三年打井及鋪設管道的花費5萬打井的號數(shù)6號井鋪設管道的長度0公里管道每年可提供的水量五、模型的評價、改進和推廣對于模型一，它的不可行性是由于該地區(qū)地下水資源相當匱乏導致的。為了該村村民能有充足的

23、水資源，為了他們能夠正常生活，地質(zhì)專家應該多多勘察該村的地質(zhì)構造，幫助他們多找到幾個可供打井的位置。對于模型二和模型三，由得到的結果可知：3 年總共鋪設的管道均為20.11公里（略大于20 公里），這是因為每年投入的費用必須為整數(shù)的原因得到的結果，其實管道可以只鋪設20 公里，則由管道鋪設費用的計算公式可以知道，只鋪設20 公里均可以減小費用0.76 萬元，因此，最小目標值均變?yōu)?71.24萬元，略優(yōu)于172萬元。分析模型二和模型三，我們可以這樣說，在費用可以拖欠的情況下，模型三的建立，有其現(xiàn)實的緊迫性。一方面，隨著中國經(jīng)濟的騰飛，在未來三年內(nèi)物價上漲的幅度很大。另一方面，據(jù)科學家和氣象專家預

24、測表明，2012年即將是一個非常時期，各種自然災害將頻繁發(fā)生。這告訴我們，抗旱方案要盡早實施。然而模型二更符合政府的要求，它體現(xiàn)了循序漸進，統(tǒng)籌規(guī)劃的理念。由表六和表七（見附錄），我們可以得出以下結論：不管是模型二還是模型三，隨著時間的推移，計劃所打的井的水資源會慢慢減少并將逐漸接近枯竭，該村的供水量將由管道直接提供。顯然管道的供水量是不足以提供該村的生活用水及農(nóng)業(yè)用水，所以我們應該重新選擇合理的位置進行打井，配合管道提供該村的用水量。六、幾點建議建議一：有可能的話建議進行人工降雨。建議二：建議重新選擇合理的位置進行打井。建議三：由于該村水資源相當匱乏，在不遠的將來旱情即將繼續(xù)，所以建議政府建

25、立大型調(diào)水工程。參考文獻1 謝金星、薛毅.優(yōu)化建模與lindo/lingo軟件M.北京： 清華大學出版社2孫祥.matlab7.0基礎入門M.北京:清華大學出版社3劉洋、胡恩召.組合優(yōu)化的城市供水量預測模型J.中國期刊全文數(shù)據(jù)庫,2008年08期4葛芬莉、韓廣鈞.陜西干旱趨勢與水資源利用J.水資源與水工程學報,2004年02期5李志剛、楊繼榮.基于LINGO9.0軟件的混合離散優(yōu)化方法研究J.機床與液壓,2007年1月第33卷第1期附錄：附錄一： 20012009年原有四口井的年產(chǎn)水量隨時間變化的散點圖x=1:9;x1=2:9;y1=32.2 31.3 29.7 28.6 27.5 26.1

26、25.3 23.7 22.7;y2=21.5 15.9 11.8 8.7 6.5 4.8 3.5 2.6 2.0;y3=27.9 25.8 23.8 21.6 19.5 17.4 15.3 13.3 11.2;y4=46.2 32.6 26.7 23.0 20.0 18.9 17.5 16.3;plot(x,y1,'or',x,y2,'+b',x,y3,'xm',x1,y4,'*g')附錄二： 模型一：“只打井”模型的lingo源代碼model:min=5*(r11+r21+r31)+7*(r12+r22+r32)+5*(r13+

27、r23+r33)+4*(r14+r24+r34)+6*(r15+r25+r35)+5*(r16+r26+r36)+5*(r17+r27+r37)+3*(r18+r28+r38);r11+r21+r31<=1;r12+r22+r32<=1;r13+r23+r33<=1;r14+r24+r34<=1;r15+r25+r35<=1;r16+r26+r36<=1;r17+r27+r37<=1;r18+r28+r38<=1;41.8+r11*25+r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12>=15

28、0;29.51+(r11*25+r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12)*0.9+r21*25+r22*36+r23*32+r24*15+r25*31+r26*28+r27*22+r28*12>=160;24.61+(r11*25+r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12)*0.9+r21*25+r22*36+r23*32+r24*15+r25*31+r26*28+r27*22+r28*12)*0.9+r31*25+r32*36+r33*32+r34*15+r35*31+r3

29、6*28+r37*22+r38*12>=170;20.68+(r11*25+r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12)*0.9+r21*25+r22*36+r23*32+r24*15+r25*31+r26*28+r27*22+r28*12)*0.9+r31*25+r32*36+r33*32+r34*15+r35*31+r36*28+r37*22+r38*12)*0.9>=180;17.39+(r11*25+r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12)*0.9+r21*25+

30、r22*36+r23*32+r24*15+r25*31+r26*28+r27*22+r28*12)*0.9+r31*25+r32*36+r33*32+r34*15+r35*31+r36*28+r37*22+r38*12)*0.9*0.9>=190;gin(r11) ;gin(r12) ;gin(r13) ;gin(r14) ;gin(r15) ;gin(r16) ;gin(r17) ;gin(r18) ;gin(r21) ;gin(r22) ;gin(r23) ;gin(r24) ;gin(r25) ;gin(r26) ;gin(r27) ;gin(r28) ;gin(r31) ;gin

31、(r32) ;gin(r33) ;gin(r34) ;gin(r35) ;gin(r36) ;gin(r37) ;gin(r38) ;EndNo feasible solution found. Extended solver steps: 0 Total solver iterations: 0 Variable Value R11 1.234568 R21 1.234568 R31 1.234568 R12 1.234568 R22 1.234568 R32 1.234568 R13 1.234568 R23 1.234568 R33 1.234568 R14 1.234568 R24 1

32、.234568 R34 1.234568 R15 1.234568 R25 1.234568 R35 1.234568 R16 1.234568 R26 1.234568 R36 1.234568 R17 1.234568 R27 1.234568 R37 1.234568 R18 1.234568 R28 1.234568 R38 1.234568 Row Slack or Surplus 1 0.000000 2 1.234568 3 1.234568 4 1.234568 5 1.234568 6 1.234568 7 1.234568 8 1.234568 9 1.234568 10

33、-1.234568 11 -1.234568 12 -1.234568 13 -1.234568 14 -1.234568附錄三： 模型二：“既打井又鋪設管道引水”模型的lingo源代碼model:min=g1+g2+g3 ;r11+r21+r31<=1;r12+r22+r32<=1;r13+r23+r33<=1;r14+r24+r34<=1;r15+r25+r35<=1;r16+r26+r36<=1;r17+r27+r37<=1;r18+r28+r38<=1;g1=0.66*Q0.51*L1+r11*5+r12*7+r13*5+r14*4+r1

34、5*6+r16*5+r17*5+r18*3 ;g2=0.66*Q0.51*L2+r21*5+r22*7+r23*5+r24*4+r25*6+r26*5+r27*5+r28*3 ;g3=0.66*Q0.51*L3+r31*5+r32*7+r33*5+r34*4+r35*6+r36*5+r37*5+r38*3 ;Q>=100 ;L1+L2+L3>=20 ;g1<=60 ;g2<=60 ;g3<=60 ;41.8+r11*25+r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12>=150 ;29.51+(r11*25+

35、r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12)*0.9+r21*25+r22*36+r23*32+r24*15+r25*31+r26*28+r27*22+r28*12>=160 ;24.61+(r11*25+r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12)*0.9+r21*25+r22*36+r23*32+r24*15+r25*31+r26*28+r27*22+r28*12)*0.9+r31*25+r32*36+r33*32+r34*15+r35*31+r36*28+r37*22+r38

36、*12>=170 ;20.68+(r11*25+r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12)*0.9+r21*25+r22*36+r23*32+r24*15+r25*31+r26*28+r27*22+r28*12)*0.9+r31*25+r32*36+r33*32+r34*15+r35*31+r36*28+r37*22+r38*12)*0.9+Q>=180 ;17.39+(r11*25+r12*36+r13*32+r14*15+r15*31+r16*28+r17*22+r18*12)*0.9+r21*25+r22*36+r23*

37、32+r24*15+r25*31+r26*28+r27*22+r28*12)*0.9+r31*25+r32*36+r33*32+r34*15+r35*31+r36*28+r37*22+r38*12)*0.9*0.9+Q>=190 ;gin(r11) ;gin(r12) ;gin(r13) ;gin(r14) ;gin(r15) ;gin(r16) ;gin(r17) ;gin(r18) ;gin(r21) ;gin(r22) ;gin(r23) ;gin(r24) ;gin(r25) ;gin(r26) ;gin(r27) ;gin(r28) ;gin(r31) ;gin(r32) ;g

38、in(r33) ;gin(r34) ;gin(r35) ;gin(r36) ;gin(r37) ;gin(r38) ;gin(g1) ;gin(g2) ;gin(g3) ;EndLocal optimal solution found. Objective value: 172.0000 Extended solver steps: 134 Total solver iterations: 5497 Variable Value Reduced Cost G1 52.00000 1.000000 G2 60.00000 1.000000 G3 60.00000 1.000000 R11 0.0

39、00000 0.000000 R21 0.000000 0.000000 R31 1.000000 0.000000 R12 1.000000 0.000000 R22 0.000000 0.000000 R32 0.000000 0.000000 R13 1.000000 0.000000 R23 0.000000 0.000000 R33 0.000000 0.000000 R14 0.000000 0.000000 R24 0.000000 0.000000 R34 0.000000 0.000000 R15 0.000000 0.000000 R25 1.000000 0.000000 R35 0.000000 0.000000 R16 1.000000 0.000000 R26 0.000000 0.000000 R36 0.000000 0.000000 R17 1.000000 0.000000 R27 0.000000 0.000000 R37 0.000000 0.000000 R18 0.000000 0.000000 R28 0.000000 0.000000 R38 0.000000 0.000000 Q 100.0000 0.000000 L1 4.340875 0.000000