建模論文 最優(yōu)化捕魚措施

1、最佳捕魚策略摘要漁業(yè)作為一種再生資源產(chǎn)業(yè)，在可持續(xù)發(fā)展的時代主題下，保證其持續(xù)穩(wěn)產(chǎn)是形勢所趨。本文利用微分方程和非線性規(guī)劃理論，探討在可持續(xù)收獲的條件下，如何通過調(diào)整捕撈強度系數(shù)，實現(xiàn)捕魚量的最大化。針對問題一，首先推導出魚群產(chǎn)卵、自然死亡、年齡隨時間變化等諸因素與各年齡組魚群數(shù)量的數(shù)學表達式，結合可持續(xù)捕撈，形成一組約束條件，以年捕獲量最大作為目標函數(shù)，建立非線性規(guī)劃模型。用Lingo編程求解得到：當捕撈強度系數(shù)取17.36時，年捕獲量最大，為3.88×1011克。然后利用Matlab畫出了在保證可持續(xù)捕獲的前提下，年度捕獲量隨捕撈強度系數(shù)變化的圖象，并經(jīng)過多次計算，驗證了結果的準

2、確性和穩(wěn)定性。針對問題二，在問題一模型的基礎之上，修改約束條件。首先采用每年的捕撈努力量固定，但各年彼此之間的捕撈努力量不盡相同的方式，然后采用每年的捕撈努力量都保持不變的方式，并將兩個模型比較得出采用模型二收益更大。鑒于此問是多元非線性規(guī)劃問題，且數(shù)據(jù)較大，為了得到全局最優(yōu)解，我們采用Matlab進行求解，最終得到結果為：13.8815.8818.3633.095.52得到最大的捕獲量為1.721012克，從而制定出最佳捕魚策略。此外，在模型的推廣中，改變模型一的假設，在認為4齡魚一年后仍為4齡魚的基礎上，對問題一進行了改進，得出的結果雖相差甚微，但是思路更具邏輯性。關鍵詞：微分方程 多元非

3、線性規(guī)劃 馬爾薩斯人口增長模型一、 問題重述為了保護人類賴以生存的自然環(huán)境，可再生資源（如漁業(yè)，林業(yè)資源）的開發(fā)必須適度。一種合理、簡化的策略是，在實現(xiàn)可持續(xù)收獲的前提下，追求最大產(chǎn)量或最佳效益?？紤]對鳀魚的最優(yōu)捕撈策略，該種魚的基本信息如表1所示；表1. 鳀魚的基本信息1齡魚2齡魚3齡魚4齡魚平均重量5.0711.5517.8622.99自然死亡率 0.8產(chǎn)卵量000.5545×1051.109×105這種魚為季節(jié)性集中產(chǎn)卵繁殖，產(chǎn)卵和孵化期為每年的最后4個月，卵孵化并成活為1齡魚，成活率（1齡魚條數(shù)與產(chǎn)卵量n之比）為1.221011/（1.221011 + n）.漁業(yè)管

4、理部門規(guī)定，每年只允許在產(chǎn)卵孵化期前的8個月進行捕撈作業(yè)。如果每年投入的捕撈能力（如漁船數(shù)，下網(wǎng)次數(shù)等）固定不變，這時單位時間捕撈量將與各年齡組魚群條數(shù)成正比，比例系數(shù)不妨稱捕撈強度系數(shù)。通常使用13mm網(wǎng)眼的拉網(wǎng)，這種網(wǎng)只能捕撈3齡魚和4齡魚，其兩個捕撈強度系數(shù)之比為0.421。漁業(yè)上稱這種方式為固定努力量捕撈。1分析如何實現(xiàn)可持續(xù)捕撈并且在此前提下得到最高的年收獲量2某漁業(yè)公司承包這種魚的捕撈業(yè)務5年，合同要求5年后魚群的生產(chǎn)能力不能受到太大破壞。已知承包時各年齡組魚群的數(shù)量分別為：122，29.7，10.1，3.29（×109條），如果仍用固定努力量的捕撈方式，該公司應采取怎樣

