Ing Mathematik: Differentialgleichungen

Einleitung und Definition

Gewöhnliche Differenzialgleichungen sind Gleichungen der Form $F(x,y,y^{\prime },\dots ,y^{(n)})=0$. Die Ordnung der Differentialgleichung ist $n$.

Eine gewöhnliche Differenzialgleichung 1. Ordnung ist somit eine Gleichung der Form $F(x,y,y^{\prime })=0$.

Der Begriff gewöhnliche Differentialgleichung wird im Englischen auch mit ODE (Ordinary Differential Equation) abgekürzt. Wenn hier von Dgl. die Rede ist, dann ist damit immer eine gewöhnliche Differentialgleichung gemeint. Es gibt nämlich auch partielle Differentialgleichungen (PDE). Die werden aber erst später in dieser Buchreihe behandelt.

Beispiele für Dgln. 1. Ordnung:

• $y^{\prime }=x$
• $y^{\prime }=x^2+y^2$

Beispiel für eine Dgl. 2. Ordnung:

• $A{y}^{″}(x)+B{y}^{\prime }(x)+Cy(x)=f(x)$

Richtungsfeld

Nachfolgend sei das Richtungsfeld der Dgl. $y^{\prime }=y-x$ dargestellt (siehe auch Vorlage:W):

Existenz- und Eindeutigkeitssatz von Picard-Lindelöf

Siehe vorerst Vorlage:W

• 
  Ernst Leonard Lindelöf (1870-1946) finnischer Mathematiker
• 
  Charles Émile Picard (1856-1941) französischer Mathematiker

Separation der Variablen

Nun kommen wir zu den analytischen Lösungsmöglichkeiten für Dgln. Das Problem dabei ist, dass es keine geschlossene Lösungstheorie für Dgln. gibt. Stattdessen gibt es eine Vielzahl von Methoden, die jeweils auf eine bestimmte Klasse von Dgln. zugeschnitten ist. Die wahrscheinlich einfachste Möglichkeit ist die Separation (Trennung) der Variablen.

Beispiel: Löse die Gleichung $y^{\prime }=\frac{\mathrm{d}y}{\mathrm{d}x}={\text{e}}^x$.

Diese Dgl. ist elementar zu lösen. Man kann $\mathrm{d}x$ auf die rechte Gleichungsseite bringen und dann integrieren. Man separiert sozusagen die Gleichungsbestandteile $x$ und $y$.

$\int \mathrm{d}y=\int {\text{e}}^x\mathrm{d}x\Rightarrow y={\text{e}}^x+C$

Die Konstante C ist z.B. aus einer Anfangsbedingung zu ermitteln.

Beispiel: $y^{\prime }=xy$

$\frac{\mathrm{d}y}{y}=x\mathrm{d}x$

$\int \frac{\mathrm{d}y}{y}=\int x\mathrm{d}x$

$\log |y|-\stackrel{-D}{\overbrace{\log C}}=\frac{x^2}{2}$

$\log |\frac{y}{C}|=\frac{x^2}{2}$

$y=C{\text{e}}^{\frac{x^2}{2}}$

Übung: Löse folgende Dgl.

• $(1+y^2)\mathrm{d}x+xy\mathrm{d}y=0$

Manchmal muss man zuerst substituieren, um anschließend eine Separation der Variablen durchführen zu können.

Beispiel: $y^{\prime }=1+\frac{y}{x}$

Substitution: $z=\frac{y}{x}\Rightarrow$ $y=xz;$ $y^{\prime }=z+x{z}^{\prime }$

$z+x{z}^{\prime }=1+\frac{xz}{x}=1+z$

$x\frac{\mathrm{d}z}{\mathrm{d}x}=1$

$\int \mathrm{d}z=\int \frac{\mathrm{d}x}{x}$

$z=\log |x|+C$

$y=x(\log |x|+C)$

Beispiel: $2y^{\prime }+y^2+\frac{1}{x^2}=0$

Das lässt sich so umformen: $2y^{\prime }{x}^2+(xy{)}^2+1=0$

Probieren wir nun folgende Substitution aus: $z=xy\Rightarrow y=\frac{z}{x};y^{\prime }=\frac{z^{\prime }x-z}{x^2}$

$2\frac{z^{\prime }x-z}{x^2}{x}^2+z^2+1=0$

$2\frac{\mathrm{d}z}{\mathrm{d}x}x=2z-z^2-1$

$-2\frac{\mathrm{d}z}{\mathrm{d}x}x=z^2-2z+1=(z-1{)}^2$

