Kurs:Analysis II/Kapitel VI: Gewöhnliche Differentialgleichungen/Lineare Differentialgleichungen m-ter Ordnung (§9)

Auf dem kompakten Intervall $I:=\{x\in ℝ:|x-\xi |\le a\}$ um den Punkt $\xi \in ℝ$ von der halben Länge $a\in (0,+\mathrm{\infty })$ wählen wir die reellwertigen, stetigen Koeffizientenfunktionen

$p_k(x)\in C^0(I,ℝ)$ für $k=0,1,\dots ,m$.

Dabei setzen wir $p_m(x)\ne 0$ für alle $x\in I$ voraus. Dann betrachten wir den reellen Vektorraum $C^m(I):=C^m(I,ℝ)$ mit den Verknüpfungen

$c\in ℝ,u,v\in C^m(I)\Rightarrow cu,u+v\in C^m(I)$.

Wir erklären nun den linearen Differentialoperator $m$-ter Ordnung

(1) $\begin{array}{c}ℒ:C^m(I)\to C^0(I)\mathrm{v}\mathrm{e}\mathrm{r}\mathrm{m}\ddot{\mathrm{o}}\mathrm{g}\mathrm{e}\\ ℒ(y){|}_x:=p_m(x)y^{(m)}(x)+p_{m-1}(x)y^{(m-1)}(x)+\dots +p_0(x)y(x),x\in I.\end{array}$

Offenbar gilt die Linearitätsregel

(2) $ℒ(c\cdot u+d\cdot v)=c\cdot ℒ(u)+d\cdot ℒ(v)$ für alle $c,d\in ℝ$ und $u,v\in C^m(I)$.

Zu einer gegebenen rechten Seite $f(x)\in C^0(I)$ wollen wir nun alle Lösungen $y(x)\in C^m(I)$ von $ℒ(y)=f$ bestimmen. Zunächst berechnen wir alle Lösungen der homogenen Gleichung $ℒ(y)=0$. Diese bilden – gemäß unserer nachfolgenden Untersuchungen – einen $m$-dimensionalen Vektorraum

$𝒰:=\{y\in C^m(I):ℒ(y)=0\}$.

Wir bestimmen dann eine Lösung $y_0(x)\in C^m(I)$ der inhomogenen Gleichung $ℒ(y)=f$ mittels Variation der Konstanten. Wir erhalten schließlich die Lösungsgesamtheit der inhomogenen Gleichung in der Form $y_0(x)+y(x)$ mit $y(x)\in 𝒰$.

Satz 1

Sei der Vektor $\eta =({\eta }_1,\dots ,{\eta }_m{)}^{*}\in {ℝ}^m$gewählt, so sind die folgenden Aussagen äquivalent:

1. Die Funktion $y(x)\in C^m(I)$ genügt dem Anfangswertproblem

$ℒ(y)=f(x),x\in I,y^{(k-1)}(\xi )={\eta }_k$ für $k=1,\dots ,m$.

2. Die Funktion $\hat{y}(x)={(y(x),y^{\prime }(x),\dots ,y^{m-1}(x))}^{*}\in C^1(I,{ℝ}^m)$ genügt dem System

${\hat{y}}^{\prime }(x)=P(x)\circ \hat{y}(x)+q(x),\hat{y}(\xi )=\eta$

mit der Matrixfunktion

(3) $P(x)=(p_{jk}(x){)}_{j,k=1,\dots ,m}=\left(\begin{array}{ccccc}0& 1& 0& \dots & 0\\ ⋮& \ddots & \ddots & \ddots & ⋮\\ ⋮& & \ddots & \ddots & 0\\ 0& \dots & \dots & 0& 1\\ \frac{-p_0(x)}{p_m(x)}& \dots & \dots & \dots & \frac{-p_{m-1}(x)}{p_m(x)}\end{array}\right)$

und der Vektorfunktion

(4) $q(x):={(0,\dots ,0,\frac{f(x)}{p_m(x)})}^{*},x\in I$.

Beweis

Wir beachten zunächst

(5) $\begin{array}{c}ℒ(y)=f(x)\iff p_m(x)y^{(m)}(x)+\dots +p_0(x)y(x)=f(x)\iff \\ y^{(m)}(x)=\frac{f(x)}{p_m(x)}-\frac{p_0(x)}{p_m(x)}y(x)-\dots -\frac{p_{m-1}(x)}{p_m(x)}{y}^{(m-1)}(x)=:f(x,y,\dots ,y^{(m-1)}).\end{array}$

