Lotka-Volterramodellen - Home Page of Lars Holm Jensen
Lotka-Volterramodellen - Home Page of Lars Holm Jensen
Lotka-Volterramodellen - Home Page of Lars Holm Jensen
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
4.1 Linearisering 59<br />
f2(x1, x2) omkring ligevægtspunktet opn˚as en god approksimation af f1(x1, x2)<br />
og f2(x1, x2) nær (x1 ∗ , x2 ∗ ) jf. sætning A.8. Vi bestemmer Taylorpolynomiet 1.<br />
orden, da vi ellers vil f˚a ikke-lineære led med. S˚aledes vil vi ikke være kommet<br />
nærmere en løsning til den oprindelige ikke-lineære differentialligning.<br />
f1(x1, x2) = f1(x1 ∗ , x2 ∗ ) + ∂f1(x1 ∗ , x2 ∗ )<br />
+ ∂f1(x1 ∗ , x2 ∗ )<br />
(x2 − x2<br />
∂x2<br />
∗ )<br />
∂x1<br />
f2(x1, x2) = f2(x1 ∗ , x2 ∗ ) + ∂f2(x1 ∗ , x2 ∗ )<br />
+ ∂f2(x1 ∗ , x2 ∗ )<br />
(x2 − x2<br />
∂x2<br />
∗ )<br />
∂x1<br />
(x1 − x1 ∗ )<br />
(x1 − x1 ∗ )<br />
(4.7)<br />
For at vurdere hvorvidt det er en god approksimation, er det nødvendigt at se<br />
p˚a restleddet:<br />
R f1,x∗<br />
2 (x) = 1<br />
2! D(2) f(c, x − x ∗ )<br />
R f2,x∗<br />
2 (x) = 1<br />
2! D(2) g(c, x − x ∗ )<br />
Eftersom vi krævede, at f er C 2 , kan vi tillade os at vurdere p˚a restleddet.<br />
Punktet (x1, x2) ligger tæt p˚a ligevægtspunktet, derfor vil (x1 − x1 ∗ ) 2 være<br />
endnu mindre, n˚ar x1 − x1 ∗ < 1. Da punktet c er valgt mellem x og x ∗ , vil de<br />
forskellige partielle afledede af 2. orden i c ogs˚a blive meget sm˚a. Dermed bliver<br />
restleddet meget lille.<br />
Eftersom (x1 ∗ , x2 ∗ ) er ligevægtspunktet, giver leddene f(x1 ∗ , x2 ∗ ) og g(x1 ∗ , x2 ∗ )<br />
i ligningerne (4.7) begge nul.<br />
Vi opn˚ar s˚aledes følgende approksimation af f(x1, x2) og g(x1, x2):<br />
f(x1, x2) = ∂f(x1 ∗ , x2 ∗ )<br />
∂x1<br />
g(x1, x2) = ∂g(x1 ∗ , x2 ∗ )<br />
∂x1<br />
Dette kan ogs˚a skrives som<br />
⎡<br />
x ′ ⎢<br />
= ⎣<br />
∂f(x1 ∗ , x2 ∗ )<br />
∂x1<br />
∂g(x1 ∗ , x2 ∗ )<br />
∂x1<br />
(x1 − x1 ∗ ) + ∂f(x1 ∗ , x2 ∗ )<br />
(x2 − x2<br />
∂x2<br />
∗ )<br />
(x1 − x1 ∗ ) + ∂g(x1 ∗ , x2 ∗ )<br />
(x2 − x2<br />
∂x2<br />
∗ )<br />
∂f(x1 ∗ , x2 ∗ )<br />
∂x2<br />
∂g(x1 ∗ , x2 ∗ )<br />
∂x2<br />
⎤<br />
⎥<br />
⎦ (x − x ∗ )<br />
Da x ∗ er en konstant vektor, er det muligt at erstatte x ′ med [x − x ∗ ] ′ .<br />
Ved at sætte y = x − x ∗ opn˚as følgende linearisering af systemet x ′ = f(x):<br />
⎡<br />
y ′ ⎢<br />
= ⎣<br />
∂f(x1 ∗ , x2 ∗ )<br />
∂x1<br />
∂g(x1 ∗ , x2 ∗ )<br />
∂x1<br />
∂f(x1 ∗ , x2 ∗ )<br />
∂x2<br />
∂f(x1 ∗ , x2 ∗ )<br />
∂x2<br />
⎤<br />
⎥<br />
⎦ y