🐇 Lotka-Volterra — dynamika populacji drapieżnik-ofiara
Ustaw współczynniki narodzin, śmierci i drapieżnictwa dla królików i lisów, a potem obserwuj oscylacje populacji, orbity zamknięte i nullkliny na portrecie fazowym.
O modelu drapieżnik-ofiara Lotki-Volterry
Ta symulacja animuje klasyczne równania Lotki-Volterry, parę sprzężonych równań różniczkowych opisujących, jak populacje ofiar (króliki) i drapieżników (lisy) rosną i maleją w czasie. Ofiary rosną w tempie r i są konsumowane proporcjonalnie do liczby spotkań; drapieżnicy rosną, jedząc ofiary, i maleją w wyniku naturalnej śmierci. Trajektoria jest całkowana numerycznie metodą Rungego-Kutty czwartego rzędu (RK4) z krokiem 0,02 roku.
Suwaki ustawiają tempo wzrostu ofiar r, tempo drapieżnictwa α, efektywność konwersji drapieżnika β, współczynnik śmiertelności drapieżnika m oraz populacje początkowe H₀ i P₀. Opcje widoku przełączają między krzywymi populacji w funkcji czasu, portretem fazowym z nullklinami lub oboma widokami jednocześnie. Takie modele leżą u podstaw ekologii teoretycznej, zarządzania rybołówstwem i historycznych badań zapisów odłowów futrzarskich, takich jak cykle rysia i zająca.
Najczęściej zadawane pytania
Co pokazuje ta symulacja?
Pokazuje model drapieżnik-ofiara Lotki-Volterry, w którym populacja ofiar i populacja drapieżników oscylują w powtarzających się cyklach. Można obserwować populacje w czasie, zobaczyć zamkniętą orbitę na płaszczyźnie fazowej lub wyświetlić oba widoki jednocześnie.
Jakie są dwa równania stojące za modelem?
Ofiary zmieniają się według dH/dt = r·H − α·H·P, a drapieżnicy według dP/dt = β·α·H·P − m·P. Ofiary rosną wykładniczo, ale są zjadane przy kontakcie z drapieżnikami; drapieżnicy rosną, zamieniając zjedzone ofiary w potomstwo, i giną w stałym tempie.
Co robią elementy sterujące?
Suwaki ustawiają tempo wzrostu ofiar r, tempo drapieżnictwa α, efektywność drapieżnika β, współczynnik śmiertelności drapieżnika m oraz początkowe liczby ofiar i drapieżników. Regulacja dowolnej wartości zmienia kształt cyklu i przesuwa punkt równowagi pokazany w panelu.
Gdzie znajduje się punkt równowagi?
Nietrywialna równowaga to H* = m/(β·α) dla ofiar i P* = r/α dla drapieżników. Przy tych wartościach obie populacje przestają się zmieniać. Symulacja zaznacza ten punkt na portrecie fazowym i podaje H* oraz P* w panelu równowagi.
Dlaczego populacje wciąż oscylują, zamiast się ustabilizować?
Równowaga jest neutralnym centrum, a nie atraktorem, więc trajektorie tworzą zamknięte pętle, które ani nie rosną, ani nie zanikają. Każdy warunek początkowy wyznacza własną orbitę, a amplituda cyklu zależy od tego, jak daleko populacje początkowe znajdują się od równowagi. W podstawowym modelu nie ma tłumienia.
Dlaczego szczyty drapieżników opóźniają się względem szczytów ofiar?
Ofiary muszą stać się liczne, zanim drapieżnicy będą mieli wystarczająco pożywienia, by się rozmnażać, więc fala drapieżników rośnie po fali ofiar. Rosnąca liczba drapieżników następnie załamuje populację ofiar, co z kolei głodzi drapieżników. To opóźnienie fazowe o ćwierć cyklu jest charakterystycznym kształtem oscylacji.
Czym jest portret fazowy i nullkliny?
Portret fazowy przedstawia ofiary na osi poziomej względem drapieżników na osi pionowej, więc cykl pojawia się jako zamknięta pętla. Nullkliny to linie, na których każda populacja jest chwilowo stała: nullklina ofiar to pozioma linia P = r/α, a nullklina drapieżników to pionowa linia H = m/(β·α).
Jak dokładny biologicznie jest ten model?
Oddaje jakościowe wzajemne oddziaływanie drapieżników i ofiar, ale wprowadza istotne uproszczenia. Pomija pojemność środowiska, ruch przestrzenny, strukturę wiekową, pokarm alternatywny i zdarzenia losowe, a jego idealnie powtarzające się cykle są nierealistycznie kruche. Ekolodzy traktują go jako koncepcyjny punkt wyjścia, a nie precyzyjną prognozę.
Jaka metoda numeryczna jest używana i dlaczego?
Symulacja całkuje równania metodą Rungego-Kutty czwartego rzędu (RK4) z krokiem czasowym 0,02 roku. RK4 jest znacznie dokładniejsza niż prosty krok Eulera przy tym samym nakładzie, co zapobiega sztucznemu dryfowaniu zamkniętych orbit do wewnątrz lub na zewnątrz w ciągu osiemdziesięcioletniego przebiegu.
Jak szacowany jest okres cyklu?
Kod wykrywa kolejne szczyty krzywej ofiar i uśrednia czas między nimi, aby podać okres cyklu w latach. W pobliżu równowagi okres zbliża się do 2π podzielonego przez pierwiastek kwadratowy z r razy m, ale wyświetlana wartość jest mierzona bezpośrednio z symulowanej trajektorii.
Gdzie model Lotki-Volterry jest stosowany w prawdziwym świecie?
Pojawia się w ekologii i planowaniu ochrony przyrody, rybołówstwie i zarządzaniu szkodnikami, a także w klasycznej analizie zapisów odłowów rysia i zająca znad Zatoki Hudsona. Matematycznie podobne równania opisują też oscylujące reakcje chemiczne, cykle koniunkturalne w gospodarce oraz niektóre systemy epidemiologiczne.