Stránky: 1
↑ noacoked:
2. Newton method
3. Jednou z podmienok (ako som sa dočítal na Wikipédii) je, že matica musí byť diagonálne dominantná, čo tvoja nie je, druhá podmienka je, že spektrálny polomer
musí byť menší ako 1.
4. Neviem čo je NIP, LIP.
Offline
2. 
potřebujeme derivaci.. 
potřebujeme odhad výsledku, např.
(můžeš zvolit cokoli..)
první iterace 
druhá iterace
třetí iterace
čtvrtá iterace
výsledek už se nemění, takže jeden z kořenů bude přibližně 1,16.
Offline
Najít LU rozklad matice 3x3 by ti nemělo dát moc práce
1 0 0
L=3/5 1 0
4/5 8/21 1
5 3 4
U= 0 21/5 8/5
0 0 25/21
dál se řeší soustava
Ly=b (b je pravá strana, tedy (12,13,13))
z toho snadno 
a nakonec se řeší soustava
Ux=y
výsledek x=(1,1,1)
A samozřejmě, že to jde řešit jak Jakobim tak Gauss-Seidlem.
LU má výhodu pokud máme více pravých stran a pro ně jednu společnou matici A. Pak je výpočet přes LU rychlejší než počítat každou soustavu zvlášť.
Offline
Offline
↑ noacoked:Já nevím, jak jste to brali, ale já nikdy žádný mezikrok pro výpočet LU rozkladu neměla. Vše probíhá v hlavě.
L je dolní trojúhelník s jedničkami na diagonále. U je horní trojúhelník.
Násobí se mezi sebou, takže se dopočítávají postupně jejich prvky tak, aby vyšla původní matice. Algoritmus najdeš kdekoli, ten ti sem fakt psát nebudu.
A řešení těch dvou malých soustav si určitě zvládneš sepsat sám.
Offline
tak si to napiš, jednu matici s nulami nad diagonálou a jedničkami na diagonále a druhou s nualmi pod diagonálou a zkus je takhle prázdné násobit. Hned zjištíš, co kam máš napsat, aby ti po násobení vznikla původní matice.
Offline
Offline
↑ noacoked:Pronásob je s ypsilony a sečti..viz postup např. zde
Offline
Offline
Stránky: 1