Differentiaalvergelijking

Een differentiaalvergelijking (afk.: DV) is een wiskundige vergelijking voor een functie waarin, naast eventueel de functie zelf, een of meer van de afgeleiden van die functie voorkomen. Betreft het een functie van meer dan één onafhankelijke veranderlijke, dan zijn het de partiële afgeleiden die in de vergelijking voorkomen en spreekt men van een partiële differentiaalvergelijking. Is er slechts één onafhankelijke veranderlijke, dan spreekt men van een gewone differentiaalvergelijking.

Voorbeeld

Een eenvoudig voorbeeld is een vat gevuld met een vloeistof dat langzaam leegstroomt door een opening onder in het vat ter grootte $a$. De hoogte van het vloeistofniveau in het vat $h$ neemt af met de tijd. Het is dus een functie $h(t)$ van de tijd $t$. Voor een vat met constante doorsnede $D$ is de uitgestroomde hoeveelheid vloeistof in het tijdsinterval $(t,t+\mathrm {d} t)$ gelijk aan:

${\big (}h(t)-h(t+\mathrm {d} t){\big )}D=a\cdot v(t)\mathrm {d} t$,

zodat

$v(t)=-{\frac {D}{a}}h'(t)$

De uitstroomsnelheid van de vloeistof hangt af van de hoogte. In een eenvoudig geval is die snelheid evenredig met de hoogte en is de hoogteverandering $h'(t)$ dus evenredig met de hoogte. Dat levert de differentiaalvergelijking:

$h'(t)=-c\cdot h(t)$

met $c$ een positieve constante.

Omdat het eigenlijk een vergelijking is voor de functie $h$, wordt vaak geschreven:

$h'=-c\cdot h$

De algemene oplossing is:

$h(t)=A\,e^{-ct}$,

waarin de constante bepaald wordt door de voorwaarde $A=h(0)$.

Stel namelijk dat ook ${\tilde {h}}$ een oplossing is, dan volgt:

$\left({\frac {\tilde {h}}{h}}\right)'={\frac {{\tilde {h}}'h-{\tilde {h}}h'}{h^{2}}}={\frac {-c{\tilde {h}}h+c{\tilde {h}}h}{h^{2}}}=0$

Dus is ${\tilde {h}}$ een veelvoud van $h$.

Toepassingen

Allerlei verschijnselen in de natuurkunde en de toepassingen daarvan in de techniek worden door differentiaalvergelijkingen beschreven. Het voorbeeld bij uitstek vormen trillingen en golven. Bij de eenvoudigste trillingsvergelijking is er een evenredig verband tussen de tweede afgeleide en de functie zelf. De oplossing is een sinusfunctie. Golfvoortplanting in de ruimte wordt door een partiële differentiaalvergelijking beschreven met als veranderlijken de drie ruimtelijke coördinaten en de tijd.

Ook bevolkingsgroei kan door een differentiaalvergelijking worden beschreven. Gaat men uit van de (sterk vereenvoudigde) veronderstelling van een constante vruchtbaarheid, dat wil zeggen dat de bevolking van een land groeit met een snelheid die evenredig is met het aantal inwoners, dan is de eerste afgeleide van het aantal inwoners als functie van de tijd, evenredig met dat aantal zelf. De oplossing van deze differentiaalvergelijking is een exponentiële functie; vandaar de welhaast spreekwoordelijke exponentiële groei.

Differentiaalvergelijkingen die in de natuur voorkomende verschijnselen beschrijven, zijn vaak niet wiskundig oplosbaar. Soms kan dan met numerieke methoden een benaderende oplossing gevonden worden. Maar dan is er alsnog een hoop wiskunde nodig voordat zo'n vergelijking kan worden 'opgelost' met een computer die eigenlijk slechts elementaire rekenkundige bewerkingen kan uitvoeren. Sommige differentiaalvergelijkingen kunnen zelfs met computers niet of niet nauwkeurig worden opgelost, bijvoorbeeld vergelijkingen die turbulente stromingen beschrijven. Dat is een van de redenen dat er geen betrouwbare weersvoorspellingen mogelijk zijn op langere termijn.

Definitie

