最優(yōu)捕魚策略

最優(yōu)捕魚策略第1頁，課件共38頁，創(chuàng)作于2023年2月

種群生態(tài)學模型

用數(shù)學模型去定量或定性地描述人口、生物種群等一些不易量化的復雜現(xiàn)象的變化規(guī)律,并將較抽象的數(shù)學概念與一些生態(tài)意義結(jié)合起來。第2頁，課件共38頁，創(chuàng)作于2023年2月一人口數(shù)量的變化規(guī)律

1.Malthus模型模型假設(shè)記x(t)為t時刻該國家或地區(qū)的人口總數(shù),1)忽略遷移對人口變化的影響,2)假設(shè)人口變化與出生率與死亡率有關(guān)，且每一個社會成員的死亡與生育水平相同,人口的出生率與死亡率之差與人口總數(shù)成正比,B-D=rx,

這個比例常數(shù)r稱為自然增長率,第3頁，課件共38頁，創(chuàng)作于2023年2月結(jié)論當r>0時,人口將以指數(shù)規(guī)律增長。當r<0時,人口將以指數(shù)規(guī)律減少。當r=0時,人口將將保持常數(shù)。模型建立與求解考慮t時刻到t+Δt時刻人口的變化Malthus模型

x(t+Δt)-x(t)=(B-D)Δt=rx(t)Δt得其解為：第4頁，課件共38頁，創(chuàng)作于2023年2月2Logistic模型當人口的數(shù)量特別大時，每個社會成員之間為生存的食物、空間和自然資源的競爭就不能忽略。必須在模型中增加一個競爭項，從統(tǒng)計學的觀點來看，社會成員在單位時間內(nèi)相遇的概率與x^2成比例，故選取競爭項為x^2。于是，Malthus模型修改為Logistic模型其中r稱為內(nèi)稟增長率，k為環(huán)境容納量。第5頁，課件共38頁，創(chuàng)作于2023年2月Logistic模型的解為：結(jié)論

人口總數(shù)趨于其環(huán)境容納量。2.當x(t)>k時,人口的數(shù)量將減少;

當x(t)<k時,人口的數(shù)量將增加。第6頁，課件共38頁，創(chuàng)作于2023年2月

3Leslie模型(具年齡結(jié)構(gòu)的模型)前述模型的不足(1)僅有人口總數(shù)，不能滿足需要；(2)沒有考慮到社會成員之間的個體差異,即不同年齡、不同體質(zhì)的人在死亡、生育方面存在的差異。1)模型假設(shè)同一年齡的人有相同的死亡機會和生育能力。這樣建立的模型不但使我們能夠更細致的預測人口總數(shù),而且能夠預測老年人口、學齡人口等不同年齡組的人口信息。2)模型建立設(shè)xk(t)為第t年年齡為k的人口數(shù)量,k=0,1,2,...,100,(忽略百歲以上的人口)。記

bk

是k歲人口的年生育率；

pk=1-dk

是k歲人口的年存活率，

dk為k歲人口的年死亡率。第7頁，課件共38頁，創(chuàng)作于2023年2月根據(jù)人口發(fā)展變化的特點：時間和年齡同步增長得模型如下：根據(jù)人的生理特征和人口學中的習慣,育齡區(qū)間一般取為15歲至49歲,即當k<15或k>49時,bk=0。此模型稱為Leslie模型利用此模型遞推計算，就可以得到每年各年齡組的人口數(shù)。第8頁，課件共38頁，創(chuàng)作于2023年2月模型的分析模型的矩陣形式其中利用Leslie矩陣模型遞推得第9頁，課件共38頁，創(chuàng)作于2023年2月解的漸近性態(tài)L是一個非負矩陣，它有主特征值0和特征向量v0。設(shè)L的特征值和特征向量為：0,v0;1,v1;…100,v100;如果所有的特征值是單根時,將x0表示為:x0=c0v0+c1v1+…+c100v100;則

x(t)=Ltx0=c0t0v0+c1

t1v1+…+c100