Mit den Überlegungen zu Beginn von §8 erhalten wir für $\hat{y}(x)$ das folgende System

(6) ${\left(\begin{array}{c}y(x)\\ y^{\prime }(x)\\ ⋮\\ y^{(m-1)}\end{array}\right)}^{\prime }=\left(\begin{array}{ccccc}0& 1& 0& \dots & 0\\ ⋮& \ddots & \ddots & \ddots & ⋮\\ ⋮& & \ddots & \ddots & 0\\ 0& \dots & \dots & 0& 1\\ \frac{-p_0}{p_m}& \dots & \dots & \dots & \frac{-p_{m-1}}{p_m}\end{array}\right)\circ \left(\begin{array}{c}y(x)\\ y^{\prime }(x)\\ ⋮\\ ⋮\\ y^{(m-1)}\end{array}\right)+\left(\begin{array}{c}0\\ ⋮\\ ⋮\\ 0\\ \frac{f(x)}{p_m}\end{array}\right)$.

Hieraus ersehen wir sofort die Behauptung des Satzes.

q.e.d.

Definition 1

Ein System von $m$ Lösungen $\{y_1,\dots ,y_m\}$ der Differentialgleichung $ℒ(y)=0$ heißt Fundamentalsystem von $ℒ(y)=0$, wenn die Funktionen $y_1(x),\dots ,y_m(x)$ im Intervall $I$ linear unabhängig sind. Letzteres bedeutet, dass aus der Identität

$c_1{y}_1(x)+\dots +c_m{y}_m(x)\equiv 0$ in $I$ mit Konstanten $c_1,\dots ,c_m\in ℝ$

die Aussage $c_1=\dots =c_m=0$ folgt.

Definition 2

Wir erklären die Wronskische Determinante des Systems $\{y_1,\dots ,y_m\}$ durch

(7) $W(x)=W(y_1,\dots ,y_m){|}_x:=\left|\begin{array}{ccc}{y}_1(x)& \dots & y_m(x)\\ {y_1}^{\prime }(x)& \dots & {y_m}^{\prime }(x)\\ ⋮& & ⋮\\ y_1^{(m-1)}(x)& \dots & y_m^{(m-1)}(x)\end{array}\right|,x\in I$.

Bemerkung

Ist $\{y_1,\dots ,y_m\}$ ein System von Lösungen von $ℒ(y)=0$ im Intervall $I$ mit der Wronskischen Determinante $W(x)=W(y_1,\dots ,y_m){|}_x$. Dann genügt diese der Differentialgleichung

(8) $\frac{d}{dx}W(x)=\left|\begin{array}{ccc}{y}_1(x)& \dots & y_m(x)\\ {y_1}^{\prime }(x)& \dots & {y_m}^{\prime }(x)\\ ⋮& & ⋮\\ y_1^{(m-2)}(x)& \dots & y_m^{(m-2)}(x)\\ y_1^{(m)}(x)& \dots & y_m^{(m)}(x)\end{array}\right|=-\frac{p_{m-1}(x)}{p_m(x)}\cdot W(x),x\in I$.

Somit ist $W(x)\equiv 0$ in $I$ genau dann erfüllt, wenn in einem Punkt $x_0\in I$ die Aussage $W(x_0)=0$ richtig ist.

Satz 2 (Wronskische Determinante)

Die Funktionen $y_1,\dots ,y_m$ der Klasse $C^m(I)$ seien Lösungen der Differentialgleichung $ℒ(y)=0$ in $I$ und sei $x_0\in I$ beliebig gewählt.

Dann ist $\{y_1,\dots ,y_m\}$ ein Fundamentalsystem von \mathcal{L}(y) = 0</math> genau dann, wenn $W(y_1,\dots ,y_m){|}_{x=x_0}\ne 0$ gilt.

Beweis

„$\Rightarrow$“ Wäre $W(x_0)=0$ erfüllt, so existiert ein Vektor

$c=(c_1,\dots ,c_m{)}^{*}\in {ℝ}^m\setminus \{0\}$ mit der Eigenschaft $Y(x_0)\circ c=0$.

Dabei haben wir die Matrixfunktion

(9) $Y(x):={(y_k^{(j-1)}(x))}_{j,k=1,\dots ,m},x\in I$

erklärt. Also löst $\mathrm{\Phi }(x):=Y(x)\circ c:I\to {ℝ}^m$ das lineare Anfangswertproblem

