aboutsummaryrefslogtreecommitdiffstats
path: root/buch/papers
diff options
context:
space:
mode:
authorhaddoucher <reda.haddouche@ost.ch>2022-08-18 12:31:12 +0200
committerhaddoucher <reda.haddouche@ost.ch>2022-08-18 12:31:12 +0200
commitb3611aa8b6f2c56c8940c18c582de0fd3dd205f2 (patch)
tree4aa9eae2436e8742f5445f2a7406acf57b4cc805 /buch/papers
parenteinleitung (diff)
downloadSeminarSpezielleFunktionen-b3611aa8b6f2c56c8940c18c582de0fd3dd205f2.tar.gz
SeminarSpezielleFunktionen-b3611aa8b6f2c56c8940c18c582de0fd3dd205f2.zip
tschebyscheff kapitel
fertig geschrieben
Diffstat (limited to '')
-rw-r--r--buch/papers/sturmliouville/einleitung.tex13
-rw-r--r--buch/papers/sturmliouville/tschebyscheff_beispiel.tex31
2 files changed, 28 insertions, 16 deletions
diff --git a/buch/papers/sturmliouville/einleitung.tex b/buch/papers/sturmliouville/einleitung.tex
index 700ea1d..d497622 100644
--- a/buch/papers/sturmliouville/einleitung.tex
+++ b/buch/papers/sturmliouville/einleitung.tex
@@ -27,9 +27,9 @@ Alle homogene 2. Ordnung lineare gewöhnliche Differentialgleichungen können in
\subsection{Randbedingungen\label{sub:was-ist-das-slp-randbedingungen}}
Wenn von der Funktion $y(x)$ die Werte $x$ des jeweiligen Randes des Definitionsbereiches anzunehmen sind, also
\begin{equation}
- y(a) = y(b) = 0
+ y(a) = y(b) = 0,
\end{equation}
-, so spricht man von einer Dirichlet-Randbedingung\footnote{Die Dirichlet-Randbedingung oder auch Randbedingung des ersten Typs genannt ist nach dem deutschen Mathematiker Peter Gstav Lejeune Dirichlet benannt. Sie findet Anwendung auf gewöhnliche oder patielle Differentialgleichungen und gibt mit der Bedingung die Werte an, die für die abgeleitete Lösung innerhalb der Domänengrenze gelten.}, und von einer Neumann-Randbedingung\footnote{Die Neumann-Randbedingung oder auch Randbedingung des zweiten Typs genannt, ist nach dem deutschen Mathematiker Carl Neumann benannt. Sie legt die Werte fest, die eine Lösung entlang der Domänengrenze annehmen muss, wenn eine gewöhnliche oder partielle Differentialgleichung gestellt wird.} spricht man, wenn
+so spricht man von einer Dirichlet-Randbedingung\footnote{Die Dirichlet-Randbedingung oder auch Randbedingung des ersten Typs genannt ist nach dem deutschen Mathematiker Peter Gstav Lejeune Dirichlet benannt. Sie findet Anwendung auf gewöhnliche oder patielle Differentialgleichungen und gibt mit der Bedingung die Werte an, die für die abgeleitete Lösung innerhalb der Domänengrenze gelten.}, und von einer Neumann-Randbedingung\footnote{Die Neumann-Randbedingung oder auch Randbedingung des zweiten Typs genannt, ist nach dem deutschen Mathematiker Carl Neumann benannt. Sie legt die Werte fest, die eine Lösung entlang der Domänengrenze annehmen muss, wenn eine gewöhnliche oder partielle Differentialgleichung gestellt wird.} spricht man, wenn
\begin{equation}
y'(a) = y'(b) = 0
\end{equation}
@@ -53,15 +53,16 @@ Es ist genau das gleiche Prinzip wie bei den Matrizen, andere Eigenwerte ergeben
Es besteht eine Korrespondenz zwischen den Eigenwerten und den Eigenvektoren.
Das gleiche gilt auch beim Sturm-Liouville-Problem, und zwar
\begin{equation}
- \lambda \overset{Korrespondenz}\leftrightarrow y
-\end{equation}.
+ \lambda \overset{Korrespondenz}\leftrightarrow y.
+\end{equation}
Die Theorie besagt, wenn $y_m$, $y_n$ Eigenfuktionen des Sturm-Liouville-Problems sind, die verschiedene Eigenwerte $\lambda_m$, $\lambda_n$ ($\lambda_m \neq \lambda_n$) entsprechen, so sind $y_m$, $y_n$ orthogonal zu y -
dies gilt für das Intervall (a,b).
Somit ergibt die Gleichung
\begin{equation}
- \int_{a}^{b} w(x)y_m y_n = 0
-\end{equation}.
+ \label{eq:skalar-sturm-liouville}
+ \int_{a}^{b} w(x)y_m y_n = 0.
+\end{equation}
\subsection{Koeffizientenfunktionen}
Die Funktionen $p(x)$, $q(x)$ und $w(x)$ werden als Koeffizientenfunktionen mit ihren freien Variablen $x$ bezeichnet.
diff --git a/buch/papers/sturmliouville/tschebyscheff_beispiel.tex b/buch/papers/sturmliouville/tschebyscheff_beispiel.tex
index a18684f..3817dc0 100644
--- a/buch/papers/sturmliouville/tschebyscheff_beispiel.tex
+++ b/buch/papers/sturmliouville/tschebyscheff_beispiel.tex
@@ -4,8 +4,9 @@
% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil
%
-\subsection{Tschebyscheff-Polynome\label{sub:tschebyscheff-polynome}}
-Im Kapitel \ref{sub:beispiele_sturm_liouville_problem} sind die Koeffizientenfunktionen die man braucht schon aufgeliste, und zwar mit
+\subsection{Sind Tschebyscheff-Polynome orthogonal zueinander?\label{sub:tschebyscheff-polynome}}
+\subsubsection*{Definition der Koeffizientenfunktion}
+Im Kapitel \ref{sub:beispiele_sturm_liouville_problem} sind die Koeffizientenfunktionen, die man braucht, schon aufgeliste, und zwar mit
\begin{align*}
w(x) &= \frac{1}{\sqrt{1-x^2}} \\
p(x) &= \sqrt{1-x^2} \\
@@ -14,10 +15,12 @@ Im Kapitel \ref{sub:beispiele_sturm_liouville_problem} sind die Koeffizientenfun
Da die Sturm-Liouville-Gleichung
\begin{equation}
\label{eq:sturm-liouville-equation-tscheby}
- \frac{d}{dx}\lbrack \sqrt{1-x^2} \frac{dy}{dx} \rbrack + \lbrack 0 + \lambda \frac{1}{\sqrt{1-x^2}} \rbrack y = 0
+ \frac{d}{dx} (\sqrt{1-x^2} \frac{dy}{dx}) + (0 + \lambda \frac{1}{\sqrt{1-x^2}}) y = 0
\end{equation}
nun mit den Koeffizientenfunktionen aufgestellt werden kann, bleibt die Frage, ob es sich um ein reguläres oder singuläres Sturm-Liouville-Problem handelt.
-Für das reguläre Problem laut der Definition \ref{def:reguläres_sturm-liouville-problem} muss die funktion $p(x) = \sqrt{1-x^2}$, $p'(x) = -2x$, $q(x) = 0$ und $w(x) = \frac{1}{\sqrt{1-x^2}}$ stetig und reell sein - und sie sind es auch.
+
+\subsubsection*{regulär oder singulär?}
+Für das reguläre Problem laut der Definition \ref{def:reguläres_sturm-liouville-problem} muss die funktion $p(x) = \sqrt{1-x^2}$, $p'(x) = -2x$, $q(x) = 0$ und $w(x) = \frac{1}{\sqrt{1-x^2}}$ stetig und reell sein --- und sie sind es auch.
Auf dem Intervall $(-1,1)$ sind die Tschebyscheff-Polynome erster Art mit Hilfe von Hyperbelfunktionen
\begin{equation}
T_n(x) = \cos n (\arccos x)
@@ -28,22 +31,23 @@ Für $x>1$ und $x<-1$ sehen die Polynome wie folgt aus:
(-1)^n \cosh (n \arccos (-x)), & x<-1 \end{array}\right.
\end{equation},
jedoch ist die Orthogonalität nur auf dem Intervall $[ -1, 1]$ sichergestellt.
-Die nächste Bedingung beinhaltet, dass die Funktion $p(x)^{-1}$ und $w(x)>0$ sein müssen.
+Die nächste Bedingung beinhaltet, dass die Funktion $p(x)$ und $w(x)>0$ sein müssen.
Die Funktion
\begin{equation*}
p(x)^{-1} = \frac{1}{\sqrt{1-x^2}}
\end{equation*}
-ist die gleiche wie $w(x)$.
+ist die gleiche wie $w(x)$ und erfüllt die Bedingung.
+\subsubsection*{Randwertproblem}
Für die Verifizierung der Randbedingungen benötigt man erneut $p(x)$.
Da sich die Polynome nur auf dem Intervall $[ -1,1 ]$ orthogonal verhalten, sind $a = -1$ und $b = 1$ gesetzt.
Beim einsetzen in die Randbedingung \ref{eq:randbedingungen}, erhält man
\begin{equation}
\begin{aligned}
k_a y(-1) + h_a y'(-1) &= 0
- k_b y(-1) + h_b y'(-1) &= 0
+ k_b y(-1) + h_b y'(-1) &= 0.
\end{aligned}
-\end{equation}.
+\end{equation}
Die Funktion $y(x)$ und $y'(x)$ sind in diesem Fall die Tschebyscheff Polynome (siehe \label{sub:definiton_der_tschebyscheff-Polynome}).
Es gibt zwei Arten von Tschebyscheff Polynome: die erste Art $T_n(x)$ und die zweite Art $U_n(x)$.
Jedoch beachtet man in diesem Kapitel nur die Tschebyscheff Polynome erster Art (\ref{eq:tschebyscheff-polynome}).
@@ -52,12 +56,19 @@ Somit erhält man
\begin{equation}
\begin{aligned}
k_a T_2(-1) + h_a T_{2}'(-1) &= k_a = 0\\
- k_b T_2(1) + h_b T_{2}'(1) &= k_b = 0
+ k_b T_2(1) + h_b T_{2}'(1) &= k_b = 0.
\end{aligned}
-\end{equation}.
+\end{equation}
Ähnlich wie beim Beispiel der Wärmeleitung in einem homogenen Stab kann man, damit die Bedingung $|k_i|^2 + |h_i|^2\ne 0$ erfüllt ist, können beliebige $h_a \ne 0$ und $h_b \ne 0$ gewählt werden.
Somit ist erneut gezeigt, dass die Randbedingungen der Tschebyscheff-Polynome auf die Sturm-Liouville-Randbedingungen erfüllt und alle daraus resultierenden Lösungen orthogonal sind.
+\begin{beispiel}
+ Die Gleichung \ref{eq:skalar-sturm-liouville} mit $y_m = T_1(x)$ und $y_n(x) = T_2(x)$ eingesetzt sowie $a=-1$ und $b = 1$ ergibt
+ \[
+ \int_{-1}^{1} w(x) x (2x^2-1) dx = 0.
+ \]
+\end{beispiel}
+