數(shù)值分析課件第四章_線性方程組的迭代解法

1、第四章第四章 線性方程組的迭代解法線性方程組的迭代解法第四章第四章 線性方程組的迭代解法線性方程組的迭代解法引言引言4.1 向量和矩陣的范數(shù)向量和矩陣的范數(shù)向量的范數(shù)向量的范數(shù)222212|nxxxx證明（2）: 22211max |max |kkk nk nxxnx 2|xxnx所以|211nxxxx 證明（1）: |max|max111knkknkxnxx |1xnxx所以222212|nxxxx212|xxnx所以|211nxxxx 證明（3）: 222121|xxxx2222221121232221|()()()()()0nnnnn xxxxxxxxxxxx矩陣的范數(shù)矩陣的范數(shù)基礎(chǔ)知識(shí)

2、回顧基礎(chǔ)知識(shí)回顧Axx()0AI x0AI11AAA迭代法的基本思想：迭代法的基本思想： 3612363311420238321321321xxxxxxxxx )3636(121)334(111)2023(81213312321xxxxxxxxx 01231261110114828300B 12361133820f任取初始值，如取任取初始值，如取x(0)=(0,0,0)T，代入，代入x=B0 x+f右邊，右邊，若等式成立則求得方程組的解，否則，得新值若等式成立則求得方程組的解，否則，得新值x(1)=(x1(1),x2(1),x3(1)T=(2.5,3,3)T，再將，再將x(1)代入，反代入，反

3、復(fù)計(jì)算，得一向量序列復(fù)計(jì)算，得一向量序列x(k)和一般的計(jì)算公式和一般的計(jì)算公式（迭代公式）：（迭代公式）： )3636(121)334(111)2023(81)(2)(1)1(3)(3)(1)1(2)(3)(2)1(1kkkkkkkkkxxxxxxxxx4.2 Jacobi迭代法迭代法4.3 GaussSeidel迭代法迭代法4.4 迭代法的收斂性迭代法的收斂性fBxxkk )()1(fBxx *)(*)()1(xxBxxkk *)()(xxkk 令令,2, 1 ,0 k)()1(kkB )1(2kB )0(1 kB為為非非零零常常數(shù)數(shù)向向量量注注意意*)0()0(xx 1lim0kkB(0

4、)(0)kkkBB*limkkxxB*kxx|1|*|)0()1()(xxBBxxkk |1|*|)1()()( kkkxxBBxx（1）（2）*lim)(xxkk |11|*|)()1()(kkkxxBxx |1|*|)1()()( kkkxxBBxx|1|*|)0()1()(xxBBxxkk |1|*|)0()1()(xxBBxxkk |1|*|)1()()( kkkxxBBxxu（2）式）式說(shuō)明，說(shuō)明，|B|越小，越小， x(k) 收斂越快，可作誤差收斂越快，可作誤差估計(jì)式。估計(jì)式。fBxxkk )()1(nnrJJJJ 21iinniiiJ 110lim kkBlim0kkJlim0k

5、ik1| i 1|max1 iri 111122111221321xxx)(1ULDBJ 100010001 022101220 022101220)det(JBI 221122det3 0 0 10|)max(|)( JBULDBG1)( 1122011001 000100220 GB 2003202200 2 12|)max(|)( GB nijjijiiaa0niaaijijii, 3 , 2 , 1| niaaijijii, 3 , 2 , 11|1 0002122222211111112nnnnnnnnJaaaaaaaaaaaaB JB ijijiiiaa|1max1 )(1ULDB

6、J 0)det( GBI 0)(det1 ULDI 0)(det)det(1 ULDLD 0)(det ULD ULDBG1)( ijijiiaa| nijijijijiiaaa111| nijijnijijijijaaa1111|)1|(| 則則有有如如果果, 1| nijijijijiiaaa111| 為為嚴(yán)嚴(yán)格格對(duì)對(duì)角角占占優(yōu)優(yōu)矩矩陣陣則則)(ULD 0)(det ULD 從從而而, 1| 所以所以, 1)( GB 即即4.5 松弛迭代法松弛迭代法bLDUxLDxkk1)(1)1()()( GS迭代法GS法為SOR法的特例, SOR法為G-S法的加速。45. 1 取取 321242124 321xxx 320ULDBSG1)( 1321042004 000200120 5 . 03/10625. 025. 0025. 05 . 00bLDfSG1)( 1321042004 320 3/25 . 00Tx)1 , 1 , 1()0( 取取初初值值bLDxUDLDxkk1)(1)1()()1()( 45. 1 取取4.6 誤差分析誤差分析 220001. 111121xx 0001. 220001. 111121xxbbAAxx 1AAAAxxx 1A-1A-1 A 1211111131211211nnnnnHn 5/