12.09.2013 Views

Lotka-Volterramodellen - Home Page of Lars Holm Jensen

Lotka-Volterramodellen - Home Page of Lars Holm Jensen

Lotka-Volterramodellen - Home Page of Lars Holm Jensen

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!