Getb:Vorgehen zur Lösung linearer Differenzialgleichungen 2. Ordnung

Gegeben sei eine inhomogene, lineare Differenzialgleichung 2. Ordnung mit konstanten Koeffizienten. Diese DGL lässt sich allgemein schreiben als:

$y''(t)+a \cdot y'(t)+b \cdot y(t)=g(t)$

Wobei $a$ und $b$ Konstanten sind. Die Lösung erfolgt genauso wie die Lösung einer DGL 1. Ordnung in vier Schritten:

1. Lösung der dazugehörigen homogenen DGL
2. Finden einer partikulären Lösung der inhomogenen DGL
3. Addieren der beiden Lösungen liefert die allgemeine Lösung der inhomogenen DGL
4. Berücksichtigung von Anfangsbedingungen

Lösung der homogenen linearen Differenzialgleichung

Die zur inhomogenen DGL dazugehörige homogene DGL lautet:

$y''(t)+a \cdot y'(t)+b \cdot y(t)=0$

Die allgemeine Lösung einer Differenzialgleichung 2. Ordnung enthält zwei Parameter. Im Falle einer linearen DGL 2. Ordnung mit konstanten Koeffizienten ist die allgemeine Lösung eine Linearkombination von zwei linear unabhängigen partikulären Lösungen. Die Koeffizienten der partikulären Lösungen sind die beiden Parameter. Formal gilt dann für die allgemeine Lösung:

$y(t)=K_1 \cdot y_1(t)+K_2 \cdot y_2(t)$ mit $K_1,K_2 \in \R$, wobei weiterhin gelten muss:

• $y_1''(t)+a\cdot y_1'(t)+b \cdot y_1(t)=0$      (y1 ist partikuläre Lösung der homogenen DGL)
• $y_2''(t)+a\cdot y_2'(t)+b \cdot y_2(t)=0$      (y2 ist partikuläre Lösung der homogenen DGL)
• $y_1$ und $y_2$ sind linear unabhängig

Beide Lösungen müssen also die DGL erfüllen und linear unabhängig sein. Allgemein gilt, dass die Lösungen $y_1$ und $y_2$ auch komplexwertig sein können.

Zwei linear unabhängige partikuläre Lösungen lassen sich mit Hilfe des komplexen Exponentialansatzes $y(t)=\underline{K} \cdot \operatorname{e}^{s \cdot t}$ finden. Hierbei ist zu beachten, dass mit $\underline{K}$ schon eine Konstante (bzw. ein Parameter) in den Ansatz eingebracht wurde. Der Unterstrich unter dem $\underline{K}$ deutet an, dass $\underline{K}$ komplexe Werte annehmen kann. Auch die Variable $s$ kann komplexe Werte annehmen, wird üblicherweise aber nicht unterstrichen. Für den Fall, dass $s$ komplexe Werte annimmt, wird der Ausdruck der e-Funktion selbst auch eine komplexe Zahl (vgl. z.B. Papula (2011): „Mathematik für Ingenieure und Naturwissenschaftler“, Bd. 1, Kapitel VII p. 652 ff.). In den meisten Fällen (und auch in der Vorlesung GET B) ist man aber nur an reell wertigen Lösungen für die DGL interessiert, da die gesuchten Ströme oder Spannungen reelle Größen sind. Deshalb ist es notwendig, komplexwertige Konstanten $\underline{K}$ anzunehmen. Im weiteren Verlauf wird deutlich, dass dadurch, mit bestimmten Anforderungen an die komplexen Konstanten, eine allgemeine komplexe Lösung auf eine allgemeine reelle Lösung der DGL zurückgeführt werden kann.

Der Exponentialansatz und dessen Ableitungen $y'(t)=s \cdot \underline{K} \cdot \operatorname{e}^{s \cdot t}$ und $y''(t)=s^2 \cdot \underline{K} \cdot \operatorname{e}^{s \cdot t}$ können nun in die homogene DGL eingesetzt werden:

$y''(t)+a \cdot y'(t)+b \cdot y(t) = s^2 \cdot \underline{K} \cdot \operatorname{e}^{s \cdot t} + a \cdot s \cdot \underline{K} \cdot \operatorname{e}^{s \cdot t} + b \cdot \underline{K} \cdot \operatorname{e}^{s \cdot t} = 0$

$\Leftrightarrow (s^2 + a \cdot s + b) \cdot \underline{K} \cdot \operatorname{e}^{s \cdot t} = 0$

