aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorAdrian Kummerlaender2017-02-12 22:37:40 +0100
committerAdrian Kummerlaender2017-02-12 22:37:40 +0100
commit548cc449bc37bf8609da1415b2430358b5733fe0 (patch)
tree910f2ad70b82c72773f74b04fb47210a2b3f823b
parent18e369cbf05462e7659bc1512b9b29ae13172e46 (diff)
downloadmath_reference_sheets-548cc449bc37bf8609da1415b2430358b5733fe0.tar
math_reference_sheets-548cc449bc37bf8609da1415b2430358b5733fe0.tar.gz
math_reference_sheets-548cc449bc37bf8609da1415b2430358b5733fe0.tar.bz2
math_reference_sheets-548cc449bc37bf8609da1415b2430358b5733fe0.tar.lz
math_reference_sheets-548cc449bc37bf8609da1415b2430358b5733fe0.tar.xz
math_reference_sheets-548cc449bc37bf8609da1415b2430358b5733fe0.tar.zst
math_reference_sheets-548cc449bc37bf8609da1415b2430358b5733fe0.zip
Add QR decomposition, Householder reflexion section
-rw-r--r--numerik_1.tex33
-rw-r--r--zusammenfassung.tex2
2 files changed, 34 insertions, 1 deletions
diff --git a/numerik_1.tex b/numerik_1.tex
index 53736f7..473259d 100644
--- a/numerik_1.tex
+++ b/numerik_1.tex
@@ -110,7 +110,7 @@ $\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}
+\subsection*{Positiv definite Matrizen}
$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$.
@@ -161,6 +161,37 @@ Für $A > 0$ existiert untere Dreiecksmatrix $L$ mit positiver Diagonale, so das
\subsection*{QR-Faktorisierung}
+Für alle $A \in \R^{m \times n}$ mit $m \geq n$ und $Rang(A)=n$ existiert $A=QR$ mit unitärem $Q \in \R^{m \times m}$ und injektiver oberer Dreiecksmatrix $R$.
+
+\subsubsection*{Householder-Reflexionen}
+
+$$H(v) := Id_m - 2 \frac{vv^T}{v^Tv} = Id_m - 2 \frac{vv^T}{\|v\|_2^2} \text{ für } \forall v \in \R^m \setminus \{0\}$$
+
+Solche Householder-Reflexionen $H(v)$ sind orthogonal, d.h. $H(v)^T=H(v)$ und $H(v)^2=Id_m$.
+
+Wegen $H(v)v=v-2v=-v$ und $\forall w \in spann\{v\}^\perp : H(w)w=w$ ist $H(v)$ Spiegelung an der Hyperebene $spann\{v\}^\perp$.
+
+Solche Reflexionen können durch wiederholte Anwendung Matrizen in obere Dreiecksgestalt überführen:
+
+\vspace{1mm}
+
+Gesucht ist $v \in \R^m$ für $y \in \R^m$ s.d.:
+
+\vspace{-2mm}
+$$H(v)y=y - 2 \frac{\skp{v,y}}{\|v\|_2^2}v \overset{!}{=} \alpha e_1$$
+
+Vermeidung von Auslöschung: $\alpha := -sign(y_1)\|y\|_2$
+
+Es ergibt sich mit $v:=y-\alpha e_1$: $H(y-\alpha e_1)y=\alpha e_1$.
+
+\vspace{1mm}
+
+Seien $Q_k$ die sukzessiven, auf $m \times m$ erweiterten, Householder-Reflexionen. Dann gilt:
+
+\vspace{1mm}
+
+$R:=Q_p \cdots Q_1 A$, $Q:=Q_1^T \cdots Q_p^T$ s.d. $A=QR$.
+
\section*{Lineare Ausgleichsprobleme}
\section*{Iterative Verfahren zur LGS Lösung}
diff --git a/zusammenfassung.tex b/zusammenfassung.tex
index d6c4ae9..23f0203 100644
--- a/zusammenfassung.tex
+++ b/zusammenfassung.tex
@@ -20,6 +20,8 @@
\newcommand{\N}{\mathbb{N}}
\newcommand{\K}{\mathbb{K}}
+\newcommand{\skp}[1]{\langle #1 \rangle}
+
\begin{document}
\small{
\begin{multicols*}{4}