Mehrgittermethoden

Übungsaufgaben

Robert Speck & Dieter Moser, Sommersemester 2016

</center>


Über Systemmatrizen

In der Vorlesung haben wir das Gleichungssytem $Au = f$ für das Possion-Problem in 1D mit Dirichlet-Rändern durch eine Diskretisierung mit Hilfe von zentrierten Finiten Differenzen aufgestellt.

  1. Wie sieht das System für periodische Randbedingungen $u(0) = u(1)$ und $u'(0) = u'(1)$ aus? Nutzen Sie die Gitterpunkte $x_0$, $x_1$, ..., $x_{N-1}$. Wie viele Unbekannte hat das Gleichungsssystem?
  2. Wie sieht das System mit Dirichlet-Rändern in zwei Dimensionen aus? Nutzen Sie Blockmatrizen zur Darstellung von $A$. Wie viele Unbekannte hat das System nun und wie viele Einträge ungleich Null hat die Systemmatrix?

Über das Kroneckerprodukt (Teil 1)

Ist $A$ eine $m\times n$-Matrix und $B$ eine $p\times r$-Matrix, so ist das Kronecker-Produkt $C = A \otimes B$ definiert als $$C = (a_{ij} \cdot B) =\begin{pmatrix} a_{11} B & \cdots & a_{1n} B \\ \vdots & \ddots & \vdots \\ a_{m1} B & \cdots & a_{mn} B \end{pmatrix}.$$

  1. Ist das Kronecker-Produkt kommutativ?

  2. Zeigen Sie

    a.) $A\otimes (B+C)=A\otimes B+A\otimes C $

    b.) $AC\otimes BD=(A\otimes B)(C\otimes D)$

    c.) $\mathrm{Spur}(A \otimes B) = \mathrm{Spur}(A) \cdot \mathrm{Spur}(B)$

Über das Kroneckerprodukt (Teil 2)

Ist $A$ eine $m\times n$-Matrix und $B$ eine $p\times r$-Matrix, so ist das Kronecker-Produkt $C = A \otimes B$ definiert als $$C = (a_{ij} \cdot B) =\begin{pmatrix} a_{11} B & \cdots & a_{1n} B \\ \vdots & \ddots & \vdots \\ a_{m1} B & \cdots & a_{mn} B \end{pmatrix}.$$

Zeigen Sie

a.) Sind $\{\lambda_i\}_{i=1..n}$ die Eigenwerte von $A$ und $\{\mu_j\}_{j=1..m}$ die Eigenwerte von $B$ dann gilt $\{\lambda_i \, \mu_j\}_{i=1..n \atop j=1..m} \mbox{sind die Eigenwerte von }\; A \otimes B$.

b.) Für die Spektralnorm gilt $\| A \otimes B \|_2 = \| A \|_2 \cdot \| B \|_2$.

c.) Sind $A,B$ invertierbar, so ist die Inverse $(A\otimes B)^{-1}=A^{-1} \otimes B^{-1}$.

Über den Rayleigh-quotienten

Sei $H$ eine hermitsche Matrix, so definiert sich der zugehörige Rayleigh quotient als $$R_H(x) = \frac{x^*Hx}{x^*\cdot x}.$$

  1. Zeigen Sie, dass für den größten bzw. kleinsten Eigenwert $\lambda_M,\lambda_m$ von $H$ $$ \lambda_m = \min_x R_H(x) \; \mbox{und} \; \lambda_M = \max_x R_H(x) $$ gilt.
  2. Sei $H = A^*A$. Zeigen sie, dass $\max_x R_H(x)^{1/2}$ eine Norm von $A$ ist.

Über persymmetrische Matrizen

Sei $\mathcal{P}_1 = \{ A \in K^{N \times N} | a_{i,j} = a_{n-j+1,n-i+1} \}$ und $\mathcal{P}_2 = \{ A \in K^{N \times N} | JA = A^TJ \}$, mit $ J = (\delta_{i,n-j+1})_{ij} = \begin{pmatrix} 0 & & 1 \\ & \scriptstyle\cdot^{\,\scriptstyle\cdot^{\,\scriptstyle\cdot}} & \\ 1 & & 0 \end{pmatrix}$.

