Variation der Konstanten

Variation der Konstanten

Studium
yp(x)=C(x) e−A(x),C′(x)=b(x) eA(x)y_p(x) = C(x)\,e^{-A(x)}, \qquad C'(x) = b(x)\,e^{A(x)}

Voraussetzung: y′+a(x) y=b(x),y' + a(x)\,y = b(x), A′(x)=a(x)A'(x) = a(x)

ypy_p
partikuläre (spezielle) Lösung der inhomogenen Gleichung
C(x)C(x)
Funktion, die an die Stelle der Konstanten C der homogenen Lösung tritt
A(x)A(x)
eine Stammfunktion von a(x)
a(x)a(x)
Koeffizient vor y in der Normalform y′ + a(x)·y = b(x)
b(x)b(x)
Störfunktion (Inhomogenität)

Man ersetzt die Konstante der homogenen Lösung durch eine Funktion C(x) und erhält so eine Lösung der inhomogenen Gleichung.

Erklärung

Die Variation der Konstanten liefert eine partikuläre Lösung einer inhomogenen linearen Differentialgleichung, wenn die homogene Lösung bekannt ist. Für y′+a(x) y=b(x)y' + a(x)\,y = b(x) geht man in drei Schritten vor:

  1. Die homogene Gleichung y′+a(x) y=0y' + a(x)\,y = 0 lösen: yh=C e−A(x)y_h = C\,e^{-A(x)}.
  2. Die Konstante durch eine Funktion ersetzen: Ansatz yp=C(x) e−A(x)y_p = C(x)\,e^{-A(x)}.
  3. Einsetzen. Die Terme mit C(x)C(x) heben sich weg, übrig bleibt C′(x) e−A(x)=b(x)C'(x)\,e^{-A(x)} = b(x), also C′(x)=b(x) eA(x)C'(x) = b(x)\,e^{A(x)}. Integrieren liefert C(x)C(x).

Die allgemeine Lösung ist dann y=yh+ypy = y_h + y_p. Das Ergebnis ist dasselbe wie mit der Lösungsformel für lineare DGL erster Ordnung; die Variation der Konstanten zeigt, woher diese Formel kommt.

Das Verfahren funktioniert für jede stetige Störfunktion, auch wenn ein Ansatz vom Typ der rechten Seite versagt. Es lässt sich auf lineare Gleichungen höherer Ordnung übertragen; dort variiert man alle Konstanten der homogenen Lösung und erhält ein lineares Gleichungssystem für die Ci′(x)C_i'(x).

Beispiele

y′−y=xy' - y = x. Hier ist a=−1a = -1, A(x)=−xA(x) = -x, die homogene Lösung yh=C exy_h = C\,e^{x}. Ansatz yp=C(x) exy_p = C(x)\,e^{x}:

C′(x) ex=x⇒C′(x)=x e−x⇒C(x)=−(x+1) e−xC'(x)\,e^{x} = x \quad\Rightarrow\quad C'(x) = x\,e^{-x} \quad\Rightarrow\quad C(x) = -(x + 1)\,e^{-x}

(partielle Integration; Probe durch Ableiten). Damit ist yp=−(x+1)y_p = -(x + 1) und die allgemeine Lösung

y=K ex−x−1y = K\,e^{x} - x - 1

Probe: yp′−yp=−1+x+1=xy_p' - y_p = -1 + x + 1 = x.

y′+2y=e−2xy' + 2y = e^{-2x}. Die homogene Lösung ist C e−2xC\,e^{-2x}. Ein Ansatz yp=c e−2xy_p = c\,e^{-2x} scheitert, weil er selbst homogene Lösung ist. Mit der Variation der Konstanten:

C′(x)=e−2x⋅e2x=1⇒C(x)=x⇒yp=x e−2xC'(x) = e^{-2x} \cdot e^{2x} = 1 \quad\Rightarrow\quad C(x) = x \quad\Rightarrow\quad y_p = x\,e^{-2x}

Probe: yp′+2yp=e−2x−2x e−2x+2x e−2x=e−2xy_p' + 2y_p = e^{-2x} - 2x\,e^{-2x} + 2x\,e^{-2x} = e^{-2x}.

Typische Fehler

  • Beim Einsetzen nicht beachten, dass C(x)C(x) jetzt abgeleitet werden muss: yp′=C′(x) e−A−a(x) C(x) e−Ay_p' = C'(x)\,e^{-A} - a(x)\,C(x)\,e^{-A}.
  • C′(x)=b(x) e−A(x)C'(x) = b(x)\,e^{-A(x)} statt b(x) e+A(x)b(x)\,e^{+A(x)} schreiben.
  • Am Ende die homogene Lösung vergessen: Die allgemeine Lösung ist yh+ypy_h + y_p, nicht nur ypy_p.
  • Die Gleichung nicht auf Normalform mit Faktor 11 vor y′y' bringen, bevor man b(x)b(x) abliest.

Verwandte Formeln