Matematické Fórum

↑ frenkiss:

 
L=tril(A);% L .... dolni trojuhelniková
U=triu(A);% U .... horni trojuhelniková
[L,U] = lu(A);
GeaA=[u]

Ještě ti to můžu dát v C-čku.

↑ gladiator01:
takze nebude fungovat to asi ,vsak?

Mně to teda počítá správně (ještě jsem to kontrolovala maplem jestli jsou kořeny opravdu správně a jsou):

format compact
%02_J: definujte matici 4.řádu a algoritmem GEA ji převeďte na schodovitý tvar
Kom1=['02_J: definujte matici 4.řádu a algoritmem GEA ji převeďte na schodovitý tvar';
      '-----------------------------------------------------------------------------']
Kom2=['----------------------';
      'Matice A              ';
      '----------------------']
A=[5 2 3 4 1;  0 14 15 16 5; 3 6 7 8 3; 1 10 5 12 7;];

Kom3=['---------------------------------------------';
      'Matice A převedená pomocí GEA                ';
      '---------------------------------------------']
L=tril(A);% L .... dolni trojuhelniková
U=triu(A);% U .... horni trojuhelniková
[L,U] = lu(A);
GeaA=[ U]

x=GeaA(4,5)/GeaA(4,4)
y=(GeaA(3,5)-x*GeaA(3,4))/GeaA(3,3)
z=(GeaA(2,5)-y*GeaA(2,3)-x*GeaA(2,4))/GeaA(2,2)
a=(GeaA(1,5)-z*GeaA(1,2)-y*GeaA(1,3)-x*GeaA(1,4))/GeaA(1,1)

↑ frenkiss:
Tam má být Velké U ne malé. Edit. ne mod Code, ale hranát závorka, to udělala.

A máš to přesně jako já? Zkopíruj to tam znova: lu (dohormady) je malýma, Písmena L a U (ty zvlášť) jsou velkýma.

format compact
A=[5 2 3 4 1;  0 14 15 16 5; 3 6 7 8 3; 1 10 5 12 7;];

Kom3=['---------------------------------------------';
      'Matice A převedená pomocí GEA                ';
      '---------------------------------------------']
L=tril(A);% L .... dolni trojuhelniková
U=triu(A);% U .... horni trojuhelniková
[L,U] = lu(A);
GeaA=[ U]

x=GeaA(4,5)/GeaA(4,4)
y=(GeaA(3,5)-x*GeaA(3,4))/GeaA(3,3)
z=(GeaA(2,5)-y*GeaA(2,3)-x*GeaA(2,4))/GeaA(2,2)
a=(GeaA(1,5)-z*GeaA(1,2)-y*GeaA(1,3)-x*GeaA(1,4))/GeaA(1,1)

↑ frenkiss:
Nevím jak to počítá matlab - myslím s částečným výběrem pivota (postupuješ přesně po řádcích a vždycky první (nenulové) číslo v řádku si označíš jako pivot a vynuluješ čísla pod ním, atd. (Prostě jdeš systematicky a ne na přeskáčku.)

Když to budeš počítat ručně, tak ti asi nevyjdou přesně tatáž čísla, protože způsobů jak nakombinovat řádky (aby ti vyšly nuly pod diagonálou) je mnoho. Kořeny budou samozřejmě pořád stejné. (třeba prohodíš dva řádku a hned máš trochu jiná čísla).