Aufgabe2.tex 2.6 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879
  1. \section*{Aufgabe 2}
  2. Zeige die Aussage für $2\times2$ Matrizen durch Gauß-en mit
  3. Spaltenpivotwahl.
  4. \subsection*{Lösung}
  5. \subsubsection*{Behauptung:}
  6. Für alle tridiagonalen Matrizen gilt:
  7. \begin{enumerate}
  8. \item[(i)] Die Gauß-Elimination erhält die tridiagonale Struktur
  9. \item[(ii)] $\rho_n(A) := \frac{\alpha_\text{max}}{\max_{i,j} |a_{ij}|} \leq 2$
  10. \end{enumerate}
  11. \subsubsection*{Beweis:}
  12. \paragraph{Teil 1: (i)}
  13. \begin{align}
  14. A &= \begin{gmatrix}[p]
  15. * & * & & \\
  16. * & \ddots & \ddots & \\
  17. & \ddots & \ddots & * \\
  18. & & * & *
  19. \rowops
  20. \add[\cdot \frac{-a_{21}}{a_{11}}]{0}{1}
  21. \end{gmatrix}
  22. \end{align}
  23. Offensichtlich ändert diese Operation nur Zeile 2. $a_{21}$ wird zu 0,
  24. $a_{22}$ ändert sich irgendwie, alles andere bleibt unverändert.
  25. Die gesammte Matrix ist keine tridiagonale Matrix mehr, aber die
  26. um Submatrix in $R^{(n-1) \times (n-1)}$ ist noch eine.
  27. Muss man zuvor Zeile 1 und 2 tauschen (andere Zeilen kommen nicht in
  28. Frage), so ist später die Stelle $a_{21} = 0$, $a_{22}$ ändert sich
  29. wieder irgendwie und $a_{23}$ ändert sich auch. Dies ändert aber nichts
  30. an der tridiagonalen Struktur der Submatrix.
  31. \paragraph{Teil 2: (ii) für $A \in \mathbb{R}^{2 \times 2}$}
  32. Sei $\begin{pmatrix}a_{11} & a_{12}\\a_{21} & a_{22} \end{pmatrix} \in \mathbb{R}^{2 \times 2}$
  33. beliebig.
  34. O.B.d.A sei die Spaltenpivotwahl bereits durchgeführt, also $|a_{11}| \geq |a_{21}|$.
  35. Nun folgt:
  36. \begin{align}
  37. \begin{gmatrix}[p]
  38. a_{11} & a_{12}\\
  39. a_{21} & a_{22}
  40. \rowops
  41. \add[\cdot \frac{-a_{21}}{a_{11}}]{0}{1}
  42. \end{gmatrix}\\
  43. \leadsto
  44. \begin{gmatrix}[p]
  45. a_{11} & a_{12}\\
  46. 0 & a_{22} - \frac{a_{12} \cdot a_{21}}{a_{11}}
  47. \end{gmatrix}
  48. \end{align}
  49. Wegen $|a_{11}| \geq |a_{21}|$ gilt:
  50. \begin{align}
  51. \|\frac{a_{21}}{a_{11}}\| \leq 1
  52. \end{align}
  53. Also insbesondere
  54. \begin{align}
  55. \underbrace{a_{22} - a_{12} \cdot \frac{a_{21}}{a_{11}}}_{\leq \alpha_\text{max}} \leq 2 \cdot \max_{i,j}|a_{ij}|
  56. \end{align}
  57. Damit ist Aussage (ii) für $A \in \mathbb{R}^{2 \times 2}$ gezeigt.
  58. \paragraph{Teil 3: (ii) für allgemeinen Fall}
  59. Aus Teil 2 folgt die Aussage auch direkt für größere Matrizen.
  60. Der worst case ist, wenn man beim Addieren einer Zeile auf eine
  61. andere mit $\max_{i,j}|a_{ij}|$ multiplizieren muss um das erste nicht-0-Element
  62. der Zeile zu entfernen und das zweite auch $\max_{i,j}|a_{ij}|$ ist.
  63. Dann muss man aber im nächsten schritt mit einem Faktor $\leq \frac{1}{2}$
  64. multiplizieren, erhält also nicht einmal mehr $2 \cdot \max_{i,j}|a_{ij}|$.