5、的策略才能使總收獲量最高。二、 問題分析由題設條件得出各年齡組魚群的轉化規(guī)律如1所示：1齡魚 死亡4齡魚2齡魚3齡魚 產(chǎn)卵捕獲 收獲 9-12月 1-8月 圖1.各年齡組魚群的轉化規(guī)律由上圖可知，3、4齡魚在第年產(chǎn)的卵所孵化成活的在第年初全部為1齡魚；第年末的（= 1，2,）齡魚在第年初全部變?yōu)榈邶g魚；為了建模的簡便和易于求解，不妨設4齡魚在第年初全部死亡。對于問題一，我們可以列出一組微分方程，求出各年齡組魚群上一年與今年之間的相互關系，作為一組約束條件，以年捕獲量最大作為目標函數(shù)，建立非線性規(guī)劃模型。對于問題二，考慮五年的最優(yōu)捕魚策略，對初始的魚群數(shù)量進行分析，發(fā)現(xiàn)即使不捕撈也不可能恢復到原

6、來的魚群數(shù)量，只能盡量減小對魚群的生產(chǎn)能力的破壞。利用各個年齡組魚群數(shù)量之間的遞推關系，根據(jù)初始魚群數(shù)量，遞推出五年后魚群的數(shù)量，從而計算得到五年內(nèi)的總收獲量。魚群的生產(chǎn)能力不造成太大的破壞可以理解為五年后魚群數(shù)量不少于初始魚群的一定比例，根據(jù)種群增長規(guī)律，選擇一個既能保證增長速度，又能保證種群數(shù)量優(yōu)勢的比例，即70。由此，該問題變?yōu)橐粋€多變量約束非線性優(yōu)化問題，由于模型較為復雜，采用智能算法進行優(yōu)化較為快速。三、 模型假設1) 漁場是非開放式漁場，不與其他水域發(fā)生關系，從而構成獨立的生態(tài)群落；2) 魚群產(chǎn)卵在九月初短時間內(nèi)完成，產(chǎn)卵魚的自然死亡發(fā)生在此之后；3) 13年齡組的魚群都在翌年年初

7、進入下一個年齡組，而原先的4齡魚魚群由于捕撈或者自然死亡等原因全部消失；4) 捕撈上來的魚全為活魚，即死亡的與被捕撈的魚分開計算；5) 3，4齡魚的平均產(chǎn)卵量掩蓋了性別差異四、 符號說明符號意義齡魚在一年中某時刻的數(shù)量（）齡魚魚群在第年年初的條數(shù)魚群的自然死亡率魚群的自然死亡率（）4齡魚的捕撈強度系數(shù)（4齡魚第年的捕撈強度系數(shù)）n（）一年內(nèi)總產(chǎn)卵量（第年的總產(chǎn)卵量）（）齡魚的產(chǎn)卵量（齡魚第年的產(chǎn)卵量）（）齡魚的年收獲量（齡魚在第年的收獲量）齡魚每條魚的平均重量年度總捕獲量（重量）承包期5年的總捕撈量五、 模型的建立與求解5.1問題一5.1.1模型的建立模型一（1）首先找到各個年齡組魚群在一年內(nèi)

8、的數(shù)量變化規(guī)律。由于1、2齡魚不會被捕撈，因此它們在全年內(nèi)所以只受自然死亡率的影響，即滿足以下關系： 對于3齡魚和4齡魚，它們在產(chǎn)卵開始前經(jīng)過了捕撈期間的自然死亡和被捕撈的雙重淘汰，在18月它們數(shù)量變化的關系為： 在912月，3齡魚和4齡魚數(shù)量變化的關系則為： 根據(jù)上述微分方程，在8個月的捕撈期過后，3齡魚的數(shù)量為，4齡魚的數(shù)量為。由此前的假設，可知第年年初1齡魚的數(shù)量應該為第年3齡魚和4齡魚產(chǎn)下的卵中成活下來的數(shù)量之和，即 （1）其中，若記3齡魚魚群產(chǎn)下的卵數(shù)量為，4齡魚產(chǎn)下的卵數(shù)量為，那么 且 同時，第年年初齡魚的數(shù)量為第年年底齡魚的數(shù)量。其中2齡魚的數(shù)量為 （2）3齡魚的數(shù)量為 （3）

