CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH
Transcript of CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH
![Page 1: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/1.jpg)
CHÖÔNG 3
HEÄ PHÖÔNG TRÌNHTUYEÁN TÍNH
![Page 2: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/2.jpg)
I. ÑAËT BAØI TOAÙN :
Heä phöông trình tuyeán tính n pt vaø n aån coù daïng
Ax = b
11 12 1 1 1
21 22 2 2 2
1 2
...
...( )
... ... ... ... ... ...
...
n
n
ij
n n nn n n
a a a x b
a a a x bA a x b
a a a x b
vôùi
![Page 3: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/3.jpg)
Caùc phöông phaùp giaûi
Phöông phaùp giaûi chính xaùc
Phöông phaùp Gauss
Phöông phaùp Gauss-Jordan
Phöông phaùp nhaân töû LU
Phöông phaùp Cholesky
Phöông phaùp giaûi gaàn ñuùng
Phöông phaùp laëp Jacobi
Phöông phaùp laëp Gauss-Seidel
![Page 4: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/4.jpg)
II. PHÖÔNG PHAÙP GAUSS
1. Caùc daïng ma traän ñaëc bieät :
a. Ma traän cheùo :
11
22
0 ... 0
0 ... 0
... ... ... ...
0 0 ... nn
a
aA
a
detA = a11a22 . . . ann 0 aii 0, i
Nghieäm xi = bi/aii
![Page 5: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/5.jpg)
b. Ma traän tam giaùc döôùi
11
21 22
1 2
0 ... 0
... 0
... ... ... ...
...n n nn
a
a aA
a a a
detA = a11a22 . . . ann 0 aii 0, i
Phöông trình coù nghieäm
11
11
1
1
1[ ] , 2,
k
k k kj jjkk
bx
a
x b a x k na
![Page 6: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/6.jpg)
c. Ma traän tam giaùc treân :
11 12 1
22 2
...
0 ...
... ... ... ...
0 0 ...
n
n
nn
a a a
a aA
a
detA = a11a22 . . . ann 0 aii 0, i
Phöông trình coù nghieäm
1
1[ ] , 1, 1
nn
nn
n
k k kj jj kkk
bx
a
x b a x k na
![Page 7: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/7.jpg)
2. Phöông phaùp Gauss :
Ta söû duïng caùc pheùp bieán ñoåi sô caáp theo doøng ñeå chuyeån ma traän A veà ma traân tam giaùc treân
Caùc pheùp bieán ñoåi sô caáp theo doøng
hoaùn chuyeån 2 doøng
nhaân 1 doøng vôùi 1 soá khaùc 0
coäng 1 doøng vôùi doøng khaùc
![Page 8: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/8.jpg)
Ví duï : Giaûi heä phöông trình 1 2 3 4
1 2 3 4
1 2 3
1 2 3 4
2 8
2 2 3 3 20
2
4 3 4
x x x x
x x x x
x x x
x x x x
Giaûi1 1 2 1 8
2 2 3 3 20[ / ]
1 1 1 0 2
1 1 4 3 4
A b
2 3
4 4 /2
1 1 2 1 8
0 2 1 1 6
0 0 1 1 4
0 0 1 2 6
h h
h h
Giaûi pt ma traän tam giaùc treân, ta ñöôïc nghieäm
x = (-7, 3, 2, 2)t
2 2 1
3 3 1
4 4 1
2 1 1 2 1 8
0 0 1 1 4
0 2 1 1 6
0 0 2 4 12
h h h
h h h
h h h
4 4 3
1 1 2 1 8
0 2 1 1 6
0 0 1 1 4
0 0 0 1 2
h h h
![Page 9: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/9.jpg)
III. PHÖÔNG PHAÙP NHAÂN TÖÛ LU
Phaân tích ma traän A thaønh tích 2 ma traän L vaø UA = LU
L : ma traän tam giaùc döôùiU : ma traän tam giaùc treân
Phöông trình Ax = b L(Ux) = b
Ta ñöa veà giaûi 2 heä phöông trình
Ly b
Ux y
![Page 10: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/10.jpg)
Phöông phaùp Doolittle :
Giaû söû A ma traän khoâng suy bieán vaø a11 0 Ta coù theå phaân tích A thaønh
A = LU
21
1 2
1 0 ... 0
1 ... 0
... ... ... ...
... 1n n
lL
l l
11 12 1
22 2
...
0 ...
... ... ... ...
0 0 ...
n
n
nn
u u u
u uU
u
Ma traân döôùi
Ma traân treân
![Page 11: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/11.jpg)
Caùc phaàn töû cuûa L vaø U ñöôïc xaùc ñònh theo coâng thöùc
1 1
11
11
1
1
1
1
, 1
, 2
, 1
1[ ], 1
j j
ii
i
ij ij ik kjk
j
ij ij ik kjkjj
u a j n
al i n
u
u a l u i j
l a l u j iu
![Page 12: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/12.jpg)
Ví duï : Giaûi heä phöông trình
1 2 3
1 2 3
1 2 3
2 2 3 9
4 3 4 15
2 2 3
x x x
x x x
x x x
Giaûi Ta phaân tích
22 23
32 33
2 2 3 1 0 0 2 2 3
4 3 4 2 1 0 0
2 1 2 1 1 0 0
A u u
l u
22 22 21 12
23 23 21 13
32 32 31 12
22
33 33 31 13 32 23
1
2
1( ) 1
3
u a l u
u a l u
l a l uu
u a l u l u
Ví duï : Giaûi heä phöông trình
1 2 3
1 2 3
1 2 3
2 2 3 9
4 3 4 15
2 2 3
x x x
x x x
x x x
![Page 13: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/13.jpg)
Giaûi heä Ly = b
1
2
3
1 0 0 9 9
2 1 0 15 3
1 1 1 3 3
y
y y
y
Giaûi heä Ux = y
1
2
3
2 2 3 9 2
0 1 2 3 1
0 0 3 3 1
x
x x
x
Nghieäm x1 = 2, x2 = 1, x3= -1
![Page 14: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/14.jpg)
TH ñaëc bieät : A ma traän 3 ñöôøng cheùo
11 12
21 22 23
32 33
1
0 ... 0 0
... 0 0
0 ... 0 0
... ... ... ... ... ...
0 0 0 ... nn nn
a a
a a a
A a a
a a
Ta phaân tích A thaønh LU vôùi
11 12
22 2321
3332
0 ... 01 0 0 ... 0
0 ... 01 0 ... 0
0 0 ... 00 1 ... 0
... ... ... ... ...... ... ... ... ...
0 0 0 ...0 0 0 ... 1 nn
u u
u ul
L U ul
u
![Page 15: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/15.jpg)
Caùc phaàn töû cuûa L vaø U ñöôïc xaùc ñònh theo coâng thöùc
2111 11 12 12 21
11
1 1
1 1
1
1
, ,
, 2,
, 2, 1
, 2, 1
ii ii i i i i
i i i i
i i
i i
ii
au a u a l
u
u a l u i n
u a i n
al i n
u
![Page 16: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/16.jpg)
Ví duï : Giaûi heä phöông trình Ax = b
2 1 0 2
1 2 1 1
0 1 2 2
A b
Giaûi Ta phaân tích
22 23
32 33
1 0 0 2 1 0
1 / 2 1 0 0
0 1 0 0
A u u
l u
22 22 21 12
3223 23 32
22
33 33 32 23
3 / 2
1, 2 / 3
4 / 3
u a l u
au a l
u
u a l u
![Page 17: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/17.jpg)
Giaûi heä Ly = b
1
2
3
1 0 0 2 2
1 / 2 1 0 1 2
0 2 / 3 1 2 10 / 3
y
y y
y
Giaûi heä Ux = y
1
2
3
2 1 0 2 5 / 2
0 3 / 2 1 2 3
0 0 4 / 3 10 / 3 5 / 2
x
x x
x
Nghieäm x1 = 5/2, x2 = 3, x3= 5/2
![Page 18: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/18.jpg)
III. PHÖÔNG PHAÙP CHOLESKY
Ñònh nghóa :
Ma traân A goïi laø ñoái xöùng neáuA = At
Ma traân A goïi laø xaùc ñònh döông neáu
1 21 1
0, ( , ,..., ) , 0
n n
t t nij i j n
i j
x Ax a x x x x x x R x
![Page 19: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/19.jpg)
Ñeå kieåm tra xaùc ñònh döông, ta duøng ñình lyù sau:
Ñònh lyù :
Ma traän A laø xaùc ñònh döông khi vaø chæ khi taát caû caùc ñònh thöùc con chính cuûa noù ñeàu döông
Ví duï : Kieåm tra tính xaùc ñònh döông cuûa ma traän1 1 1
1 2 0
1 0 4
A
Giaûi
Caùc ñònh thöùc con chính: 1 2
1 11 0, 1 0
1 2
3
1 1 11 2 1 1 1 1
1 2 0 1 0 4 2 01 0 1 0 1 2
1 0 4
Vaäy A laø xaùc ñònh döông
![Page 20: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/20.jpg)
Ñònh lyù (Cholesky) :
Neáu A ma traän ñoái xöùng vaø xaùc ñònh döông, thì toàn taïi ma traän döôùi, khaû ñaûo B sao cho
A = BBt
Ma traän B = (bij) tìm theo coâng thöùc sau :
11 11
11
11
12
1
1
1
, 2
, 2
1[ ], 2
ii
i
ii ii ikk
j
ij ij ik jkkjj
b a
ab i n
b
b a b i n
b a b b j ib
![Page 21: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/21.jpg)
Ví duï : Giaûi heä phöông trình Ax = b
1 1 1 1
1 2 0 2
1 0 4 3
A b
Giaûi
Ta coù A ma traän ñoái xöùng vaø xaùc ñònh döôngPhaân tích A = BBt
22
32 33
1 0 0
1 0
1
B b
b b
Caùc heä soá
222 22 21
32 32 31 21
22
2 233 33 31 32
1
1[ ] 1
2
b a b
b a b bb
b a b b
![Page 22: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/22.jpg)
Giaûi heä By = b
1
2
3
1 0 0 1 1
1 1 0 2 1
31 1 2 3 / 2
y
y y
y
Giaûi heä Bt x = y
1
2
3
1 1 1 1 3
0 1 1 1 1 / 2
3 / 20 0 2 3 / 2
x
x x
x
Nghieäm x1 = 3 , x2 = -1/2, x3= 3/2
![Page 23: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/23.jpg)
IV. PHÖÔNG PHAÙP LAËP
1. Chuaån :
a. Chuaån vector :
Ñònh nghóa :Chuaån cuûa vector xRn laø haøm soá thöïc kyù hieäu laø ||x||, thoûa 3 ñieàu kieän sau :
(i) ||x||≥0, xRn vaø ||x|| = 0 x=0
(ii) ||x|| = || ||x||, xRn, R
(iii) ||x+y|| ≤||x|| + ||y||, x,yRn
![Page 24: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/24.jpg)
Coù nhieàu coâng thöùc chuaån khaùc nhau, xeùt 2 coâng thöùc
1
11
|| || max {| |}
|| || | |
ix n
n
ii
x x
x x
x= (x1,x2,…, xn)t
Deã daøng kieåm tra ||x||, ||x||1 laø caùc chuaån goïi laø chuaån vaø chuaån 1
![Page 25: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/25.jpg)
b. Chuaån ma traän :
Ñònh nghóa :Chuaån cuûa ma traân A ñöôïc xaùc ñònh theo coâng thöùc
0 || || 1
|| |||| || max max || ||
|| ||
x x
AxA Ax
x
Ñònh lyù : Cho ma traän A = (aij), ta coù
11
|| || max{ | |}
n
iji n
j
A a
11
1
|| || max{ | |}
n
ijj n
i
A a
![Page 26: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/26.jpg)
c. Hoäi tuï theo chuaån :
Ñònh nghóa :
Daõy caùc vector {x(m)}Rn hoäi tuï veà x theo chuaån neáu ||x(m) –x|| 0 khi m
Ñònh lyù :Daõy {x(m)=(x1
(m), x2(m),…, xn
(m) )}Rn hoäi tuï veà x = (x1, x2, …, xn) theo chuaån neáu vaø chæ neáu daõy {xk
(m)}hoäi tuï veà xk khi m, k=1,n
![Page 27: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/27.jpg)
2. Phöông phaùp laëp :
Ta chuyeån heä pt veà daïngx = Tx + c
Vôùi T laø ma traän vuoâng caáp n vaø c laø 1 vector
Ñeå tìm nghieäm gaàn ñuùng, vôùi vector ban ñaàu x(0), ta xaây döïng daõy laëp theo coâng thöùc
x(m) = Tx(m-1)+ c, m=1,2,…
Ta caàn khaûo saùt söï hoäi tuï cuûa daõy {x(m)}
![Page 28: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/28.jpg)
Ta coù ñònh lyù sau
Ñònh lyù :
Neáu ||T|| < 1 thì daõy laëp x(m) seõ hoäi tuï veà nghieäm x cuûa heä pt, vôùi moïi vector ban ñaàu x(0).
Ta coù coâng thöùc ñaùnh giaù sai soá :
( ) ( ) ( 1)|| ||(2) || || || ||
1 || ||m m mTx x x x haäu nghieäm
T
( ) (1) (0)|| ||(1) || || || ||
1 || ||
mm Tx x x x tieàn nghieäm
T
hoaëc
![Page 29: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/29.jpg)
V. PHÖÔNG PHAÙP LAËP JACOBI
Ta phaân tích
A = D - L - U
11
22
0 ... 0
0 ... 0
... ... ... ...
0 0 0nn
a
aD ma traän cheùo
a
21
1 2
0 0 ... 0
0 ... 0
... ... ... ...
... 0n n
aL ma traän döôùi
a a
12 1
2
0 ...
0 0 ...
... ... ... ...
0 0 0 0
n
n
a a
aU ma traän treân
trong ñoù
![Page 30: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/30.jpg)
Phöông trình Ax = b
(D-L-U)x = b
Dx = (L+U)x + b
x = D-1(L+U)x + D-1b
x = Tx + c
vôùi T = D-1(L+U) vaø c = D-1b
pp laëp theo phaân tích treân goïi laø pp laëp Jacobi
Baây giôø ta tìm ñieàu kieän ñeå pp laëp Jacobi HT
![Page 31: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/31.jpg)
Ñònh nghóa :
Ma traän A goïi laø ma traän ñöôøng cheùo troäi nghieâm
ngaët neáu noù thoûa ñieàu kieän sau :
1,
| | | |, 1,n
ij iij j i
a a i n
Nhaän xeùt :
Neáu A laø ma traän ñöôøng cheùo troäi nghieâm ngaët thì
detA 0 vaø aii 0 i=1,n
![Page 32: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/32.jpg)
Ñònh lyù :
Neáu A laø ma traän ñöôøng cheùo troäi nghieâm ngaët, thì pp laëp Jacobi hoäi tuï vôùi moïi giaù trò ban ñaàu x(0)
Ta coù coâng thöùc laëp Jacobi
x(m) = Tx(m-1)+ c, m=1,2,…
1( ) ( 1) ( 1)
1 1
1[ ], 1,
i nm m m
i ij i ij i ij j iii
x a x a x b i na
![Page 33: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/33.jpg)
11
1 22
1 / 0 ... 0
0 1 / ... 0
... ... ... ...
0 0 0 1 /nn
a
aD
a
112 1
11 11 11
2 221
2222 22
1 2
0 ...
0 ...
..... ... ... ...
... 0
n
n
nn n
nnnn nn
aa b
a a a
a ba
aa aT C
ba a
aa a
CM
A ma traân ñöôøng cheùo troäi nghieâm ngaët neân
Ta coù x(m) = Tx(m-1)+ c, T = D-1(L+U) vaø c = D-1b
11,
| ||| || max{ } 1
| |
nij
i nj j i ii
aT
a
pp laëp hoäi tuï
![Page 34: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/34.jpg)
Ví duï : Cho heä phöông trình
1 2 3
1 2 3
1 2 3
10 7
10 8
10 9
x x x
x x x
x x x
a. Tìm nghieäm gaàn ñuùng x(5) vôùi vector ban ñaàu x(0) = 0
b. Tính ma traän T vaø c
c. Tính sai soá cuûa nghieäm x(5)
Ví duï : Cho heä phöông trình
1 2 3
1 2 3
1 2 3
10 7
10 8
10 9
x x x
x x x
x x x
![Page 35: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/35.jpg)
a.10 1 1
1 10 1
1 1 10
A ma traän ñöôøng cheùo troäi nghieâm ngaët
Coâng thöùc laëp Jacobi( ) ( 1) ( 1)
1 2 3
( ) ( 1) ( 1)
2 1 3
( ) ( 1) ( 1)
3 1 2
1( 7)
10
1( 8)
10
1( 9)
10
m m m
m m m
m m m
x x x
x x x
x x x
m 0 1 2 3 4 5
x1(m) 0 0.7 0.71 0.725 0.7267 0.72717
x2(m) 0 0.8 0.64 0.640 0.6368 0.63648
x3(m) 0 0.9 0.89 0.907 0.9085 0.90899
![Page 36: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/36.jpg)
b. Ta coù 0 0.1 0.1 0.7
0.1 0 0.1 0.8
0.1 0.1 0 0.9
T c
c. Coâng thöùc sai soá
(5) (5) (4)|| |||| || || ||
1 || ||
Tx x x x
T
Ta coù ||T||=0.2, neân
(5) 4 30.2|| || 4.9*10 0.1225*10
0.8 x x
![Page 37: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/37.jpg)
VI. Phöông phaùp laëp Gauss-Seidel : Ta phaân tích
A = D - L – Unhö trong phaàn tröôùc
Phöông trình Ax = b
(D-L-U)x = b
(D-L)x = Ux + b
x = (D-L)-1Ux + (D-L)-1b
x = Tx + c
vôùi T = (D-L)-1U vaø c = (D-L)-1b
pp laëp theo phaân tích naøy goïi laø pp laëp Gauss-Seidel
![Page 38: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/38.jpg)
Ñònh lyù :
Neáu A laø ma traän ñöôøng cheùo troäi nghieâm ngaët, thì pp laëp Gauss-Seidel hoäi tuï vôùi moïi giaù trò ban ñaàu x(0)
Ta coù coâng thöùc laëp Gauss-Seidel
x(m) = Tx(m-1)+ c, m=1,2,…1
( ) ( ) ( 1)
1 1
1[ ], 1,m
i nm m
i ij i ij i ij j iii
x a x a x b i na
![Page 39: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/39.jpg)
Ví duï : Cho heä phöông trình
a. Tìm nghieäm gaàn ñuùng x(4) vôùi vector ban ñaàu x(0) = 0
b. Tính ma traän T vaø c
c. Tính sai soá cuûa nghieäm x(4)
Ví duï : Cho heä phöông trình
1 2 3
1 2 3
1 2 3
20 2 12
20 13
2 20 14
x x x
x x x
x x x
![Page 40: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/40.jpg)
a.20 1 2
1 20 1
2 1 20
A ma traän ñöôøng cheùo troäi nghieâm ngaët
Coâng thöùc laëp Gauss-Seidel( ) ( 1) ( 1)
1 2 3
( ) ( ) ( 1)
2 1 3
( ) ( ) ( )
3 1 2
1( 2 12)
20
1( 13)
20
1(2 14)
20
m m m
m m m
m m m
x x x
x x x
x x x
m 0 1 2 3 4
x1(m) 0 0.6 0.5519 0.554268975 0.554233852
x2(m) 0 0.62 0.661955 0.661700938 0.661713904
x3(m) 0 0.791 0.78828775 0.788511944 0.788509080
![Page 41: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/41.jpg)
b. Ta coù
1
20 0 0 0.05 0 0
1 20 0 ( ) 0.0025 0.05 0
2 1 20 0.004875 0.0025 0.05
D L D L
1
0 1 2 0 0.05 0.1
0 0 1 ( ) 0 0.0025 0.055
0 0 0 0 0.004875 0.00725
U T D L U
1
0.6
( ) 0.62
0.791
c D L b
![Page 42: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/42.jpg)
(4) 5 50.15|| || 3.5123*10 0.6199*10
0.85 x x
Ta coù ||T||=0.15, neân
c. Coâng thöùc sai soá
(4) (4) (3)|| |||| || || ||
1 || ||
Tx x x x
T
![Page 43: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/43.jpg)
VII. Heä pt oån ñònh vaø soá ñieàu kieän :
Xeùt heä phöông trình Ax = b
Ñònh nghóa :
Heä phöông trình goïi laø oån ñònh neáu moïi thay ñoåi nhoû cuûa A hay b thì nghieäm cuûa heä chæ thay ñoåi nhoû
1. Heä pt oån ñònh :
![Page 44: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/44.jpg)
Ví duï : Xeùt heä phöông trình Ax = b vôùi
1 2 3
1 2.01 3.01A b
Heä phöông trình coù nghieäm x = (1, 1)T
Thay ñoåi 3
3.1b
Nghieäm cuûa heä : x=(-17, 10)T
Ta thaáy nghieäm cuûa heä khaùc raát xa khi b thay ñoåi nhoû. Vaäy heä khoâng oån ñònh
![Page 45: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/45.jpg)
Ví duï : Xeùt heä phöông trình Ax = b vôùi10 7 8 7 32
7 5 6 5 23
8 6 10 9 33
7 5 9 10 31
A b
Heä coù nghieäm x = (1, 1, 1, 1)T
Nghieäm cuûa heä : x=(-81, 137, -34, 22)T
Ta thaáy nghieäm cuûa heä khaùc raát xa khi A thay ñoåi nhoû. Vaäy heä khoâng oån ñònh
Thay ñoåi A moät ít
10 7 8.1 7.2
7.08 5.04 6 5
8 5.98 9.98 9
6.99 4.99 9 9.98
A
![Page 46: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/46.jpg)
2. Soá ñieàu kieän :
Ta tìm ñieàu kieän ñeå heä oån ñònh
Ñònh nghóa : Soá
k(A) = ||A|| ||A||-1
Goïi laø soá ñieàu kieän cuûa ma traän A
Ta coù caùc tính chaát :
( ) 1 ( )
|| || || ||( ) ( )
|| || || ||
|| || || ||( ) ( )
|| || || ||
i k A
x bii k A
x b
x Aiii k A
x x A A
![Page 47: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/47.jpg)
Nhaän xeùt :
Soá ñieàu kieän cuûa ma traän ñaëc tröng cho tính oån ñònh cuûa heä phöông trình
k(A) caøng gaàn 1 thì heä caøng oån ñònh
k(A) caøng xa 1 thì heä caøng khoâng oån ñònh
Ví duï : 1 2
1 2.01A
Ta coù 1 201 200
100 100A
k(A) = 3.01 x 401 = 1207.01 >> 1
Vaäy heä khoâng oån ñònh
![Page 48: CHÖÔNG 3 HEÄ PHÖÔNG TRÌNH TUYEÁN TÍNH](https://reader035.fdocument.pub/reader035/viewer/2022081503/62a22b00b5588a440408c095/html5/thumbnails/48.jpg)
Ví duï : 3 1 1
1 2 1
1 1 4
A
Ta coù1
7 /13 5 /13 3 /13
5 /13 11/13 4 /13
3 /13 4 /13 5 /13
A
k(A) = 6 x 20/13 = 9.2308 >> 1
Vaäy heä khoâng oån ñònh