(10) $\begin{array}{c}{\mathrm{\Phi }}^{\prime }(x)=Y^{\prime }(x)\circ c=P(x)\circ Y(x)\circ c=P(x)\circ \mathrm{\Phi }(x),x\in I\\ \mathrm{\Phi }(x_0):=Y(x_0)\circ c=0.\end{array}$

Somit liefert der Eindeutigkeitssatz für Differentialgleichungssysteme

$\mathrm{\Phi }(x)=0$ bzw. $c_1{y}_1(x)+\dots +c_m{y}_m(x)=0$ für alle $x\in I$.

Also ist $\{y_1,\dots ,y_m\}$ kein Fundamentalsystem – im Widerspruch zur Voraussetzung!

„$\Leftarrow$“ Wäre $\{y_1,\dots ,y_m\}$ kein Fundamentalsystem, dann existiert ein Vektor $c=(c_1,\dots ,c_m{)}^{*}\in {ℝ}^m\setminus \{0\}$ mit der Eigenschaft

$0=c_1{y}_1(x)+\dots +c_m{y}_m(x),x\in I$.

Eine $(m-1)$-fache Differentiation liefert für alle $x\in I$ die $m-1$ Gleichungen

$0=c_1{y_1}^{\prime }(x)+\dots +c_m{y_m}^{\prime }(x),0=c_1{y}_1^{m-1}(x)+\dots +c_m{y}_m^{m-1}(x)$.

Wir erhalten die Matrixidentität $Y(x_0)\circ c=0$ mit einem Vektor $c\in {ℝ}^m\setminus \{0\}$ und es folgt $W(x_0)=\det Y(x_0)=0$ – im Widerspruch zur Voraussetzung!

q.e.d.

Satz 3 (Fundamentalsystem)

1. Es gibt ein Fundamentalsystem $\{y_1,\dots ,y_m\}$ der homogenen, linearen Differentialgleichung $m$-ter Ordnung $ℒ(y)=0$.

2. Jede Lösung $y=y(x)\in C^m(I)$ von $ℒ(y)=0$ lässt sich in der Form

$y(x)=c_1{y}_1(x)+\dots +c_m{y}_m(x),x\in I$

mit gewissen Konstanten $c_1,\dots ,c_m\in ℝ$ darstellen. Somit folgt

$𝒰=\{y(x)=c_1{y}_1(x)+\dots +c_m{y}_m(x):c_1,\dots ,c_m\in ℝ\}$.

Beweis

1. Mit der Matrixfunktion (3) lösen wir zu den Einheitsvektoren

$\stackrel{j}{e}:=({\delta }_{1j},\dots ,{\delta }_{mj}{)}^{*}\in {ℝ}^m$ für $j=1,\dots ,m$

die folgenden Anfangswertprobleme

(11) $\stackrel{j}{y}(x)\in C^1(I,{ℝ}^m)$ mit $\stackrel{j}{y^{\prime }}(x)=P(x)\circ \stackrel{j}{y}(x),x\in I$ und $\stackrel{j}{y}(\xi )=\stackrel{j}{e}$

über den Satz 1 aus §7. Wir fassen dann diese Lösungen zu einer Matrixfunktion zusammen, welche das folgende Anfangswertproblem löst:

(12) $\begin{array}{c}Y(x):=(\stackrel{1}{y}(x),\dots ,\stackrel{m}{y}(x))\in C^1(I,{ℝ}^{m\times m})\\ \text{mit }{Y}^{\prime }(x)=P(x)\circ Y(x),x\in I\text{ und }Y(\xi )=E:=({\delta }_{ij}{)}_{i,j=1,\dots ,m}.\end{array}$

Wir definieren die ersten Komponentenfunktionen

(13) $y_j(x):={\stackrel{j}{y}}_1(x),x\in I$ für $j=1,\dots ,m$,

welche gemäß Satz 1 die homogene Differentialgleichung lösen. Weiter bilden $\{y_1,\dots ,y_m\}$ ein Fundamentalsystem von $ℒ(y)=0$ wegen Satz 2 und $W(\xi )=\det Y(\xi )=\det E=1$.

2. Wir hatten $Y(x),x\in I$ aus (12) die Fundamentallösung des homogenen Differentialgleichungssystems genannt, welches den Lösungsraum