9、4齡魚的數(shù)量為 （4）（2）考慮可持續(xù)發(fā)展要達到可持續(xù)發(fā)展的目的，即為每一年各個年齡組魚的數(shù)量都相等，就必須有 （5）（3） 目標函數(shù)即年收獲量的表達式3齡魚的捕撈條數(shù) （6）4齡魚的捕撈條數(shù) （7）因此年捕撈量（重量）為 （8）將（6）、（7）式帶入（8）式即可得到目標函數(shù)。綜上，以（8）為目標函數(shù)，以（1）、（2）、（3）、（4）、（5）為約束條件，建立一個非線性規(guī)劃模型為其中，5.1.2模型的求解用Lingo11.0軟件對以上模型求解（源碼見附錄1），輸出結果如下：Local optimal solution found. Objective value: 0.3887076E+12 E

10、xtended solver steps: 5 Total solver iterations: 209 Variable Value Reduced Cost K 17.36293 0.000000 N 0.6078058E+13 0.000000 Row Slack or Surplus Dual Price 1 0.3887076E+12 1.000000 2 0.000000 0.1258410E-02即可知當取17.36時，年捕獲量最大，為3.88×1011克。為了更好的描述捕撈強度系數(shù)與年捕獲量之間的對應變化關系，我們在Matlab中編寫程序，得到在可持續(xù)捕撈的前提下，年

11、捕獲量隨捕撈強度系數(shù)變化的圖象圖2 可持續(xù)捕撈下的年捕撈量與捕撈強度系數(shù)的關系圖由圖可知，年捕撈量G關于捕撈強度系數(shù)的關系為單峰函數(shù)，當捕撈強度系數(shù)k=17.36時，年最大捕撈量G=3.88×1012克。當捕撈強度系數(shù)為0時，捕撈量為0，符合實際；而當k趨于很大的時候，即過渡捕撈，魚群一定不能維持可持續(xù)發(fā)展狀態(tài)，所以捕撈強度系數(shù)k有極限，利用Lingo求得k最大為31.39，即為圖中曲線與橫軸的右交點。另外由于曲線靠近中間的斜率較小，在k=17.36附近對應的年捕撈量相對變化不大，實際最優(yōu)捕魚策略只需維持在k17.36的附近。 5.2問題二5.2.1模型的建立模型二在問題二中，需要考

12、慮五年的捕撈策略，對于每一年都采用固定努力量捕撈，但各年捕撈強度系數(shù)都可以不同，由不可對魚群造成太大破壞，定義五年后各魚群數(shù)量不少于初始魚群的70，即利用遞推公式(1) (2) (3) (4)，根據(jù)初始魚群數(shù)量，得到之后每年的魚群數(shù)量。由年捕撈量G 式（8），將5年的年捕撈量累加得到總捕撈量，即為目標函數(shù)其中，這是一個復雜的多變量非線性優(yōu)化問題，應用pattern search模式搜索算法可以快速進行模型的求解。因為該算法用于計算最小值，所以目標函數(shù)的負值作為算法的評價函數(shù)，求得最優(yōu)解，下圖為算法迭代過程，算法迭代40次后快速收斂到最優(yōu)值。圖3 算法迭代過程求得最后結果為 =1.72×

13、;1012克，具體如下表：表2 模型二的計算結果13.8815.8818.3633.095.52為了獲得最大的總收獲量，五年中各年的捕撈強度系數(shù)分別為13.88 ,15.88,18.36,33.09,5.52，在開始的三年要維持一定的捕撈強度，為之后的捕撈提供更多的儲備，在第4年應有比較大的捕撈強度，以達到最大的收獲量，第五年則應減少捕撈。模型三：然后我們采用5年的捕撈努力量不變的捕撈方式。只需令模型二中的（=1,2,3,4,5,），編制程序計算得結果為表3 模型三的結果5.601.49×1012最后結果為=1.491012克比較模型二與模型三，分別繪制4齡魚數(shù)量變化圖：模型二 ：捕

