Variation der Konstanten

Die Variation der Konstanten ist ein Verfahren aus der Theorie linearer gewöhnlicher Differentialgleichungen zur Bestimmung einer speziellen Lösung eines inhomogenen linearen Differentialgleichungssystems erster Ordnung bzw. einer inhomogenen linearen Differentialgleichung beliebiger Ordnung. Vorausgesetzt wird hierfür eine vollständige Lösung (Fundamentalsystem) der zugehörigen homogenen Differentialgleichung.

Leonhard Euler benutzte einen Vorläufer dieser Methode bereits 1748 im Zusammenhang mit astronomischen Problemen. In seiner heutigen Form wurde das Verfahren von dem Mathematiker Joseph-Louis Lagrange entwickelt.

Motivation

Lineare Differentialgleichung erster Ordnung

Seien {\displaystyle a:\mathbb {R} \rightarrow \mathbb {R} } und {\displaystyle b:\mathbb {R} \rightarrow \mathbb {R} } stetige Funktionen, dann lautet die lineare Differentialgleichung erster Ordnung

{\displaystyle y'(x)=a(x)y(x)+b(x)\ .}

Sei weiter A eine Stammfunktion von a, so gilt

{\displaystyle A(x):=\int _{x_{0}}^{x}a(t){\rm {d}}t\ ,}

wobei x_{0} geeigneten Randbedingungen genügen muss. Dann ist

{\displaystyle \left\{y(x)=ce^{A(x)}\ |\ c\in \mathbb {R} \right\}}

die Menge aller Lösungen der homogenen Differentialgleichung y'(x)=a(x)y(x).

Zur Lösung der inhomogenen Differentialgleichung wird nun die Funktion {\displaystyle c(x)} eingeführt und der Ansatz der Variation der Konstanten gewählt

{\displaystyle y(x)=c(x)e^{A(x)}}.

Dies ergibt eine eindeutige Zuordnung zwischen den Funktionen y und c, denn {\displaystyle \exp(A(x))} ist eine stets positive, stetig differenzierbare Funktion. Die Ableitung dieser Ansatzfunktion ist

{\displaystyle y'(x)=c(x)a(x)e^{A(x)}+c'(x)e^{A(x)}=a(x)y(x)+c'(x)e^{A(x)}\ .}

Also löst y die inhomogene Differentialgleichung

{\displaystyle y'(x)=a(x)y(x)+b(x)\ }

genau dann, wenn

{\displaystyle \ c'(x)=b(x)e^{-A(x)}}

gilt. Beispielsweise ist

{\displaystyle c(x):=\int _{x_{0}}^{x}b(t)e^{-A(t)}{\rm {d}}t}

eine solche Funktion und somit

{\displaystyle y_{sp}(x):=e^{A(x)}\cdot \int _{x_{0}}^{x}b(t)e^{-A(t)}{\rm {d}}t}

die spezielle Lösung mit y_{{sp}}(x_{0})=0. Also ist

{\displaystyle \left\{y(x)=e^{A(x)}\cdot \left[\int _{x_{0}}^{x}b(t)e^{-A(t)}{\rm {d}}t+c\right]\ |\ c\in \mathbb {R} \right\}}

die Menge aller Lösungen der inhomogenen Differentialgleichung {\displaystyle y'(x)=a(x)y(x)+b(x)}.

Beispiel

Liegt an einer Spule mit der Induktivität L und dem elektrischen Widerstand R eine Gleichspannung U_{0} an, so gilt für die Spannung an dem Widerstand

U(t)=U_{0}-L{\dot  I}(t).\!\,

Nach dem ohmschen Gesetz gilt zudem

{\displaystyle I(t)={\frac {U(t)}{R}}={\frac {U_{0}}{R}}-{\frac {L}{R}}{\dot {I}}(t)}.

Es handelt sich also um eine inhomogene lineare Differentialgleichung erster Ordnung mit konstanten Koeffizienten, die nun mithilfe des Verfahrens der Variation der Konstanten gelöst werden soll.

Für die zugehörige homogene Differentialgleichung I_{h}

{\displaystyle {\dot {I}}_{h}(t)=-{\frac {R}{L}}I_{h}(t)}

lautet die Lösung für jede Konstante c

{\displaystyle I_{h}(t)=ce^{-{\frac {R}{L}}t}}.

Als Ansatz für die Lösung der inhomogenen Differentialgleichung ersetze man die Konstante c durch einen variablen Ausdruck c(t). Man setzt also

I(t):=c(t)e^{{-{\frac  {R}{L}}t}}

und versucht, eine differenzierbare Funktion c so zu bestimmen, dass I die inhomogene Differentialgleichung erfüllt. Es folgt

{\displaystyle {\begin{aligned}{\frac {U_{0}}{R}}-{\frac {L}{R}}{\dot {I}}(t)&={\frac {U_{0}}{R}}-{\frac {L}{R}}{\dot {c}}(t)e^{-{\frac {R}{L}}t}+{\frac {L}{R}}c(t){\frac {R}{L}}e^{-{\frac {R}{L}}t}\\&={\frac {U_{0}}{R}}-{\frac {L}{R}}{\dot {c}}(t)e^{-{\frac {R}{L}}t}+c(t)e^{-{\frac {R}{L}}t}\\&={\frac {U_{0}}{R}}-{\frac {L}{R}}{\dot {c}}(t)e^{-{\frac {R}{L}}t}+I(t).\end{aligned}}}

Demnach ist die inhomogene Differentialgleichung genau dann gelöst, wenn gilt

{\frac  {U_{0}}{R}}-{\frac  {L}{R}}{\dot  c}(t)e^{{-{\frac  {R}{L}}t}}=0.

Diese Randwertbedingung ist gleichbedeutend mit {\displaystyle \textstyle {\dot {c}}(t)={\frac {U_{0}}{L}}e^{{\frac {R}{L}}t}} oder nach Integration mit \textstyle c(t)={\frac  {U_{0}}{R}}e^{{{\frac  {R}{L}}t}}+d. Somit lautet die Lösung der inhomogenen Differentialgleichung

I(t)={\frac  {U_{0}}{R}}+de^{{-{\frac  {R}{L}}t}}.

Die Konstante d lässt sich aus der Anfangsbedingung bestimmen und ergibt für I(0)=0 die Lösung

{\displaystyle I(t)={\frac {U_{0}}{R}}-{\frac {U_{0}}{R}}e^{-{\frac {R}{L}}t}}.

Inhomogene lineare Differentialgleichungssysteme erster Ordnung

Das obige Verfahren lässt sich auf folgende Weise verallgemeinern:

Formulierung

Seien A:{\mathbb  {R}}\rightarrow {\mathbb  {R}}^{{n\times n}} und b:{\mathbb  {R}}\rightarrow {\mathbb  {R}}^{n} stetige Funktionen und \Phi (x)=(y_{1}(x)\ |\ \cdots \ |\ y_{n}(x)) eine Fundamentalmatrix des homogenen Problems y'(x)=A(x)y(x) sowie \Phi _{k}(x) diejenige Matrix, die aus \Phi(x) entsteht, indem man die k-te Spalte durch b(x) ersetzt. Dann ist

y_{{sp}}(x):=\sum _{{k=1}}^{n}c_{k}(x)y_{k}(x)

mit

c_{k}(x):=\int _{{x_{0}}}^{x}{\frac  {\det \Phi _{k}(s)}{\det \Phi (s)}}{{\rm {d}}}s

die Lösung des inhomogenen Anfangswertproblems y'(x)=A(x)y(x)+b(x) und y(x_{0})=0.

Beweis

Setze

y_{{sp}}(x):=\Phi (x)\int _{{x_{0}}}^{x}\Phi (s)^{{-1}}b(s){{\rm {d}}}s\ .

Es ist y_{{sp}}(x_{0})=0, und wegen \Phi '(x)=A(x)\Phi (x) sieht man durch Differenzieren, dass y_{{sp}} die Differentialgleichung y_{{sp}}'(x)=A(x)y_{{sp}}(x)+b(x) erfüllt. Nun löst

a(s):=\Phi ^{{-1}}(s)b(s)\in {\mathbb  {R}}^{n}

für festes s das lineare Gleichungssystem

\Phi (s)\cdot a(s)=b(s)\ .

Nach der cramerschen Regel ist somit

a_{k}(s)={\frac  {\det \Phi _{k}(s)}{\det \Phi (s)}}\ ,\ k=1,\ldots ,n\ .

Also gilt

y_{{sp}}(x)=\int _{{x_{0}}}^{x}\Phi (x)a(s){{\rm {d}}}s=\sum _{{k=1}}^{n}\left[\int _{{x_{0}}}^{x}{\frac  {\det \Phi _{k}(s)}{\det \Phi (s)}}{{\rm {d}}}s\right]y_{k}(x)\ .

Spezialfall: Resonanzfall

Falls die Inhomogenität b selber Lösung des homogenen Problems ist, d. h. b'(x)=A(x)b(x), so bezeichnet man dies als Resonanzfall. In diesem Fall ist

\ y_{{sp}}(x):=(x-x_{0})b(x)

die Lösung des inhomogenen Anfangswertproblems y'(x)=A(x)y(x)+b(x) und y(x_{0})=0.

Inhomogene lineare Differentialgleichungen höherer Ordnung

Das Lösen einer Differentialgleichung höherer Ordnung ist äquivalent zum Lösen eines geeigneten Differentialgleichungssystems erster Ordnung. Auf diese Weise kann man obiges Verfahren nutzen, um eine spezielle Lösung für eine Differentialgleichung höherer Ordnung zu konstruieren.

Formulierung

Seien a_{0},\ldots ,a_{{n-1}},b:{\mathbb  {R}}\rightarrow {\mathbb  {R}} stetige Funktionen und \Phi(x) eine Fundamentalmatrix des homogenen Problems \textstyle y^{{(n)}}(x)=\sum _{{k=0}}^{{n-1}}a_{k}(x)y^{{(k)}}(x), deren erste Zeile (y_{1}(x)\ |\ \cdots \ |\ y_{n}(x)) lautet, sowie \Phi _{k}(x) diejenige Matrix, die aus \Phi(x) entsteht, indem man die k-te Spalte durch \textstyle {\begin{pmatrix}0\\\vdots \\0\\b(x)\\\end{pmatrix}} ersetzt. Dann ist

y_{{sp}}(x):=\sum _{{k=1}}^{n}c_{k}(x)y_{k}(x)

mit

c_{k}(x):=\int _{{x_{0}}}^{x}{\frac  {\det \Phi _{k}(s)}{\det \Phi (s)}}{{\rm {d}}}s

die Lösung des inhomogenen Anfangswertproblems \textstyle y^{{(n)}}(x)=\sum _{{k=0}}^{{n-1}}a_{k}(x)y^{{(k)}}(x)+b(x) und y(x_{0})=0.

Beweis

Man betrachte zunächst das hierzu korrespondierende Differentialgleichungssystem erster Ordnung, bestehend aus n Gleichungen

A(x):={\begin{pmatrix}0&1&&0\\&\ddots &\ddots &\\&&\ddots &1\\a_{0}(x)&a_{1}(x)&\cdots &a_{{n-1}}(x)\\\end{pmatrix}}\ ,\ B(x):={\begin{pmatrix}0\\\vdots \\0\\b(x)\\\end{pmatrix}}\ .

Es gilt: y(x) löst die skalare Gleichung n-ter Ordnung genau dann, wenn \textstyle Y(x):={\begin{pmatrix}y(x)\\y'(x)\\\vdots \\y^{{(n-1)}}(x)\\\end{pmatrix}} Lösung obigen Systems erster Ordnung ist. Per definitionem ist \Phi eine Fundamentalmatrix für dieses System erster Ordnung. Darauf wende man schließlich das oben bewiesene Verfahren der Variation der Konstanten an.

Alternative: Grundlösungsverfahren

Im Fall konstanter Koeffizienten ist es gelegentlich von Vorteil, das Grundlösungsverfahren zur Konstruktion einer speziellen Lösung zu verwenden: Ist y_{h} diejenige homogene Lösung von \textstyle y^{{(n)}}(x)=\sum _{{k=0}}^{{n-1}}a_{k}y^{{(k)}}(x), welche

y_{h}^{{(k)}}(x_{0})=0\ ,\ k=0,\ldots ,n-2\ ,\ y_{h}^{{(n-1)}}(x_{0})=1

erfüllt, dann ist

y_{{sp}}(x):=\int _{{x_{0}}}^{x}y_{h}(x_{0}+x-t)b(t){{\rm {d}}}t

diejenige spezielle Lösung von \textstyle y^{{(n)}}(x)=\sum _{{k=0}}^{{n-1}}a_{k}y^{{(k)}}(x)+b(x) mit y_{{sp}}(x_{0})=0.

Beweis

Durch Differenzieren überprüft man

y_{{sp}}^{{(k)}}(x)=\int _{{x_{0}}}^{x}y_{h}^{{(k)}}(x_{0}+x-t)b(t){{\rm {d}}}t\ ,\ k=0,\ldots ,n-1

und

y_{{sp}}^{{(n)}}(x)=b(x)+\int _{{x_{0}}}^{x}y_{h}^{{(n)}}(x_{0}+x-t)b(t){{\rm {d}}}t\ .

Es ergibt sich

y_{{sp}}^{{(n)}}(x)-\sum _{{k=0}}^{{n-1}}a_{k}y_{{sp}}^{{(k)}}(x)=b(x)+\int _{{x_{0}}}^{x}\left[y_{h}^{{(n)}}-\sum _{{k=0}}^{{n-1}}a_{k}y_{h}^{{(k)}}\right](x_{0}+x-t)b(t){{\rm {d}}}t=b(x)\ .
Trenner
Basierend auf einem Artikel in: Wikipedia.de
Seitenende
Seite zurück
©  biancahoegel.de
Datum der letzten Änderung:  Jena, den: 24.09. 2022