$𝒱:=\{\tilde{y}(x)\in C^1(I,{ℝ}^m):P(x)\circ \tilde{y}(x)=0,x\in I\}$

besitzt. Ist nun $y(x)$ eine beliebige Lösung von $ℒ(y)=0$, so liegt die Funktion $\hat{y}(x):=(y(x),y^{\prime }(x),\dots ,y^{(m-1)}(x){)}^{*}$ im Lösungsraum $𝒱$. Also gibt es nach Satz 2 aus §7 einen Vektor $(c_1,\dots ,c_m)\in {ℝ}^m$, so dass

(14) $\hat{y}(x)=c_1\stackrel{1}{y}(x)+\dots +c_m\stackrel{m}{y}(x),x\in I$

und folglich

(15) $y(x)=c_1{y}_1(x)+\dots +c_m{y}_m(x),x\in I$

richtig ist.

q.e.d.

Satz 4 (Variation der Konstanten)

Sei $\{y_1,\dots ,y_m\}$ ein Fundamentalsystem der Differentialgleichung $ℒ(y)=0$. Dann lässt sich die allgemeine Lösung der Differentialgleichung $ℒ(y)=f(x)$ in der Form

(16) $y(x)=\sum _{k=1}^m{c}_k{y}_k(x)+\sum _{k=1}^m[y_k(x)\left(\underset{\xi }{\overset{x}{\int }}\left\{\frac{z_k(t)}{W(t)}\frac{f(t)}{p_m(t)}\right\}dx\right)],x\in I$

mit Konstanten $c_1,\dots ,c_m\in ℝ$ darstellen. Dabei ist der Faktor

(17) $z_k(t):=(-1{)}^{m+k}W(y_1,\dots ,y_{k-1},y_{k+1},\dots ,y_m){|}_t$ für $k=1,\dots ,m$

durch die reduzierte Wronski-Determinante erklärt. Der erste Summand in (16) stellt die allgemeine Lösung der homogenen Gleichung dar und der zweite gibt eine partikuläre Lösung der inhomogenen Gleichung an.

Beweis

Aus dem Fundamentalsystem bilden wir die Fundamentallösung

(18) $Y(x):={\left(y_k^{(j-1)}\right)}_{j,k=1,\dots ,m}$ mit $Y^{\prime }(x)=P(x)\circ Y(x),x\in I$.

Gemäß Satz 1 gilt $ℒ(y)=f(x)$ genau dann, wenn die Funktion $\mathrm{\Phi }(x)={(y(x),y^{\prime }(x),\dots ,y^{(m-1)}(x))}^{*}$ das System ${\mathrm{\Phi }}^{\prime }(x)=P(x)\circ \mathrm{\Phi }(x)+q(x),x\in I$ löst. Letzteres können wir vollständig mittels Satz 3 aus §7 über inhomogene Differentialgleichungssysteme lösen. Wir ermitteln jedoch direkt eine Lösung der inhomogenen Gleichung durch den Ansatz der Variation der Konstanten

(19) $\mathrm{\Phi }(x)=Y(x)\circ c(x)$ mit $c(x)=(c_1(x),\dots ,c_m(x){)}^{*},x\in I$.

Damit erhalten wir

(20) $\begin{array}{c}P(x)\circ \mathrm{\Phi }(x)+q(x)={\mathrm{\Phi }}^{\prime }(x)=(Y(x)\circ c(x){)}^{\prime }=Y^{\prime }(x)\circ c(x)+Y(x)\circ c^{\prime }(x)\\ =P(x)\circ Y(x)\circ c(x)+Y(x)\circ c(x)\text{ bzw. }Y(x)\circ c^{\prime }(x)=q(x)\mathrm{f}\ddot{\mathrm{u}}\mathrm{r}\mathrm{a}\mathrm{l}\mathrm{l}\mathrm{e}x\in I.\end{array}$

Dieses lineare Gleichungssystem lösen wir mit der Cramerschen Regel und erhalten

(21) ${c_k}^{\prime }(x)=(-1{)}^{k+m}\frac{f(x)}{p_m(x)}\frac{W(y_1,\dots ,y_{k-1},y_{k+1},\dots ,y_m)}{W(y_1,\dots ,y_m)}$ in $I$

für $k=1,\dots ,m$. Hieraus folgt durch Integration die Behauptung.

q.e.d.

Bemerkungen