Zeigen Sie

  1. $\mathcal{P}_1 = \mathcal{P}_2$
  2. $\mathcal{P}_1$ is ein Untervektorraum des $K^{N \times N}$.

Wie muss $\mathcal{P}_1$ eingeschränkt werden, damit

  1. $(\mathcal{P}_1,+,0,-)$ eine abelsche Gruppe ist
  2. $(\mathcal{P}_1,\cdot)$ eine Gruppe ist
  3. $(\mathcal{P}_1,+,0,-,\cdot)$ ein Ring ist

Über zyklische Matrizen (Teil 1)

Sei $A$ von der Gestalt $$A:=\begin{pmatrix} a_0&a_{n-1}&a_{n-2}&\ldots&a_1\\ a_1&a_0&a_{n-1}&\ldots&a_2\\ a_2&a_1&a_0&\ldots&a_3\\ &\ddots&\ddots&\ddots\\ a_{n-1}&a_{n-2}&a_{n-3}&\ldots&a_0\end{pmatrix}. $$

  1. Zeigen Sie, dass A persymmetrisch ist.
  2. Berechnen Sie das charakteristische Polynom, die Eigenwerte, Eigenvektoren und die Determinante von $A$.
  3. Seien $B$ und $C$ zyklische miteinander kommutierende Matrizen mit den Eigenwerten $\{\beta_k\},\{\zeta_k\}$. Was sind die Eigenwerte von $BC$.
  4. Zeigen Sie, dass der Raum der zyklischen Matrizen ein Untervektorraum des $K^{N \times N}$

Fourier Reihen

Sei $\{t_k\}_{-\infty}^{\infty}$ eine absolut summierbare Folge, d.h. $\sum_{k=-\infty}^{\infty}|t_k| < \infty$. Sei desweiteren $f(\lambda) = \lim_{n \to \infty} \sum_{k=-n}^{n}t_k e^{-ik\lambda}$.

  1. Zeigen Sie, dass $S_n(\lambda) = \sum_{k=-n}^{n}t_k e^{-ik\lambda}$ gleichmässig gegen $f(\lambda)$ konvergiert.
  2. Folgern Sie, dass $f(\lambda)$ Riemann-integrierbar und beschränkt auf $[0,2\pi]$ ist.
  3. Finden Sie mithilfe der inversen Fouriertransformation eine Darstellung von $t_k$ unter Verwendung von $f(\lambda)$.

Wir nennen $f(\lambda)$ eine Funktion der Wiener Klasse.

Über Toeplitz-Matrizen (Teil 1)

Eine Matrix nennt man eine Toeplitz-Matrix falls die Werte auf der Hauptdiagonale und allen Nebendiagonalen konstant sind. Mithilfe einer Funktion der Wiener Klasse $f(\lambda)$ lässt sich die Klasse der Toeplitz-Matrizen $$ T_n(f) = \{\frac{1}{2 \pi} \int_0^{2\pi} f(\lambda) e^{-i(k-j)\lambda} \mathrm{d}\lambda ; k,j = 0,1,\ldots,n-1\}$$ konstruiren.

  1. Welches $f(\lambda)$ und $n$ generiert den Finite Differenzen Poisson Operator.
  2. Zeigen Sie, dass $T_n(f(\lambda))$ genau dann hermitesch ist wenn $f(\lambda)$ eine reelwertige Funktion ist.
  3. Sei $f(\lambda)$ reelwertig und begrenzt durch $m_f \leq f(\lambda) \leq M_f$. Beweisen Sie mithilfe des Raleigh Koeffizienten, dass für die Eigenwerte $\tau_{n,k}$ von $T_n(f)$ gilt $m_f\leq\tau_{n,k}\leq M_f$.
  4. Zeigen Sie für nicht hermitesche $T_n(f)$, dass die Ungleichung $\| T_n(f) \| \leq M_{|f_r|}+M_{|f_i|} \leq 2 M_{|f|}$ gilt, wobei $f(\lambda) = f_r(\lambda) + i f_i(\lambda)$ mit den realwertigen Funktionen $f_i,f_r$ aus der Wiener Klasse.

