Hoppa till innehållet

Linjära ordinära differentialekvationer med konstanta koefficienter

Från Wikipedia

Linjära ordinära differentialekvationer med konstanta koefficienter är av formen

Den allmänna lösningen till ekvationen är

lösning = partikulärlösning + homogen lösning.

Homogena ekvationer

[redigera | redigera wikitext]

Differentialekvationen ovan sägs vara homogen när högerledet är 0. För att få den homogena lösningen till en ekvation vars högerled inte är 0, sätter man högerledet till 0.

Den första lösningsmetoden för ordinära differentialekvationer med konstanta koefficienter gavs av Euler.

Som ett exempel kan vi ta

,

där är reella konstanter. Betrakta den algebraiska ekvationen

(byt ut mot , för alla k från 0 till n).
Om nu denna ekvation (kallad den karakteristiska ekvationen) har lösningarna med multiplicitet : respektive, så ges lösningen till differentialekvationen av

,

där är ett polynom av grad i x.

Det exakta utseendet av polynomen kan bestämmas utifrån begynnelsevärden eller randvillkor.

Om lösningarna till den karakterisitiska ekvationen inte är reella, så kommer lösningen innehålla termer på formen

vilket kan överföras på trigonometrisk form.

Partikulärlösningar

[redigera | redigera wikitext]

Partikulärlösningarna beror på hur högerledet i differentialekvationen ser ut, det vill säga vad är.

Högerledet är ett polynom

[redigera | redigera wikitext]

Antag att högerledet är ett polynom av grad k. Knepet för att hitta partikulärlösningen går ut på att göra en ansats, det vill säga att man låter , deriverar så många gånger som behövs, sätter in i ekvationen, och konstaterar att man får ett polynom på var sida likhetstecknet; och att två polynom är lika om alla koefficienter stämmer överens.

Antag att ekvationen är

.

Den homogena ekvationen har lösningarna , där A och B är två konstanter. För partikulärlösningen skriver vi

Nu är a, b c och d konstanter som vi försöker bestämma. Andraderivatan, , blir då

.

Insättning i ekvationen ger

Om polynomen till höger och till vänster ska vara lika, så måste:

Detta linjära ekvationssystem har lösningen a=1/4, b=0, c=-7/8 och d=1/4. Alltså är partikulärlösningen

och totala lösningen

Konstanterna A och B kan bestämmas ur randvillkor.

Högerledet är en exponentialfunktion

[redigera | redigera wikitext]

Om högerledet är en exponentialfunktion så måste man titta lite grand på den homogena lösningen innan man kan göra ansättningen. Om högerledet ges av , och k är en lösning till den karakteristiska ekvationen, med multiplicitet p, så är det lämpligt att ansätta

och sedan bestämma a. Om å andra sidan k inte är en lösning till den karakteristiska ekvationen, räcker det att ansätta yp(x)=bekx.

Antag att vi vill lösa

Eftersom den homogena lösningen är

så kan vi sätta yp(x)=be3x. Derivatan y′p(x)=3be3x, så om vi sätter in i ekvationen får vi:

Då ser vi att 4b=2, eller b=1/2 och att allmänna lösningen är

Betrakta ekvationen

De homogena lösningarna är på formen

Den karakteristiska ekvationen har en dubbelrot, varvid partikulärlösningen bör ansättas som

Efter deriveringar:

Insättning i ekvationen ger då att

För att högerled och vänsterled ska stämma överens, måste alltså a=1/2, och totala lösningen till ekvationen blir

där A och B bestäms ur randvillkor

Högerledet är en trigonometrisk funktion

[redigera | redigera wikitext]

Om högerledet är en trigonometrisk funktion med sinus- och/eller cosinus-termer antar man alltid att

och löser som i tidigare exempel.
Observera att den antagna partikulärlösningen alltid ska ha både cosinus- och sinusterm, även om den ena termen saknas i differentialekvationen.

Högerledet är en diskontinuerlig funktion

[redigera | redigera wikitext]

Antag att högerledet f(x) är kontinuerligt i intervallen (a,b) och (b,c), samt att problemet är ett begynnelsevärdesproblem. Då löses ekvationen som vanligt i intervallet (a,b). (Där beter sig ekvationen precis som vanligt). Därefter beräknas lösningens värde, liksom så många derivator som behövs, då x=b. Dessa värden får sedan fungera som begynnelsevärden för den nya ekvationen som definieras på intervallet (b,c).

