Linjära ordinära differentialekvationer med konstanta koefficienter
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
redigeraDifferentialekvationen 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
redigeraPartikulärlösningarna beror på hur högerledet i differentialekvationen ser ut, det vill säga vad är.
Högerledet är ett polynom
redigeraAntag 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.
Exempel
redigeraAntag 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
redigeraOm 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.
Exempel 1
redigeraAntag 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
Exempel 2
redigeraBetrakta 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
redigeraOm 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
redigeraAntag 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.
Exempel
redigeraLå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:
- då
- då
Högerledet är en godtycklig funktion
redigeraDenna 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
redigeraEtt 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.