Matematické Fórum

KateřinaDardová · 29. 10. 2012 14:57

Máme za úkol napsat algoritmus LU rozkladu v programu wxMaxima a nevím si s tím rady. Na internetu jsem našla:

matrix([1,2],[3,4])
lu_factor(%);
matrix([1,2],[3,-2]),[1,2],generalring
get_lu_factors(%);
[matrix([1,0],[0,1]),matrix([1,0],[3,1]),matrix([1,2],[0,-2])]

, ale moc se mi to nezdá, sama jsem si načrtla toto:

VSTUP: matice A (respektive LU), stopkrok(kolikrát se bude algoritmus opakovat)
$L_{0}=L$
$U_{0}=U$
cyklus i=1,...,stopkrok
$A_{i}=U_{i-1}\cdot L_{i-1}$
rozklad $A_{i}$ na $L_{i}\cdot U_{i}$
$A_{stopkrok}$
VÝSTUP: vlastní čísla

Jen nevím jak to přepsat do Maximy...

Matematické Fórum

#1 29. 10. 2012 14:57

algoritmus LU rozkladu wxMaxima