LU rozklad s pivotaci
Napsal: 8. 3. 2006 00:05
mel bych jeden technicky...
kdyz pocitam LU rozklad se sloupcovou pivotaci pro matici A, tak mi plati takova ta rovnost P*A = L*U, kde P je matice permutaci, ktera mi v kazdem eliminacnim kroku, prehodi radky tak jak potrebuju... U je horni trojuhelnikova a L je dolni trojuhelnikova...
jde mne o to, jak spocitat tu matici L
bud je takovy ten zpusob nadrzo, kdy spocitam permutacni matici P a matici U, coz je jednoduche a pak teda prenasobim celou rovnici zprava inverzem U
nebo postup algoritmicky, kdy pocitam postupne pomaloucku L1 az Lk a ruzne to prenasobuju jednotlivyma permutacnima maticema a pak je teda L soucin tech Li - tento postup je ovsem o hodne zdlouhavejsi (tedy nachylnejsi na chyby) ale asi spravnejsi
Janovsky, kdyz nam ukazoval ty zkusebni pisemky, tam presne tenhle priklad mel - otazka teda je, jestli staci pouzit ten prvni postup a clovek je hotovy do tri minut nebo se s tim alespon 20 minut patlat pres ruzne inverni matice atd...
co myslite?
kdyz pocitam LU rozklad se sloupcovou pivotaci pro matici A, tak mi plati takova ta rovnost P*A = L*U, kde P je matice permutaci, ktera mi v kazdem eliminacnim kroku, prehodi radky tak jak potrebuju... U je horni trojuhelnikova a L je dolni trojuhelnikova...
jde mne o to, jak spocitat tu matici L
bud je takovy ten zpusob nadrzo, kdy spocitam permutacni matici P a matici U, coz je jednoduche a pak teda prenasobim celou rovnici zprava inverzem U
nebo postup algoritmicky, kdy pocitam postupne pomaloucku L1 az Lk a ruzne to prenasobuju jednotlivyma permutacnima maticema a pak je teda L soucin tech Li - tento postup je ovsem o hodne zdlouhavejsi (tedy nachylnejsi na chyby) ale asi spravnejsi
Janovsky, kdyz nam ukazoval ty zkusebni pisemky, tam presne tenhle priklad mel - otazka teda je, jestli staci pouzit ten prvni postup a clovek je hotovy do tri minut nebo se s tim alespon 20 minut patlat pres ruzne inverni matice atd...
co myslite?