Een gewone differentiaalvergelijking is een vergelijking waarin als onbekende een functie $f$ van één veranderlijke voorkomt en een of meer van zijn afgeleiden, met steeds hetzelfde argument.

De algemene impliciete vorm van een gewone differentiaalvergelijking is:

$F(x,f(x),f^{\prime }(x),f^{\prime \prime }(x),\ldots ,f^{(n)}(x))=0$

Daarnaast zijn er varianten zoals $f^{\prime }(x)-f(x-1)=0$ (een delay differential equation), die buiten deze definitie vallen.

De orde $n$ van de hoogste voorkomende afgeleide $f^{(n)}$ heet de orde van de differentiaalvergelijking.

Als de vergelijking opgelost is naar de hoogste afgeleide, dus geschreven kan worden in de vorm:

$f^{(n)}(x)=\varphi (x,f(x),f^{\prime }(x),f^{\prime \prime }(x),\ldots ,f^{(n-1)}(x))$

noemt men deze de differentiaalvergelijking in expliciete vorm.

Oplossing

Er zijn maar weinig differentiaalvergelijkingen met een "gesloten" oplossing, een oplossing in een formule. Dat stelt een grens aan de mogelijkheden van een wiskundige benadering van natuurverschijnselen. Nog de meeste mogelijkheden van een puur wiskundige oplossing hebben lineaire differentiaalvergelijkingen: vergelijkingen waarbij de som van twee oplossingen ook een oplossing is. Trillingen en golven worden vaak door zulke differentiaalvergelijkingen beheerst. Transformaties vanuit het domein van tijd en ruimte naar frequenties (waaraan de naam van Fourier is verbonden) kunnen het differentiëren tot een vermenigvuldiging reduceren.

Veel verschijnselen in de natuur dragen echter een niet-lineair karakter. Dat geldt met name voor turbulente stromingen. Natuurlijk zijn er benaderingsmethoden om aan zulke verschijnselen te rekenen, maar de mogelijkheden zijn toch zo beperkt dat experimenten hier vooralsnog onontbeerlijk blijven. Het zal de leek verbazen dat ook met moderne supercomputers nog steeds niet alles uit te rekenen is.

Lineaire differentiaalvergelijking

Een lineaire $n$-de-orde-differentiaalvergelijking is een gewone differentiaalvergelijking van de vorm:

$c_{n}(x)f^{(n)}(x)+c_{n-1}(x)f^{(n-1)}(x)+\ldots +c_{1}(x)f'(x)+c_{0}(x)f(x)=g(x)$

Als $g(x)\equiv 0$, heet de differentiaalvergelijking homogeen of gereduceerd, anders inhomogeen of compleet. De verzameling functies die oplossing zijn van een bepaalde homogene lineaire differentiaalvergelijking vormt een vectorruimte (elke lineaire combinatie van oplossingen is ook een oplossing).

De algemene oplossing van een inhomogene vergelijking is te schrijven als een particuliere oplossing plus de algemene oplossing van de bijbehorende homogene vergelijking. De verzameling functies die oplossing zijn van een bepaalde inhomogene lineaire differentiaalvergelijking vormt een affiene ruimte.

Constante coëfficiënten

Een speciaal geval is een lineaire differentiaalvergelijking met constante coëfficiënten, waarin de coëfficiënten $c_{i}(x)$ constanten $c_{i}$ zijn. Voor dergelijke differentiaalvergelijkingen bestaat een relatief simpele algemene oplossingsmethode.

Algemene oplossingsmethode

Een oplossing van de homogene vergelijking:

$c_{n}f^{(n)}+c_{n-1}f^{(n-1)}+\ldots +c_{1}f'+c_{0}f=0$

is de functie

$f(x)=e^{ax}$,

mits de parameter $a$ voldoet aan de zogenaamde karakteristieke vergelijking:

$c_{n}a^{n}+c_{n-1}a^{n-1}+\ldots +c_{1}a+c_{0}=0$,