t100v100=t0(c0v0+c1

t1/

t0v1+…+c100

t100

/

t0v100)由上式得我們希望(c)發(fā)生，這可以通過適當?shù)挠媱澤邅韺崿F(xiàn)。第10頁，課件共38頁，創(chuàng)作于2023年2月預測與控制1三種模型預測比較(見書p210）2控制調(diào)節(jié)生育率bi

考慮Leslie模型，設(shè)bi(i=0,1,…100)是1982年的生育率。給他們乘以一個常數(shù)r,使得Leslie矩陣的主特征值為1。求得第11頁，課件共38頁，創(chuàng)作于2023年2月二多種群模型研究在同一環(huán)境中兩種或兩種以上的生物種群數(shù)量的變化規(guī)律.1兩種群模型用x(t),y(t)分別表示兩種群在t時刻的數(shù)量或密度,

考察x(t),y(t)各自的相對增長率,兩種群模型常用的形式是第12頁，課件共38頁，創(chuàng)作于2023年2月伏特拉模型假設(shè)函數(shù)f1(x),f(x),g1(ｙ),g2(ｙ)都是線性的.互惠共存型:

每一種群的存在,都對對方有利,對對方的數(shù)量增長起促進作用,這時c１≥０，b2≥０．(2)

捕食與被捕食型:

種群y以種群x為食物來源(或相反),這時種群x的存在對種群y的增長有利,而種群y的存在對種群x不利.因而c1≤0,b2≥0．(3)相互競爭型:

兩種群或者互相殘殺,或者競爭同一食物資源,各自的存在對對方的增長都是不利的,因而ｃ1≤０，ｂ2≤０．第13頁，課件共38頁，創(chuàng)作于2023年2月定性分析相互競爭型:記根據(jù)x,y的不同位置，可以得到兩種群變化趨勢。第14頁，課件共38頁，創(chuàng)作于2023年2月2三種群模型例1

一個捕食者，兩個食餌:

種群z是捕食者，不受密度制約；

x，y為食餌，受密度制約，且相互競爭。第15頁，課件共38頁，創(chuàng)作于2023年2月例2

捕食鏈種群z是y的捕食者，種群y是x的捕食者.

均受密度制約。第16頁，課件共38頁，創(chuàng)作于2023年2月三地中海鯊魚問題意大利生物學家Ancona曾致力于魚類種群相互制約關(guān)系的研究，他從第一次世界大戰(zhàn)期間,地中海各港口捕獲的幾種魚類捕獲量百分比的資料中，發(fā)現(xiàn)鯊魚等的比例有明顯增加（見下表），而供其捕食的食用魚的百分比卻明顯下降.顯然戰(zhàn)爭使捕魚量下降，食用魚增加，鯊魚等也隨之增加，但為何鯊魚的比例大幅增加呢？他無法解釋這個現(xiàn)象，于是求助于著名的意大利數(shù)學家V.Volterra，希望建立一個食餌—捕食系統(tǒng)的數(shù)學模型，定量地回答這個問題.第17頁，課件共38頁，創(chuàng)作于2023年2月記x(t)——食餌在t時刻的數(shù)量；y(t)——捕食者在t時刻的數(shù)量；a——食餌獨立生存時的增長率；b——食餌對捕食者的供養(yǎng)能力;c——捕食者獨自存在時的死亡率；d——捕食者掠取食餌的能力；h——捕獲強度系數(shù)2．基本假設(shè)：

1)

食餌由于捕食者的存在使增長率降低，假設(shè)降低的程度與捕食者數(shù)量成正比；

2)捕食者由于食餌為它提供食物的作用使其死亡率降低或使之增長，假定增長的程度與食餌數(shù)量成正比。第18頁，課件共38頁，創(chuàng)作于2023年2月

該模型反映了在沒有人工捕獲的自然環(huán)境中食餌與捕食者之間的制約關(guān)系，沒有考慮食餌和捕食者自身的阻滯作用，是Volterra提出的最簡單的模型.3．模型建立與求解

（一）不考慮人工捕獲第19頁，課件共38頁，創(chuàng)作于2023年2月

得首次積分求解方程有兩個平衡點兩方程相除得根據(jù)微分方程的理論,是閉曲線.第20頁，課件共38頁，創(chuàng)作于2023年2月故x(t)與y(t)都是周期函數(shù).該閉曲線就給出了生態(tài)學解釋.第21頁，課件共38頁，創(chuàng)作于2023年2月所以同樣得求平均值由方程解得第22頁，課件共38頁，創(chuàng)作于2023年2月（二）考慮人工捕獲設(shè)捕獲強度系數(shù)為h相當于食餌的自然增長率由a

降為a-h，捕食者的死亡率由c

增為c+h.此時,同樣得所以捕獲強度h減小,將使食餌的平均值減少,捕食者的平均值增加第23頁，課件共38頁，創(chuàng)作于2023年2月數(shù)值求解a=1,b=0.1,c=0.5,d=0.02x(0)=25,y(0)=2無人工捕獲時,對方程組求解得解曲線與相曲線為

x為實線，y為“*”線.

第24頁，課件共38頁，創(chuàng)作于2023年2月考慮人工捕獲設(shè)戰(zhàn)前捕獲強度系數(shù)為h=0.3,戰(zhàn)爭中降為h=0.1a=1,b=0.1,c=0.5,d=0.02x(0)=25,y(0)=2對方程組求數(shù)值解得解曲線為

實線為戰(zhàn)前的鯊魚比例，“*”線為戰(zhàn)爭中的鯊魚比例第25頁，課件共38頁，創(chuàng)作于2023年2月最優(yōu)捕魚策略

為了保護人類賴以生存的自然環(huán)境,可再生資源(如漁業(yè)、林業(yè)資源)的開發(fā)必須適度.一種合理、簡化的策略是,在實現(xiàn)可持續(xù)收獲的前提下,追求最大產(chǎn)量或最佳效益.考慮對某種魚的最優(yōu)捕撈策略.