$\int \frac{\mathrm{d}z}{(z-1{)}^2}=-\int \frac{\mathrm{d}x}{2x}$

$\frac{1}{1-z}=-\frac{1}{2}\log |x|$

$\frac{1}{1-xy}=-\frac{1}{2}\log |x|+C$

$2=(1-xy)(C-\log |x|)$

Beispiel: $y^{\prime }=ax+by+c$

Hier substituieren wir direkt $z=ax+by+c$

$z^{\prime }=a+b{y}^{\prime };y^{\prime }=\frac{z^{\prime }-a}{b}$

$\frac{z^{\prime }-a}{b}=z$

Und nun separieren wir wieder die Variablen: $\int \frac{\mathrm{d}z}{a+bz}=\int \mathrm{d}x$

Integriert führt dies auf folgende Formel: $\log |a+bz|-\log D=bx$

Rücksubstituieren: $\log \frac{a+b(ax+by+c)}{D}=bx$

$a+b(ax+by+c)=D{\text{e}}^{bx}$

Homogene Dgl.

Homogene Funktion: $f(\lambda x_1,\lambda x_2,\dots \lambda x_n)={\lambda }^{m}f(x_1,x_2,\dots x_n)$ mit $m$ als Homogenitätsgrad.

Seien $M(x,y)$ und $N(x,y)$ homogene Funktionen vom gleichen Grad. Die Differentialgleichung laute $M(x,y)\mathrm{d}x+N(x,y)\mathrm{d}y=0$. Dann führt die Substitution $z=\frac{y}{x}$ zum Ziel.

Beispiel: $(x+y)\mathrm{d}x+(x-y)\mathrm{d}y=0$

Dies ist eine homogene Gleichung vom Homogenitätsgrad 1.

Substitution: $y=zx;y^{\prime }=z^{\prime }x+z$

$x+zx\mathrm{+}(x-zx)(z^{\prime }x+z)=0$

$x+zx+x^2{z}^{\prime }-z{x}^2{z}^{\prime }+zx-z^{2}x=0$

$1+z+x{z}^{\prime }-zx{z}^{\prime }+z-z^2=0$

$z^{\prime }x(1-z)+2z-z^2+1=0$

$\frac{\mathrm{d}z}{\mathrm{d}x}x=\frac{z^2-2z-1}{1-z}$

$\int \frac{1-z}{z^2-2z-1}\mathrm{d}z=\int \frac{\mathrm{d}x}{x}$

Diese Integrale können gelöst werden und dann wird wieder rücksubstituiert. Dies machen wir hier aber nicht und beenden dieses Beispiel somit. Die Berechnung mit Maxima ergibt:

Exakte Dgl. und integrierender Faktor

Siehe z.B. Vorlage:W

Lineare Dgl. 1. Ordnung

Homogene lineare Dgl.

Obwohl sie ähnlich heißen, haben homogene lineare Dgln. nichts mit den vorher behandelten homogenen Dgln. zu tun. Es liegen hier Dgln. folgender Form vor:

$y^{\prime }-g(x)y=0$

Eine Separation der Variablen führt zum Ziel:

$\frac{\mathrm{d}y}{y}=g(x)\mathrm{d}x$

$\log |y|-\log C=\int g(x)\mathrm{d}x$

$|y|=C{\text{e}}^{\int g(x)\mathrm{d}x}$

Beispiel: $y^{\prime }-2\sin (x)y=0$

$\int \frac{\mathrm{d}y}{y}=2\int \sin (x)\mathrm{d}x$

$\log |y|-\log C=-2\cos (x)$

$|y|=C{\text{e}}^{-2\cos (x)}$

Inhomogene lineare Dgl.

Es liegen Dgln. folgender Form vor:

$y^{\prime }-g(x)y=h(x)$

Wir können die homogene lineare Gleichung wie oben lösen:

$y^{\prime }-g(x)y=0$

$y=C{\text{e}}^{G(x)}$ mit $G(x)=\underset{x_0}{\overset{x}{\int }}(g(x))\mathrm{d}x$.

Variation der Konstanten

Wir setzen bei dieser Methode

$y(x)=C(x){\text{e}}^{G(x)}$

$y^{\prime }(x)=C^{\prime }(x){\text{e}}^{G(x)}+C(x)G^{\prime }(x){\text{e}}^{G(x)}$

$C^{\prime }(x){\text{e}}^{G(x)}+\cancel{C(x)\stackrel{g(x)}{\overbrace{G^{\prime }(x)}}{\text{e}}^{G(x)}}=\cancel{g(x)C(x){\text{e}}^{G(x)}}+h(x)$

