eine gewöhnliche Differentialgleichungn-ter Ordnung für die Funktion y, die in y, y^′, …, y⁽ⁿ⁾ linear ist.

Mit 𝕂 = ℝ oder 𝕂 = ℂ, einem offenen Intervall I ⊂ ℝ, und stetigen Funktionen b, a_i : I → 𝕂 hat sie also die Form

\begin{eqnarray}{y}^{(n)}+{a}_{n-1}(x){y}^{(n-1)}+\cdots +{a}_{1}(x){y}^{\prime}+{a}_{0}(x)y=b(x).\end{eqnarray}

b heißt Inhomogenität der Differentialgleichung. Falls b(x) = 0 für alle x ∈ I, so heißt die lineare Differentialgleichung homogene Differentialgleichung, sonst inhomogene Differentialgleichung. Eine Differentialgleichung für y, die nicht in allen y, y^′, …, y⁽ⁿ⁾ linear ist, heißt nichtlineare Differentialgleichung.

Die 𝕂-wertigen Lösungen (Lösung einer Differentialgleichung) der Gleichung (1) bilden einen n-dimensionalen affinen Raum über 𝕂. Die Lösungen der zugehörigen homogenen Gleichung

\begin{eqnarray}{y}^{(n)}+{a}_{n-1}(x){y}^{(n-1)}+\ldots +{a}_{1}(x){y}^{\prime}+{a}_{0}(x)y=0(2)\end{eqnarray}

bilden einen n-dimensionalen Vektorraum.

Jede Lösung von (1) ist von der Form y = y_p + y_h. Dabei ist y_p eine partikuläre (spezielle) Lösung der inhomogenen Gleichung (1) und y_h eine geeignete Lösung der homogenen Gleichung (2), d. h., mit einem Fundamentalsystemy₁, …, y_n von (2) ist

\begin{eqnarray}{y}_{h}={c}_{1}{y}_{1}+\ldots +{c}_{n}{y}_{n}\end{eqnarray}

mit geeigneten c_i ∈ 𝕂.

Sind also y₁, y₂ zwei Lösungen der linearen Gleichung (1), so ist y₁ − y₂ eine Lösung der zugehörigen homogenen Gleichung (2). Zum Lösen einer linearen Differentialgleichung benötigt man also ein Fundamentalsystem der homogenen Gleichung und eine partikuläre Lösung der inhomogenen Gleichung. Ein allgemeines Verfahren zur Bestimmung eines Fundamentalsystems existiert nur für den Spezialfall der linearen Differentialgleichung mit konstanten Koeffizienten. Ist allerdings ein Fundamentalsystem der homogenen Gleichung bekannt, so kann man eine partikuläre Lösung der inhomogenen Gleichung durch Variation der Konstanten erzeugen.

Mit y₁ := y, y₂ := y^′, …, y_n := y⁽ⁿ⁻¹⁾, y := (y₁, …, y_n)^T, b(t) := (0, …, 0, b(t))^T und

\begin{eqnarray}A(t):=\left(\begin{array}{ccccc}0 & 1 & 0 & \cdots & 0\\ \vdots & \ddots & \ddots & & \vdots \\ \vdots & & \ddots & \ddots & \vdots \\ 0 & 0 & \cdots & 0 & 1\\ -{a}_{0}(t) & -{a}_{1}(t) & \cdots & -{a}_{n-2}(t) & -{a}_{n-1}(t)\end{array}\right)\end{eqnarray}

ist die lineare Differentialgleichung n-ter Ordnung (1) äquivalent zu dem folgenden System von n linearen Gleichungen 1. Ordnung (lineares Differentialgleichungssystem):

\begin{eqnarray}{\bf{y}}^{\prime}=A(t){\bf{y}}+{\bf{b}}(t).\end{eqnarray}

Viele Sätze über lineare Differentialgleichungssysteme können so direkt auf lineare Differentialgleichungen n-ter Ordnung übertragen werden.

Von besonderem Interesse ist der Spezialfall der linearen Differentialgleichung zweiter Ordnung, der in vielen physikalischen Anwendungen auftritt. Mit Funktionen p, q, f hat sie die allgemeine Form

\begin{eqnarray}{y}^{\prime\prime}(x)+p(x){y}^{\prime}(x)+q(x)y(x)=f(x).\end{eqnarray}

Die Inhomogenität f wird bisweilen als Störung der homogenen Differentialgleichung bezeichnet, da sie in vielen Anwendungen als eine von außen wirkende Störung aufgefaßt werden kann.

Falls p und q Konstanten sind (lineare Differentialgleichung mit konstanten Koeffizienten), gibt es geschlossene Formeln für die allgemeine Lösung der homogenen linearen Differentialgleichung

\begin{eqnarray}{y}^{\prime\prime}(x)+p{y}^{\prime}(x)+qy(x)=0.\end{eqnarray}

Das zugehörige charakteristische Polynom

\begin{eqnarray}{\lambda }^{2}+p\lambda +q=0\end{eqnarray}

hat die Nullstellen

\begin{eqnarray}{\lambda }_{1,2}=-\frac{p}{2}\pm \sqrt{\frac{{p}^{2}}{4}-q}.\end{eqnarray}

Abhängig vom Radikanden in (4) ergeben sich drei Fälle:

1.) p²/4−q > 0 : Die beiden reellen verschiedenen Nullstellen λ_1,2 führen zu den linear unabhängigen Lösungen

\begin{eqnarray}\begin{array}{cc}{y}_{1}(x)={e}^{{\lambda }_{1}\cdot x}, & {y}_{2}(x)={e}^{{\lambda }_{2}\cdot x}\end{array}.\end{eqnarray}

2.) p²/4 − q = 0 : Die reelle doppelte Nullstelle λ führt zu den beiden linear unabhängigen Lösungen

\begin{eqnarray}\begin{array}{cc}{y}_{1}(x)={e}^{\lambda \cdot x}, & {y}_{2}(x)=x\cdot {e}^{\lambda \cdot x}\end{array}.\end{eqnarray}

3.) p²/4 − q< 0 : Die beiden konjugiert komplexen Nullstellen λ_1,2 = α ± i β führen zu den beiden linear unabhängigen reellen Lösungen

\begin{eqnarray}\begin{array}{cc}{y}_{1}(x)={e}^{\alpha \cdot x}\cdot \cos (\beta \cdot x), & {y}_{2}(x)={e}^{\alpha \cdot x}\cdot \sin (\beta \cdot x)\end{array}.\end{eqnarray}

Damit ergibt sich für die allgemeine Lösung von (4) die Formel

\begin{eqnarray}y(x)={c}_{1}\cdot {y}_{1}(x)+{c}_{2}\cdot {y}_{2}(x)\end{eqnarray}

mit Konstanten c₁, c₂ ∈ ℝ, die man für eine konkrete Lösung aus Anfangs- bzw. Randbedingungen erhält.

Für die systematische Berechnung einer partikulären Lösung der inhomogenen Differentialgleichung kann das Verfahren der Variation der Konstanten verwendet werden.

[1] Kamke, E.: Differentialgleichungen, Lösungsmethoden und Lösungen I. B. G. Teubner Stuttgart, 1977.
[2] Walter, W.: Gewöhnliche Differentialgleichungen. Springer-Verlag Berlin, 1972.