Tipp: Nutzen Sie den Satz von Parseval für $x^* \cdot x$.

Über die Hilbert-Schmidt Norm

Wir definieren die Hilbert-Schmidt Norm einer Matrx $A \in K^{n \times n}$ als $$ |A| = \left( \frac{1}{n}\sum_{i = 0}^{n-1}\sum_{i = 0}^{n-1} |a_{i,j}|^2 \right)^{1/2}.$$ Zeigen Sie

  1. $|A| = \left( \frac{1}{n}\mbox{Spur}(A^*A) \right)^{1/2}$
  2. $|A| = \left( \frac{1}{n}\sum_{k=0}^{n-1}\lambda_k\right)^{1/2}$, wobei $\lambda_k$ die Eigenwerte von $A^*A$ sind
  3. $|A| \leq \|A\|$

Überlegen Sie sich warum diese Norm, auch schwache Norm genannt wird.

Über asymptotisch äquivalente Folgen von Matrizen

Seien $\{A_n\}$ und $\{B_n\}$ Folgen von $n\times n$ Matrizen, welche beschränkt bzgl. der starken Norm sind: $$ \|A_n\|,\|B_n\| \leq M \le \infty, n=1,2,\ldots $$ und bzgl. der schwachen Norm konvergieren $$\lim_{n \to \infty} |A_n -B_n| = 0$$. Wir nennen diese Folgen asymptotisch äquivalent und notieren dies als $A_n \sim B_n$. Zeigen Sie nun für $\{A_n\}$ , $\{B_n\}$ und $\{C_n\}$, welche jeweils die Eigenwerte $\{\alpha_{n,i}\}$,$\{\beta_{n,i}\}$ und $\{\zeta_{n,i}\}$ haben, folgende Zusammenhänge.

  1. Wenn $A_n \sim B_n$, dann $\lim_{n \to \infty} |A_n| = \lim_{n \to \infty} |B_n| $
  2. Wenn $A_n \sim B_n$ und $B_n \sim C_n$, dann $A_n \sim C_n$
  3. Wenn $A_nB_n \sim C_n$ und $\|A_n^{-1}\|\leq K \le \infty$, dann gilt $B_n \sim A_n^{-1}C_n$
  4. Wenn $A_n \sim B_n$, dann $\exists -\infty \le m,M\le \infty$, s.d. $m\leq \alpha_{n,i}, \beta_{n,i}\leq M \; \forall n\geq 1 \mbox{und}\; k\geq 0$

Tipp: Nutzen Sie ohne Beweis, dass gilt $|GH|\leq \|G\| \cdot |H|$. Wer Lust hat kann es aber auch Beweisen.

Über zyklische Matrizen (Teil 2)

Sei $f(\lambda):[0,2\pi]\mapsto C$ eine Funktion der Wiener Klasse und $\{t_k\}_{-s}^{s}$ die zugehörigen Fourier Koeffizienten. D.h. $$f(\lambda) = \sum_{k=-s}^{s} t_k \exp(-i\lambda)^k.$$
Wir nutzen diese Koeffizienten zur Konstruktion von

$$C_{n}(f)=\begin{pmatrix} t_0 & \ldots & t_{-s}& \mathbf{0}_m &t_{s} &\ldots & t_{1} \\ \vdots & & & \ddots& &\ddots & \vdots \\ t_{s} & & & & & & t_{s} \\ \mathbf{0}_m^T &\ddots & & t_0 & & \ddots& \mathbf{0}_m^T \\ t_{-s} & & & & & & t_{-s} \\ \vdots & \ddots & & \ddots& & & \vdots \\ t_{-1} & \ldots & t_{-s}& \mathbf{0}_m & t_s & \ldots & t_0 \end{pmatrix} \in K^{n \times n}, \mbox{mit}\; n > 2s+1.$$

