From 7e8f10448d910fdc938383ce4ce7904a60be51c5 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Tue, 30 Nov 2021 17:25:41 +0100 Subject: add new bessel images --- buch/chapters/060-integral/Makefile.inc | 1 + buch/chapters/060-integral/legendredgl.tex | 369 ++++++++++++++++++++ buch/chapters/060-integral/orthogonal.tex | 499 +-------------------------- buch/chapters/090-pde/Makefile.inc | 12 + buch/chapters/090-pde/bessel/Makefile | 15 + buch/chapters/090-pde/bessel/besselnodes.m | 106 ++++++ buch/chapters/090-pde/bessel/besselzeros.m | 70 ++++ buch/chapters/090-pde/bessel/besselzeros.tex | 17 + buch/chapters/090-pde/bessel/pauke.pdf | Bin 0 -> 24066 bytes buch/chapters/090-pde/bessel/pauke.tex | 21 ++ buch/chapters/090-pde/chapter.tex | 30 ++ buch/chapters/090-pde/gleichung.tex | 11 + buch/chapters/090-pde/kreis.tex | 219 ++++++++++++ buch/chapters/090-pde/rechteck.tex | 7 + buch/chapters/090-pde/separation.tex | 391 +++++++++++++++++++++ buch/chapters/part1.tex | 2 +- buch/common/macros.tex | 1 + 17 files changed, 1280 insertions(+), 491 deletions(-) create mode 100644 buch/chapters/060-integral/legendredgl.tex create mode 100644 buch/chapters/090-pde/Makefile.inc create mode 100644 buch/chapters/090-pde/bessel/Makefile create mode 100644 buch/chapters/090-pde/bessel/besselnodes.m create mode 100644 buch/chapters/090-pde/bessel/besselzeros.m create mode 100644 buch/chapters/090-pde/bessel/besselzeros.tex create mode 100644 buch/chapters/090-pde/bessel/pauke.pdf create mode 100644 buch/chapters/090-pde/bessel/pauke.tex create mode 100644 buch/chapters/090-pde/chapter.tex create mode 100644 buch/chapters/090-pde/gleichung.tex create mode 100644 buch/chapters/090-pde/kreis.tex create mode 100644 buch/chapters/090-pde/rechteck.tex create mode 100644 buch/chapters/090-pde/separation.tex diff --git a/buch/chapters/060-integral/Makefile.inc b/buch/chapters/060-integral/Makefile.inc index 9cc5356..09be355 100644 --- a/buch/chapters/060-integral/Makefile.inc +++ b/buch/chapters/060-integral/Makefile.inc @@ -9,5 +9,6 @@ CHAPTERFILES = $(CHAPTERFILES) \ chapters/060-integral/differentialkoerper.tex \ chapters/060-integral/risch.tex \ chapters/060-integral/orthogonal.tex \ + chapters/060-integral/legendredgl.tex \ chapters/060-integral/gaussquadratur.tex \ chapters/060-integral/chapter.tex diff --git a/buch/chapters/060-integral/legendredgl.tex b/buch/chapters/060-integral/legendredgl.tex new file mode 100644 index 0000000..9aeac40 --- /dev/null +++ b/buch/chapters/060-integral/legendredgl.tex @@ -0,0 +1,369 @@ +% +% legendredgl.tex +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\subsection{Orthogonale Polynome und Differentialgleichungen} +Legendre hat einen ganz anderen Zugang zu den nach ihm benannten +Polynomen gefunden. +Er hat sie gefunden als die Lösungen einer speziellen Differentialgleichungen. +In diesem Abschnitt sollen diese Funktionen mit der Potenzreihen-Methode +wiedergefunden werden. +Dabei stellt sich heraus, dass diese Polynome auch Eigenfunktionen eines +selbstadjungierten Differentialgoperator sind. +Die Orthogonalität wird dann aus einer Verallgemeinerung der bekannten +Eingeschaft folgen, dass Eigenvektoren einer symmetrischen Matrix zu +verschiedenen Eigenwerten orthogonal sind. + +\subsubsection{Legendre-Differentialgleichung} +Die {\em Legendre-Differentialgleichung} ist die Differentialgleichung +\begin{equation} +(1-x^2) y'' - 2x y' + n(n+1) y = 0 +\label{buch:integral:eqn:legendre-differentialgleichung} +\end{equation} +für eine Funktion $y(x)$ auf dem Intervall $[-1,1]$. + +Sei $y(x)$ eine Lösung der Differentialgleichung +\eqref{buch:integral:eqn:legendre-differentialgleichung}. +Setzt man $y_s(x)=y(-x)$ in die Differentialgleichung ein, erhält +man +\[ +(1-x^2)y_s''(x) - 2x y'_s(x) + n(n+1)y_s(x) += +(1-x^2)y''(-x) +2x y(-x) +n(n+1)y(-x). +\] +Ersetzt man $t=-x$, dann wird daraus +\[ +(1-x^2)y''(t) -2t y(t) + n(n+1) y(t) = 0 +\] +aus der Differentialgleichung +\eqref{buch:integral:eqn:legendre-differentialgleichung}. +Insbesondere ist die gespiegelte Funktion $y_s(x)$ ebenfalls +eine Lösung der Differentialgleichung. + +Ist $y(x)$ eine Lösung der Differentialgleichung ist, dann lässt +sie sich in die Summe einer geraden und einer ungeraden Funktion +\[ +\left. +\begin{aligned} +y_g(x) &= \frac{y(x)+y(-x)}{2}\\ +y_u(x) &= \frac{y(x)-y(-x)}{2} +\end{aligned} +\quad +\right\} +\quad +\Rightarrow +\quad +y(x) = y_g(x) + y_u(x) +\] +zerlegen, die als Linearkombinationen der beiden Lösungen +$y(x)$ und $y_s(x)$ ebenfalls Lösungen der Differentialgleichung +sind. + +\subsubsection{Potenzreihenlösung} +Wir suchen eine Lösung in Form einer Potenzreihe um $x=0$ und +verwenden dazu den Ansatz +\[ +y(x) = a_0+a_1x+a_2x^2+ \dots = \sum_{k=0}^\infty a_kx^k. +\] +\begin{align*} +(1-x^2) \sum_{k=2}^\infty k(k-1)a_kx^{k-2} +-2x\sum_{k=0}^\infty ka_kx^{k-1} ++ +n(n+1)\sum_{k=0}^\infty a_kx^k +&= +0 +\\ +\sum_{k=0}^\infty (k+2)(k+1)a_{k+2}x^k +- +\sum_{k=2}^\infty k(k-1)a_kx^k +- +2\sum_{k=1}^\infty ka_kx^k ++ +n(n+1)\sum_{k=0}^\infty a_kx^k +&= +0 +\end{align*} +Die Koeffizienten zur Potenz $k$ sind daher +\begin{align} +k&=0: +& +0&= +2a_2+n(n+1)a_0 +\notag +\\ +&& +a_2&=-\frac{n(n+1)}{2}a_0 +\notag +\\ +k&=1: +& +0&= +6a_3-2a_1+n(n+1)a_1 +\notag +\\ +&& +a_3&= \frac{2-n(n+1)}{6}a_1 +\notag +\\ +k&>1: +& +0&= +(k+2)(k+1)a_{k+2} -k(k-1)a_k -2ka_k +n(n+1) a_k +\notag +\\ +&& +a_{k+2} +&= +\frac{ k(k+1)-n(n+1) }{(k+2)(k+1)} +a_k +\label{buch:integral:legendre-dgl:eqn:akrek} +\end{align} +Wenn $a_1=0$ und $a_0\ne 1$ ist, dann ist die Funktion $y(x)$ gerade, +alle ungeraden Koeffizienten verschwinden. +Ebenso verschwinden alle geraden Koeffizienten, wenn $a_0=0$ und $a_1\ne 0$. +Für jede Lösung $y(x)$ der Differentialgleichung ist +$y_g(x)$ ein Lösung mit $a_1=0$ und $y_u(x)$ eine Lösung mit $a_0=0$. +Wir können die Diskussion der Lösungen daher auf gerade oder ungerade +Lösungen einschränken. + +Gesucht ist jetzt eine Lösung in Form eines Polynoms. +In diesem Fall müssen die Koeffizienten $a_k$ ab einem +gewissen Index verschwinden. +Dies tritt nach \eqref{buch:integral:legendre-dgl:eqn:akrek} genau +dann auf, wenn der Zähler für ein $k$ verschwindet. +Folglich gibt es genau dann Polynomlösungen der Differentialgleichungen, +wenn $n$ eine natürlich Zahl ist. +Ausserdem ist die Lösung ein Polynom $\bar{P}_n(x)$ vom Grad $n$. +Das Polynom soll wieder so normiert sein, dass $\bar{P}_n(1)=1$ ist. + +Die Lösungen der Differentialgleichungen können jetzt explizit +berechnet werden. +Zunächst ist $\bar{P}_0(x)=1$ und $\bar{P}_1(x)=x$. +Für $n=2$ setzen wir zunächst $a_0=1$ und $a_1=0$ und erhalten +\[ +y(x) += +1 + \frac{0(0+1) - 2(2+1)}{(0+2)(0+1)}a_0 x^2 += +1 +-3x^2 +\qquad\text{oder}\qquad +\bar{P}_3(x) = \frac12(3x^2-1). +\] +Für $n=3$ starten wir von $a_1=1$ und $a_0=0$, was zunächst $a_2=0$ +impliziert. +Für $a_3$ finden wir +\[ +a_3=\frac{1(1+1)-3(3+1)}{(1+2)(1+1)} = -\frac53 +\qquad\Rightarrow\qquad +y(x) = x-\frac53x^3 +\qquad\Rightarrow\qquad +\bar{P}_3(x) = \frac12(5x^3-3x). +\] +Dies stimmt überein mit den früher gefundenen Ausdrücken für +die Legendre-Polynome. + +Die Potenzreihenlösung zeigt zwar, dass es für jedes $n\in\mathbb{N}$ +eine Polynomlösung $\bar{P}_n(x)$ vom Grad $n$ gibt. +Dies kann aber nicht erklären, warum die so gefundenen Polynome +orthogonal sind. + +\subsubsection{Eigenfunktionen} +Die Differentialgleichung +\eqref{buch:integral:eqn:legendre-differentialgleichung} +Kann mit dem Differentialoperator +\[ +D = \frac{d}{dx}(1-x^2)\frac{d}{dx} +\] +als +\[ +Dy + n(n+1)y = 0 +\] +geschrieben werden. +Tatsächlich ist +\[ +Dy += +\frac{d}{dx} (1-x^2) \frac{d}{dy} += +\frac{d}{dx} (1-x^2)y' += +(1-x^2)y'' -2x y'. +\] +Dies bedeutet, dass die Lösungen $\bar{P}_n(x)$ Eigenfunktionen +des Operators $D$ zum Eigenwert $n(n+1)$ sind: +\[ +D\bar{P}_n = -n(n+1) \bar{P}_n. +\] + +\subsubsection{Orthogonalität von $\bar{P}_n$ als Eigenfunktionen} +Ein Operator $A$ auf Funktionen heisst {\em selbstadjungiert}, wenn +für zwei beliebige Funktionen $f$ und $g$ gilt +\[ +\langle Af,g\rangle = \langle f,Ag\rangle +\] +gilt. +Im vorliegenden Zusammenhang möchten wir die Eigenschaft nutzen, +dass Eigenfunktionen eines selbstadjungierten Operatores zu verschiedenen +Eigenwerten orthogonal sind. +Dazu seien $Df = \lambda f$ und $Dg=\mu g$ und wir rechnen +\begin{equation*} +\begin{array}{rcccl} +\langle Df,g\rangle &=& \langle \lambda f,g\rangle &=& \lambda\langle f,g\rangle +\\ +=\langle f,Dg\rangle &=& \langle f,\mu g\rangle &=& \mu\langle f,g\rangle +\\ + & & 0 &=& (\lambda-\mu)\langle f,g\rangle +\end{array} +\end{equation*} +Da $\lambda-\mu\ne 0$ ist, muss $\langle f,g\rangle=0$ sein. + +Der Operator $D$ ist selbstadjungiert, d.~h. +für zwei beliebige zweimal stetig differenzierbare Funktion $f$ und $g$ +auf dem Intervall $[-1,1]$ gilt +\begin{align*} +\langle Df,g\rangle +&= +\int_{-1}^1 (Df)(x) g(x) \,dx +\\ +&= +\int_{-1}^1 +\biggl(\frac{d}{dx} (1-x^2)\frac{d}{dx}f(x)\biggr) g(x) +\,dx +\\ +&= +\underbrace{ +\biggl[ +\biggl((1-x^2)\frac{d}{dx}f(x)\biggr) g(x) +\biggr]_{-1}^1 +}_{\displaystyle = 0} +- +\int_{-1}^1 +\biggl((1-x^2)\frac{d}{dx}f(x)\biggr) \frac{d}{dx}g(x) +\,dx +\\ +&= +- +\int_{-1}^1 +\biggl(\frac{d}{dx}f(x)\biggr) \biggl((1-x^2)\frac{d}{dx}g(x)\biggr) +\,dx +\\ +&= +- +\underbrace{ +\biggl[ +f(x) \biggl((1-x^2)\frac{d}{dx}g(x)\biggr) +\biggr]_{-1}^1}_{\displaystyle = 0} ++ +\int_{-1}^1 +f(x) \biggl(\frac{d}{dx}(1-x^2)\frac{d}{dx}g(x)\biggr) +\,dx +\\ +&= +\langle f,Dg\rangle. +\end{align*} +Dies beweist, dass $D$ selbstadjungiert ist. +Da $\bar{P}_n$ Eigenwerte des selbstadjungierten Operators $D$ zu +den verschiedenen Eigenwerten $-n(n+1)$ sind, folgt auch, dass +die $\bar{P}_n$ orthogonale Polynome vom Grad $n$ sind, die die +gleiche Standardierdisierungsbedingung wie die Legendre-Polyonome +erfüllen, also ist $\bar{P}_n(x)=P_n(x)$. + +\subsubsection{Legendre-Funktionen zweiter Art} +Siehe Wikipedia-Artikel \url{https://de.wikipedia.org/wiki/Legendre-Polynom} + +Die Potenzreihenmethode liefert natürlich auch Lösungen der +Legendreschen Differentialgleichung, die sich nicht als Polynome +darstellen lassen. +Ist $n$ gerade, dann liefern die Anfangswerte $a_0=0$ und $a_1=1$ +eine ungerade Funktion, die Folge der Koeffizienten bricht +aber nicht ab, vielmehr ist +\begin{align*} +a_{k+2} +&= +\frac{k(k+1)}{(k+1)(k+2)}a_k += +\frac{k}{k+2}a_k. +\end{align*} +Durch wiederholte Anwendung dieser Rekursionsformel findet man +\[ +a_{k} += +\frac{k-2}{k}a_{k-2} += +\frac{k-2}{k}\frac{k-4}{k-2}a_{k-4} += +\frac{k-2}{k}\frac{k-4}{k-2}\frac{k-6}{k-4}a_{k-6} += +\dots += +\frac{1}{k}a_1. +\] +Die Lösung hat daher die Reihenentwicklung +\[ +Q_0(x) = x+\frac13x^3 + \frac15x^5 + \frac17x^7+\dots += +\frac12\log \frac{1+x}{1-x} += +\operatorname{artanh}x. +\] +Die Funktion $Q_0(x)$ heisst {\em Legendre-Funktion zweiter Art}. + +Für $n=1$ wird die Reihenentwicklung $a_0=1$ und $a_1=0$ etwas +interessanter. +Die Rekursionsformel für die Koeffizienten ist +\[ +a_{k+2} += +\frac{k(k+1)-2}{(k+1)(k+2)} a_k. +\qquad\text{oder}\qquad +a_k += +\frac{(k-1)(k-2)-2}{k(k-1)} +a_{k-2} +\] +Man erhält der Reihe nach +\begin{align*} +a_2 &= \frac{-2}{2\cdot 1} a_0 = -1 +\\ +a_3 &= 0 +\\ +a_4 &= \frac{3\cdot 2-2}{4\cdot 3} a_2 = \frac{4}{4\cdot 3}a_2 = \frac13a_2 = -\frac13 +\\ +a_5 &= 0 +\\ +a_6 &= \frac{5\cdot 4-2}{6\cdot 5}a_4 = \frac{18}{6\cdot 5}a_4 = -\frac15 +\\ +a_7 &= 0 +\\ +a_8 &= \frac{7\cdot 6-2}{8\cdot 7}a_6 = \frac{40}{8\cdot 7} = -\frac17 +\\ +a_9 &= 0 +a_{10} &= \frac{9\cdot 8-2}{10\cdot 9}a_8 = \frac{70}{10\cdot 9} = -\frac19, +\end{align*} +woraus sich die Reihenentwicklung +\begin{align*} +y(x) +&= +-x^2 -\frac13x^4 -\frac15x^6 - \frac17x^8 -\frac19x^{10}-\dots +\\ +&= +-x\biggl(x+\frac13x^3 + \frac15x^5 + \frac17x^7 + \frac19x^9+\dots\biggr) += +-x\operatorname{artanh}x. +\end{align*} +Die {\em Legendre-Funktionen zweiter Art} $Q_n(x)$ werden allerdings +so definiert, dass gewisse Rekursionsformeln für die Legendre-Polynome, +die wir hier nicht hergeleitet haben, auch für die $Q_n(x)$ gelten. +In dieser Normierung muss statt des eben berechneten $y(x)$ die Funktion +\[ +Q_1(x) = x \operatorname{artanh}x-1 +\] +verwendet werden. + + + + + + diff --git a/buch/chapters/060-integral/orthogonal.tex b/buch/chapters/060-integral/orthogonal.tex index ceba53a..109cd61 100644 --- a/buch/chapters/060-integral/orthogonal.tex +++ b/buch/chapters/060-integral/orthogonal.tex @@ -508,496 +508,15 @@ Durch weitere Durchführung des Verfahrens liefert die Polynome in Tabelle~\ref{buch:integral:table:legendre-polynome}. -% -% Differentialgleichungen -% -\subsection{Orthogonale Polynome und Differentialgleichungen} -\subsubsection{Legendre-Differentialgleichung} -\subsubsection{Legendre-Polyome} -\subsubsection{Legendre-Funktionen zweiter Art} -Siehe Wikipedia-Artikel \url{https://de.wikipedia.org/wiki/Legendre-Polynom} - %% -%% Anwendung: Gauss-Quadratur +%% Differentialgleichungen %% -%\subsection{Anwendung: Gauss-Quadratur} -%Orthogonale Polynome haben eine etwas unerwartet Anwendung in einem -%von Gauss erdachten numerischen Integrationsverfahren. -%Es basiert auf der Beobachtung, dass viele Funktionen sich sehr -%gut durch Polynome approximieren lassen. -%Wenn man also sicherstellt, dass ein Verfahren für Polynome -%sehr gut funktioniert, darf man auch davon ausgehen, dass es für -%andere Funktionen nicht allzu schlecht sein wird. -% -%\subsubsection{Interpolationspolynome} -%Zu einer stetigen Funktion $f(x)$ auf dem Intervall $[-1,1]$ -%ist ein Polynome vom Grad $n$ gesucht, welches in den Punkten -%$x_0{$}c<{$}|>{$}r<{$}|>{$}r<{$}|} -%\hline -% n & \text{Gauss-Quadratur} & \text{Trapezregel} \\ -%\hline -%\phantom{0}2 & 0.\u{95}74271077563381 & 0.\u{95}63709682242596 \\ -%\phantom{0}4 & 0.\u{95661}28333449730 & 0.\u{956}5513401768598 \\ -%\phantom{0}6 & 0.\u{9566114}812034364 & 0.\u{956}5847489712136 \\ -%\phantom{0}8 & 0.\u{956611477}5028123 & 0.\u{956}5964425360520 \\ -% 10 & 0.\u{9566114774905}637 & 0.\u{9566}018550715587 \\ -% 12 & 0.\u{956611477490518}7 & 0.\u{9566}047952369826 \\ -% 14 & 0.\u{95661147749051}72 & 0.\u{9566}065680717177 \\ -% 16 & 0.\u{956611477490518}7 & 0.\u{9566}077187127541 \\ -% 18 & 0.\u{956611477490518}3 & 0.\u{9566}085075898731 \\ -% 20 & 0.\u{956611477490518}4 & 0.\u{9566}090718697414 \\ -%\hline -% \infty & 0.9566114774905183 & 0.9566114774905183 \\ -%\hline -%\end{tabular} -%\caption{Integral von $\sqrt{1-x^2}$ zwischen $-\frac12$ und $\frac12$ -%berechnet mit Gauss-Quadratur und der Trapezregel, aber mit zehnmal -%so vielen Stützstellen. -%Bereits mit 12 Stützstellen erreicht die Gauss-Quadratur -%Maschinengenauigkeit, die Trapezregel liefert auch mit 200 Stützstellen -%nicht mehr als 4 korrekte Nachkommastellen. -%\label{buch:integral:gaussquadratur:table0.5}} -%\end{table} -% -%%\begin{table} -%%\def\u#1{\underline{#1}} -%%\centering -%%\begin{tabular}{|>{$}c<{$}|>{$}r<{$}|>{$}r<{$}|} -%%\hline -%% n & \text{Gauss-Quadratur} & \text{Trapezregel} \\ -%%\hline -%%\phantom{0}2 & 1.\u{5}379206741571556 & 1.\u{5}093105464758343 \\ -%%\phantom{0}4 & 1.\u{51}32373472933831 & 1.\u{51}13754509594814 \\ -%%\phantom{0}6 & 1.\u{512}1624557410367 & 1.\u{51}17610879524799 \\ -%%\phantom{0}8 & 1.\u{51207}93479994321 & 1.\u{51}18963282632112 \\ -%% 10 & 1.\u{51207}13859966004 & 1.\u{51}19589735776959 \\ -%% 12 & 1.\u{512070}5317779943 & 1.\u{51}19930161260693 \\ -%% 14 & 1.\u{5120704}334802813 & 1.\u{5120}135471596636 \\ -%% 16 & 1.\u{5120704}216176006 & 1.\u{5120}268743889558 \\ -%% 18 & 1.\u{5120704}201359081 & 1.\u{5120}360123137213 \\ -%% 20 & 1.\u{5120704199}459651 & 1.\u{5120}425490275837 \\ -%%\hline -%% \infty & 1.5120704199172947 & 1.5120704199172947 \\ -%%\hline -%%\end{tabular} -%%\end{table} -% -%%\begin{table} -%%\def\u#1{\underline{#1}} -%%\centering -%%\begin{tabular}{|>{$}c<{$}|>{$}r<{$}|>{$}r<{$}|} -%%\hline -%% n & \text{Gauss-Quadratur} & \text{Trapezregel} \\ -%%\hline -%%\phantom{0}2 & 1.\u{}6246862220133462 & 1.\u{5}597986803933712 \\ -%%\phantom{0}4 & 1.\u{5}759105515463101 & 1.\u{56}63563456168101 \\ -%%\phantom{0}6 & 1.\u{5}706630058381434 & 1.\u{56}77252866190838 \\ -%%\phantom{0}8 & 1.\u{56}94851106536780 & 1.\u{568}2298707696152 \\ -%% 10 & 1.\u{56}91283195332679 & 1.\u{568}4701957758742 \\ -%% 12 & 1.\u{56}90013806299465 & 1.\u{568}6030805941198 \\ -%% 14 & 1.\u{5689}515434853885 & 1.\u{568}6841603070025 \\ -%% 16 & 1.\u{5689}306507843050 & 1.\u{568}7372230731711 \\ -%% 18 & 1.\u{5689}214761291217 & 1.\u{568}7738235496322 \\ -%% 20 & 1.\u{56891}73062385982 & 1.\u{568}8001228530786 \\ -%%\hline -%% \infty & 1.5689135396691616 & 1.5689135396691616 \\ -%%\hline -%%\end{tabular} -%%\end{table} -% -%\begin{table} -%\def\u#1{\underline{#1}} -%\centering -%\begin{tabular}{|>{$}c<{$}|>{$}r<{$}|>{$}r<{$}|} -%\hline -% n & \text{Gauss-Quadratur} & \text{Trapezregel} \\ -%\hline -%\phantom{0}2 & 1.\u{}6321752373234928 & 1.\u{5}561048774629949 \\ -%\phantom{0}4 & 1.\u{57}98691557134743 & 1.\u{5}660124134617943 \\ -%\phantom{0}6 & 1.\u{57}35853681692993 & 1.\u{5}683353001877542 \\ -%\phantom{0}8 & 1.\u{57}19413565928206 & 1.\u{5}692627503425400 \\ -% 10 & 1.\u{57}13388119633434 & 1.\u{5}697323578543481 \\ -% 12 & 1.\u{57}10710489948883 & 1.\u{570}0051217458713 \\ -% 14 & 1.\u{570}9362135398341 & 1.\u{570}1784766276063 \\ -% 16 & 1.\u{570}8621102742815 & 1.\u{570}2959121005231 \\ -% 18 & 1.\u{570}8186779483588 & 1.\u{570}3793521168343 \\ -% 20 & 1.\u{5707}919411931615 & 1.\u{570}4408749735932 \\ -%\hline -% \infty & 1.5707367072605671 & 1.5707367072605671 \\ -%\hline -%\end{tabular} -%\caption{Integral von $\sqrt{1-x^2}$ zwischen $-0.999$ und $0.999$ -%berechnet mit Gauss-Quadratur und der Trapezregel, aber mit zehnmal -%so vielen Stützstellen. -%Wegen der divergierenden Steigung des Integranden bei $\pm 1$ tun -%sich beide Verfahren sehr schwer. -%Trotzdem erreich die Gauss-Quadrator 4 korrekte Nachkommastellen -%mit 20 Stütztstellen, während die Trapezregel auch mit 200 Stützstellen -%nur 3 korrekte Nachkommastellen findet. -%\label{buch:integral:gaussquadratur:table0.999}} -%\end{table} -% -%\begin{figure} -%\centering -%\includegraphics{chapters/060-integral/gq/gq.pdf} -%\caption{Approximationsfehler des -%Integrals~\eqref{buch:integral:gaussquadratur:bspintegral} -%in Abhängigkeit von $a$. -%Die Divergenz der Ableitung des Integranden an den Intervallenden -%$\pm 1$ führt zu schlechter Konvergenz des Verfahrens, wenn $a$ -%nahe an $1$ ist. -%\label{buch:integral:gaussquadratur:fehler}} -%\end{figure} -% -%Zur Illustration der Genauigkeit der Gauss-Quadratur berechnen wir -%das Integral -%\begin{equation} -%\int_{-a}^a \sqrt{1-x^2}\,dx -%= -%\arcsin a + a \sqrt{1-a^2} -%\label{buch:integral:gaussquadratur:bspintegral} -%\end{equation} -%mit Gauss-Quadratur einerseits und dem Trapezverfahren -%andererseits. -%Da Gauss-Quadratur mit sehr viel weniger Sützstellen auskommt, -%berechnen wir die Trapeznäherung mit zehnmal so vielen Stützstelln. -%In den Tabellen~\ref{buch:integral:gaussquadratur:table0.5} -%und -%\ref{buch:integral:gaussquadratur:table0.999} -%sind die Resultate zusammengestellt. -%Für $a =\frac12$ zeigt -%Tabelle~\ref{buch:integral:gaussquadratur:table0.5} -%die sehr schnelle Konvergenz der Gauss-Quadratur, schon mit -%12 Stützstellen wird Maschinengenauigkeit erreicht. -%Das Trapezverfahren dagegen erreicht auch mit 200 Stützstellen nur -%4 korrekte Nachkommastellen. -% -%An den Stellen $x=\pm 1$ divergiert die Ableitung des Integranden -%des Integrals \eqref{buch:integral:gaussquadratur:bspintegral}. -%Da grösste und kleinste Stützstelle der Gauss-Quadratur immer -%deutlich vom Rand des Intervalls entfernt ist, kann das Verfahren -%diese ``schwierigen'' Stellen nicht erkennen. -%Tabelle~\ref{buch:integral:gaussquadratur:table0.999} zeigt, wie -%die Konvergenz des Verfahrens in diesem Fall sehr viel schlechter ist. -%Dies zeigt auch der Graph in -%Abbildung~\ref{buch:integral:gaussquadratur:fehler}. -% -%\subsubsection{Skalarprodukte mit Gewichtsfunktion} +%\subsection{Orthogonale Polynome und Differentialgleichungen} +%\subsubsection{Legendre-Differentialgleichung} +%\subsubsection{Legendre-Polyome} +%\subsubsection{Legendre-Funktionen zweiter Art} +%Siehe Wikipedia-Artikel \url{https://de.wikipedia.org/wiki/Legendre-Polynom} + +\input{chapters/060-integral/legendredgl.tex} + \input{chapters/060-integral/gaussquadratur.tex} diff --git a/buch/chapters/090-pde/Makefile.inc b/buch/chapters/090-pde/Makefile.inc new file mode 100644 index 0000000..191bad6 --- /dev/null +++ b/buch/chapters/090-pde/Makefile.inc @@ -0,0 +1,12 @@ +# +# Makefile.inc -- Makefile dependencies for chapter 9 +# +# (c) 2022 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +# + +CHAPTERFILES = $(CHAPTERFILES) \ + chapters/090-pde/gleichung.tex \ + chapters/090-pde/separation.tex \ + chapters/090-pde/rechteck.tex \ + chapters/090-pde/kreis.tex \ + chapters/090-pde/chapter.tex diff --git a/buch/chapters/090-pde/bessel/Makefile b/buch/chapters/090-pde/bessel/Makefile new file mode 100644 index 0000000..c189517 --- /dev/null +++ b/buch/chapters/090-pde/bessel/Makefile @@ -0,0 +1,15 @@ +# +# Makefile +# +# (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +# +all: besselzeros.tex besselnodes.tex pauke.pdf + +besselzeros.tex: besselzeros.m + octave besselzeros.m + +besselnodes.tex: besselnodes.m + octave besselnodes.m + +pauke.pdf: pauke.tex besselnodes.tex + pdflatex pauke.tex diff --git a/buch/chapters/090-pde/bessel/besselnodes.m b/buch/chapters/090-pde/bessel/besselnodes.m new file mode 100644 index 0000000..0dcba3e --- /dev/null +++ b/buch/chapters/090-pde/bessel/besselnodes.m @@ -0,0 +1,106 @@ +# +# besselnodes.m +# +# (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +# +global maxmu; +maxmu = 3; +global maxk; +maxk = 4; +global mu; + +nachkommastellen = 4; + +function retval = f(x) + global mu; + retval = besselj(mu, x); +endfunction + +global kzeros; +kzeros = zeros(maxk+1, maxmu+1); +for mu = (0:maxmu) + k = 0; + x = 0.0001; + while (k <= maxk) + bracket = [ x, x+1 ]; + if (f(bracket(1)) * f(bracket(2)) < 0) + kzeros(k+1,mu+1) = fzero("f", bracket); + k = k + 1; + endif + x = x + 1; + endwhile +endfor + +xshift = 4; +yshift = 4; +global r; +r = 1.8; + +function retval = anderefarbe(f) + if (1 == strcmp("red", f)) + retval = "blue"; + else + retval = "red"; + endif +endfunction + +function sektor(fn, mu, k, w0, w1, startfarbe) + global kzeros; + global r; + fprintf(fn, "\\begin{scope}\n"); + fprintf(fn, "\\clip (0,0)--(%.4f:%.4f) arc (%.4f:%.4f:%.4f)--cycle;\n", + w0, r, w0, w1, r); + faktor = kzeros(k+1,mu+1); + + K = k + 1; + farbe = startfarbe; + while (K > 0) + R = r * kzeros(K, mu+1) / faktor; + fprintf(fn, "\\fill[color=%s!20] ", farbe); + fprintf(fn, "(0,0) circle[radius=%.4f];\n", R); + farbe = anderefarbe(farbe); + K = K-1; + end + fprintf(fn, "\\end{scope}\n"); +endfunction + +fn = fopen("besselnodes.tex", "w"); + +#fprintf(fn, "\\begin{tikzpicture}[>=latex,thick]\n"); + +for mu = (0:maxmu) + if (mu > 0) + winkel = 180 / mu; + else + winkel = 360; + endif + for k = (0:maxk) + fprintf(fn, "\\begin{scope}[xshift=%.3fcm,yshift=-%.3fcm]\n", + mu * xshift, k * yshift); + for w0 = (0:2*winkel:360) + sektor(fn, mu, k, w0, w0 + winkel, "red"); + if (winkel < 270) + sektor(fn, mu, k, w0 + winkel, w0 + 2 * winkel, "blue"); + endif + endfor + + fprintf(fn, "\\draw (0,0) circle[radius=%.4f];\n", r); + + fprintf(fn, "\\end{scope}\n\n"); + endfor +endfor + +for mu = (0:maxmu) + fprintf(fn, "\\node at (%.4f,%.4f) [above] {$\\mu=%d$};\n", + mu * xshift, 0.5 * yshift, mu); +endfor + +for k = (0:maxk) + fprintf(fn, "\\node at (%.4f,%.4f) [above,rotate=90] {$k=%d$};\n", + -0.5 * xshift, -k * yshift, k); +endfor + +#fprintf(fn, "\\end{tikzpicture}\n"); + +fclose(fn); + diff --git a/buch/chapters/090-pde/bessel/besselzeros.m b/buch/chapters/090-pde/bessel/besselzeros.m new file mode 100644 index 0000000..9c8fa9d --- /dev/null +++ b/buch/chapters/090-pde/bessel/besselzeros.m @@ -0,0 +1,70 @@ +# +# besselzeros.m -- find zeros of bessel functions +# +# (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +#  +global maxmu; +maxmu = 7; +global maxk; +maxk = 10; +global mu; + +nachkommastellen = 4; + +function retval = f(x) + global mu; + retval = besselj(mu, x); +endfunction + +kzeros = zeros(maxk+1, maxmu+1); +for mu = (0:maxmu) + k = 0; + if (mu > 0) + kzeros(1, mu+1) = 0; + k = k+1; + endif + x = 0.0001; + while (k <= maxk) + bracket = [ x, x+1 ]; + if (f(bracket(1)) * f(bracket(2)) < 0) + kzeros(k+1,mu+1) = fzero("f", bracket); + k = k + 1; + endif + x = x + 1; + endwhile +endfor + +# kzeros + +fn = fopen("besselzeros.tex", "w"); + +fprintf(fn, "\\begin{tabular}{|>{$}c<{$}"); +for mu = (0:maxmu) + fprintf(fn, "|>{$}r<{$}"); +endfor +fprintf(fn, "|}\n"); + +fprintf(fn, "\\hline\n"); +fprintf(fn, " k "); +for mu = (0:maxmu) + fprintf(fn, "& \\mu = %d ", mu); +endfor +fprintf(fn, "\\\\\n"); +fprintf(fn, "\\hline\n"); + +for k = (0:maxk) + fprintf(fn, " %d ", k); + for mu = (0:maxmu) + value = kzeros(k+1, mu+1); + if (value == 0) + fprintf(fn, "& 0\\phantom{.%0*d}", nachkommastellen, 0); + else + fprintf(fn, "& %*.*f", nachkommastellen+4, nachkommastellen, kzeros(k+1, mu+1)); + endif + endfor + fprintf(fn, "\\\\\n"); +endfor +fprintf(fn, "\\hline\n"); +fprintf(fn, "\\end{tabular}\n"); + +fclose(fn); diff --git a/buch/chapters/090-pde/bessel/besselzeros.tex b/buch/chapters/090-pde/bessel/besselzeros.tex new file mode 100644 index 0000000..1b8d33b --- /dev/null +++ b/buch/chapters/090-pde/bessel/besselzeros.tex @@ -0,0 +1,17 @@ +\begin{tabular}{|>{$}c<{$}|>{$}r<{$}|>{$}r<{$}|>{$}r<{$}|>{$}r<{$}|>{$}r<{$}|>{$}r<{$}|>{$}r<{$}|>{$}r<{$}|} +\hline + k & \mu = 0 & \mu = 1 & \mu = 2 & \mu = 3 & \mu = 4 & \mu = 5 & \mu = 6 & \mu = 7 \\ +\hline + 0 & 2.4048& 0\phantom{.0000}& 0\phantom{.0000}& 0\phantom{.0000}& 0\phantom{.0000}& 0\phantom{.0000}& 0\phantom{.0000}& 0\phantom{.0000}\\ + 1 & 5.5201& 3.8317& 5.1356& 6.3802& 7.5883& 8.7715& 9.9361& 11.0864\\ + 2 & 8.6537& 7.0156& 8.4172& 9.7610& 11.0647& 12.3386& 13.5893& 14.8213\\ + 3 & 11.7915& 10.1735& 11.6198& 13.0152& 14.3725& 15.7002& 17.0038& 18.2876\\ + 4 & 14.9309& 13.3237& 14.7960& 16.2235& 17.6160& 18.9801& 20.3208& 21.6415\\ + 5 & 18.0711& 16.4706& 17.9598& 19.4094& 20.8269& 22.2178& 23.5861& 24.9349\\ + 6 & 21.2116& 19.6159& 21.1170& 22.5827& 24.0190& 25.4303& 26.8202& 28.1912\\ + 7 & 24.3525& 22.7601& 24.2701& 25.7482& 27.1991& 28.6266& 30.0337& 31.4228\\ + 8 & 27.4935& 25.9037& 27.4206& 28.9084& 30.3710& 31.8117& 33.2330& 34.6371\\ + 9 & 30.6346& 29.0468& 30.5692& 32.0649& 33.5371& 34.9888& 36.4220& 37.8387\\ + 10 & 33.7758& 32.1897& 33.7165& 35.2187& 36.6990& 38.1599& 39.6032& 41.0308\\ +\hline +\end{tabular} diff --git a/buch/chapters/090-pde/bessel/pauke.pdf b/buch/chapters/090-pde/bessel/pauke.pdf new file mode 100644 index 0000000..54edc20 Binary files /dev/null and b/buch/chapters/090-pde/bessel/pauke.pdf differ diff --git a/buch/chapters/090-pde/bessel/pauke.tex b/buch/chapters/090-pde/bessel/pauke.tex new file mode 100644 index 0000000..bba092e --- /dev/null +++ b/buch/chapters/090-pde/bessel/pauke.tex @@ -0,0 +1,21 @@ +% +% pauke.tex -- template for standalon tikz images +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\documentclass[tikz]{standalone} +\usepackage{amsmath} +\usepackage{times} +\usepackage{txfonts} +\usepackage{pgfplots} +\usepackage{csvsimple} +\usetikzlibrary{arrows,intersections,math} +\begin{document} +\def\skala{1} +\begin{tikzpicture}[>=latex,thick,scale=\skala] + +\input{besselnodes.tex} + +\end{tikzpicture} +\end{document} + diff --git a/buch/chapters/090-pde/chapter.tex b/buch/chapters/090-pde/chapter.tex new file mode 100644 index 0000000..543a92d --- /dev/null +++ b/buch/chapters/090-pde/chapter.tex @@ -0,0 +1,30 @@ +% +% chapter.tex -- Kapitel zu partiellen Differentialgleichungen +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +% !TeX spellcheck = de_CH +\chapter{Partielle Differentialgleichungen +\label{buch:chapter:pde}} +\lhead{Partielle Differentialgleichungen} +\rhead{} +Partielle Differentialgleichungen sind eine besonders ergiebige +Quelle für Anwendungen spezieller Funktionen. +Die Separationsmethode zum Beispiel für die Wellengleichung +auf gewissen, besonders einfachen Gebieten wie Rechtecken, +Kreisscheiben oder Kugel führt auf gewöhnliche Differentialgleichungen, +deren Lösungen spezielle Funktionen sind. + +\input{chapters/090-pde/gleichung.tex} +\input{chapters/090-pde/separation.tex} +\input{chapters/090-pde/rechteck.tex} +\input{chapters/090-pde/kreis.tex} + +%\section*{Übungsaufgaben} +%\rhead{Übungsaufgaben} +%\aufgabetoplevel{chapters/020-exponential/uebungsaufgaben} +%\begin{uebungsaufgaben} +%\uebungsaufgabe{0} +%\uebungsaufgabe{1} +%\end{uebungsaufgaben} + diff --git a/buch/chapters/090-pde/gleichung.tex b/buch/chapters/090-pde/gleichung.tex new file mode 100644 index 0000000..07dd2ff --- /dev/null +++ b/buch/chapters/090-pde/gleichung.tex @@ -0,0 +1,11 @@ +% +% gleichung.tex +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\section{Gleichungen und Randbedingungen +\label{buch:pde:section:gleichungen-und-randbedingungen}} + +\subsection{Laplace-Operator} + +\subsection{Orthogonalität} diff --git a/buch/chapters/090-pde/kreis.tex b/buch/chapters/090-pde/kreis.tex new file mode 100644 index 0000000..a54ce38 --- /dev/null +++ b/buch/chapters/090-pde/kreis.tex @@ -0,0 +1,219 @@ +% +% kreis.tex +% +% (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +% +\section{Kreisförmige Membran +\label{buch:pde:section:kreis}} +In diesem Abschnitt soll die Differentialgleichung einer kreisförmigen +Membran mit Hilfe der Separationsmethode gelöst werden. +Dabei werden die Bessel-Funktionen als Lösungsfunktionen +auftreten und die Eigenfrequenzen werden durch ihre Nullstellen +berechnet. + +\subsection{Differentialgleichung und Randbedingung} +Die Wellengleichung auf einem Kreisgebiet mit Radius $r_0$ +lässt sich am besten mit Hilfe von Polarkoordinaten $(r,\varphi)$ +ausdrücken. +Gesucht ist also eine Funktion $u(t,r,\varphi)$ gesucht, wobei +$0\le r0$. +Dies bedeutet für die Funktion $U(r,\varphi)$, dass +$U(r_0,\varphi)=0$ sein muss für alle $\varphi$. + +Die Bedingungen an $U$ reichen aber nicht ganz. +Alle Koordinaten $(0,\varphi)$ bezeichnen ja gleichermassen +den Nullpunkt des Koordinatensystems, es muss also auch sichergestellt +sein, dass $U(0,\varphi)$ für alle $\varphi$ den gleichen Wert gibt. + +\subsection{Separation} +Das Eigenwertproblem $\Delta U=-\lambda^2 U$ soll jetzt in Polarkoordinaten +separiert werden. +Dazu schreiben wir die Lösung als +\[ +U(r,\varphi) += +R(r)\cdot \Phi(\varphi). +\] +Die Randbedingungen an $U$ werden zu $R(r_0)=0$. + +Im Ursprung des Koordinatensystems ist die Randbedingung etwas +komplizierter. +Wenn $R(0)=0$ ist, dann ist sichergestellt, dass +$U(0,\varphi)=R(0)\Phi(\varphi)0$ ist, dass also der Wert unabhängig +ist von $\varphi$. +Wenn aber $R(0)\ne 0$ ist, dann kann die geforderte Unabhängigkeit +von $\varphi$ nur erfüllt werden, wenn $\Phi(\varphi)$ konstant ist. +Da die Funktion aber auch noch differenzierbar sein soll, darf es +an der Stelle $r=0$ keine ``Spitze'' geben, die Ableitung $R'(0)$ +muss also auch $=0$ sein. +% XXX Evtl Bild zur Illustration dieses Problems + +Die Differntialgleichungen wird mit der Form~\eqref{buch:pde:kreis:laplace} +des Laplace-Operators +\[ +\Delta U += +R''(r) \Phi(\varphi) ++ +\frac1r R'(r)\Phi(\varphi) ++ +\frac{1}{r^2} R(r)\Phi''(\varphi) += +-\lambda^2 +R(r)\Phi(\varphi) +\] +Nach Division durch die rechte Seite erhalten wir +\[ +\frac{R''(r)}{R(r)} ++ +\frac1r \frac{R'(r)}{R(r)} ++ +\frac{1}{r^2} \frac{\Phi''(\varphi)}{\Phi(\varphi)} += +-\lambda^2 +\] +Im letzten Term auf der linken Seite kommen die Variablen $r$ und $\varphi$ +gemischt vor, man muss also die Gleichung erst mit $r^2$ multiplizieren, +bevor man sie in +\[ +\frac{r^2R''(r)+rR'(r)+\lambda^2 r^2R(r)}{R(r)} += +-\frac{\Phi''(\varphi)}{\Phi(\varphi)} +\] +separieren kann. +Die beiden Seiten sind also konstant, wir nennen die gemeinsame +Konstanten $\mu^2$, das vereinfacht die Lösung der Gleichung +für $\Phi(\varphi)$. + +Die Gleichung für $\Phi$ hat für $\mu\ne 0$ die Lösungen +\begin{align*} +\Phi(\varphi) &= \cos\mu\varphi +\text{und}\qquad +\Phi(\varphi) &= \sin\mu\varphi. +\end{align*} +Die Lösung muss aber auch stetig sein, d.~h.~es muss $\Phi(0)=\Phi(2\pi)$ +gelten. +Dies ist nur möglich, wenn $\mu$ eine ganze Zahl ist. + +Für $\mu=0$ hat das charakteristische Polynome eine doppelte Nullstelle, +die allgemeine Lösung lautet daher +\[ +\Phi(\varphi)= C \varphi + D. +\] +Die Funktion $\Phi$ muss aber auch stetig sein, d.~h.~$\Phi(0)=\Phi(2\pi)$, +das ist mit $C\ne 0$ nicht möglich, somit kommt für $\mu=0$ nur die +Lösung $\Phi(\varphi)=D$ in Frage. + +Die Gleichung für $R(r)$ wird jetzt +\begin{equation} +r^2R''(r) + rR'(r)+(\lambda^2 r^2-\mu^2)R(r) += +0. +\label{buch:pde:kreis:Rdgl} +\end{equation} +Bis auf den Faktor $\lambda^2$ ist dies eine Besselsche Differentialgleichung. + +\subsection{Umformung in eine Besselsche Differentialgleichung} +Die Funktion $y(x) = J_\mu(sx)$ hat die Ableitungen +\begin{align*} +y'(x) &= sJ'_mu(sx) +\\ +y''(x) &= s^2J''_\mu(sx) +\end{align*} +Setzt man dies in die Besselsche Differentialgleichung für $J_\mu$ an +der Stelle $sx$ ein, erhält man +\[ +s^2x^2 J''_\mu(sx) + sx J'_\mu(sx) + (s^2x^2 -\mu^2) J_\mu(sx) = 0. +\] +Die Differentialgleichung \eqref{buch:pde:kreis:Rdgl} der Funktion $R(r)$ +wird also gelöst von den Funktionen $R(r) = J_\mu(\lambda r)$. + +\subsection{Eigenfrequenzen} +Im vorangegangenen Abschnitt haben wir gefunden, dass die Lösungen +für $R(r)$ die Funktionen $J_\mu(\lambda r)$ sind. +Bis jetzt haben wir aber nicht nachgeprüft, dass die Randbedingung +eingehalten wird. +Diese ist erfüllt, $R(r_0)=0$ ist. +Es muss also +$J_\mu(\lambda r_0)=0$ sein, oder $\lambda r_0$ muss eine +Nullstelle von $J_{\mu}$ sein. +Bezeichnen wir die Nullstellen von $J_\mu$ mit $j_{\mu k}$, wobei $k$ +eine natürliche Zahl ist, dann muss +\[ +\lambda = \frac{j_{\mu k}}{r_0} +\] +sein. +Die Eigenfrequenzen der kreisförmigen Membran werden also im Wesentlichen +durch die Nullstellen der Bessel-Funktionen gegeben. + +Zu jedem ganzzahligen $\mu$ gibt es also eine Folge $j_{\mu k}/r_0$ von +Eigenfrequenzen. +Die Lösungen mit Index $k$ der Differentialgleichung mit Index $k$ hat +die Form +\[ +U_{\mu k}(r,\varphi) += +C \cos(\mu \varphi+\delta) +J_{\mu}\biggl( +\frac{j_{\mu k}}{r_0}r +\biggr) +\] +Der Faktor $J_{\mu}$ hat $k$ weitere Nullstellen für Radien $r0$ in einem +Punkt $x\in G$ des Gebietes $G$ ist. +Die Randbedingungen zerfallen in zwei Teile: +\begin{itemize} +\item +Bedingungen, die wiedergeben, dass die Membran in einen +Rahmen eingespannt und damit unbeweglich ist. +Dies bedeutet, dass $u(t,x)=0$ für alle Zeiten $t>0$ und für +Randpunkte $x\in\partial G$ von $G$ ist. +\item +Bedingungen, die Auslenkung und Geschwindigkeit der Membran zur +Zeit $t=0$ beschreiben, typischerweise ind er Form +\begin{align*} +u(0,x) = f(x), +\frac{\partial u}{\partial t}(0,x) = g(x) +\end{align*} +wobei $f(x)$ und $g(x)$ Funktionen auf dem Gebiet $G$ sind. +\end{itemize} + +In der Zeitableitung auf der linken Seite +von~\eqref{buch:pde:separation:wellengleichung} +kommen die Ortskoordinaten nicht vor und im Laplace-Operator +auf der rechten Seite tritt die Zeit nicht auf. +Es ist daher naheliegend zu versuchen, die Lösung der Differntialgleichung +als Produkt +\[ +u(t,x) = T(t) \cdot U(x) +\] +zu schreiben. +Wendet man die Differentialgleichung darauf an, wird daraus die Gleichung +\[ +\frac{1}{c^2} +T''(t)\cdot U(x) += +T(t) \cdot \Delta U(x). +\] +Indem man druch $T(t)$ und $U(x)$ teilt, entsteht die separierte Gleichung +\[ +\frac{1}{c^2} \frac{T''(t)}{T(t)} += +\frac{\Delta U(x)}{U(x)}. +\] +Die linke Seite ist nur von der Zeit abhängig, die rechte nur von den +Ortskoordinaten. +Damit ist die Differentialgleichung separiert und das Problem darauf +reduziert, die gewöhnliche Differentialgleichung +\[ +T''(t) = \lambda T(t) +\] +und die partielle Differentialgleichung +\[ +\Delta U(x) = \lambda U(x) +\] +niedrigerer Dimension zu lösen. + +\subsubsection{Allgemeine Situation} +Das Definitionsgebiet der partiellen Differentialgleichung ist +also von der Form $\mathbb{R}^+\times G$, wobei $G\subset\mathbb{R}^n$ +ein räumliches Gebiet ist und $\mathbb{R}^+$ die Zeitachse. +Auch die Randbedingungen zerfallen in zwei Arten: +\begin{itemize} +\item +Bedingungen über die Lösungsfunktion zur Zeit $t=0$ im inneren des +räumliche Gebietes $G$, zum Beispiel +die Anfangsauslenkung und/oder Anfangsgeschwindigkeit einer schwingenden +Saite oder Membran. +\item +Bedingungen über die Lösungsfunktion auf dem Rand $\partial G$ von +$G$ für alle Zeiten $t>0$, zum Beispiel die Bedingung, dass die +Membran fest eingespannt ist. +\end{itemize} +Oft zerfällt auch der Differentialoperator in Zeitableitungen +und einen zeitunabhängigen Teil der nur Ableitungen nach den +Ortsvariablen enthält. +Die Wellengleichung +\[ +\frac{1}{c^2} +\frac{\partial^2}{\partial t^2} u += +\Delta u +\qquad\Leftrightarrow\qquad +\biggl( +\frac{1}{c^2}\frac{\partial^2}{\partial t^2} - \Delta +\biggr) u = 0 +\] +enthält Ableitungen nach der Zeit, die nicht von Ortskoordinaten +abhängig sind. +Der Laplace-Operator $\Delta$ ist nicht von der Zeitabhängig und das +Gebiet $G$ hängt ebenfalls nicht von der Zeit ab. + +\subsubsection{Separation der Zeit} +Unter den gegeben Voraussetzungen ist es naheliegend, die Lösungsfunktion +$u(t,x)$ als Produkt +\[ +u(t,x) = T(t) \cdot U(x),\qquad t\in\mathbb{R}^+, x\in G +\] +anzusetezen. +Die Wellengleichung wird dann +\[ +\frac{1}{c^2} +T''(t)\cdot U(x) += +T(t)\cdot\Delta U(x) +\] +und nach Separation +\[ +\frac{1}{c^2} \frac{T''(t)}{T(t)} += +\frac{\Delta U(x)}{U(x)}. +\] +Es gibt also eine gemeinsame Konstante. +Da wir Schwingungslösungen erwarten, für die $T''(t) = -\omega^2 T(t)$ +ist, schreiben wir die gemeinsame Konstante als $-\lambda^2$, was +später die Formeln vereinfachen wird. +Die separierten Differentialgleichungen werden jetzt +\begin{align*} +\frac{1}{c^2} +\frac{T''(t)}{T(t)} +&= +-\lambda^2 +&&\Rightarrow& +T''(t)-c^2\lambda T(t)&=0 +&&\Rightarrow& +T''(t) &= A \cos(c\sqrt\lambda t) + B \sin(c \lambda t) +\\ +&&&&&&&& + &= C \cos(c \lambda t+\delta) +\\ +\frac{\Delta U(x)}{U(x)}&=-\lambda^2 +&&\Rightarrow& +\Delta U &= -\lambda^2 U +\end{align*} +Die letzte Gleichung für die Funktion $U(x)$ hat die Form +eines Eigenwertproblems mit dem Eigenwert $-\lambda^2$. + +\begin{definition} +Eine Eigenfunktion eines Operators $L$ zum Eigenwert $\lambda$ +ist eine Funktion $U$ derart, dass $LU=\lambda U$. +\end{definition} + +Die Separation ermöglich also, das ursprüngliche Problem aufzuspalten +in ein Eigenwertproblem für eine nur ortsabhängige Funktion $U(x)$ +und eine Schwingungsgleichung für $T(t)$. +Die Schwingungsfrequenz $c \lambda $ hängt direkt mit dem +Eigenwert zusammen. +Die Funktion $U(x)$ beschreibt die Form der Membran, die Amplitude +in jedem Punkt, der Faktor $T(t)$ beschreibt die Schwingung. + + diff --git a/buch/chapters/part1.tex b/buch/chapters/part1.tex index b449aa2..51134ba 100644 --- a/buch/chapters/part1.tex +++ b/buch/chapters/part1.tex @@ -19,7 +19,7 @@ \input{chapters/060-integral/chapter.tex} %\input{chapters/070-reihenprodukte/chapter.tex} \input{chapters/080-funktionentheorie/chapter.tex} -%\input{chapters/090-funktional/chapter.tex} +\input{chapters/090-pde/chapter.tex} % Gamma und Pi % Eulersche Beta-Funktion diff --git a/buch/common/macros.tex b/buch/common/macros.tex index 2c6eea2..7c82180 100644 --- a/buch/common/macros.tex +++ b/buch/common/macros.tex @@ -104,6 +104,7 @@ \newtheorem{frage}[satz]{Frage} \newtheorem{problem}[satz]{Problem} \newtheorem{aufgabe}[satz]{Aufgabe} +\newtheorem{prinzip}[satz]{Prinzip} \newtheorem*{problem*}{Problem} \newtheorem{forderung}{Forderung}[chapter] \newtheorem{konsequenz}[satz]{Konsequenz} -- cgit v1.2.1