Om problemet är ett randvärdesproblem, ansätts konstanter för de begynnelsevärden som är okända. I princip räknar man sedan på samma sätt som ovan. Då kommer konstanterna följa med genom samtliga räkningar, och man kan avsluta med att jämföra värdet i den givna punkten på den lösning som räknats fram, med det värde som specificerats av randvillkoret.

Det är ofta enklare att använda Laplace-transformen för differentialekvationer där högerledet är en diskontinuerlig funktion.

Låt y′′(x) + 4π2 y(x) = f(x), där f(x)= 2-x då 0 ≤ x <1, och f(x)= x-2 då 1 ≤ x ≤ 2.

Antag även att vi har randvillkor y(0) = y(2) = 0

Den homogen lösningen i första intervallet (d.v.s. (0,1)) ges av

yh(x)= A sin(2πx) + B cos(2πx).

För att få fram partikulärlösningen i det första intervallet ser vi på att den bör se ut som yp = (2-x) / (4π2), en kontroll visar att det stämmer. Alltså har vi totala lösningen

y(x) = A sin (2πx) + Bcos(2πx) + (2-x) / (4π2).

Använder vi det begynnelsevillkor vi har, så får vi att

B + 2 / (4π2) = 0.

Därmed är

y(x) = A sin(2πx) - 2 / (4π2) cos(2πx) + (2-x) / (4π2),

så länge som x tillhör intervallet [0,1].

Nu gäller att få fram så mycket begynnelsevärden som möjligt till den andra delen av ekvationen: Insättning i den förra lösningen ger att

y(1) = -1 / (4π2)

(i just detta fall oberoende av A), och

y′(1) = 2πA - 1 / (4π2).

Den homogena lösningen i det andra intervallet skriver vi som

yh = C sin(2πx) + D cos(2πx),

och partikulärlösningen är trevligt nog

yp = (x - 2) / (4π2).

Totala lösningen alltså

y(x) = C sin(2πx) + D cos(2πx) + (x - 2) / (4π2).

Jämför med det begynnelsevärde vi känner, så fås y(1) = D - 1 / (4π2) = - 1 / (4π2), och y′(1) = 2πC + 1 / (4π2) = 2πA - 1 / (4π2).

Därmed är D = 0, och C = A - 1 / (4π3). Insättning i "dellösning 2" ger att

y(x) = (A - 1 / (4π3))sin(2πx) + (x - 2) / (4π2),

och det enda som borde återstå är att använda villkoret y(2) = 0 för att bestämma A. Dock ser vi att i just detta exempel så är y(2) = 0 oberoende av A. Kort sagt: vi har stött på ett fenomen som förekommer hos randvärdesproblem: det är inte säkert att en differentialekvation med randvillkor har en unik (entydig) lösning, det är till och med så att det inte ens är säkert att den har någon lösning över huvud taget. Betänk till exempel det ovanstående fallet där vi valt y(2) = 1 som ett av randvillkoren - den ekvationen är uppenbarligen inte möjlig att lösa.

För att sammanfatta exemplet:

Högerledet är en godtycklig funktion

[redigera | redigera wikitext]

Denna metod går ut på att ersätta konstanterna i den homogena lösningen med funktioner. Här visas metoden för differentialekvationer av andra ordningen. Då är den homogena lösningen av formen

Man antar partikulärlösningen

(I fortsättningen utlämnas (x) och högerledets funktion betecknas med H.)

Lösningen är

Integrering ger funktionerna och och integrationskonstanterna sätts till noll, eftersom endast en partikulärlösning behövs. Slutligen sätts och in i och lösningen är färdig.

System av ekvationer

[redigera | redigera wikitext]

Ett system av linjära ODEer med konstanta koefficienter kan skrivas

där är en vektorvärd funktion, en vektorvärd funktion och en konstant matris.

Den allmänna lösningen till detta system är summan av en partikulärlösning och den allmänna homogena lösningen, precis som ovan. Den allmänna homogena lösningen kan skrivas där är matrisexponentialen och en godtycklig konstant vektor.

Eulers formel

[redigera | redigera wikitext]
Huvudartikel: Eulers formel