Wir dividieren durch ${\text{e}}^{G(x)}$

$C^{\prime }(x)=h(x){\text{e}}^{-G(x)}$.

und integrieren

$C(x)=\underset{x_0}{\overset{x}{\int }}h(t){\text{e}}^{-G(t)}\mathrm{d}t+C_1$

Jetzt setzen wir ein

$y(x)=C(x){\text{e}}^{G(x)}=(\underset{x_0}{\overset{x}{\int }}h(t){\text{e}}^{-G(t)}\mathrm{d}t+C_1){\text{e}}^{G(x)};C_1=y_0$

Beispiel: $y^{\prime }-\frac{y}{x}=x;y(1)=0$

$g(x)=\frac{1}{x};h(x)=x;(x_0,y_0)=(1,0)$

$G(x)=\underset{x_0}{\overset{x}{\int }}g(t)\mathrm{d}t=\underset{1}{\overset{x}{\int }}\frac{dt}{t}=\log |x|$

$y(x)=(\underset{x_0}{\overset{x}{\int }}h(t){\text{e}}^{-G(t)}\mathrm{d}t+y_0){\text{e}}^{G(x)}=(\underset{1}{\overset{x}{\int }}t{\mathrm{e}}^{-\log |t|}\mathrm{d}t+0){\mathrm{e}}^{\log |x|}=\underset{1}{\overset{x}{\int }}t\frac{1}{t}\mathrm{d}tx=x(x-1)$

Beispiel: $y^{\prime }+ay=x;$ mit $a=\text{konstant}$

Lösung der homogenen Dgl.: $y^{\prime }=-ay\Rightarrow y=C{\mathrm{e}}^{-ax}$

Variation der Konstanten: $y^{\prime }=C^{\prime }{\mathrm{e}}^{-ax}-Ca{\mathrm{e}}^{-ax}$

$C^{\prime }{\mathrm{e}}^{-ax}-\cancel{Ca{\mathrm{e}}^{-ax}}+\cancel{aC{\mathrm{e}}^{-ax}}=x$

$C^{\prime }=x{\mathrm{e}}^{ax}$

Aus einer Integraltafel (oder mittels partieller Integration) findet man damit $C=\frac{{\mathrm{e}}^{ax}}{a^2}(ax-1)+D$

und somit $y=C{e}^{-ax}=(\frac{{\mathrm{e}}^{ax}}{a^2}(ax-1)+D){\mathrm{e}}^{-ax}$

Die Konstante D kann man wieder aus einer Anfangsbedingung ermitteln. Auch hier sei wieder auf die Möglichkeit der Berechnung mittels Computeralgebrasystemen hingewiesen. So ergibt die Berechnung mit dem CAS Maxima natürlich obige Formel:

Bernoullische Dgl.

Jakob Bernoulli (1654–1705) schweizer Mathematiker und Physiker

Die bernoullische Dgl. $y^{\prime }=f(x)y+g(x)y^{\alpha };\alpha \ne 0,\alpha \ne 1$ lässt sich durch Transformation von $z=y^{-\alpha +1}$ auf eine lineare Dgl. 1. Ordnung zurückführen.

Siehe auch Vorlage:W.

Riccatische Dgl.

Jacopo Francesco Riccati (1676-1754) italienischer Mathematiker

Siehe vorerst Vorlage:W.

Lineare Systeme 1. Ordnung

Homogene

${𝐲}^{\prime }-𝐀(x)𝐲=\mathrm{𝟎}$

Vorlage:W, Vorlage:W

Josef Hoëné-Wronski (1776-1853) polnischer Mathematiker

Sei ${𝐲}_1,\dots ,{𝐲}_n$ ein Fundamentalsystem von ${𝐲}^{\prime }$. Dann lässt sich $𝐲$ in der Form $𝐲=C_1{𝐲}_1+\dots +C_n{𝐲}_n$ darstellen ($C_i=\text{konst.}$).

Inhomogene

${𝐲}^{\prime }-𝐀(x)𝐲=𝐁(x)$

Sei ${𝐲}_p$ eine (partikuläre) Lösung des inhomogenen linearen Systems und bilden ${𝐲}_1,\dots ,{𝐲}_n$ ein Fundamentalsystem des zugehörigen homogenen linearen Systems. Dann ist ${𝐲}_p(x)+C_1{𝐲}_1+\dots +C_n{𝐲}_n$ eine Lösung des Systems ($C_i=\text{konst.}$).