Iki Tür Modeli
· Iki Tür Modeli
P
1ve P
2nüfuslu, iki türden olu¸san küçük bir ekosistem göz önüne alal¬m.
Her bir türün de¼ gi¸sim oran¬n¬n sadece türlerin nüfusuna ba¼ g¬ml¬oldu¼ gunu kabul edelim. O halde, matematiksel modelimiz
dP
1dt = g ( P
1, P
2) (1)
dP
2dt = f ( P
1, P
2) (2)
¸seklindedir.
Genel olarak bir türün etkisi, di¼ ger türün nüfusunu art¬rmak veya azaltmak
¸seklindedir. Böylece iki tür aras¬nda üç tip etkile¸sim olu¸sur.
+ ( +) av-avc¬; Bir tür di¼ gerinin nüfusunu azalt¬yor.
++ (ortakl¬k); Her iki tür birbirinin nüfusunu art¬r¬yor.
· Iki tür için verilen (1)-(2) modelinde, e¼ ger hiç bir türün d¬¸s göçü yoksa bu durumda
g ( 0, P
2) = 0 ve f ( P
1, 0 ) = 0 olur. (1)-(2) sisteminin faz düzlem denklemi,
dP
2dP
1= f ( P
1, P
2)
g ( P
1, P
2) (3)
iki-nüfus ekosisteminin fazlar¬ olarak adland¬r¬l¬r. Bu denklemin çözümü nüfuslar¬n yörüngelerini verir.
g ( P
1, P
2) ve f ( P
1, P
2) zamana aç¬k olarak ba¼ gl¬olmad¬klar¬ndan, verilen
model otonomdur. Nüfuslar¬n etkile¸simini çal¬¸smak için ya (1)-(2) sistemini
ya da (3) denklemini analiz etmek gerekir. Ço¼ gu zaman her ikisinin de aç¬k
bir çözümünü bulmak mümkün de¼ gildir.
Iki Tür Modeli Kararl¬l¬k
· Iki Tür Modeli - Kararl¬l¬k
Denge nüfusu-her iki nüfusun da zamanla de¼ gi¸smedi¼ gi nüfus.
(do¼ gum say¬s¬=ölüm say¬s¬) ( P
1e, P
2e) bir denge nüfusu ise
f ( P
1e, P
2e) = 0 (4)
g ( P
1e, P
2e) = 0 (5)
demektir. Denge durumunda, dP
2/dP
1belirsizdir. Böylece, ( P
1e, P
2e) noktas¬faz düzlem denkleminin bir ayk¬r¬(tekil) noktas¬d¬r.
Denge noktas¬(nüfusu) kararl¬m¬d¬r?
E¼ ger her iki nüfus da denge nüfuslar¬na yak¬n ise, zaman içinde ne
olur?
Kararl¬l¬k analizi: 0 < j ε j 1 için
P
1( t ) = P
1e+ εP
11( t ) P
2( t ) = P
2e+ εP
21( t )
yi (1)-(2) sisteminde yerlerine yaz¬p, Taylor aç¬l¬m¬nda e
2li terimleri atarsak
dP
11( t )
dt = ∂g ( P
1e, P
2e)
∂P
1P
11( t ) + ∂g ( P
1e, P
2e)
∂P
2P
21( t ) (6) dP
21( t )
dt = ∂f ( P
1e, P
2e)
∂P
1P
11( t ) + ∂f ( P
1e, P
2e)
∂P
2P
21( t ) (7)
Iki Tür Modeli Kararl¬l¬k
(6)-(7) sistemi birinci basamaktan bir lineer denklem sistemidir ve çözümü denge noktas¬kom¸sulu¼ gunda türlerin davran¬¸s¬n¬verir. Bu denklem
sistemini matris formunda
dP dt = JP
¸seklinde de yazabiliriz. Burada J
J = 2 6 4
∂g(P1e,P2e)
∂P1
∂g(P1e,P2e)
∂P2
∂f(P1e,P2e)
∂P1
∂f(P1e,P2e)
∂P2
3 7 5
olup, sistemin Jakobiyen matrisidir.
Sistem (6)-(7) daha basit bir gösterimle dx
dt = ax + by (8)
dy
dt = cx + dy (9)
veya matris formda
dx dydt dt
= dX
dt = JX = a b c d
x y
¸seklinde yaz¬labilir. Burada, a, b, c, d sabitler x ve y de denge noktas¬ndan
sapma olarak al¬nabilir.
Iki Tür Modeli Kararl¬l¬k
c = 0 ise (9) denklemi çözülebilir olup, çözümü y ( t ) = k
1e
dt, ( k
1sabit ) ve
x ( t ) = 8 >
<
> :
k
2e
at+ bk
1d a e
dt, d 6= a k
2e
at+ k
1bte
at, d = a
( k
2sabit )
c 6= 0 ise, (9) denkleminde türev alarak, d
2y
dt
2( a + d ) dy
dt + ( ad bc ) y = 0 (10)
sabit katsay¬l¬denklemi elde edilir ki karakteristik denklemi r
2( a + d ) r + ( ad bc ) = 0 veya
r
2iz ( J ) r + det ( J ) = 0 olup, dikkat edilirse bu denklem dX
dt = JX sisteminde J matrisinin karakteristik denklemidir. Kökleri ise
r
1,2= a + d p
( a + d )
24 ( ad bc ) 2
olup J nin özde¼ gerlerini verir.
Iki Tür Modeli Kararl¬l¬k
Durum 1. r
1ve r
2reel ve r
16= r
2. (10) denkleminin genel çözümü
y ( t ) = k
3e
r1t+ k
4e
r2t, ( k
3, k
4isteksel sabitler ) (11) olup, (9) e¸sitli¼ ginden,
x ( t ) = k
3( r
1d )
c e
r1t+ k
4( r
2d )
c e
r2t(12)
elde edilir. O halde (11) ve (12), t ! ∞ için ( x ( t ) , y ( t )) ! ( 0, 0 ) , ( r
1, r
2< 0 )
xy –de bir nokta, ( bir kök < 0, di¼ ger kök = 0 )
verir.
Durum 2. r
1= r
2.
y ( t ) = k
3e
r1t+ k
4te
r1t, ( k
3, k
4isteksel sabitler ) (13)
x ( t ) = k
3( r
1d )
c e
r1t+ k
4( 1 + r
1t dt )
c e
r1t(14)
elde edilir. Böylece, t ! ∞ için
( x ( t ) , y ( t )) ! ( 0, 0 ) , ( r
1< 0 )
∞ da bir nokta, ( r
1= 0 )
olur.
Iki Tür Modeli Kararl¬l¬k
Durum 3. Kompleks e¸slenik kökler: r
1= α + i β, r
2= α i β.
y ( t ) = k
3e
αtcos βt + k
4e
αtsin βt (15)
x ( t ) = k
3e
αt
c ( α cos βt β sin βt d cos βt ) + k
4e
αt
c ( α sin βt + β cos βt d sin βt ) (16) bulunur. t ! ∞ için
( x ( t ) , y ( t )) ! ( 0, 0 ) , ( α < 0 )
elde edilir.
Yukar¬daki durumlar¬göz önüne alarak a¸sa¼ g¬daki sonuçlar¬elde ederiz:
Bir denge çözümü,
(i) kararl¬d¬r e¼ ger her iki kökün reel k¬s¬mlar¬negatif ise, (ii) etkisiz kararl¬d¬r (yani çözüm t ! ∞ için bir denge
noktas¬na yakla¸smaz), e¼ ger her iki kök s¬f¬r de¼ gilse, ve bir kökün reel k¬sm¬s¬f¬r ve di¼ geri küçük veya e¸sit s¬f¬r ise (Durum 1 ve 3),
(iii) cebirsel karars¬zd¬r (yani çözüm t nin bir cebirsel kuvveti ile büyür) e¼ ger her iki kök de s¬f¬r ise (Durum 2, dengenin etkisiz kararl¬oldu¼ gu a = b = c = d = 0 hariç),
(iv) karars¬zd¬r e¼ ger bir kökün reel k¬sm¬pozitif ise.
Iki Tür Modeli Kararl¬l¬k
Bu veriler ¬¸s¬¼ g¬nda a¸sa¼ g¬daki tabloyu olu¸sturabiliriz:
r = a + d p
( a + d )
24 ( ad bc ) 2
ifadesinde p = a + d , q = ad bc, 4 = p
24q olsun.
DURUM KÖKLER KARARLILIK
I. 4 > 0 Reel ve farkl¬
A. q > 0 ayn¬i¸saret
( a ) p > 0 her ikisi de + Karars¬z ( b ) p < 0 her ikisi de Kararl¬
B. q = 0 bir kök s¬f¬r
( a ) p > 0 di¼ ger kök + Karars¬z ( b ) p < 0 di¼ ger kök Etkisiz kararl¬
C. q < 0 z¬t i¸saretler Karars¬z
II. 4 = 0 Reel ve e¸sit
A. p > 0 her ikisi de + Karars¬z
B. p = 0
( a ) en az bir
a, b, c, d 6= 0 her ikisi de s¬f¬r Cebirsel karars¬z ( b ) a, b, c, d = 0 her ikisi de s¬f¬r Etkisiz kararl¬
C. p < 0 her ikisi de Kararl¬
III. 4 < 0 Karma¸s¬k e¸slenik
A. p > 0 reel k¬s¬m + Karars¬z
B. p = 0 reel k¬s¬m s¬f¬r Etkisiz kararl¬
C. p < 0 reel k¬s¬m Kararl¬
Iki Tür Modeli Kararl¬l¬k