Diese Gleichung ist erfüllt, wenn einer der Faktoren 0 ist. Die e-Funktion kann niemals 0 sein, und der Fall $\underline{K}=0$ führt zu der trivialen Lösung $y=0$. Somit muss der Ausdruck in der Klammer zu Null werden. Die daraus resultierende Bedingung $s^2+a \cdot s+b=0$ heißt charakteristische Gleichung und kann zum Beispiel mit der p-q-Formel aufgelöst werden:

$s_{1/2} = - \dfrac{a}{2} \pm \sqrt{\left( \dfrac{a}{2} \right)^2 - b} = -\dfrac{a}{2} \pm \dfrac{\sqrt{a^2-4b}}{2}$

Die Lösungen $s_{1/2}$ heißen Eigenwerte. Es ergeben sich drei verschiedene Lösungsfälle:

• $a^2-4b > 0$ → Die charakteristische Gleichung hat zwei verschiedene, reelle Lösungen

• $a^2-4b = 0$ → Die charakteristische Gleichung hat zwei gleiche, reelle Lösungen

• $a^2-4b < 0$ → Die charakteristische Gleichung hat zwei zueinander konjugiert komplexe Lösungen (komplexe Lösungen treten stets als konjugiert komplexe Paare auf)

Fall 1 $(a^2-4b > 0)$

Die beiden voneinander verschiedenen, reellen Eigenwerte $\textstyle s_{1/2}=-\frac{a}{2} \pm \frac{\sqrt{a^2-4b}}{2}$ führen mit dem Exponentialansatz auch zu zwei verschiedenen reellen Lösungen der homogenen DGL:

$y_1 (t)=K_1 \cdot \operatorname{e}^{s_1 \cdot t}$ und $y_2 (t)=K_2 \cdot \operatorname{e}^{s_2 \cdot t}$ mit $K_1,K_2 \in \R$. Da die Exponentialfunktionen reell sind, müssen auch die Konstanten $K_1$ und $K_2$ reell sein, damit die resultierende allgemeine Lösung reell ist. Die beiden Lösungen $y_1 (t)$ und $y_2 (t)$ sind linear unabhängig. Deshalb gilt für die allgemeine Lösung der homogenen DGL 2. Ordnung mit konstanten Koeffizienten, bei der zwei unterschiedliche reelle Eigenwerte auftreten:

$y_h (t)=y_1 (t)+y_2 (t)=K_1 \cdot \operatorname{e}^{s_1 \cdot t}+K_2 \cdot \operatorname{e}^{s_2 \cdot t}$

Fall 2 $(a^2-4b = 0)$

Da die beiden Eigenwerte $\textstyle s_{1/2}=-\frac{a}{2}$ gleich sind, ergibt sich aus dem Exponentialansatz nur eine Lösung:

$y_1 (t)=K_1 \cdot \operatorname{e}^{s_1 \cdot t}=K_1 \cdot \operatorname{e}^{s_2 \cdot t}=K_1 \cdot \operatorname{e}^{-\frac{a \cdot t}{2}} \text{ mit } K_1 \in \R$

Die Exponentialfunktion ist reell, deshalb muss auch der Parameter $K_1$ reell sein, da für $y_1$ nur reelle Lösungen in Frage kommen. Der Exponentialansatz führt in diesem Fall aber nur zu einer Lösungsfunktion, es werden jedoch zwei linear unabhängige Lösungsfunktionen für die allgemeine Lösung benötigt. Mit dem Verfahren „Variation der Konstanten“ (Idee: Einsetzen der Lösung aus dem Exponentialansatz und Annahme einer veränderlichen Konstante $K_1 (t)$) lässt sich eine weitere unabhängige Lösung $y_2 (t)=K_2 \cdot t \cdot \operatorname{e}^{-\frac{a \cdot t}{2}}$ mit $K_2 \in \R$ finden. Dass dies tatsächlich eine Lösung der DGL ist, lässt sich durch Einsetzen in die DGL verifizieren. Die allgemeine Lösung für eine homogene DGL 2. Ordnung mit konstanten Koeffizienten, bei der zwei gleiche Eigenwerte auftreten, ergibt sich damit durch:

$y_h (t)=y_1 (t)+y_2 (t)=K_1 \cdot \operatorname{e}^{s \cdot t}+K_2 \cdot t \cdot \operatorname{e}^{s \cdot t}=K_1 \cdot \operatorname{e}^{-\frac{a \cdot t}{2}}+K_2 \cdot t \cdot \operatorname{e}^{-\frac{a \cdot t}{2}}$