petsaal.nb

Peto- ja saaliskanta

Peto- ja saaliskantojen keskenäistä vuorovaikutusta voidaan mallintaa toisistaan riippuvien differentiaaliyhtälöiden avulla. Tässä tarkastellaan yksinkertaista mallia, joka perustuu ns. Volterran–Lotkan yhtälöihin. Yhtälöt konstruoidaan seuraavalla logiikalla. Saaliseläinkantaa kuvaavan funktion $s\left(t\right)$ muutokset (derivaatan arvot) riippuvat kahdesta tekijästä. Toisaalta kantaan vaikuttaa  ympäristötekijöiden muodostama kasvutermi $\alpha s\left(t\right)$, joka kuvaa populaation (eksponentiaalista) kasvua olosuhteissa, jossa petokannan vaikutusta ei ole. Tämän lisäksi kantaan vaikuttaa petokannan $p\left(t\right)$ aiheuttama kantaa vähentävä tekijä. Perusyhtälö saaliskannalle on siis

$\frac{ds\left(t\right)}{dt}=\left(\alpha -\beta p\left(t\right)\right)s\left(t\right),$

jossa $\alpha$ on saalispopulaation kasvukerroin ja $\beta$ on petopopulaation syömiskerroin. Vastaavasti petokannan yhtälössä tulee olla saaliskannan puuttumisesta aiheutuva kantaa heikentävä (eksponentiaalinen) vähenemistekijä $\gamma p\left(t\right)$ sekä saaliskannan runsaudesta aiheutuva kasvutekijä. Niinpä perusyhtälö petokannalle on muotoa  

$\frac{dp\left(t\right)}{dt}=\left(-\gamma +\delta \beta s\left(t\right)\right)p\left(t\right),$

missä $\gamma$ on saalistajien vähenemistekijä ja $\delta$ on saalistajien tehokerroin saaliseläinten ollessa saatavilla.

Tarkastellaan seuraavaksi, millaisen mallin kantojen käyttäytymiselle yllä esitetyt mallit antavat. Olkoot tarkasteltavina lajeina ilves- ja lumikenkäjäniskannat Kanadassa.   

Laskujen aluksi on syytä hävittää mahdollisista aiemmista laskuista jääneet muuttujat.