假設(shè)這種魚分4個年齡組,稱1齡魚,…,4齡魚.各年齡組每條魚的平均重量分別為5.07克,11.55克,17.86克,22.99克,各年齡組魚的自然死亡率均為0.8/年,這種魚為季節(jié)性集中產(chǎn)卵繁殖,平均每條4齡魚的產(chǎn)卵量為1.109×105個,3齡魚的產(chǎn)卵量為這個數(shù)的一半,2齡魚和1齡魚不產(chǎn)卵,產(chǎn)卵和孵化期為每年的最后4個月,卵孵化并成活為1齡魚,成活率(1齡魚條數(shù)與產(chǎn)卵總量n之比)為1.22×１011／(１.2２×1011＋n).

漁業(yè)管理部門規(guī)定,每年只允許在產(chǎn)卵孵化期前的8個月內(nèi)進行捕撈作業(yè).如果每年投入的捕撈能力(如漁船數(shù)、下網(wǎng)次數(shù)等)固定不變,這時單位時間捕撈量將與各年齡組魚群條數(shù)成正比,比例系數(shù)不妨稱捕撈強度系數(shù).通常使用13ｍｍ網(wǎng)眼的拉網(wǎng),第26頁，課件共38頁，創(chuàng)作于2023年2月這種網(wǎng)只能捕撈3齡魚和4齡魚,其兩個捕撈強度系數(shù)之比042:1.漁業(yè)上稱這種方式為固定努力量捕撈.(1)建立數(shù)學模型分析如何實現(xiàn)可持續(xù)捕獲(即每年開始捕撈時漁場中各年齡組魚群條數(shù)不變),并且在此前提下得到最高的年收獲量(捕撈總重量).(2)某漁業(yè)公司承包這種魚的捕撈業(yè)務(wù)5年,合同要求5年后魚群的生產(chǎn)能力不能受到太大破壞.已知承包時各年齡組魚群的數(shù)量分別為:122×109條,29.7×109條,10.1×109條,3.29×109條,如果仍用固定努力量的捕撈方式,該公司應(yīng)采取怎樣的策略才能使總收獲量最高.第27頁，課件共38頁，創(chuàng)作于2023年2月這是一個典型的可再生資源開發(fā)的問題參考解答模型假設(shè)(1)魚分為1、2、3、4齡魚,4齡魚存活一年后仍劃為4齡魚.(2)各年齡組魚的自然死亡率為0.8／年,且死亡是一連續(xù)過程,

不是在某一時刻突發(fā).(3)3、4齡魚在一年的最后4個月集中產(chǎn)卵,且在該4個月的開始時刻進行;3齡魚產(chǎn)卵量為0.5×1.109×105個／條,4齡魚產(chǎn)卵量為1.109×105個／條;卵孵化成活為1齡魚,成活率為

1.22×１011／(１.2２×1011＋n).(4)捕撈在產(chǎn)卵孵化前8個月進行,且捕撈是一連續(xù)過程,不是在某一時刻發(fā)生;捕撈強度系數(shù)固定,只能捕撈3、4齡魚,它們的捕撈強度系數(shù)之比為0.42∶1(5)經(jīng)濟效益以捕撈總量來衡量.第28頁，課件共38頁，創(chuàng)作于2023年2月2.模型的建立及求解引入記號

T---年份(0,1,2,…).t---時間.

ｔ---間隔時間.Ni---ｉ齡魚的數(shù)量.r---自然死亡率(0.8／年).n---年產(chǎn)卵數(shù)量.E3---3齡魚捕撈強度系數(shù).E4---4齡魚捕撈強度系數(shù).3---3齡魚每年產(chǎn)卵量.4---4齡魚每年產(chǎn)卵量.第29頁，課件共38頁，創(chuàng)作于2023年2月(1).模型的建立考察1、2齡魚的生長過程,則有其解為為每年年初的ｉ齡魚總數(shù)。T年初的ｉ齡魚在T＋1年初變?yōu)椋椋饼g魚,所以有第30頁，課件共38頁，創(chuàng)作于2023年2月考察3、4齡魚的生長過程,在前8個月,由于捕撈與死亡均影響魚的變化,因而模型為其解為為每年年初的ｉ齡魚總數(shù)。故第8個月末時的i(i=3，4)齡魚總數(shù)為在后4個月,只有死亡率起作用,因而模型為第31頁，課件共38頁，創(chuàng)作于2023年2月其解為為第8個月末時的i齡魚總數(shù).故第T年末(T+1年初)i齡魚總數(shù)為第T+1年初的1齡魚由3、4齡魚產(chǎn)卵而來,產(chǎn)卵總數(shù)為故第T+1年初的1齡魚總數(shù)為第32頁，課件共38頁，創(chuàng)作于2023年2月總結(jié)以上結(jié)果,得出整個生存過程中滿足的關(guān)系為其中（1）第33頁，課件共38頁，創(chuàng)作于2023年2月(2).模型的求解在可持續(xù)捕撈情況下,Ni(T)趨于平衡,因而Ni(T)與T無關(guān),可得以下方程組:（2）第34頁，課件共38頁，創(chuàng)作于2023年2月求解得其中（3）（4）第35頁，課件共38頁，創(chuàng)作于2023年2月方程（4）的解為:由方程組（1）求得所以當