DIFFERENTIALLIGNINGER NOTER TIL CALCULUS 2003 AARHUS UNIVERSITET H.A. NIELSEN
I NDHOLD 1. Lineær ligning 2. Lineært system 3. Generel ligning 4. Stabilitet Litteratur
2 8 16 18 21
Noterne er til 4 timers forelæsninger i Calculus 2003. De supplerer og præciserer lærebogen "J.Stewart - Calculus, concepts and contexts, 2nd ed." Fremstillingen er bevidst kortfattet og præcis, fremfor beskrivende og detaljeret. Dette kræver en lidt større indsats ved første læsning, men vil betale sig i overskuelighed senere. Det forudsættes, at noterne suppleres med både lærebogen, forelæsninger og mange øvelser. I afsnit 1 behandles den generelle 1. ordens lineære differentialligning. Løsningsmængdens lineære struktur præciseres. Løsningen omskrives til et stamfunktionsproblem og en samlet formel angives. Videre i afsnit 2 beskrives det lineære 1. ordens differentialligningssystem i vilkålig mange variable. Kun ligningssystemet med konstante koefficienter behandles eksplicit. Løsningen kan ved inddragelse af matrixmetoder og egenvektorer reduces til én variabel metoder. Ikke-lineære differentialligninger er vanskelige at løse, men i afsnit 3 omtaler vi en entydighedsog eksistenssætning for 1. ordens differentialligninger, som i mange tilfælde sikrer en løsning. I sidste afsnit 4 introduces ligevægt og stabilitet. I omegnen af sådanne punkter kan en ikke-lineær differentialligning i nogle tilfælde med fordel tilnærmes med en lineær differentialligning. Dette gælder også for systemer af differentialligninger.
Version: 1.0. 15. november, 2003. 1
2
1. L INEÆR LIGNING
1. L INEÆR
LIGNING
Den lineære 1. ordens differentialligning er den simpleste. Men den spiller en stor rolle, da den i princippet nemt kan løses og løsningen kan bruges til at tilnærme løsningen af en mere vanskelig differentialligning. Den lineære differentialligning defineres rimeligt præcist og den lineære struktur af løsningsmængden angives. Ligningen med konstante koefficienter er separabel og løses ved integration. Den generelle ligning reduceres på analog måde og en samlet formel angives. Resultaterne anvendes på en populær opgavetype. En rimelig præcis definition af den lineære ligning og lidt almindelig sprogbrug, der bruges i de fleste fremstillinger, følger her. Definition 1.1. Den lineœre 1. ordens differentialligning er dy = a(x)y + b(x) dx En partikulær løsning er en differentiabel funktion y(x) som opfylder y 0 (x) = a(x)y(x) + b(x) Den fuldstœndige løsning er en angivelse af alle løsninger, også kaldet løsningsrumdy met. Ligningen dx = a(x)y kaldes homogen og er den homogene part af den inhomogene, b 6= 0, ligning ovenfor. Den lineære differentiallignings form har en afgørende betydning for strukturen af løsningsrummet. I det homogene tilfælde er linearkombinationer af løsninger igen løsninger. Det kaldes i anvendelsessammenhænge ofte for ”superpositionsprincippet”. Formuleringen af følgende sætning er inspireret af lineær algebra. Det ses, at løsningen af det inhomogene problem reduceres til det tilsvarende homogene problem samt angivelse af bare én partikulær løsning. Sætning 1.2. Hvis z1 (x), z2 (x) er løsninger til den homogene lineœre differentialligning dy = a(x)y dx så er enhver linearkombination z(x) = C1 z1 (x) + C2 z2 (x) også en løsning. Hvis z0 (x) er en løsning til den inhomogene lineœre differentialligning dy = a(x)y + b(x) dx så er enhver løsning af formen y(x) = z(x) + z0 (x) hvor z(x) er en løsning til den homogene part af systemet.
1. L INEÆR LIGNING
3
Bevis. z 0 = C1 z10 + C2 z20 = C1 az1 + C2 az2 = az giver første del. For anden del ses, at y − z0 er løsning til den homogene part.
Ligningen med konstante koefficienter er særlig nem. Den er separabel og løses ved integration. Ved brug af den lineære struktur kan løsningen opdeles i det homogene problem og angivelse af bare én partikulær løsning. Sætning 1.3. Den lineœre ligning med konstante koefficienter dy = ay + b dx har fuldstœndig løsning givet ved a = 0: y(x) = C + bx a 6= 0: b y(x) = Ceax − a hvor C er arbitrœr. Specielt er den konstante funktion y(x) = −b/a en løsning. Bevis. Den homogene part dy = ay dx
er separabel med løsninger
Z dy = a(x)dx y ln |y| = ax + K
Z
y(x) = Ceax Afslut ved Sætning 1.2.
y
1
0
1
x
Grafer af løsninger Samme metode som anvendt på ligningen med konstante koefficienter kan bruges på den generelle homogene ligning. Denne er igen separabel og løses ved integration i følgende sætning.
4
1. L INEÆR LIGNING
Sætning 1.4. Den homogene lineœre ligning dy = a(x)y dx har fuldstœndig løsning y(x) = CeA(x) hvor C er arbitrœr og A(x) = Bevis.
Z
ea(x) dx
dy = a(x)y dx
er separabel med løsninger Z dy = a(x)dx y ln |y| = A(x) + K
Z
y(x) = CeA(x) På snedig vis reduceres den inhomogene ligning til et stamfunktionsproblem. Der opnås en færdig formel for den fuldstændige løsning. Det er hovedresultatet i dette afsnit. Efterfølgende samles fremgangsmåden i en klar metode. Sætning 1.5. Den generelle lineœre ligning dy = a(x)y + b(x) dx har fuldstœndig løsning y(x) = CeA(x) + B(x)eA(x) hvor C er arbitrœr og A(x) =
Z
a(x)
e
dx,
B(x) =
Z
e−A(x) b(x) dx
Bevis. z(x) = e−A(x) y(x) opfylder ligningen dz = e−A(x) b(x) dx som integreres til z(x) = B(x) + C og forlænges til y(x) = CeA(x) + B(x)eA(x)
1. L INEÆR LIGNING
5
Læg mærke til den efterfølgende metode, som med fordel kan bruges i mange populære opgavetyper. Bemærkning 1.6 (Metode). dy = a(x)y + b(x) dx 1. Bestem en stamfunktion Z
A(x) =
ea(x) dx
2. Bestem en stamfunktion B(x) =
Z
e−A(x) b(x) dx
3. Skriv løsningen y(x) = CeA(x) + B(x)eA(x) 4. Konstanten C fastlægges ved indsættelse i løsningen fra 3. Metoden giver altså den fuldstændige løsning samt eventuelt en partikulær løsning der tilfredsstiller yderligere betingelser. To repræsentative opgaver af en ofte stillet type løses ved brug af resultater og metoder fra dette afsnit. Opgave 1.7 (Opgave 7, Matematik Alfa 1, August 2002). Angiv den fuldstændige løsning til differentialligningen y 0 + 2y = xe−2x + 3 Angiv endvidere den partikulære løsning y(x), der opfylder y(0) = 2. Løsning. a(x) = −2, b(x) = xe−2x + 3
A(x) =
Z
a(x) dx =
Z
−2 dx = −2x Z Z B(x) = e−A(x) b(x) dx = e2x (xe−2x + 3)dx 1 3 = x2 + e2x 2 2 Som giver fuldstændig løsning y(x) = CeA(x) + B(x)eA(x) 1 3 = Ce−2x + ( x2 + e2x )e−2x 2 2 1 3 = Ce−2x + x2 e−2x + 2 2
6
1. L INEÆR LIGNING
I den partikulære løsning bestemmes C ved y(0) = 2. y(0) = Ce0 +
3 =2 2
I alt er 1 3 1 y(x) = e−2x + x2 e−2x + 2 2 2 y
1
0
x
1
Grafen af løsningen
Opgave 1.8 (Opgave 5, Matematik Alfa 1, Januar 2003). Angiv den fuldstændige løsning y(x) til differentialligningen (for x > 0) y y 0 + = 2x−1 . x Angiv endvidere den løsning, der opfylder betingelsen y(2) = 5. Løsning. 2 1 a(x) = − , b(x) = x x 1 − dx = − ln x x Z Z 2 −A(x) B(x) = e b(x) dx = eln x dx x Z Z 2 = x dx = 2 dx x = 2x A(x) =
Z
a(x) dx =
Z
Som giver fuldstændig løsning y(x) = CeA(x) + B(x)eA(x) = Ce− ln x + 2xe− ln x 1 =C +2 x
1. L INEÆR LIGNING
7
I den partikulære løsning bestemmes C ved y(0) = 2. 1 y(2) = C + 2 = 5 2 I alt er 6 y(x) = + 2 x y
1
0
x
1
Grafen af løsningen
Opgave 1.9. Angiv den fuldstændige løsning til differentialligningen dy =y+1 dx Angiv endvidere den partikulære løsning y(x), der opfylder y(0) = 1. Opgave 1.10. Angiv den fuldstændige løsning til differentialligningen y0 = y − x Angiv endvidere den partikulære løsning y(x), der opfylder y(0) = 1. Opgave 1.11. Angiv den fuldstændige løsning til differentialligningen dy 2 = 2xy + ex dx Angiv endvidere den partikulære løsning y(x), der opfylder y(1) = e + 1. Opgave 1.12. Angiv for alle a 6= 0 den fuldstændige løsning til differentialligningen y 0 = ay + ex
8
2. L INEÆRT SYSTEM
2. L INEÆRT
SYSTEM
Det lineære differentialligningssystem er en umiddelbar, men meget kraftig udvidelse af den lineære differentiallingning. Kun tilfældet med konstante koefficienter behandles. Ved indragelse af matrixmetoder og egenvektorer/værdier kan et system reduceres til lineære ligninger, som kan løses ved metoder fra afsnit 1. En præcis definition efterfølges af en sætning om løsningsrummets lineære struktur. En egenvektor giver en 1-parameter mængde af løsninger. For en diagonaliserbar matrix kan den fuldstændige løsning angives. Et par opgaveforslag behandles. Den præcise definition og gængs sprogbrug er en umiddelbar udvidelse af tilsvarende definition og sprogbrug i afsnit 1. Tilfældet med ikke-konstante koefficienter overlades til læseren. Definition 2.1. Ved et lineœrt 1. ordens differentialligningssystem med konstante koefficienter forstås dy1 dx dy2 dx
= a11 y1 + . . . = a21 y1 + . . . .. .
+ a1n yn + b1 + a2n yn + b2
dyn dx
= an1 y1 + . . .
+ ann yn + bn
En partikulær løsning er differentiable funktioner x 7→ y1 (x), . . . , x 7→ yn (x) som indsat opfylder lignningerne. Løsningsrummet, den fuldstændige løsning er angivelsen af alle løsninger. For n × n−matricen A = (aij ), koefficientmatricen, og n−søjlerne b = (bi ), y(x) = (yi (x)) skrives det lineœre differentialligningssystem dy = Ay + b dx En løsning skrives
y1 (x) x 7→ y(x) = ... yn (x)
Bemærkning 2.2. Givet n × n−matricen A = (aij ) og n−søjlerne b = (bi ), y(x) = (yi (x)) kaldes systemet dy = Ay dx homogent og er den homogene part af det inhomogene, b 6= 0, system dy = Ay + b dx Den lineære struktur af løsningsrummet går igen fra den lineære ligning. I det homogene tilfælde er linearkombinationer af løsninger igen løsninger. Det ses,
2. L INEÆRT SYSTEM
9
at løsningen af det inhomogene problem reduceres til det tilsvarende homogene problem samt angivelse af bare én partikulær løsning. Sætning 2.3. Betragt n × n−matricen A = (aij ) og n−søjlen y(x) = (yi (x)). Hvis z1 (x), z2 (x) er løsninger til det homogene lineœre differentialligningssystem dy = Ay dx så er enhver linearkombination z(x) = C1 z1 (x) + C2 z2 (x) også en løsning. Betragt yderligere n−søjlen b. Hvis z0 (x) er en løsninger til det inhomogene lineœre differentialligningssystem dy = Ay + b dx så er enhver løsning af formen y(x) = z(x) + z0 (x) hvor z(x) er en løsning til den homogene part af systemet. Bevis. Som for Sætning 1.2.
Inddragelse af teorien for egenværdier og egenvektorer forenkler teknikken betydeligt. En egenvektor giver en 1-parameter mængde af løsninger for det homogene system. Sætning 2.4. Betragt n × n−matricen A = (aij ) og n−søjlen y(x) = (yi (x)) samt det homogene lineœre differentialligningssystem dy = Ay dx Hvis u er en egenvektor for A med egenvœrdi λ, så er y(x) = Ceλx u løsninger, hvor C er arbitrær. Bevis. Gør prøve dy = λeλx u = eλx Au = Ay dx Det er ofte muligt at finde en konstant løsning. Kombineres dette med Sætning 2.3 og 2.4 kan en én parameter mængde af løsninger angives. Sætning 2.5. Betragt n×n−matricen A = (aij ) og n−søjlerne b = (bi ), y(x) = (yi (x)) samt det lineœre differentialligningssystem dy = Ay + b dx
10
2. L INEÆRT SYSTEM
En konstant funktion y(x) = v er en løsning, hvis Av = −b. Hvis yderligere u er en egenvektor for A med egenvœrdi λ, så er y(x) = Ceλx u + v løsninger, hvor C er arbitrær. Bevis. Gør prøve ved brug af Sætning 2.4 dy = eλx Au = A(y − v) = Ay + b dx For en diagonaliserbar koefficientmatrix kan man finde den fuldstændige løsning for det homogene tilfælde. Kan man samtidig finde en konstant løsning til det inhomogene problem, så kan den fuldstændige løsning også angives i dette tilfælde. Sætning 2.6. Betragt n × n−matricen A = (aij ) og n−søjlen y(x) = (yi (x)) samt det homogene lineœre differentialligningssystem dy = Ay dx Hvis matricen U med søjler u1 , . . . , un diagonaliserer A med egenvœrdier λ1 , . . . , λn , Auj = λj uj , så er den fuldstændige løsning givet ved y(x) = C1 eλ1 x u1 + · · · + Cn eλn x un hvor C1 , . . . , Cn er arbitrœre. Bevis. Fra Sætning 2.3, 2.4 følger, at linearkombinationerne er løsninger. Omvendt for en given løsning z, findes C1 , . . . , Cn så z(0) = C1 u1 + · · · + Cn un For y(x) = C1 eλ1 x u1 + · · · + Cn eλn x un gælder d(z − y) = 0, z(0) = y(0) dx hvoraf resultatet følger.
Sætning 2.7. Betragt n×n−matricen A = (aij ) og n−søjlerne b = (bi ), y(x) = (yi (x)) samt det lineœre differentialligningssystem dy = Ay + b dx En konstant funktion y(x) = v er en løsning, hvis Av = −b. Hvis matricen U med søjler u1 , . . . , un diagonaliserer A med egenvœrdier λ1 , . . . , λn , Auj = λj uj , så er den fuldstændige løsning givet ved y(x) = C1 eλ1 x u1 + · · · + Cn eλn x un + v hvor C1 , . . . , Cn er arbitrœre.
2. L INEÆRT SYSTEM
11
Følgende opgavetype overvejes som repræsentativ for resultaterne i dette afsnit. Opgave 2.8 (Opgave måske, Calculus 2, Januar 2004). Betragt differentialligningssystemet y10 = y1 + y2 y20 = 8y1 − y2 Det oplyses, at vektoren u = (1, 2) er en egenvektor for matricen 1 1 A= 8 −1 Angiv den løsning y(x) = (y1 (x), y2 (x)) der opfylder y(0) = u, altså (y1 (0), y2 (0)) = (1, 2). Løsning. Egenværdien λ = 3 fås af udregningen 1 1 1 3 Au = = = 3u 8 −1 2 6 I følge Sætning 2.4 er 1 y(x) = Ce 2 løsninger for alle valg af C. For C = 1 giver dette den ønskede løsning. Skrevet ud 3x
y1 (x) = e3x y2 (x) = 2e3x En mere komplet, men også ret omfangsrig opgave kunne være. Opgave 2.9. Angiv den fuldstændige løsning til differentialligningssystemet dy1 = y1 + 2y2 − 8 dx dy2 = 2y1 + y2 − 7 dx Løsning. Koefficientmatricen er 1 2 A= 2 1 Egenværdierne findes som rødder i det karakteristiske polynomium 1 − λ 2 |A − λI2 | = 2 1 − λ = λ2 − 2λ − 3
Egenværdierne er λ1 = −1,
λ2 = 3
12
2. L INEÆRT SYSTEM
Egenvektorer hørende til egenværdien −1: 2 2 1 1 A+I= ∼ 2 2 0 0 giver egenvektorer −x2 −1 x1 = = x2 1 x2 x2 Egenvektorer hørende til egenværdien 3: −2 2 1 −1 A − 3I = ∼ 2 −2 0 0
giver egenvektorer
x1 x2 1 = = x2 x2 x2 1 Den fuldstændige løsning til den homogene part dy1 = y1 + 2y2 dx dy1 = 2y1 + y2 dx er i følge Sætning 2.6 −x −1 3x 1 y(x) = C1 e + C2 e 1 1 Skrevet ud
y1 (x) = −C1 e−x + C2 e3x y2 (x) = C1 e−x + C2 e3x
En konstant løsning y(x) = v = (v1 , v2 ) skal opfylde 0 = v1 + 2v2 − 8 0 = 2v1 + v2 − 7 Dette løses 2 v= 3 Den fuldstændige løsning til systemet dy1 = y1 + 2y2 − 8 dx dy2 = 2y1 + y2 − 7 dx er i følge Sætning 2.7 2 −x −1 3x 1 y(x) = C1 e + C2 e + 1 1 3 Skrevet ud
y1 (x) = −C1 e−x + C2 e3x + 2 y2 (x) = C1 e−x + C2 e3x + 3
Hvor C1 , C2 er arbitrære konstanter.
2. L INEÆRT SYSTEM
13
Et hastighedsfelt giver en grafisk fornemmelse for løsningskurvernes x 7→ y(x) forløb. y2
y1
Hastighedsfelt I det ikke-diagonaliserbare tilfælde kan man lidt mere besværligt også finde løsningerne. Eksempel 2.10 (Ingen egenværdier). Betragt det lineære system y10 = y1 − y2 y20 = y1 + y2 Koefficientmatricen A=
1 −1 1 1
har karakteristisk polynomium λ2 − 2λ + 2 med diskriminant −4 og dermed ingen egenværdier. Ved brug af komplekse tal findes løsningen ved metoden med egenvektorer x cos x x − sin x y(x) = C1 e + C2 e sin x cos x Skrevet ud
y1 (x) = C1 ex cos x − C2 ex sin x y2 (x) = C1 ex sin x + C2 ex cos x y2
y1
Hastighedsfelt
14
2. L INEÆRT SYSTEM
Eksempel 2.11 (1 egenværdi). Betragt det lineære system y10 = 3y1 + y2 y20 = 3y2 Koefficientmatricen 3 1 A= 0 3 har egenværdi 3 og egenrum E3 = span(e1 ). Løsningen bestemmes ved metoden med egenvektorer og/eller substitution. 3x 1 3x x + C2 e y(x) = C1 e 0 1 Skrevet ud y1 (x) = C1 e3x + C2 e3x x y2 (x) = C2 e3x
y2
y1
Hastighedsfelt
Opgave 2.12. Betragt differentialligningssystemet dy1 = 3y1 + 2y2 dx dy1 = y1 + 4y2 dx Det oplyses, at vektoren u = (2, −1) er en egenvektor for matricen 3 2 A= 1 4 Angiv den løsning y(x) = (y1 (x), y2 (x)) der opfylder y(0) = −u, altså (y1 (0), y2 (0)) = (−2, 1)
2. L INEÆRT SYSTEM
15
Opgave 2.13. Betragt differentialligningssystemet y10 = −y1 + y2 y20 = y2 Det oplyses, at vektoren u = (1, 0) er en egenvektor for matricen −1 1 A= 0 1 Angiv den løsning y(x) = (y1 (x), y2 (x)) der opfylder y(0) = 2u, altså (y1 (0), y2 (0)) = (2, 0) Opgave 2.14. Betragt differentialligningssystemet y10 = 2y1 + 3y2 y20 = 3y1 + 2y2 Det oplyses, at vektorerne u1 = (1, 1), u2 = (1, −1) er en egenvektorer for systemets koefficientmatrix. Angiv den fuldstændige løsning. Opgave 2.15. Angiv den fuldstændige løsning til differentialligningssystemet y10 = 7y1 + 2y2 + 7 y20 = 3y1 + 8y2 − 3 Opgave 2.16. Angiv den fuldstændige løsning til det homogene differentialligningssystem dy1 = −y1 − 2y2 dx dy2 = y1 − 4y2 dx Angiv den løsning y(x) der opfylder y(0) = (1, −1).
16
3. G ENEREL LIGNING
3. G ENEREL
LIGNING
Emnet differentialligninger er meget omfattende. Det er kun i specialtilfælde muligt at angive løsninger ved elementære funktionsudtryk. For en ren matematisk behandling af differentialligninger, indføres en mere præcis definition af en ”differentialligning og en løsning” som er hensigtsmæssig for formulering og bevis af en såkaldt ”eksistens- og entydighedssætning”. Her er så en lidt mere præcis sprogbrug for differentialligninger. Formuleringen for differentialligningssystemer overlades til læseren. Definition 3.1. Lad I, J være åbne intervaller og F (x, y) : I × J → R en reel funktion. En løsning til 1. ordens differentialligningen dy = F (x, y) dx er en differentiabel funktion y(x) : I 0 → J på et åbent delinterval I 0 ⊆ I, som indsat giver y 0 (x) = F (x, y(x)), x ∈ I 0 Følgende ikke helt optimale sætning er ofte anvendelig til at sikre eksistens og entydighed af løsninger til 1. ordens differentialligninger. Sætning 3.2 (Eksistens og entydighed). Antag at F (x, y) er kontinuert og ∂F ∂y (x, y) eksisterer og er kontinuert i I × J. For et givet (x0 , y0 ) ∈ I × J findes entydigt bestemt et maximalt åbent delinterval I 0 ⊆ I om x0 og en differentiabel funktion y(x) : I 0 → J, som er en løsning til 1. ordens differentialligningen dy = F (x, y) dx og opfylder y(x0 ) = y0 Bevis. Se f.eks. [Iversen].
Bemærkning 3.3. Den udvidede ligning dy = F (x, y), y(x0 ) = y0 dx kaldes et begyndelsesværdiproblem. Eksistens- og entydighedssætningen 3.2 for begyndelsesværdiproblemer har en vigtig udvidelse til differentialligningssystemer, som det overlades til læseren at formulere. Eksempel 3.4. Differentialligningen dy = x3 y + exy dx
3. G ENEREL LIGNING
17
har løsningskurver igennem ethvert (x0 , y0 ) ∈ R × R. Løsninger kan ikke umiddelbart udtrykkes ved kendte elementære funktioner. y
1
0
x
1
Retningsfelt Som en anvendelse af Eksistens- og entydighedssætningen 3.2 kan de elementære funktioner genfindes som løsninger til simple differentialligninger. Eksempel 3.5 (Elementære funktioner). Den entydigt bestemte løsning til begyndelsesværdiproblemet dy = y, y(0) = 1 dx er eksponentialfunktionen y(x) = ex Den entydigt bestemte løsning til begyndelsesværdiproblemet dy1 dx dy2 dx
= = y1
−y2
y1 (0) = 1, y2 (0) = 0 er de trigonometriske funktioner y1 (x) = cos x y2 (x) = sin x Den entydigt bestemte løsning til begyndelsesværdiproblemet dy1 dx dy2 dx
= = y1
y2
y1 (0) = 1, y2 (0) = 0 er de hyperbolske funktioner (se Stewart, p. 251.) y1 (x) = cosh x = y2 (x) = sinh x =
ex +e−x 2 ex −e−x 2
18
4. S TABILITET
4. S TABILITET Da det normalt ikke er muligt at løse en differentialligning ved eksplicitte funktionsudtryk, er det vigtigt at kunne beskrive en løsnings egenskaber på anden vis. Her kommer begreberne ligevægt og stabilitet til deres ret. I sådanne punkter er en tilnærmelse med en lineær differentialligning ofte meningsfuld. Den følgende opremsning er ultra kort og bør opfattes som en smagsprøve. Eksemplerne refererer til [Stewart]. Definition 4.1. En differentialligning dy = F (y) dx kaldes autonom. En konstant løsning y(x) = b, F (b) = 0 kaldes en ligevægt. En ligevægt kaldes (lokal) stabil, hvis enhver løsning y(x) som kommer tilstrækkelig tæt på b, vil konvergere mod b for x gående mod uendelig. I modsat fald kaldes ligevægten ustabil. Grafen for funktionen F (y) kaldes fasediagrammet. Der er en oplagt udvidelse til differentialligningsystemer, som det overlades læseren at formulere. Bemærkning 4.2. I en omegn af en ligevægt y(x) = b, F (b) = 0 kan det autonome begyndelsesværdiproblem dy = F (y), y(x0 ) = b + dx tilnærmes med den lineære ligning dz = F 0 (b)z, z(x0 ) = dx hvor y(x) = b + z(x). Bemærkning 4.3. For en ligevægt y(x) = b, F (b) = 0 for det autonome system dy = F (y) dx gælder (forudsat F (y) er tilstrækkelig ’pæn’) F 0 (b) < 0: Stabil ligevægt. F 0 (b) > 0: Ustabil ligevægt. F 0 (b) = 0: Ingen konklussion.
4. S TABILITET
19
y
x
Fasediagram Eksempel 4.4 (Logistisk ligning). Den logistiske ligning dP P = kP (1 − ) = F (P ) dt K har ligevægts løsninger P (t) = 0, P (t) = K og 2k P +k K F 0 (0) = k > 0: P = 0 er en ustabil ligevœgt. F 0 (K) = −k: P = K er en Stabil ligevœgt. F 0 (P ) = −
y
x
Fasediagram
Eksempel 4.5 (Lotka-Volterra). For Lotka-Volterra systemet, a, b, k, r > 0, dR = kR − aRW dt dW = −rW + bRW dt er der to ligevægtsløsninger (R, W ) = (0, 0), (R, W ) = (r/b, k/a)
20
4. S TABILITET
I (R, W ) = (0, 0) er den lineære approximation dR = kR dt dW = −rW dt som giver en ustabil ligevægt. ¯ W ¯ ) = (R−r/b, W − I (R, W ) = (r/b, k/a) er den lineære approximation for (R, k/a) ¯ dR ar ¯ =− W dt b ¯ dW bk ¯ = R dt a som ifølge Definition 4.1 giver en ustabil ligevægt. Man kan vise, at løsningskurverne t 7→ (R(t), W (t) for det oprindelig system er deformationer af cirkler omkring ligevægtspunktet. Der er altså en cyklisk udvikling i modellen. W
100
0
1000
Hastighedsfelt
R
L ITTERATUR
21
L ITTERATUR Arnold, Ordinary differential equations, MIT press 1980. (Frisk velskrevet og moderne tilgang). Birkhoff & Rota, Ordinary differential equations, Wiley and Sons 1978. (Fine figurer, meget populær). Braun, Differential equations and their applications, Springer Verlag 1975. (Udførlige eksempler og gennemregninger af konkrete modeller). Coddington & Levinson, Theory of ordinary differential equations, McGraw-Hill 1955. (Omfattende, men alligevel tilgængelig). Edwards & Penney, Differential equations, computing and modeling, Prentice Hall 1996. (Elementær, supplerer forfatternes ”calculus” bog). Hartmann, Ordinary differential equations, Wiley & Sons 1964. (Den mest omfattende nyere bog om emnet). Hirsh & Smale, Differential equations, dynamical systems and linear algebra, Academic Press 1974. (Bogen om lineære systemer). Ince, Ordinary differential equations, Dover Publ. 1917. (Klassisk bog med mange ikke-lineære typer). Iversen, Reelle funktioner af flere variable, Aarhus Universitetsforlag 1988. (Har et godt afsnit om lineære systemer med udførlige beviser). Simmons, Differential equations with applications and historical notes, McGraw-Hill 1972. (Fornøjelig læsning med mange oplysninger). Stewart, Calculus - concepts and contexts, Brooks/Cole 2001. (Lærebogen, lidt tyndt om differentialligninger og modelering).