14、撈強度K不同，五年后各魚群數(shù)目不少于初始數(shù)目的70，圖為4齡魚數(shù)量的變化趨勢，第9年恢復平衡模型三：捕撈強度K五年都相同，五年后各魚群數(shù)目不少于初始數(shù)目的70，圖為4齡魚數(shù)量的變化趨勢，第7年恢復平衡圖4.4齡魚數(shù)量變化 由圖可知，從魚群恢復的角度來看，五年的捕撈強度系數(shù)相同時魚群更快恢復。顯然，此種情況下得到的捕撈總量與每年捕撈努力量不同的情況向比較要少的多，此方法不可取。那么可知，該公司要想獲得最大的總收獲量，應該采取5年內(nèi)的固定努力捕撈量互不相同的策略，其具體情況如表2所示。六模型的評價本文采用了非線性規(guī)劃的思想建立模型，通過求解有約束的非線性最大值問題，求的最優(yōu)解。文在建模過程中對魚的

15、產(chǎn)卵孵化過程進行了適當合理的簡化，建模簡單、清晰，具有一定的普遍意義。（1）模型的優(yōu)點問題一，在實現(xiàn)可持續(xù)不老的前提下，采用固定努力量捕撈，確定捕撈策略以得到最大捕撈總量。用微積分的方式分析各年各年齡魚的數(shù)量關系，并建立了非線性規(guī)劃的模型，運用Lingo和Matlab兩種軟件分別求解，得到的結果誤差很小，確保了結果的正確性問題二，在文章主體部分，采用了每年的捕撈努力量固定，但是各年之間不一定相同的方法建模，運用軟件進行多元搜索，得到最優(yōu)解，并在進一步討論中計算了各年之間捕撈努力量相同的總捕撈量，與原結果進行對比。而且在各種模型下，承包期結束后，我們對魚的恢復情況進行了模擬，證明了模型的可行性。

16、（2）模型的缺點與發(fā)展方向本模型對于產(chǎn)卵、孵化的簡化使得計算結果稍微偏離實際。未考慮自然資源的影響，使得結果與實際有所偏離。七模型的推廣與改進7.1改進一在以上模型，我們假設的是這種魚只活到4齡，以后它就死掉了。但是有實際情況仍可以假設上一年存活下來的4齡魚仍是4齡魚，這對模型沒有太大的差別，只是我們所做的假設的分析計算稍有變化。相比模型一，只有4齡魚的數(shù)量發(fā)生了變化。目標函數(shù)：約束條件：用lingo軟件進行求解，算法見附錄2：lingo2算法直接運行得:Local optimal solution found at iteration: 92 Objective value: 0.38870

17、76E+12 Variable Value Reduced Cost K 17.36292 -1.034723 N 0.6078067E+13 0.000000 Row Slack or Surplus Dual Price 1 0.3887076E+12 1.000000 2 0.000000 0.1258406E-02可知而模型一中算得出的結果為即，這兩個結果相差甚小。相比這兩個模型，推廣后的模型在邏輯方面更加符合實際，但是兩模型的計算結果相同，所以使用模型一計算更為簡便。7.2改進二在模型二中，我們單純的以捕撈總量代表收獲量，并未考慮到經(jīng)濟因素，以及捕魚的成本。為了是模型更具可靠性，可以

18、根據(jù)實際情況引入此種魚的價格、市場需求量等因素，合理改變各年的捕撈強度，以得到最大的經(jīng)濟效益。參考文獻1韓春生, 楊黎明, 茍林. 最優(yōu)捕魚策略J. 杭州電子工業(yè)學院學報, 1997, 4: 003. 2楊麗霞, 楊桂山, 苑韶峰. 數(shù)學模型在人口預測中的應用以江蘇省為例J. 長江流域資源與環(huán)境, 2006, 15(3): 287-291.3王文波. 數(shù)學建模及其基礎知識詳解M. 2006.附錄：附錄1：模型一的求解max=17.86*0.42*k/(0.8+0.42*k)*1.22*1011/(1.22*1011+n)*n*exp(-1.6)*(1-exp(-2/3*(0.8+0.42*k)+22.99*k/(0.8+k)*1.22*1011/(1.22*1011+n)*n*exp(-0.28*k-2.4) *(1-exp(-2/3*(0.8+k);n=1.22*1011*(1.109*105*(0.5*exp(-0.28*k-6.4/3)+exp(-(0.28+2/3)*k-8.8/3)-1);附錄2：模型改進的求解l