$C_n(f)$ ist zyklisch und besitzt die Eigenwerte $\{\psi_{n,j}\}$.

  1. Zeigen Sie, dass für die Eigenwerte von $C_n(f)$ gilt $\psi_{n,j} = f(\frac{2 \pi j}{n})$.
  2. Folgern Sie daraus, dass $$ \lim_{n \to \infty} \frac{1}{n} \sum_{j=0}^{n-1} \psi_{n,j}^s = \frac{1}{2\pi}\int_{0}^{2\pi} f(\lambda)^s \mathrm{d}\lambda$$
  3. Sei $C_n(f)$ hermitesch und $F$ eine stetige Funktion, definiert auf dem Bild von $f(\lambda)$. Zeigen Sie $$ \lim_{n \to \infty} \frac{1}{n} \sum_{j=0}^{n-1} F(\psi_{n,j}) = \frac{1}{2\pi}\int_{0}^{2\pi} F(f(\lambda)) \mathrm{d}\lambda$$

Tipp: Nutzen Sie für 3. den Satz von Weierstraß. Dieser besagt, dass für jede stetige Funktion $F:[a,b] \mapsto K$ eine Folge von Polynomen $\{p_n\}$ existiert, s.d. $$ \lim_{n \to \infty} p_n(x) =F(x) $$ gleichmäßig auf $[a,b]$ konvergiert.

Über Toeplitz Matrizen (Teil 2)

Sein $f(\lambda):[0,2\pi]\mapsto C$ eine Funktion der Wiener Klasse und $\{t_k\}_{-s}^{s}$ die zugehörigen Fourier Koeffizienten. Wir nutzen diese Koeffizienten zur Konstruktion von

$$T_{n}(f)=\begin{pmatrix} t_0 & \ldots & t_{-s}& & & & \\ \vdots & & & \ddots & & & \\ t_{s} & & & & & & \\ &\ddots & & t_0 & & \ddots& \\ & & & & & & t_{-s} \\ & & & \ddots & & & \vdots \\ & & & & t_s & \ldots & t_0 \end{pmatrix} \in K^{n \times n}, \mbox{mit}\; n > 2s+1.$$

Diese nennt man eine Band-Toeplitz-Matrix. Die Eigenwerte von $T_n(f)$ nennen wir $\{\tau_{n,k}\}$.

  1. Zeigen Sie: $C_n(f) \sim T_n(f)$
  2. Sei $f(\lambda)$ zusätzlich reel-wertig. Beweisen Sie, dass für ein stetiges $F:\mbox{Bild}(f)\mapsto K$ $$ \lim_{n \to \infty} \frac{1}{n} \sum_{j=0}^{n-1} F(\tau_{n,j}) = \frac{1}{2\pi}\int_{0}^{2\pi} F(f(\lambda)) \mathrm{d}\lambda$$gilt.
  3. Es gelte nun zusätzlich $f(\lambda) >0$, zeigen Sie mithilfe eines geeigneten $F$, dass $$ \lim_{n \to \infty} \mbox{det}(T_n(f))^{1/n} =\lim_{n \to \infty} \mbox{det}(C_n(f))^{1/n}$$ gilt.

Tipp: Für zwei asymptotisch äquivalente Folgen von Matrizen $\{A_n\}$ und $\{B_n\}$, mit den Eigenwerten $\{\alpha_{n,i}\}$ und $\{\beta_{n,i}\}$, kann man $$ \lim_{n \to \infty} \frac{1}{n} \sum_{k=0}^{n-1} (\alpha_{n,k}^s - \beta_{n,k}^s) = 0$$ zeigen, solange $s$ ein positiver integer ist. Diese Aussagen kann ohne Beweis verwendet werden.


In [ ]: