From 18e369cbf05462e7659bc1512b9b29ae13172e46 Mon Sep 17 00:00:00 2001 From: Adrian Kummerlaender Date: Sun, 12 Feb 2017 21:54:45 +0100 Subject: Complete LU-decomposition section, add matrix types, Cholesky-decomposition --- numerik_1.tex | 45 ++++++++++++++++++++++++++++++++++++++++++++- 1 file changed, 44 insertions(+), 1 deletion(-) diff --git a/numerik_1.tex b/numerik_1.tex index 02c8ca5..53736f7 100644 --- a/numerik_1.tex +++ b/numerik_1.tex @@ -100,6 +100,20 @@ Für $A \in \K^{n \times n} \in GL_n{\R}$, $\|\cdot\|$ induzierte Matrixnorm: 1 = \|Id\| = \|AA^{-1}\| &\leq \|A\| \cdot \|A^{-1}\| = \kappa(A) \end{align*} +\section*{Besondere Matrizen} + +\subsection*{Diagonaldominante Matrizen} + +$A \in \R^{n \times n}$ ist diagonaldominant, falls: + +$\forall i \in \{1,\cdots,n\} : |a_{i,i}| > \sum_{j=1,j\neq i}^n |a_{i,j}|$ + +Insbesondere sind solche $A$ regulär. + +\subsection*{Positive Definitheit} + +$A \in \R^{n \times n}$ ist positiv definit d.h. $A > 0$ falls $A=A^T$ und $\forall x \in \R^n \setminus \{0\} : x^TAx > 0$. + \section*{Direkte Verfahren zur LGS Lösung} \subsection*{Cramersche Regel} @@ -114,9 +128,38 @@ Obere Dreicksmatrizen können mittels Rückwärtssubstitution, untere Dreiecksma \subsection*{LR-Zerlegung} +$A = LR$ wobei $L$ untere Dreiecksmatrix mit $1$-Diagonale und $R$ obere Dreicksmatrix. + +\subsubsection*{Berechnung LR-Zerlegung} + +Die LR-Zerlegung existiert insofern die Diagonaleinträge nicht verschwinden. Insbesondere gilt dies für diagonaldominante Matrizen. + +\begin{enumerate} + \item Spaltenweises nullen der der unteren Einträge mittels \emph{Gauß}, Matrizen $L_1, \cdots, L_{n-1}$ + \item $L = L_1^{-1} \cdots L_{n-1}^{-1}$, $R=L_{n-1} \cdots L_1 A$ +\end{enumerate} + +\subsubsection*{Lösung $Ax=b$ mittels LR-Zerlegung} + +\begin{enumerate} + \item $A=LR$ berechnen + \item $Lz=b$ Vorwärtssubstitution + \item $Rx=z$ Rückwärtssubstitution +\end{enumerate} + +\subsubsection*{Spaltenpivotsuche} + +Die normale LR-Zerlegung ist nur Vorwärts- und nicht Rückwärtsstabil. Dies kann durch Spaltenpivotsuche verbessert werden. Hier wird in jedem Schritt mittels einer Permutationsmatrix immer mit der größten verbleibenden Zeile eliminiert. + +\vspace{1mm} + +Für alle regulären Matrizen existiert eine Spaltenpivotsuchen LR-Zerlegung so, dass $PA=LR$. + \subsection*{Cholesky-Zerlegung} -\subsection*{QR-Zerlegung} +Für $A > 0$ existiert untere Dreiecksmatrix $L$ mit positiver Diagonale, so dass $A = LL^T$ + +\subsection*{QR-Faktorisierung} \section*{Lineare Ausgleichsprobleme} -- cgit v1.2.3