$\mathrm{Remove}\left[''Global`*''\right]$

Määritellään ilves- ja jäniskantoja kuvaava differentiaaliyhtälöryhmä (yksikkönä tuhat yksilöä). $\left(\alpha =1.2,\beta =0.04,\gamma =0.5,\delta =0.3\right)$

In[1]:=

$\mathrm{ryhma}=\left\{s\text{'}\left[t\right]==s\left[t\right]\left(1.2-0.04p\left[t\right]\right),p\text{'}\left[t\right]==p\left[t\right]\left(-0.5+0.012s\left[t\right]\right)\right\}$

Out[1]=

$\left\{s^{\prime }\left[t\right]⩵\left(1.2\text{}-0.04p\left[t\right]\right)s\left[t\right],p^{\prime }\left[t\right]⩵p\left[t\right]\left(-0.5+0.012s\left[t\right]\right)\right\}$

Probleeman tuntemattomat funktiot ovat jäniskantaa kuvaava $s\left(t\right)$ sekä ilveskantaa kuvaava $p\left(t\right)$. Aika $t$ kuvaa kuluneita vuosia.

In[2]:=

$\mathrm{tuntemattomat}=\left\{s\left[t\right],p\left[t\right]\right\}$

Out[2]=

$\left\{s\left[t\right],p\left[t\right]\right\}$

Valitaan kantojen kokoa kuvaavat alkuehdot tarkastelun alkaessa. Ilveksiä olkoon 15 000 ja lumikenkäjäniksiä 60 000 yksilöä.

In[3]:=

$\mathrm{alkuehto}=\left\{s\left[0\right]==60,p\left[0\right]==15\right\}$

Out[3]=

$\left\{s\left[0\right]⩵60,p\left[0\right]⩵15\right\}$

Ratkaistaan differentiaaliyhtälöryhmä ja muutetaan tulos käyttökelpoiseen muotoon:

In[4]:=

$\mathrm{rtk}=\mathrm{NDSolve}\left[\mathrm{Flatten}\left[\left\{\mathrm{ryhma},\mathrm{alkuehto}\right\}\right],\mathrm{tuntemattomat},\left\{t,0,80\right\}\right]$

Out[4]=

$\left\{\left\{s\left[t\right]\to \mathrm{InterpolatingFunction}\left[\left\{\left\{0.,80.\right\}\right\},<>\right]\left[t\right],p\left[t\right]\to \mathrm{InterpolatingFunction}\left[\left\{\left\{0.,80.\right\}\right\},<>\right]\left[t\right]\right\}\right\}$

In[5]:=

$\mathrm{maarat}=\mathrm{tuntemattomat}/.\mathrm{First}\left[\mathrm{rtk}\right]$

Out[5]=

$\left\{\mathrm{InterpolatingFunction}\left[\left\{\left\{0.,80.\right\}\right\},<>\right]\left[t\right],\mathrm{InterpolatingFunction}\left[\left\{\left\{0.,80.\right\}\right\},<>\right]\left[t\right]\right\}$

Piirretään kuvat saaliseläin- ja petokantojen vaihtelusta ajan suhteen.

In[6]:=

Out[6]=

$⁃\mathrm{Graphics}⁃$

Saamme tulokseksi erittäin säännöllisesti vaihtelevat populaatiot, jossa jäniskanta on merkitty punaisella ja ilveskanta sinisellä käyrällä.

Todellisuudessa vain harvoin ympäristöolot ovat niin tasaiset, että yllä lasketun kaltaista kantojen vaihtelua esiintyisi. Yleensä sää ja muut kannoista riippumattomat ympäristötekijät vaikuttavat merkittävästi populaatioihin. Verrataan kuitenkin saatua tulosta Kanadan lumikenkäjänis- ja ilveskantoihin vuosina 1855–1930 (datan alkuperänä on Hudson's Bay Companyn tilastot).

In[7]:=

$<<\text{Graphics`MultipleListPlot`}$

$\mathrm{janis}=\mathrm{Table}\left[\left\{\left\{1855,76\right\},\left\{1856,88\right\},\left\{1859,19\right\},\left\{1861.5,26\right\},\left\{1862.5,10\right\},\left\{1863.5,155\right\},\left\{1865,151\right\},\left\{1867,20\right\},\left\{1868,5\right\},\left\{1871.5,11\right\},\left\{1872.5,69\right\},\left\{1873,51\right\},\left\{1875,52\right\},\left\{1876,100\right\},\left\{1877.5,73\right\},\left\{1881.5,7\right\},\left\{1882.5,9\right\},\left\{1884,43\right\},\left\{1885,130\right\},\left\{1887,133\right\},\left\{1889.5,19\right\},\left\{1892,62\right\},\left\{1893,65\right\},\left\{1895,86\right\},\left\{1897.5,11\right\},\left\{1899.5,6\right\},\left\{1901,12\right\},\left\{1902,4\right\},\left\{1904,67\right\},\left\{1906.5,24\right\},\left\{1907.5,23\right\},\left\{1909,24\right\},\left\{1910,50\right\},\left\{1912.5,72\right\},\left\{1916.5,3\right\},\left\{1919,9\right\},\left\{1921,57\right\},\left\{1923,77\right\},\left\{1925,72\right\},\left\{1926,10\right\},\left\{1928,4\right\},\left\{1930,7\right\}\right\}\right];$

$\mathrm{ilves}=\mathrm{Table}\left[\left\{\left\{1855,29\right\},\left\{1857,28\right\},\left\{1861,4\right\},\left\{1862.5,4\right\},\left\{1866,68\right\},\left\{1867.5,72\right\},\left\{1869,12\right\},\left\{1872,7\right\},\left\{1873,33\right\},\left\{1875,35\right\},\left\{1876,48\right\},\left\{1879,44\right\},\left\{1881,11\right\},\left\{1882,11\right\},\left\{1884.5,50\right\},\left\{1885.5,78\right\},\left\{1887,80\right\},\left\{1889,39\right\},\left\{1892,21\right\},\left\{1893,46\right\},\left\{1895,52\right\},\left\{1901,9\right\},\left\{1903,20\right\},\left\{1905.5,62\right\},\left\{1907,55\right\},\left\{1908,5\right\},\left\{1911,14\right\},\left\{1912,36\right\},\left\{1915,43\right\},\left\{1917.5,6\right\},\left\{1920,6\right\},\left\{1925,45\right\},\left\{1927,52\right\},\left\{1930,19\right\}\right\}\right];$

$\mathrm{MultipleListPlot}\left[\mathrm{janis},\mathrm{ilves},\mathrm{PlotJoined}->\mathrm{True},\mathrm{GridLines}->\mathrm{Automatic},\mathrm{Frame}->\mathrm{True},\mathrm{PlotLegend}->\left\{''Jänis'',''Ilves''\right\},\mathrm{PlotStyle}->\left\{\mathrm{Hue}\left[\text{.95}\right],\mathrm{Hue}\left[\text{.6}\right]\right\}\right]$

Out[10]=

$⁃\mathrm{Graphics}⁃$

Havaitaan, että vaikka luonnollisia kannanvaihteluja aiheuttavat myös muut tekijät, on yhdeksän vuoden välein esiintyvä maksimi molemmille kannoille varsin säännöllinen. Volterran–Lotkan yhtälö näyttää siis kuvaavan jotain olennaista saalis-peto-suhteessa olevien kantojen keskinäisestä vuorovaikutuksesta. Toisaalta se ei kerro kovinkaan monesta kannan suuruuteen vaikuttavasta tekijästä.

Tarkastellaan vielä peto- ja saaliskantojen käyttäytymistä faasitasossa. Ensin määritellään kantojen kasvunopeudet ajan funktioina:

In[11]:=

$\mathrm{kasvut}=\left\{D\left[\mathrm{maarat}\left[\left[1\right]\right],t\right],D\left[\mathrm{maarat}\left[\left[2\right]\right],t\right]\right\}$

Out[11]=

$\left\{\mathrm{InterpolatingFunction}\left[\left\{\left\{0.,80.\right\}\right\},<>\right]\left[t\right],\mathrm{InterpolatingFunction}\left[\left\{\left\{0.,80.\right\}\right\},<>\right]\left[t\right]\right\}$

Piirretään kuvaajat faasitasoon. x-akseli kuvastaa kannan kasvunopeutta ja y-akseli kannan suuruutta. Jäniskanta on piirretty jälleen punaisella ja ilveskanta sinisellä käyrällä.  

In[12]:=

$\mathrm{kuva1}=\mathrm{ParametricPlot}\left[\left\{\mathrm{kasvut}\left[\left[1\right]\right],\mathrm{maarat}\left[\left[1\right]\right]\right\},\left\{t,0,10\right\},\left\{\mathrm{DisplayFunction}->\mathrm{Identity},\mathrm{PlotStyle}->\left\{\mathrm{Hue}\left[\text{.95}\right]\right\}\right\}\right];$

$\mathrm{kuva2}=\mathrm{ParametricPlot}\left[\left\{\mathrm{kasvut}\left[\left[2\right]\right],\mathrm{maarat}\left[\left[2\right]\right]\right\},\left\{t,0,10\right\},\left\{\mathrm{DisplayFunction}->\mathrm{Identity},\mathrm{PlotStyle}->\left\{\mathrm{Hue}\left[\text{.6}\right]\right\}\right\}\right];$

$\mathrm{Show}\left[\mathrm{kuva1},\mathrm{kuva2},\mathrm{AspectRatio}->\mathrm{Automatic},\mathrm{DisplayFunction}->\mathrm{\$DisplayFunction}\right]$

Out[14]=

$⁃\mathrm{Graphics}⁃$

Tehtäviä

Mitä tämä faasitason kuvaaja kertoo kantojen käyttäytymisestä? Tarkastele, miten eri parametrien muuttaminen vaikuttaa faasitason kuvaajaan.

Tutki myös, onko olemassa alkuehtoja, joiden avulla saavutetaan tasapaino peto- ja saaliskannan välillä. Vihje: Tarkastele alkuehtoa, jossa $s\text{'}\left(t\right)=p\text{'}\left(t\right)=0$.

Linkkejä

differentiaaliyhtälöryhmä
faasitaso

JP & SKK 09.05.2001

Created by Mathematica  (July 21, 2004)