1. Das d'Alambertsche Verfahren der Reduktion der Ordnung
Haben wir bereits eine nullstellenfreie Lösung $u=u(x)\in C^m(I)$ der homogenen Differentialgleichung $ℒ(u)=0$ mit $u(x)\ne 0,x\in I$ gefunden, so leiten wir mit einem Produktansatz $y(x)=u(x)v(x),x\in I$ eine Differentialgleichung $(m-1)$-ter Ordnung für deren Ableitung $v^{\prime }(x)$ her. Hierzu berechnen wir zunächst

$y^{(k)}(x)=\sum _{l=0}^k(\genfrac{}{}{0ex}{}{k}{l})v^{(l)}(x)u^{(k-l)}(x)$ für $k=0,\dots ,m$

und ermitteln dann

(22) $\begin{array}{c}0=ℒ(y)=\sum _{k=0}^m{p}_k(x)y^{(k)}(x)=\sum _{k=0}^m{p}_k(x)(\sum _{l=0}^k(\genfrac{}{}{0ex}{}{k}{l})v^{(l)}(x)u^{(k-l)}(x))\\ =\sum _{k=0}^m{p}_k(x)(\sum _{l=0}^k(\genfrac{}{}{0ex}{}{k}{l})u^{(k-l)}(x)v^{(l)}(x))=\widetilde{ℒ}(v^{\prime }(x))\end{array}$

Hierbei stellt $\widetilde{ℒ}$ einen linearen Differentialoperator $(m-1)$-ter Ordnung angewandt auf die Funktion $v^{\prime }(x)$ dar.
2. Wenn die Koeffizienten der Differentialgleichung lokal in Potenzreihen entwickelbar sind, ist auch die Lösung der Differentialgleichung als Potenzreihe darstellbar. Wir haben dann nur die Koeffizienten dieser Potenzreihe zu bestimmen, wenn wir einen Potenzreihenansatz machen.

Beispiel: Die Besselsche Differentialgleichung

Die Untersuchung der Eigenschwingungen einer kreisförmigen Membran führt für die radiale Komponente der Schwingung auf die folgende Differentialgleichung

(23) $x^2{y}^{″}(x)+x{y}^{\prime }(x)+(x^2-n^2)y(x)=0,0<x<+\mathrm{\infty }$

mit $n\in \{0,1,2,\dots \}$. Wir wollen eine Lösung für diese Differentialgleichung ermitteln. Hierzu setzen wir die Potenzreihe $y(x)=\sum _{\nu =0}^{\mathrm{\infty }}{a}_{\nu }{x}^{\nu }$ mit unbestimmten $a_{\nu }\in ℝ$ in (23) ein und erhalten

(24) $0=\sum _{\nu =2}^{\mathrm{\infty }}\nu (\nu -1)a_{\nu }{x}^{\nu }+\sum _{\nu =1}^{\mathrm{\infty }}\nu a_{\nu }{x}^{\nu }-n^2\sum _{\nu =0}^{\mathrm{\infty }}{a}_{\nu }{x}^{\nu }+\sum _{\nu =2}^{\mathrm{\infty }}{a}_{\nu -2}{x}^{\nu }$.

Mittels Koeffizientenvergleich ist diese Gleichung äquivalent zu

(25) $-n^2{a}_0=0,(1-n^2)a_1=0,({\nu }^2-n^2)a_{\nu }+a_{\nu -2}=0$ für $\nu =2,3,4,\dots$.

Wir erhalten die Rekursionsformel

(26) $a_{\nu }=-\frac{1}{{\nu }^2-n^2}{a}_{\nu -2}$ für $\nu =2,3,4,\dots ,n-1,n+1,n+2,\dots$.

Setzen wir $a_n=a\in ℝ$, so folgt

(27) $a_{n+2\mu }=\frac{(-1{)}^{\mu }a}{\prod _{h=1}^{\mu }((n+2h{)}^2-n^2)}=\frac{(-1{)}^{\mu }a}{2^{2\mu }\mu !\prod _{h=1}^{\mu }(n+h)}$ für $\mu =0,1,2,\dots$,

während alle übrigen Koeffizienten der Potenzreihe verschwinden. Somit genügt die Funktion

(28) $y(x)=a{x}^n\sum _{\mu =1}^{\mathrm{\infty }}\frac{(-1{)}^{\mu }{x}^{2\mu }}{2^{2\mu }\mu !(n+1)\cdot (n+2)\cdot \dots \cdot (n+\mu )}$

der Besselschen Differentialgleichung (23), wobei $a\in ℝ$ beliebig ist.