Tridiagonale Matrix - LU < Matrizen < Lineare Algebra < Hochschule < Mathe < Vorhilfe
|
| Status: |
(Frage) beantwortet  | | Datum: | 21:52 Do 02.11.2017 | | Autor: | rem |
| Aufgabe | Let A [mm] \in \IR^{n\times n} [/mm] be a tridiagonal matrix of the form
A = [mm] \pmat{ a_1 & b_2 & 0 & ...& 0 \\ c_2 & a_2 & b_2 & ... & ...\\ 0 & ... & ... & ... & 0 \\ ... & ... & c_{n-1} & a_{n-1} & b_{n-1} \\ 0 & ... & 0 & c_n & a_n },
[/mm]
with [mm] a_k, b_k, c_k \ne [/mm] 0, for all k = 1, ..., n.
If the matrix is strictly diagonally dominant, it admits a LU factorization where the matrices L and U have the special form
L = [mm] \pmat{ d_1 & 0 & ... & ...& 0 \\ c_2 & d_2 & ... & ... & ...\\ 0 & ... & ... & ... & ... \\ ... & ... & c_{n-1} & d_{n-1} & 0 \\ 0 & ... & 0 & c_n & a_n },
[/mm]
U = [mm] \pmat{ 1 & e_1 & 0 & ...& 0 \\0 & 1 & e_2 & ... & ...\\ ... & ... & ... & ... & 0 \\ ... & ... & ... & 1 & e_{n-1} \\ 0 & ... & ... & 0 & 1 }.
[/mm]
Write an algorithm for computing their entries [mm] d_k [/mm] and [mm] e_k, [/mm] k = 1, ..., n. How many operations are required? |
Hallo!
Was mich an diesem Beispiel verwirrt ist zum einen die 1er Diagonal in der oberen Dreicksmatrix U sowie das Element [mm] a_n [/mm] in der unteren Dreiecksmatrix L. Im Internet und diversen Büchern gib es zu genau diesem Thema einiges an Material allerdings ist dort die 1er Diagonale in der unteren Dreiecksmatrix. Kann es sein das [mm] a_n [/mm] in Matrix L ein Fehler ist? Bin für jede Hilfe/Anregung dankbar!
LG
|
|
| |
|
| Status: |
(Antwort) fertig  | | Datum: | 01:06 Fr 03.11.2017 | | Autor: | Fulla |
> Let A [mm]\in \IR^{n\times n}[/mm] be a tridiagonal matrix of the
> form
>
> A = [mm]\pmat{ a_1 & b_2 & 0 & ...& 0 \\ c_2 & a_2 & b_2 & ... & ...\\ 0 & ... & ... & ... & 0 \\ ... & ... & c_{n-1} & a_{n-1} & b_{n-1} \\ 0 & ... & 0 & c_n & a_n },[/mm]
>
> with [mm]a_k, b_k, c_k \ne[/mm] 0, for all k = 1, ..., n.
> If the matrix is strictly diagonally dominant, it admits a
> LU factorization where the matrices L and U have the
> special form
>
> L = [mm]\pmat{ d_1 & 0 & ... & ...& 0 \\ c_2 & d_2 & ... & ... & ...\\ 0 & ... & ... & ... & ... \\ ... & ... & c_{n-1} & d_{n-1} & 0 \\ 0 & ... & 0 & c_n & a_n },[/mm]
>
> U = [mm]\pmat{ 1 & e_1 & 0 & ...& 0 \\0 & 1 & e_2 & ... & ...\\ ... & ... & ... & ... & 0 \\ ... & ... & ... & 1 & e_{n-1} \\ 0 & ... & ... & 0 & 1 }.[/mm]
>
> Write an algorithm for computing their entries [mm]d_k[/mm] and [mm]e_k,[/mm]
> k = 1, ..., n. How many operations are required?
> Hallo!
>
> Was mich an diesem Beispiel verwirrt ist zum einen die 1er
> Diagonal in der oberen Dreicksmatrix U sowie das Element
> [mm]a_n[/mm] in der unteren Dreiecksmatrix L. Im Internet und
> diversen Büchern gib es zu genau diesem Thema einiges an
> Material allerdings ist dort die 1er Diagonale in der
> unteren Dreiecksmatrix. Kann es sein das [mm]a_n[/mm] in Matrix L
> ein Fehler ist? Bin für jede Hilfe/Anregung dankbar!
Hallo rem,
ich denke auch, dass das [mm] $a_n$ [/mm] ein Fehler ist. Denn sonst berechnete sich der Eintrag [mm] $a_{n,n}$ [/mm] zu [mm] $c_n\cdot e_{n-1}+a_n$, [/mm] d.h. [mm] $c_n\cdot e_{n-1}=0$, [/mm] bzw., da [mm] $c_k\ne [/mm] 0$ bedeutet das, dass [mm] $e_{n-1}=0$. [/mm] Das ist glaube ich nicht beabsichtigt. Da sollte ein [mm] $d_n$ [/mm] stehen. (Aber frag vielleicht sicherheitshalber nochmal nach...)
Ob die 1er jetzt in L oder U stehen ist egal. Die [mm] $d_k$ [/mm] und [mm] $e_k$ [/mm] fallen je nach dem natürlich anders aus! Wichtig ist nur, dass es möglich ist, so eine Zerlegung zu finden, bei der eine der beiden Matrizen 1er auf der Diagonalen hat. Dadurch kannst du eben der Reihe nach die [mm] $d_k$ [/mm] und [mm] $e_k$ [/mm] berechnen.
Ich habe da jetzt nicht rumgerechnet, aber ich würde erstmal versuchen mit dem [mm] $a_n$ [/mm] in L zu rechnen. Wenn das geht, super! Wenn das zu Poblemen führt, rechne mit einem [mm] $d_n$ [/mm] an der Stelle.
Ich vermute aber, dass das [mm] $a_n$ [/mm] die Sache sogar einfacher macht...
Lieben Gruß,
Fulla
|
|
|
|
