具有三個種群的食餌捕食者模型

馮廷碧+陰宇輝+陳玉雪+汪忠忠+羅習(xí)賓

550018 貴州師范學(xué)院數(shù)學(xué)與計算機科學(xué)學(xué)院 貴州貴陽

【摘 要】本文主要是在考慮具有三個種群的Volterra模型基礎(chǔ)上，初步考慮哺乳動物和爬行動物自身的阻滯作用，利用Matlab軟件得到其圖像，最后進(jìn)一步考慮植物自身的阻滯作用，加入logistic項，同時利用Matlab軟件得出其運行軌跡，判斷穩(wěn)定性。

【關(guān)鍵詞】Volterra模型；Matlab軟件；logistic項；改進(jìn)模型

20世紀(jì)20年代著名意大利數(shù)學(xué)家Volterra建立了一個簡單的數(shù)學(xué)模型，回答了意大利生物學(xué)家DAncona提出的問題。食餌的繁殖力下降，會造成捕食者的數(shù)量減少，但是其掠取能力卻會增加；捕食者的死亡率上升會導(dǎo)致食餌增加，食餌供養(yǎng)捕食者的能力增加會使食餌減少。為后來數(shù)學(xué)家和生物學(xué)家建立食餌-捕食者模型系統(tǒng)打下基礎(chǔ)。我們把Volterra建立的這種簡單的只有兩個種群的模型稱作Volterra模型，這種模型雖然能解釋一些現(xiàn)象，但是Volterra模型存在描述的周期變化狀態(tài)不是穩(wěn)定結(jié)構(gòu)等缺點，為了滿足人類自身生存與發(fā)展的需要，對資源進(jìn)行合理的開發(fā)和管理。人們一方面對種群的發(fā)展變化做定量的分析與預(yù)測，另一方面是通過對模型的研究判斷出付出多大的捕獲量，即可維持生態(tài)平衡，又能獲得最大收益。所以研究具有三個種群的食餌-捕食者模型有重大意義。

1Volterra食餌-捕食者模型

一個島嶼上棲居著食肉爬行動物和哺乳動物，又長著茂盛的植物，爬行動物以哺乳動物為食，哺乳動物又依賴植物生存，建立三者之間關(guān)系的模型。

植物，哺乳動物，爬行動物在t時刻的數(shù)量分別記住x1（t），x2（t），x3（t）。假設(shè)不受環(huán)境等因素影響，（相對）增長率為r1，又因為哺乳動物的存在使得植物的增長率減小，設(shè)減小的程度與哺乳動物數(shù)量成正比。

哺乳動物離開植物就無法生存，設(shè)它獨自生存時的死亡率為r2，而植物為哺乳動物提供食物，降低哺乳動物的死亡率，并且促進(jìn)哺乳動物增長，設(shè)這種作用與植物數(shù)量成正比，又因為爬行動物的存在對哺乳動物的增長減小，設(shè)這種減小程度與哺乳動物的數(shù)量成正比，于是x2（t）滿足x2（t）=x2（-r2+bx-λx3）比例系數(shù)b反映植物對哺乳動物的供養(yǎng)能力，λ反映爬行動物對哺乳動物的掠取能力。

同理，設(shè)爬行動物獨自生存時的死亡率為r3，而哺乳動物為爬行動物提供食物，從而減低爬行動物的死亡率促進(jìn)爬行動物增長，設(shè)這種作用與哺乳動物數(shù)量成正比，于是x3（t）滿足方程x（t）=x3（-r3+cx2）

比例系數(shù)c反映哺乳動物對爬行動物的供養(yǎng)能力。

數(shù)值解：

記植物，哺乳動物，爬行動物的初始數(shù)量分別為x1（0）=x11，x2=x22，x3（0）=x33并作圖及軌線圖，編程如下：

functionxdot=shier（t，x）；

r1=1；r2=0.5；r3=0.6；a=0.1；b=0.02；c=0.05；lambda=0.004；

xdot=[（r1-a*x（2））.*x（1），（-r2+b*x（1）-lambda*x（3））.*x（2），（-r3+c*x（2））.*x（3）]；

[t，x]=ode45（‘shier，[0，20]，[100，25，5]）；

plot（t，x（：，1），*，t，x（：，2），-，t，x（：，3），：），

legend（‘x（1），x（2），x（3））

grid

plot3（x（：，1），x（：，2），x（：，3）），grid

可以猜想x1（t），x2（t），x3（t）是周期函數(shù)，從圖中可以得出x1（t），x2（t），x3（t）的最大值分別為125.6，28.2，29.1。最小值分別為2.5，1.5，2.1。容易得出x1（t），x2（t），x3（t）在一個周期的平均值為64，30，15.6?？傮w上看出模型趨于某種平衡狀態(tài)，存在穩(wěn)定的平衡點，這是Volterra模型的特點，但是在現(xiàn)實生活中多數(shù)的食餌-捕食者模型都觀察不到Volterra模型的那種周期，而且Volterra模型容易受數(shù)值的影響。所以我們需要改進(jìn)模型。

2初步模型改進(jìn)

自然世界中植物的生長由于受到人為的保護，且人們大量的植樹造林，所以假設(shè)植物的獨自生長時其數(shù)量的增長服從指數(shù)增長規(guī)律，哺乳動物，爬行動物獨自生長是有自身的阻滯作用，服從logistic變化規(guī)律。

同樣將植物，哺乳動物，爬行動物在t時刻的數(shù)量記為：x1（t），x2（t），x3（t）。

同理，植物在獨自生存時x1（t）滿足：x（t）=x1（r1-ax2）

哺乳動物服從logistic規(guī)律，于是可得到哺乳動物增長的方程類似的，哺乳動物的存在也影響了爬行動物的增長，也可以得到爬行動物的方程，它們就是是自然界中食餌和捕食之間的關(guān)系，這里考慮到了種群自身的阻滯作用。

數(shù)值解：

設(shè)初值條件，利用Matlab軟件編程，雖然植物，哺乳動物，爬行動物的數(shù)量，在某一時刻會發(fā)生波動，但是在受到干擾之后還是會趨于某一時刻趨于穩(wěn)定。

此模型具有局限性，因為達(dá)到穩(wěn)定狀態(tài)的情況下，這三個種群都滅絕了，說明以上假設(shè)不能客觀反應(yīng)自然界中植物，哺乳動物，爬行動物之間內(nèi)在的聯(lián)系。

3模型的進(jìn)一步改進(jìn)

假設(shè)植物，哺乳動物，爬行動物存在自身阻滯作用，滿足logistic規(guī)律，同樣將植物，哺乳動物，爬行動物在t時刻的數(shù)量記為：x1（t），x2（t），x3（t）。

3.1模型建立

有植物，哺乳動物，爬行動物三個種群，當(dāng)他們獨自生存時，數(shù)量滿足logistic規(guī)律，r1，r2，r3是它們的固有增長率，N1，N2，N3自然世界對它們的最大容量。利用Matlab軟件求微分方程的數(shù)值解，做出圖像，然后從理論上研究平衡點。

3.2數(shù)值解

利用Matlab編程如下：

編寫.m文件

functionxdot=shier（t，x）；

r1=0.6；r2=0.5；r3=0.6；ct（1）=0.5；ct（2）=0.6；ct（3）=0.6；ct（4）=0.6；N（1）=600；N（2）=100；N（3）=15；

xdot=[r1*（1-x（1）/N（1）-ct（1）*x（2）/N（2））.*x（1），r2*（-1+ct（2）*（x（1）/N（1））-x（2）/N（2）-ct（3）*（x（3）/N（3）））.*x（2），r3*（-1+ct（4）*（x（2）/N（2））-x（3）/N（3））.*x（3）]；

在操作窗口輸入：

[t，x]=ode45（'shier'，[0，10]，[300，100，25]）；

plot（t，x（：，1），*，t，x（：，2），-，t，x（：，3），：），

legend（‘x（1），x（2），x（3））

grid

這個模型比很好的解決了第一個模型和第二個模型的局限性，較好的反映了，自然世界中，這三個種群的內(nèi)在聯(lián)系。

3結(jié)論

自然世界中，瞬息萬變，雖然在Volterra模型的基礎(chǔ)上加入了服從logistic規(guī)律的項，較好的反映出在自然環(huán)境中食餌和捕食者之間的依存和制約關(guān)系。但是卻不能完全的說明物種的生存關(guān)系，所以還需要大量學(xué)習(xí)。

參考文獻(xiàn)：

[1]劉漢武，王榮欣.具有性別結(jié)構(gòu)的食餌-捕食者模型[J].2005.

[2]張?zhí)烊?具有兩性的捕食者-食餌模型的漸近形態(tài)[J].2005.

[3]王靜，王克.比率型-捕食者-兩競爭食餌模型的動力學(xué)行為[J].2004.

作者簡介：

馮廷碧（1992～），女，貴州人，貴州師范學(xué)院數(shù)學(xué)與應(yīng)用數(shù)學(xué)專業(yè)學(xué)生。

基金項目：貴州師范學(xué)院學(xué)生科研項目基金（2014DXS083）。