die door invullen van de mogelijke oplossing in de differentiaalvergelijking ontstaat. Dit is een gewone polynomiale vergelijking van de graad $n$ in de parameter $a$. In het algemeen heeft deze vergelijking $n$ complexe oplossingen $a_{1},\ldots ,a_{n}$, waarvan er eventueel enkele kunnen samenvallen. Als alle oplossingen verschillend zijn, wordt de algemene oplossing van de differentiaalvergelijking gegeven door een lineaire combinatie van de afzonderlijke e-machten:

$f(x)=A_{1}e^{a_{1}x}+A_{2}e^{a_{2}x}+\ldots +A_{n}e^{a_{n}x}$,

waarin de coëfficiënten $(A_{i})$ nog vrij gekozen kunnen worden. Meestal worden deze coëfficiënten vastgelegd door de beginvoorwaarden.

Voorbeeld 1

De differentiaalvergelijking:

$f''+f=0$

kan gelezen worden als: zoek een functie die als tweede afgeleide zijn eigen tegengestelde heeft. Het is bekend dat de sinus en de cosinus deze eigenschap hebben.

De karakteristieke vergelijking is:

$a^{2}+1=0$,

met de twee oplossingen: $a_{1}=i$ en $a_{2}=-i$.

De algemene oplossing van de differentiaalvergelijking wordt dus:

$f(x)=A_{1}e^{ix}+A_{2}e^{-ix}$

Door geschikte keuze van $A_{1}$ en $A_{2}$ komen inderdaad de sinus en de cosinus als oplossing tevoorschijn.

Voorbeeld 2

Het tweede voorbeeld is een inhomogene lineaire tweede-orde-differentiaalvergelijking:

$4f^{\prime \prime }(x)+f(x)=\sin(x)$

De bijbehorende homogene DV heeft de karakteristieke vergelijking:

$4a^{2}+1=0$,

met de oplossingen: $a_{1}={\tfrac {1}{2}}i$ en $a_{2}=-{\tfrac {1}{2}}i$

Dat geeft als algemene oplossing van de homogene DV:

$f_{H}(x)=A_{1}e^{{\frac {1}{2}}ix}+A_{2}e^{-{\frac {1}{2}}ix}$

Deze kan ook geschreven worden als

$f_{H}(x)=A\sin({\tfrac {1}{2}}x)+B\cos({\tfrac {1}{2}}x)$

Vervolgens gaat het erom een particuliere oplossing van de inhomogene DV te vinden. De twee meest gebruikte methoden om een dergelijke oplossing te vinden zijn de methode van de onbepaalde coëfficiënten en de variatie van de parameters.

Een eenvoudige vorm van de eerste van deze twee methoden is te proberen of

$f_{P}(x)=a\sin(x)$

voldoet als mogelijke oplossing. Invullen in de DV levert:

$-4a\sin(x)+a\sin(x)=\sin(x)$,

waaruit volgt $a=-{\tfrac {1}{3}}$.

Een particuliere oplossing is dus:

$f_{P}(x)=-{\tfrac {1}{3}}\sin(x)$

De algemene oplossing is de som van de gevonden particuliere oplossing $f_{P}(x)$ en de algemene oplossing $f_{H}(x)$ van de homogene DV :

$f(x)=f_{P}(x)+f_{H}(x)=-{\tfrac {1}{3}}\sin(x)+A\sin({\tfrac {1}{2}}x)+B\cos({\tfrac {1}{2}}x)$

Niet-lineaire differentiaalvergelijkingen

Het oplossen van niet-lineaire differentiaalvergelijkingen is een stuk moeilijker en onoverzichtelijker. Er is geen algemene oplossingsmethode voor.

Voorbeelden

De bewegingsvergelijking van een slinger wordt gegeven door de volgende differentiaalvergelijking:

$\varphi ''(t)+\omega \,\sin {\big (}\varphi (t){\big )}=0$

Daarin is $\varphi$ de hoek die de slinger maakt met de verticale richting en $\omega$ een constante. Deze vergelijking is niet oplosbaar met standaardtechnieken. Voor kleine uitwijkingen kan een benadering toegepast worden; dan is $\sin(\varphi )\approx \varphi$, zodat de vergelijking overgaat in die van de harmonische oscillator:

$\varphi ''+\omega \varphi =0$

Een ander voorbeeld is de differentiaalvergelijking voor de vorm van een hangende ketting of touw:

$f(t)f''(t)=(f'(t))^{2}+1$

Deze heeft de zogenaamde kettinglijn als algemene oplossing

$f(t)=C\,\cosh {\frac {t-t_{0}}{C}}$,

met de integratieconstanten $C$ en $t_{0}$.

Begin- en randvoorwaarden

Om een eenduidige oplossing van een differentiaalvergelijking te krijgen, moeten randvoorwaarden opgelegd worden. In het algemeen kan gesteld worden dat voor een $n$-de orde differentiaalvergelijking $n$ verschillende randvoorwaarden nodig zijn.

Bijvoorbeeld: de eerste orde differentiaalvergelijking

$f'(t)=f(t)$

heeft als algemene oplossing $f(t)=Ae^{t}$, waarbij $A$ nog onbepaald is. Door de beginvoorwaarde $f(0)=1$ op te leggen, wordt de oplossing eenduidig bepaald als $f(t)=e^{t}$.

Lineaire vergelijkingen

Men kan bewijzen, dat een lineaire differentiaalvergelijking van $n$-de orde, met randvoorwaarden

$y(x_{0})=y_{0}$

$y'(x_{0})=y_{1}$

...

$y^{(n-1)}(x_{0})=y_{n-1}$

één unieke continu oplossing heeft.

Een mooi voorbeeld om de natuurkundige betekenis van begin- en randvoorwaarden te illustreren is de trillende snaar. De differentiaalvergelijking die een trillende snaar beschrijft heeft oneindig veel oplossingen (waaronder de nuloplossing voor een snaar in rust). Pas met de randvoorwaarden dat de snaar een amplitude nul heeft in haar bevestigingspunten, en de beginvoorwaarde van de stand van de snaar op een bepaald moment is een eenduidige oplossing te berekenen.

Een partiële differentiaalvergelijking

Zie Partiële differentiaalvergelijking voor het hoofdartikel over dit onderwerp.

Veel verschijnselen in de fysica moeten beschreven worden met behulp van partiële differentiaalvergelijkingen. Bijvoorbeeld de trilling van een snaar wordt beschreven door

${\frac {\partial ^{2}u}{\partial x^{2}}}={\frac {\partial ^{2}u}{\partial t^{2}}}$,

waarbij $u(x,t)$ de uitwijking van de snaar is, $x$ de positie van de snaar (van 0 aan het ene eind tot 1 aan het andere eind) en $t$ de tijd. Deze partiële differentiaalvergelijking is van de tweede orde, zowel in de tijd als in de plaats. Er zijn dus twee randvoorwaarden in de tijd ("beginvoorwaarde") nodig, bijvoorbeeld

$u(x,0)=f(x),\ {\frac {\partial u}{\partial t}}(x,0)=g(x)$

en twee randvoorwaarden in de plaats, bijvoorbeeld

$u(0,t)=0$,

$u(1,t)=0$.

Deze randvoorwaarden houden in dat de snaar ingeklemd is. Door middel van scheiden van veranderlijken is deze DV te herleiden tot twee gewone DV's, waarin alleen tijd of plaats voorkomt als onafhankelijk veranderlijke.

Differentievergelijkingen

Het discrete analogon van een differentiaalvergelijking is een differentievergelijking.

Software

Software die differentiaalvergelijkingen kan oplossen:

• ExpressionsInBar
• Maple
• SageMath
• Xcas

Zie ook

• Lineaire differentiaalvergelijking van eerste orde
• Partiële differentiaalvergelijking
• Scheiden van veranderlijken
• Variatie van de parameters: methode voor het vinden van een particuliere oplossing.
• Methode van de onbepaalde coëfficiënten: methode voor het vinden van een particuliere oplossing.
• Orthogonale families

Externe bronnen

• (en) videocolleges over differentiaalvergelijkingen op MIT OpenCourseWare door professor Arthur Mattuck
• Oplossingsmethode voor enkele soorten differentiaalvergelijkingen.

Referenties

Mediabestanden

Zie de categorie Differential equations van Wikimedia Commons voor mediabestanden over dit onderwerp.