Bài giảng Phương pháp tính - Chương 4: Hệ phương trình tuyến tính - Nguyễn Thị Cẩm Vân

HỆ PHƯƠNG TRÌNH TUYẾN TÍNH BÀI GIẢNG ĐIỆN TỬ Nguyễn Thị Cẩm Vân Trường Đại học Bách Khoa TPHCM Khoa Khoa học ứng dụng, bộmôn Toán ứng dụng Ngày 12 tháng 2 năm 2018Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 1 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt NỘI DUNG BÀI HỌC 1 ĐẶT VẤN ĐỀ 2 PHƯƠNG PHÁP GAUSS 3 PHƯƠNG PHÁP NHÂN TỬ LU 4 PHƯƠNG PHÁP CHOLESKI 5 CHUẨN CỦA VÉCTƠ, CHUẨN CỦA MA TRẬN 6 NHỮNG PHƯƠNG PHÁP LẶP Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 2 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt NỘI DUNG BÀI HỌC 1 ĐẶT VẤN ĐỀ 2 PHƯƠNG PHÁP GAUSS 3 PHƯƠNG PHÁP NHÂN TỬ LU 4 PHƯƠNG PHÁP CHOLESKI 5 CHUẨN CỦA VÉCTƠ, CHUẨN CỦA MA TRẬN 6 NHỮNG PHƯƠNG PHÁP LẶP Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 2 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt NỘI DUNG BÀI HỌC 1 ĐẶT VẤN ĐỀ 2 PHƯƠNG PHÁP GAUSS 3 PHƯƠNG PHÁP NHÂN TỬ LU 4 PHƯƠNG PHÁP CHOLESKI 5 CHUẨN CỦA VÉCTƠ, CHUẨN CỦA MA TRẬN 6 NHỮNG PHƯƠNG PHÁP LẶP Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 2 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt NỘI DUNG BÀI HỌC 1 ĐẶT VẤN ĐỀ 2 PHƯƠNG PHÁP GAUSS 3 PHƯƠNG PHÁP NHÂN TỬ LU 4 PHƯƠNG PHÁP CHOLESKI 5 CHUẨN CỦA VÉCTƠ, CHUẨN CỦA MA TRẬN 6 NHỮNG PHƯƠNG PHÁP LẶP Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 2 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt NỘI DUNG BÀI HỌC 1 ĐẶT VẤN ĐỀ 2 PHƯƠNG PHÁP GAUSS 3 PHƯƠNG PHÁP NHÂN TỬ LU 4 PHƯƠNG PHÁP CHOLESKI 5 CHUẨN CỦA VÉCTƠ, CHUẨN CỦA MA TRẬN 6 NHỮNG PHƯƠNG PHÁP LẶP Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 2 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt NỘI DUNG BÀI HỌC 1 ĐẶT VẤN ĐỀ 2 PHƯƠNG PHÁP GAUSS 3 PHƯƠNG PHÁP NHÂN TỬ LU 4 PHƯƠNG PHÁP CHOLESKI 5 CHUẨN CỦA VÉCTƠ, CHUẨN CỦA MA TRẬN 6 NHỮNG PHƯƠNG PHÁP LẶP Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 2 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Đặt vấn đề ĐẶT VẤN ĐỀ Hệ phương trình đại số tuyến tính a11x1+a12x2+ . . .+a1ixi + . . .+a1nxn = b1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ai1x1+ai2x2+ . . .+ai ixi + . . .+ainxn = bi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . an1x1+an2x2+ . . .+anixi + . . .+annxn = bn (1) thường xuất hiện trong các bài toán kỹ thuật. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 3 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Đặt vấn đề ĐẶT VẤN ĐỀ Hệ phương trình đại số tuyến tính a11x1+a12x2+ . . .+a1ixi + . . .+a1nxn = b1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ai1x1+ai2x2+ . . .+ai ixi + . . .+ainxn = bi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . an1x1+an2x2+ . . .+anixi + . . .+annxn = bn (1) thường xuất hiện trong các bài toán kỹ thuật. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 3 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Đặt vấn đề 1 Ta chỉ xét hệ gồm n phương trình và n ẩn số, trong đó A = (ai j ) ∈Mn(K ) và det A 6= 0. Do đó hệ sẽ có nghiệm duy nhất X = A−1B. 2 Tuy nhiên, việc tìmma trận nghịch đảo A−1 đôi khi còn khó khăn gấp nhiều lần so với việc giải trực tiếp hệ phương trình (1). Do đó cần phải có phương pháp để giải hệ (1) hiệu quả. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 4 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Đặt vấn đề 1 Ta chỉ xét hệ gồm n phương trình và n ẩn số, trong đó A = (ai j ) ∈Mn(K ) và det A 6= 0. Do đó hệ sẽ có nghiệm duy nhất X = A−1B. 2 Tuy nhiên, việc tìmma trận nghịch đảo A−1 đôi khi còn khó khăn gấp nhiều lần so với việc giải trực tiếp hệ phương trình (1). Do đó cần phải có phương pháp để giải hệ (1) hiệu quả. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 4 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Hệ phương trình tương đương SỬ DỤNG PHÉP BIẾN ĐỔI SƠ CẤP TRÊN HÀNG ĐỂ GIẢI HỆ Xét hệ phương trình tuyến tính gồm n phương trình và n ẩn a11x1+a12x2+ . . .+a1 jx j + . . .+a1nxn = b1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ai1x1+ai2x2+ . . .+ai jx j + . . .+ainxn = bi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . an1x1+an2x2+ . . .+an jx j + . . .+annxn = bn Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 5 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Hệ phương trình tương đương SỬ DỤNG PHÉP BIẾN ĐỔI SƠ CẤP TRÊN HÀNG ĐỂ GIẢI HỆ Xét hệ phương trình tuyến tính gồm n phương trình và n ẩn a11x1+a12x2+ . . .+a1 jx j + . . .+a1nxn = b1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ai1x1+ai2x2+ . . .+ai jx j + . . .+ainxn = bi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . an1x1+an2x2+ . . .+an jx j + . . .+annxn = bn Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 5 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Hệ phương trình tương đương Nếu thực hiện các phép biến đổi sơ cấp sau trên hệ (1): 1 Đổi chỗ các phương trình của hệ (hi ↔ h j ) hay ci ↔ c j có đánh số lại các ẩn. 2 Nhân vào một phương trình của hệ một số λ 6= 0(hi →λhi ). 3 Cộng vào một phương trình của hệ một phương trình khác đã được nhân với một số (hi → hi +λh j ) thì ta sẽ được một hệ phương trình mới tương đương với hệ (1). Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 6 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Hệ phương trình tương đương Nếu thực hiện các phép biến đổi sơ cấp sau trên hệ (1): 1 Đổi chỗ các phương trình của hệ (hi ↔ h j ) hay ci ↔ c j có đánh số lại các ẩn. 2 Nhân vào một phương trình của hệ một số λ 6= 0(hi →λhi ). 3 Cộng vào một phương trình của hệ một phương trình khác đã được nhân với một số (hi → hi +λh j ) thì ta sẽ được một hệ phương trình mới tương đương với hệ (1). Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 6 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Hệ phương trình tương đương Nếu thực hiện các phép biến đổi sơ cấp sau trên hệ (1): 1 Đổi chỗ các phương trình của hệ (hi ↔ h j ) hay ci ↔ c j có đánh số lại các ẩn. 2 Nhân vào một phương trình của hệ một số λ 6= 0(hi →λhi ). 3 Cộng vào một phương trình của hệ một phương trình khác đã được nhân với một số (hi → hi +λh j ) thì ta sẽ được một hệ phương trình mới tương đương với hệ (1). Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 6 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Hệ phương trình tương đương Nếu thực hiện các phép biến đổi sơ cấp sau trên hệ (1): 1 Đổi chỗ các phương trình của hệ (hi ↔ h j ) hay ci ↔ c j có đánh số lại các ẩn. 2 Nhân vào một phương trình của hệ một số λ 6= 0(hi →λhi ). 3 Cộng vào một phương trình của hệ một phương trình khác đã được nhân với một số (hi → hi +λh j ) thì ta sẽ được một hệ phương trình mới tương đương với hệ (1). Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 6 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Hệ phương trình tương đương a11 a12 . . . a1n a21 a22 . . . a2n . . . . . . . . . . . . an1 an2 . . . ann ∣∣∣∣∣∣∣∣∣ b1 b2 . . . bn  BĐ sơ cấp trên hàng−−−−−−−−−−−−−−−−−−−→  c11 c12 . . . c1n 0 c22 . . . c2n . . . . . . . . . . . . 0 0 . . . cnn ∣∣∣∣∣∣∣∣∣ d1 d2 . . . dn  với ci i 6= 0, i = 1,2, . . . ,n. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 7 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss PHƯƠNG PHÁP GAUSS 1 Viết ma trận mở rộng AB = (A|B) của hệ (1). 2 Dùng các phép biến đổi sơ cấp trên hàng biến đổi ma trận mở rộng về ma trận bậc thang. 3 Viết hệ phương trình tương ứng với ma trận bậc thang. 4 Ta giải hệ phương trình ngược từ dưới lên, tìm biến xn sau đó xn−1, . . . ,x1 ta được 1 nghiệm duy nhất. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 8 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss PHƯƠNG PHÁP GAUSS 1 Viết ma trận mở rộng AB = (A|B) của hệ (1). 2 Dùng các phép biến đổi sơ cấp trên hàng biến đổi ma trận mở rộng về ma trận bậc thang. 3 Viết hệ phương trình tương ứng với ma trận bậc thang. 4 Ta giải hệ phương trình ngược từ dưới lên, tìm biến xn sau đó xn−1, . . . ,x1 ta được 1 nghiệm duy nhất. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 8 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss PHƯƠNG PHÁP GAUSS 1 Viết ma trận mở rộng AB = (A|B) của hệ (1). 2 Dùng các phép biến đổi sơ cấp trên hàng biến đổi ma trận mở rộng về ma trận bậc thang. 3 Viết hệ phương trình tương ứng với ma trận bậc thang. 4 Ta giải hệ phương trình ngược từ dưới lên, tìm biến xn sau đó xn−1, . . . ,x1 ta được 1 nghiệm duy nhất. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 8 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss PHƯƠNG PHÁP GAUSS 1 Viết ma trận mở rộng AB = (A|B) của hệ (1). 2 Dùng các phép biến đổi sơ cấp trên hàng biến đổi ma trận mở rộng về ma trận bậc thang. 3 Viết hệ phương trình tương ứng với ma trận bậc thang. 4 Ta giải hệ phương trình ngược từ dưới lên, tìm biến xn sau đó xn−1, . . . ,x1 ta được 1 nghiệm duy nhất. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 8 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss VÍ DỤ 2.1 Giải hệ phương trình x1+2x2+3x3+4x4 = 7 2x1+x2+2x3+3x4 = 6 3x1+2x2+x3+2x4 = 7 4x1+3x2+2x3+x4 = 18 Giải.  1 2 3 4 2 1 2 3 3 2 1 2 4 3 2 1 ∣∣∣∣∣∣∣∣∣ 7 6 7 18  h2→h2−2h1 h3→h3−3h1 h4→h4−4h1−−−−−−−→ Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 9 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss VÍ DỤ 2.1 Giải hệ phương trình x1+2x2+3x3+4x4 = 7 2x1+x2+2x3+3x4 = 6 3x1+2x2+x3+2x4 = 7 4x1+3x2+2x3+x4 = 18 Giải.  1 2 3 4 2 1 2 3 3 2 1 2 4 3 2 1 ∣∣∣∣∣∣∣∣∣ 7 6 7 18  h2→h2−2h1 h3→h3−3h1 h4→h4−4h1−−−−−−−→ Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 9 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss  1 2 3 4 0 −3 −4 −5 0 −4 −8 −10 0 −5 −10 −15 ∣∣∣∣∣∣∣∣∣ 7 −8 −14 −10  h2→h2−h3−−−−−−→  1 2 3 4 0 1 4 5 0 −4 −8 −10 0 −5 −10 −15 ∣∣∣∣∣∣∣∣∣ 7 6 −14 −10  h3→h3+4h2 h4→h4+5h2−−−−−−−→ Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 10 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss  1 2 3 4 0 −3 −4 −5 0 −4 −8 −10 0 −5 −10 −15 ∣∣∣∣∣∣∣∣∣ 7 −8 −14 −10  h2→h2−h3−−−−−−→  1 2 3 4 0 1 4 5 0 −4 −8 −10 0 −5 −10 −15 ∣∣∣∣∣∣∣∣∣ 7 6 −14 −10  h3→h3+4h2 h4→h4+5h2−−−−−−−→ Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 10 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss  1 2 3 4 0 1 4 5 0 0 8 10 0 0 10 10 ∣∣∣∣∣∣∣∣∣ 7 6 10 20  h3↔h4 h3→ 110h3−−−−−→  1 2 3 4 0 1 4 5 0 0 1 1 0 0 8 10 ∣∣∣∣∣∣∣∣∣ 7 6 2 10  h4→h4−8h3−−−−−−−→ Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 11 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss  1 2 3 4 0 1 4 5 0 0 8 10 0 0 10 10 ∣∣∣∣∣∣∣∣∣ 7 6 10 20  h3↔h4 h3→ 110h3−−−−−→  1 2 3 4 0 1 4 5 0 0 1 1 0 0 8 10 ∣∣∣∣∣∣∣∣∣ 7 6 2 10  h4→h4−8h3−−−−−−−→ Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 11 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss 1 2 3 4 0 1 4 5 0 0 1 1 0 0 0 2 ∣∣∣∣∣∣∣∣∣ 7 6 2 −6  . Vậy hệ đã cho tương đương với hệ sau x1+2x2+3x3+4x4 = 7 x2+4x3+5x4 = 6 x3+x4 = 2 2x4 = −6 ⇔  x1 = 2 x2 = 1 x3 = 5 x4 = −3 Suy ra hệ đã cho có 1 nghiệm duy nhất (x1,x2,x3,x4)= (2,1,5,−3) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 12 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss 1 2 3 4 0 1 4 5 0 0 1 1 0 0 0 2 ∣∣∣∣∣∣∣∣∣ 7 6 2 −6  . Vậy hệ đã cho tương đương với hệ sau x1+2x2+3x3+4x4 = 7 x2+4x3+5x4 = 6 x3+x4 = 2 2x4 = −6 ⇔  x1 = 2 x2 = 1 x3 = 5 x4 = −3 Suy ra hệ đã cho có 1 nghiệm duy nhất (x1,x2,x3,x4)= (2,1,5,−3) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 12 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss 1 2 3 4 0 1 4 5 0 0 1 1 0 0 0 2 ∣∣∣∣∣∣∣∣∣ 7 6 2 −6  . Vậy hệ đã cho tương đương với hệ sau x1+2x2+3x3+4x4 = 7 x2+4x3+5x4 = 6 x3+x4 = 2 2x4 = −6 ⇔  x1 = 2 x2 = 1 x3 = 5 x4 = −3 Suy ra hệ đã cho có 1 nghiệm duy nhất (x1,x2,x3,x4)= (2,1,5,−3) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 12 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss 1 2 3 4 0 1 4 5 0 0 1 1 0 0 0 2 ∣∣∣∣∣∣∣∣∣ 7 6 2 −6  . Vậy hệ đã cho tương đương với hệ sau x1+2x2+3x3+4x4 = 7 x2+4x3+5x4 = 6 x3+x4 = 2 2x4 = −6 ⇔  x1 = 2 x2 = 1 x3 = 5 x4 = −3 Suy ra hệ đã cho có 1 nghiệm duy nhất (x1,x2,x3,x4)= (2,1,5,−3) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 12 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan PHƯƠNG PHÁP GAUSS-JORDAN ĐỊNH NGHĨA 2.1 Phần tử trội là phần tử có trị tuyệt đối lớn nhất, sao cho không cùng hàng và cột với những phần tử đã chọn trước. Phương pháp Gauss-Jordan 1 Chọn phần tử trội để biến đổi cho tất cả các phần tử trên cùng cột của phần tử trội bằng không. 2 Qua n bước ta sẽ tìm được nghiệm. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 13 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan PHƯƠNG PHÁP GAUSS-JORDAN ĐỊNH NGHĨA 2.1 Phần tử trội là phần tử có trị tuyệt đối lớn nhất, sao cho không cùng hàng và cột với những phần tử đã chọn trước. Phương pháp Gauss-Jordan 1 Chọn phần tử trội để biến đổi cho tất cả các phần tử trên cùng cột của phần tử trội bằng không. 2 Qua n bước ta sẽ tìm được nghiệm. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 13 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan VÍ DỤ 2.2 Giải hệ phương trình x1−x2+2x3−x4 = −8 2x1−2x2+3x3−3x4 = −20 x1+x2+x3+0x4 = −2 x1−x2+4x3+3x4 = 4 Giải.  1 −1 2 −1 2 −2 3 −3 1 1 1 0 1 −1 4 3 ∣∣∣∣∣∣∣∣∣ −8 −20 −2 4  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 14 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan VÍ DỤ 2.2 Giải hệ phương trình x1−x2+2x3−x4 = −8 2x1−2x2+3x3−3x4 = −20 x1+x2+x3+0x4 = −2 x1−x2+4x3+3x4 = 4 Giải.  1 −1 2 −1 2 −2 3 −3 1 1 1 0 1 −1 4 3 ∣∣∣∣∣∣∣∣∣ −8 −20 −2 4  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 14 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan Chọn phần tử trội là a43 = 4. Thực hiện các phép biến đổi sơ cấp h3→4h3−h4 h2→4h2−3h4 h1→2h1−h4−−−−−−−−→  1 −1 0 −5 5 −5 0 −21 3 5 0 −3 1 −1 4 3 ∣∣∣∣∣∣∣∣∣ −20 −92 −12 4  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 15 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan Chọn phần tử trội là a43 = 4. Thực hiện các phép biến đổi sơ cấp h3→4h3−h4 h2→4h2−3h4 h1→2h1−h4−−−−−−−−→  1 −1 0 −5 5 −5 0 −21 3 5 0 −3 1 −1 4 3 ∣∣∣∣∣∣∣∣∣ −20 −92 −12 4  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 15 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan Chọn phần tử trội không được nằm trên hàng 4 và cột 3 là phần tử a24 =−21. Thực hiện các phép biến đổi sơ cấp h1→21h1−5h2 h3→7h3−h2 h4→7h4+h2−−−−−−−−−→  −4 4 0 0 5 −5 0 −21 16 40 0 0 12 −12 28 0 ∣∣∣∣∣∣∣∣∣ 40 −92 8 −64  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 16 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan Chọn phần tử trội không được nằm trên hàng 4,2 và cột 3,4 là phần tử a32 = 40. Thực hiện các phép biến đổi sơ cấp h1→10h1−h3 h2→8h2+h3 h4→10h4+3h3−−−−−−−−−→  −56 0 0 0 56 0 0 −168 16 40 0 0 168 0 280 0 ∣∣∣∣∣∣∣∣∣ 392 −728 8 −616  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 17 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan Chọn phần tử trội không được nằm trên hàng 4,2,3 và cột 3,4,2 là phần tử a11 =−56. Thực hiện các phép biến đổi sơ cấp h2→h2+h1 h3→7h3+2h1 h4→h4+3h1−−−−−−−−→  −56 0 0 0 0 0 0 −168 0 280 0 0 0 0 280 0 ∣∣∣∣∣∣∣∣∣ 392 −336 840 560  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 18 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan Vậy hệ đã cho tương đương với hệ sau −56x1 = 392 −168x4 = −336 280x2 = 840 280x3 = 560 ⇔  x1 = −7 x2 = 3 x3 = 2 x4 = 2 Suy ra hệ đã cho có 1 nghiệm duy nhất (x1,x2,x3,x4)= (−7,3,2,2) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 19 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan Vậy hệ đã cho tương đương với hệ sau −56x1 = 392 −168x4 = −336 280x2 = 840 280x3 = 560 ⇔  x1 = −7 x2 = 3 x3 = 2 x4 = 2 Suy ra hệ đã cho có 1 nghiệm duy nhất (x1,x2,x3,x4)= (−7,3,2,2) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 19 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Phương pháp Gauss-Jordan Vậy hệ đã cho tương đương với hệ sau −56x1 = 392 −168x4 = −336 280x2 = 840 280x3 = 560 ⇔  x1 = −7 x2 = 3 x3 = 2 x4 = 2 Suy ra hệ đã cho có 1 nghiệm duy nhất (x1,x2,x3,x4)= (−7,3,2,2) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 19 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Bài tập BÀI TẬP BÀI TẬP 4.1 Sử dụng phương pháp phần tử trội giải hệ phương trình 2x1−1.5x2+3x3 = 1 −x1+2x3 = 3 4x1−4.5x2+5x3 = 1 Đáp số (x1,x2,x3)= (−1,0,1) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 20 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Gauss Bài tập BÀI TẬP BÀI TẬP 4.1 Sử dụng phương pháp phần tử trội giải hệ phương trình 2x1−1.5x2+3x3 = 1 −x1+2x3 = 3 4x1−4.5x2+5x3 = 1 Đáp số (x1,x2,x3)= (−1,0,1) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 20 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Những khái niệm cơ bản NHỮNG KHÁI NIỆM CƠ BẢN ĐỊNH NGHĨA 3.1 Ma trận vuông A =  a11 a12 . . . a1n 0 a22 . . . a2n ... ... . . . ... 0 0 . . . ann  được gọi là ma trận tam giác trên. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 21 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Những khái niệm cơ bản ĐỊNH NGHĨA 3.2 Ma trận vuông  a11 0 0 0 a21 a22 . . . 0 ... ... . . . ... an1 an2 . . . ann  được gọi là ma trận tam giác dưới. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 22 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp NỘI DUNG PHƯƠNG PHÁP 1 Nội dung của phương pháp nhân tử LU là phân tích ma trận A thành tích của 2 ma trận L vàU , trong đó L là ma trận tam giác dưới, cònU là ma trận tam giác trên. 2 Khi đó việc giải hệ (1) sẽ trở thành giải 2 hệ phương trình LY =B vàUX = Y . 3 Có nhiều phương pháp phân tích A = LU , tuy nhiên ta thường xét trường hợp L có đường chéo chính bằng 1 và gọi là phương pháp Doolittle. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 23 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp NỘI DUNG PHƯƠNG PHÁP 1 Nội dung của phương pháp nhân tử LU là phân tích ma trận A thành tích của 2 ma trận L vàU , trong đó L là ma trận tam giác dưới, cònU là ma trận tam giác trên. 2 Khi đó việc giải hệ (1) sẽ trở thành giải 2 hệ phương trình LY =B vàUX = Y . 3 Có nhiều phương pháp phân tích A = LU , tuy nhiên ta thường xét trường hợp L có đường chéo chính bằng 1 và gọi là phương pháp Doolittle. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 23 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp NỘI DUNG PHƯƠNG PHÁP 1 Nội dung của phương pháp nhân tử LU là phân tích ma trận A thành tích của 2 ma trận L vàU , trong đó L là ma trận tam giác dưới, cònU là ma trận tam giác trên. 2 Khi đó việc giải hệ (1) sẽ trở thành giải 2 hệ phương trình LY =B vàUX = Y . 3 Có nhiều phương pháp phân tích A = LU , tuy nhiên ta thường xét trường hợp L có đường chéo chính bằng 1 và gọi là phương pháp Doolittle. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 23 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp L vàU có dạng L =  1 0 0 0 `21 1 . . . 0 ... ... . . . ... `n1 `n2 . . . 1  U =  u11 u12 . . . u1n 0 u22 . . . u2n ... ... . . . ... 0 0 . . . unn  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 24 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp Các phần tử của 2 ma trận L vàU được xác định theo công thức u1 j = a1 j (1É j É n) `i1 = ai1 u11 (2É i É n) ui j = ai j − i−1∑ k=1 `ikuk j (1< i É j ) `i j = 1 ui j ( ai j − j−1∑ k=1 `ikuk j ) (1< j < i ) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 25 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp VÍ DỤ 3.1 Giải hệ phương trình bằng phương pháp LU 2x1+2x2−3x3 = 9 −4x1−3x2+4x3 = −15 2x1+x2+2x3 = 3 Giải. 2 2 −3−4 −3 4 2 1 2 =  1 0 0`21 1 0 `31 `32 1  .  u11 u12 u130 u22 u23 0 0 u33  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 26 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp VÍ DỤ 3.1 Giải hệ phương trình bằng phương pháp LU 2x1+2x2−3x3 = 9 −4x1−3x2+4x3 = −15 2x1+x2+2x3 = 3 Giải. 2 2 −3−4 −3 4 2 1 2 =  1 0 0`21 1 0 `31 `32 1  .  u11 u12 u130 u22 u23 0 0 u33  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 26 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp Theo công thức nhân 2 ma trận L vàU ta có 1.u11+0.0+0.0= a11 = 2⇒ u11 = 2; 1.u12+0.u22+0.0= a12 = 2⇒ u12 = 2; 1.u13+0.u23+0.u33 = a13 =−3⇒ u13 =−3. `21.u11+1.0+0.0= a21⇒ `21 = a21 u11 = −4 2 =−2; `21.u12+1.u22+0.0= a22⇒ u22 = a22−`21.u12 = −3− (−2).2= 1; `21.u13+1.u23+0.u33 = a23⇒ u23 = a23−`21.u13 = 4− (−2).(−3)=−2; Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 27 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp `31.u11+`32.0+1.0= a31⇒ `31 = a31 u11 = 2 2 = 1; `31.u12+`32.u22+1.0= a32⇒ `32 = 1 u22 (a32−`31.u12)= 1 1 (1−1.2)=−1; `31.u13+`32.u23+1.u33 = a33⇒ u33 = a33−`31.u13−`32.u23 = 2−1.(−3)− (−1).(−2)= 3 Do đó LY =B⇔  1 0 0−2 1 0 1 −1 1   y1y2 y3 =  9−15 3  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 28 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Nội dung phương pháp ⇒ Y = L−1B =  93 −3  UX = Y ⇔  2 2 −30 1 −2 0 0 3   x1x2 x3 =  93 −3  ⇒ X =U−1Y =  21 −1  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 29 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Bài tập BÀI TẬP BÀI TẬP 3.1 Sử dụng phương pháp nhân tử LU giải hệ phương trình 2x1−5x2+4x3 = 1 3x1+3x2+9x3 = 0 3x1+6x2+5x3 = 4 Đáp số (x1,x2,x3)= (89/34,2/17,−31/34) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 30 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Bài tập BÀI TẬP BÀI TẬP 3.1 Sử dụng phương pháp nhân tử LU giải hệ phương trình 2x1−5x2+4x3 = 1 3x1+3x2+9x3 = 0 3x1+6x2+5x3 = 4 Đáp số (x1,x2,x3)= (89/34,2/17,−31/34) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 30 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Bài tập trắc nghiệm BÀI TẬP TRẮC NGHIỆM BÀI TẬP 4.1 Cho A =  2 3 32 2 8 6 5 2  . Phân tích A = LU theo phương pháp Doolittle, phần tử `32 ma trận L là: 1 3.0000 2 4.0000 3 5.0000 4 6.0000 5 Các câu kia sai. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 31 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Bài tập trắc nghiệm A = LU =  1 0 0`21 1 0 `31 `32 1  .  2 3 30 u22 u23 0 0 u33  1.u11+0.0+0.0= a11 = 2⇒ u11 = 2; 1.u12+0.u22+0.0= a12 = 3⇒ u12 = 3; 1.u13+0.u23+0.u33 = a13 = 3⇒ u13 = 3. `21.u11+1.0+0.0= a21 = 2⇒ `21 = a21 u11 = 2 2 = 1; `21.u12+1.u22+0.0= a22 = 2 ⇒ u22 = a22−`21.u12 = 2−1.3=−1; Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 32 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Bài tập trắc nghiệm `31.u11+`31.0+1.0= a31 = 6⇒ `31 = a31 u11 = 6 2 = 3; `31.u12+`32.u22+1.0= a32 = 5 ⇒ `32 = a32−`31.u12 u22 = 5−3∗3−1 = 4; ĐS.⇒ Câu 2 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 33 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Bài tập trắc nghiệm BÀI TẬP 4.2 Cho A =  2 1 86 5 3 1 6 9  . Phân tích A = LU theo phương pháp Doolittle, tổng các phần tử u11+u22+u33 của ma trậnU là 1 63.7500 2 64.7500 3 65.7500 4 66.7500 5 Các câu kia sai. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 34 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Bài tập trắc nghiệm 1.u11+0.0+0.0= a11 = 2⇒ u11 = 2; 1.u12+0.u22+0.0= a12 = 1⇒ u12 = 1; 1.u13+0.u23+0.u33 = a13 = 8⇒ u13 = 8. `21.u11+1.0+0.0= a21 = 6⇒ `21 = a21 u11 = 6 2 = 3; `21.u12+1.u22+0.0= a22 = 5 ⇒ u22 = a22−`21.u12 = 5−3.1= 2; `21.u13+1.u23+0.u33 = a23 = 3 ⇒ u23 = a23−`21.u13 = 3−3.8=−21; Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 35 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp nhân tử LU Bài tập trắc nghiệm `31.u11+`31.0+1.0= a31 = 1⇒ `31 = a31 u11 = 1 2 ; `31.u12+`32.u22+1.0= a32 = 6 ⇒ `32 = a32−`31.u12 u22 = 6− 1 2 ∗1 2 = 11 4 ; `31.u13+`32.u23+1.u33 = a33 = 9 ⇒ u33 = a33−`31.u13−`32.u23 = 9− 1 2 ∗8− 11 4 ∗ (−21)= 251 4 ; Vậy u11+u22+u33 = 2+2+ 251 4 = 66.75. ⇒ Câu 4 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 36 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Ma trận xác định dương ĐỊNH NGHĨA 4.1 Ma trận vuông A được gọi là đối xứng nếu AT = A. ĐỊNH NGHĨA 4.2 Ma trận vuông A được gọi là xác định dương nếu như ∀x ∈Rn,x 6= 0 : xT Ax > 0. ĐỊNH LÝ 4.1 Ma trận vuông A xác định dương khi và chỉ khi tất cả những định thức con chính của nó đều lớn hơn 0. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 37 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Ma trận xác định dương ĐỊNH NGHĨA 4.1 Ma trận vuông A được gọi là đối xứng nếu AT = A. ĐỊNH NGHĨA 4.2 Ma trận vuông A được gọi là xác định dương nếu như ∀x ∈Rn,x 6= 0 : xT Ax > 0. ĐỊNH LÝ 4.1 Ma trận vuông A xác định dương khi và chỉ khi tất cả những định thức con chính của nó đều lớn hơn 0. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 37 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Ma trận xác định dương ĐỊNH NGHĨA 4.1 Ma trận vuông A được gọi là đối xứng nếu AT = A. ĐỊNH NGHĨA 4.2 Ma trận vuông A được gọi là xác định dương nếu như ∀x ∈Rn,x 6= 0 : xT Ax > 0. ĐỊNH LÝ 4.1 Ma trận vuông A xác định dương khi và chỉ khi tất cả những định thức con chính của nó đều lớn hơn 0. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 37 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski ĐỊNH LÝ 4.2 Cho ma trận vuông A là đối xứng và xác định dương. Khi đó A =B.BT , với B là ma trận tam giác dưới (bi i > 0, i = 1..n) và được xác định như sau: b11 =pa11,bi1 = ai1 b11 , (2É i É n) bi i = √ ai i − i−1∑ k=1 b2ik , (1< i É n) bi j = 1 b j j ( ai j − j−1∑ k=1 bikb j k ) , (1< j < i ) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 38 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski VÍ DỤ 4.1 Giải hệ phương trình bằng phương pháp Choleski  x1+x2−x3 = 1 x1+2x2 = 2 −x1+4x3 = 3 A =  1 1 −11 2 0 −1 0 4  là ma trận đối xứng Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 39 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski VÍ DỤ 4.1 Giải hệ phương trình bằng phương pháp Choleski  x1+x2−x3 = 1 x1+2x2 = 2 −x1+4x3 = 3 A =  1 1 −11 2 0 −1 0 4  là ma trận đối xứng Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 39 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski ∆1 = |1| = 1> 0, ∆2 = ∣∣∣∣∣ 1 11 2 ∣∣∣∣∣= 1> 0, ∆3 = ∣∣∣∣∣∣∣ 1 1 −1 1 2 0 −1 0 4 ∣∣∣∣∣∣∣= 2> 0. Vậy A là ma trận xác định dương. A =B.BT =  b11 0 0b21 b22 0 b31 b32 b33  .  b11 b21 b310 b22 b32 0 0 b33  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 40 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski ∆1 = |1| = 1> 0, ∆2 = ∣∣∣∣∣ 1 11 2 ∣∣∣∣∣= 1> 0, ∆3 = ∣∣∣∣∣∣∣ 1 1 −1 1 2 0 −1 0 4 ∣∣∣∣∣∣∣= 2> 0. Vậy A là ma trận xác định dương. A =B.BT =  b11 0 0b21 b22 0 b31 b32 b33  .  b11 b21 b310 b22 b32 0 0 b33  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 40 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski b11b11+0×0+0×0= a11 = 1⇒ b11 = 1. b11b21+0.b22+0×0= a12 = 1⇒ b21 = 1. b11b31+0.b32+0.b33 = a13 =−1⇒ b31 =−1. b21b11+0.b22+0×0= a21 = 1⇒ thỏa b21b21+b22.b22+0×0= a22 = 2⇒ b22 = 1. b21b31+b22.b32+0.b33 = a23 = 0⇒ b32 = 1. b31b11+0.b32+0.b33 = a31 =−1⇒ thỏa b31b21+b32.b22+0.b33 = a32 = 0⇒ thỏa b31b31+b32.b32+b33.b33 = a33 = 4⇒ b33 = p 2. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 41 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski Hệ phương trình viết lại dưới dạng ma trận Ax = b⇒B.BT .x = b⇔ { By = b BT x = y By = b⇒ y =B−1.b =  1 0 01 1 0 −1 1 p2  −1 .  12 3  BT x = y⇒ x =  1 1 −10 1 1 0 0 p 2  −1 .  1 1 3p 2 =  3−1/2 3/2  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 42 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski Hệ phương trình viết lại dưới dạng ma trận Ax = b⇒B.BT .x = b⇔ { By = b BT x = y By = b⇒ y =B−1.b =  1 0 01 1 0 −1 1 p2  −1 .  12 3  BT x = y⇒ x =  1 1 −10 1 1 0 0 p 2  −1 .  1 1 3p 2 =  3−1/2 3/2  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 42 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski VÍ DỤ 4.2 Tìmma trận B trong phép phân tích Choleski của ma trận A =  4 −3 0−3 4 −2 0 −2 4  ∆1 = 4, ∆2 = 7, ∆3 = 12. B =  2 0 0−32 p72 0 0 −4 p 7 7 2 p 21 7  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 43 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski VÍ DỤ 4.2 Tìmma trận B trong phép phân tích Choleski của ma trận A =  4 −3 0−3 4 −2 0 −2 4  ∆1 = 4, ∆2 = 7, ∆3 = 12. B =  2 0 0−32 p72 0 0 −4 p 7 7 2 p 21 7  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 43 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski BÀI TẬP 2.1 Cho A =  4 4 α4 6 2 α 2 7  . Với những giá trị nguyên nào của α thì ma trận A là xác định dương 1 α=−4 2 α=−2 3 α= 0 4 α= 6 5 Các câu kia sai Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 44 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Phương pháp Choleski ∆1 = |4| = 4> 0, ∆2 = ∣∣∣∣∣ 4 44 6 ∣∣∣∣∣= 8> 0, ∆3 = ∣∣∣∣∣∣∣ 4 4 α 4 6 2 α 2 7 ∣∣∣∣∣∣∣=−6.α2+16α+40> 0⇔ −1.5725<α< 4.2393. Vậy A là ma trận xác định dương khi α= 0⇒ Câu 3. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 45 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Bài tập trắc nghiệm BÀI TẬP 3.1 Cho A = ( 4 2 2 6 ) . Phân tích A =BBT theo phương pháp Choleski, ma trận B là 1 B = ( 2.00 0 1 2.24 ) . 2 B = ( 2.00 0 −1 2.24 ) . 3 B = ( 2.00 0 0 2.28 ) . 4 B = ( 2.00 0 −1 2.28 ) . 5 Các câu kia sai Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 46 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Bài tập trắc nghiệm A =B.BT = ( b11 0 b21 b22 ) . ( b11 b21 0 b22 ) b11b11+0×0= a11 = 4⇒ b11 = 2. b11b21+0.b22 = a12 = 2⇒ b21 = 1. b21b11+0.b22 = a21 = 2⇒ thỏa b21b21+b22.b22 = a22 = 6⇒ b22 = p 5.⇒ Câu 1 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 47 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Bài tập trắc nghiệm A =B.BT = ( b11 0 b21 b22 ) . ( b11 b21 0 b22 ) b11b11+0×0= a11 = 4⇒ b11 = 2. b11b21+0.b22 = a12 = 2⇒ b21 = 1. b21b11+0.b22 = a21 = 2⇒ thỏa b21b21+b22.b22 = a22 = 6⇒ b22 = p 5.⇒ Câu 1 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 47 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Bài tập trắc nghiệm BÀI TẬP 3.2 Cho A =  4 2 −52 3 3 −5 3 25  . Phân tích A =BBT theo phương pháp Choleski, tổng các phần tử b11+b22+b33 của ma trận B là 1 5.3182 2 5.3184 3 5.3186 4 5.3188 5 Các câu kia sai. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 48 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Bài tập trắc nghiệm A =B.BT =  b11 0 0b21 b22 0 b31 b32 b33  .  b11 b21 b310 b22 b32 0 0 b33  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 49 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Bài tập trắc nghiệm b11b11+0×0+0×0= a11 = 4⇒ b11 = 2. b11b21+0.b22+0×0= a12 = 2⇒ b21 = 1. b11b31+0.b32+0.b33 = a13 =−5⇒ b31 =−5 2 . b21b11+0.b22+0×0= a21 = 2⇒ thỏa b21b21+b22.b22+0×0= a22 = 3⇒ b22 = p 2. b21b31+b22.b32+0.b33 = a23 = 3⇒ b32 = 11 2 p 2 . b31b11+0.b32+0.b33 = a31 =−5⇒ thỏa b31b21+b32.b22+0.b33 = a32 = 3⇒ thỏa b31b31+b32.b32+b33.b33 = a33 = 25⇒ b33 = p 29 2 p 2 . b11+b22+b33 ≈ 5.3182. Câu 1 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 50 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Phương pháp Choleski Bài tập trắc nghiệm b11b11+0×0+0×0= a11 = 4⇒ b11 = 2. b11b21+0.b22+0×0= a12 = 2⇒ b21 = 1. b11b31+0.b32+0.b33 = a13 =−5⇒ b31 =−5 2 . b21b11+0.b22+0×0= a21 = 2⇒ thỏa b21b21+b22.b22+0×0= a22 = 3⇒ b22 = p 2. b21b31+b22.b32+0.b33 = a23 = 3⇒ b32 = 11 2 p 2 . b31b11+0.b32+0.b33 = a31 =−5⇒ thỏa b31b21+b32.b22+0.b33 = a32 = 3⇒ thỏa b31b31+b32.b32+b33.b33 = a33 = 25⇒ b33 = p 29 2 p 2 . b11+b22+b33 ≈ 5.3182. Câu 1 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 50 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Chuẩn của véctơ ĐỊNH NGHĨA 5.1 Trong không gian tuyến tính thực Rn. Chuẩn của véctơ X ∈Rn là một số thực, ký hiệu ||X || thỏa các điều kiện sau: 1 ∀X ∈Rn, ||X || Ê 0, ||X || = 0⇔ X = 0 2 ∀X ∈Rn,∀λ ∈R, ||λX || = |λ|.||X || 3 ∀X ,Y ∈Rn, ||X +Y || É ||X ||+ ||Y ||. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 51 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Chuẩn của véctơ Trong Rn có rất nhiều chuẩn, tuy nhiên ta chỉ xét chủ yếu 2 chuẩn thường dùng sau: ∀X = (x1,x2, . . . ,xn)T ∈Rn ||X ||1 = |x1|+ |x2|+ . . .+|xn| = n∑ k=1 |xk |. ||X ||∞ =max{|x1|, |x2|, . . . , |xn|}=max k=1,n |xk |. VÍ DỤ 5.1 Cho X = (1,2,3,−5)T . ||X ||1 = 1+2+3+5= 11 và ||X ||∞ =max{1,2,3,5}= 5 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 52 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Chuẩn của véctơ Trong Rn có rất nhiều chuẩn, tuy nhiên ta chỉ xét chủ yếu 2 chuẩn thường dùng sau: ∀X = (x1,x2, . . . ,xn)T ∈Rn ||X ||1 = |x1|+ |x2|+ . . .+|xn| = n∑ k=1 |xk |. ||X ||∞ =max{|x1|, |x2|, . . . , |xn|}=max k=1,n |xk |. VÍ DỤ 5.1 Cho X = (1,2,3,−5)T . ||X ||1 = 1+2+3+5= 11 và ||X ||∞ =max{1,2,3,5}= 5 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 52 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Chuẩn của véctơ Trong Rn có rất nhiều chuẩn, tuy nhiên ta chỉ xét chủ yếu 2 chuẩn thường dùng sau: ∀X = (x1,x2, . . . ,xn)T ∈Rn ||X ||1 = |x1|+ |x2|+ . . .+|xn| = n∑ k=1 |xk |. ||X ||∞ =max{|x1|, |x2|, . . . , |xn|}=max k=1,n |xk |. VÍ DỤ 5.1 Cho X = (1,2,3,−5)T . ||X ||1 = 1+2+3+5= 11 và ||X ||∞ =max{1,2,3,5}= 5 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 52 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Chuẩn củama trận CHUẨN CỦA MA TRẬN ĐỊNH NGHĨA 5.2 Chuẩn của ma trận tương ứng với chuẩn véctơ được xác định theo công thức ||A|| = max ||X ||=1 ||AX || = max ||X ||6=0 ||AX || ||X || Từ định nghĩa chuẩn của ma trận, ta có ||AX || É ||A||.||X || Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 53 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Chuẩn củama trận VÍ DỤ 5.2 Xác định chuẩn của ma trận A = ( 1 2 3 4 ) tương ứng với chuẩn ||X ||1. Với mọi X = ( x1 x2 ) thỏa ||X ||1 = |x1|+ |x2| = 1, ta có ||AX ||1 = |x1+2x2|+ |3x1+4x2| É 4|x1|+6|x2| = 4+2|x2| É 6. Do đó ||A|| = 6. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 54 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Chuẩn củama trận ĐỊNH LÝ 5.1 Chuẩn của ma trận A = (ai j ) được xác định như sau: ||A||1 = max 1É jÉn n∑ i=1 |ai j |− chuẩn cột ||A||∞ = max 1ÉiÉn n∑ j=1 |ai j |− chuẩn hàng Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 55 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Chuẩn củama trận ĐỊNH LÝ 5.1 Chuẩn của ma trận A = (ai j ) được xác định như sau: ||A||1 = max 1É jÉn n∑ i=1 |ai j |− chuẩn cột ||A||∞ = max 1ÉiÉn n∑ j=1 |ai j |− chuẩn hàng Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 55 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Chuẩn củama trận VÍ DỤ 5.3 Cho A =  2 −1 45 3 2 6 −7 3  . Lúc này ||A||1 =max{2+5+6,1+3+7,4+2+3}= 13, ||A||∞ =max{2+1+4,5+3+2,6+7+3}= 16. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 56 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Bài tập BÀI TẬP BÀI TẬP 3.1 Tính chuẩn ||.||1 và ||.||∞ của ma trận A =  3 1 −11 2 1 −1 1 4  . Đáp số ||A||1 = 6 và ||A||∞ = 6 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 57 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Chuẩn của véctơ, chuẩn củama trận Bài tập BÀI TẬP BÀI TẬP 3.1 Tính chuẩn ||.||1 và ||.||∞ của ma trận A =  3 1 −11 2 1 −1 1 4  . Đáp số ||A||1 = 6 và ||A||∞ = 6 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 57 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Những khái niệm cơ bản NHỮNG KHÁI NIỆM CƠ BẢN ĐỊNH NGHĨA 6.1 Xét dãy các véctơ (X (m))∞m=0 với X (m) ∈Rn.Dãy các véctơ này được gọi là hội tụ về véctơ X khim→+∞ nếu và chỉ nếu ||X (m)−X ||→ 0 khim→+∞ (hội tụ theo chuẩn). Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 58 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Những khái niệm cơ bản ĐỊNH LÝ 6.1 Để dãy các véctơ (X (m))∞m=0 hội tụ về véctơ X khim→+∞ thì điều kiện cần và đủ là những dãy (x(m)k ) hội tụ về xk ,∀k = 1,2, . . . ,n. (hội tụ theo tọa độ). Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 59 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Số điều kiện củama trận Xét hệ phương trình AX =B(det (A) 6= 0) có nghiệm x = A−1.B. Cho B một số gia ∆B , khi đó nghiệm X tương ứng sẽ có số gia ∆X và A.∆X =∆B⇔∆X = A−1.∆B.Như vậy, ta có ||∆X || = ||A−1.∆B || É ||A−1||.||∆B || và ||B || = ||AX || É ||A||.||X || Từ đây ta được ||∆X || ||X || É ||A||.||A −1||.||∆B ||||B || Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 60 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Số điều kiện củama trận ĐỊNH NGHĨA 6.2 Số nhỏ nhất k(A) thỏa điều kiện k(A)ÉCond(A)= ||A||.||A−1|| được gọi là số điều kiện của ma trận A. Số điều kiện k(A) của ma trận A thỏa 1É k(A)É+∞ Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 61 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Sự không ổn định của hệ phương trình tuyến tính Trong thực hành tính toán, ta có thể gặp những hệ phương trình tuyến tính mà những thay đổi nhỏ trên các hệ số tự do của hệ sẽ gây ra những thay đổi rất lớn về nghiệm. Hệ phương trình tuyến tính như vậy được gọi là hệ phương trình không ổn định trong tính toán. Nếu ngược lại, hệ được gọi là hệ phương trình ổn định trong tính toán Chú ý. Người ta chứng minh được rằng, số điều kiện của ma trận đặc trưng cho tính ổn định của hệ phương trình tuyến tính. Giá trị k(A) càng gần với 1 thì hệ càng ổn định. Số điều kiện k(A) càng lớn thì hệ càng mất ổn định. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 62 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Sự không ổn định của hệ phương trình tuyến tính Trong thực hành tính toán, ta có thể gặp những hệ phương trình tuyến tính mà những thay đổi nhỏ trên các hệ số tự do của hệ sẽ gây ra những thay đổi rất lớn về nghiệm. Hệ phương trình tuyến tính như vậy được gọi là hệ phương trình không ổn định trong tính toán. Nếu ngược lại, hệ được gọi là hệ phương trình ổn định trong tính toán Chú ý. Người ta chứng minh được rằng, số điều kiện của ma trận đặc trưng cho tính ổn định của hệ phương trình tuyến tính. Giá trị k(A) càng gần với 1 thì hệ càng ổn định. Số điều kiện k(A) càng lớn thì hệ càng mất ổn định. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 62 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Sự không ổn định của hệ phương trình tuyến tính Trong thực hành tính toán, ta có thể gặp những hệ phương trình tuyến tính mà những thay đổi nhỏ trên các hệ số tự do của hệ sẽ gây ra những thay đổi rất lớn về nghiệm. Hệ phương trình tuyến tính như vậy được gọi là hệ phương trình không ổn định trong tính toán. Nếu ngược lại, hệ được gọi là hệ phương trình ổn định trong tính toán Chú ý. Người ta chứng minh được rằng, số điều kiện của ma trận đặc trưng cho tính ổn định của hệ phương trình tuyến tính. Giá trị k(A) càng gần với 1 thì hệ càng ổn định. Số điều kiện k(A) càng lớn thì hệ càng mất ổn định. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 62 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Sự không ổn định của hệ phương trình tuyến tính VÍ DỤ 6.1 Xét hệ phương trình AX =B với A = ( 1 2 1 2.01 ) và B = ( 3 3.01 ) .Dễ dàng thấy được hệ có nghiệm X = ( 1 1 ) . Bây giờ xét hệ AX˜ = B˜ với B˜ = ( 3 3.1 ) . Nghiệm bây giờ của hệ là X˜ = ( −17 10 ) . Ta thấy k∞(A)= 1207.01>> 1.Do đó B ≈ B˜ nhưng X và X˜ khác nhau rất xa. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 63 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Bài tập trắc nghiệm BÀI TẬP TRẮC NGHIỆM BÀI TẬP 4.1 Cho A = ( 2 −4 6 9 ) . Số điều kiện tính theo chuẩnmột của ma trận A là: 1 3.6429 2 4.6429 3 5.6429 4 6.6429 5 Các câu kia sai. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 64 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Bài tập trắc nghiệm Mat A x−1⇒ A−1 = ( 3 14 2 21 −17 121 ) k = ||A||1.||A−1||1 =max{|2|+ |6|, |−4|+ |9|} ∗max{| 314|+ |− 17|, | 221|+ | 121|}= 13∗ 514 ≈ 4.64285 ⇒ Câu 2 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 65 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Bài tập trắc nghiệm BÀI TẬP 4.2 Cho A =  −6 −4 74 −3 −8 −4 5 −4  . Số điều kiện tính theo chuẩn vô cùng của ma trận A là: 1 4.6854 2 4.6954 3 4.7054 4 4.7154 5 Các câu kia sai. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 66 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Bài tập trắc nghiệm Mat A x−1⇒ A−1 =  − 13 112 − 19448 − 53448 − 328 − 13112 5112 − 156 − 23224 − 17224  k = ||A||∞.||A−1||∞ = max{|−6|+ |−4|+ |7|, |4|+ |−3|+ |−8|, |−4|+ |5|+ |−4|}∗max{|− 13112|+ |− 19448|+ |− 53448|, |− 328|+ |− 13112|+ | 5112|, |− 156|+ |− 23224|+ |− 17224|} = 17∗ 31112 ≈ 4.70535⇒ Câu 3 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 67 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Những phương pháp lặp là những phương pháp giải gần đúng hệ phương trình tuyến tính. Để giải hệ (1) ta được nó về dạng tương đương X = T X +C , với T là ma trận vuông cấp n và C và 1 véctơ cột đã biết. Xuất phát từ véctơ ban đầu X (0) ta xây dựng dãy (X (m))∞m=0 theo công thức X (m) = T X (m−1)+C , m = 1,2, . . . . (2) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 68 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp ĐỊNH LÝ 6.2 Nếu ||T || < 1 thì dãy các véctơ (X (m))∞m=0 xác định theo công thức lặp (2) sẽ hội tụ về véctơ nghiệm X của hệ với mọi véctơ lặp ban đầu X (0). Khi đó công thức đánh giá sai số là: ||X (m)−X || É ||T || m 1−||T ||.||X (1)−X (0)||(tiên nghiệm) ||X (m)−X || É ||T || 1−||T ||.||X (m)−X (m−1)||(hậu nghiệm) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 69 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi ĐỊNH NGHĨA 6.3 Ma trận A được gọi là ma trận đường chéo trội nghiêm ngặt nếu nó thỏa mãn điều kiện n∑ j=1, j 6=i |ai j | < |ai i |, i = 1,2, . . . ,n Chú ý. Nếu A là ma trận đường chéo trội nghiêm ngặt thì det A 6= 0 và ai i 6= 0,∀i = 1,2, . . . ,n. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 70 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi Xét hệ phương trình (1) với A là ma trận đường chéo trội nghiêm ngặt. Ta phân tích ma trận A theo dạng A =  a11 a12 . . . a1n a21 a22 . . . a2n . . . . . . . . . . . . an1 an2 . . . ann =  a11 0 . . . 0 0 a22 . . . 0 . . . . . . . . . . . . 0 0 . . . ann −  0 0 . . . 0 −a21 0 . . . 0 . . . . . . . . . . . . −an1 −an2 . . . 0  −  0 −a12 . . . −a1n 0 0 . . . −a2n . . . . . . . . . . . . 0 0 . . . 0  =D−L−U . Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 71 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi Do ai i 6= 0,∀i = 1,2, . . . ,n nên detD 6= 0.Như vậy D−1 =  1 a11 0 . . . 0 0 1a22 . . . 0 . . . . . . . . . . . . 0 0 . . . 1ann  Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 72 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi NỘI DUNG PHƯƠNG PHÁP Ta có AX =B⇔ (D−L−U )X =B⇔ (D)X = (L+U )X +B ⇔ X =D−1(L+U )X +D−1B. Ký hiệu T j =D−1(L+U ) và C j =D−1B. Khi đó công thức lặp có dạng X (m) = T jX (m−1)+C j , m = 1,2, . . . Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 73 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi Dạng tường minh của công thức lặp Jacobi là x(m)i = 1 ai i ( − i−1∑ j=1 ai jx m−1 j − n∑ j=i+1 ai jx m−1 j +bi ) . Ta có ||T j ||∞ = ||D−1(L+U )||∞ =max i=1,n n∑ j=1, j 6=i ∣∣∣∣ai jai i ∣∣∣∣< 1 do A là ma trận đường chéo trội nghiệm ngặt. Vậy ||T j || < 1 nên phương pháp Jacobi luôn hội tụ với mọi véctơ lặp ban đầu X (0). Thường thì ta chọn X (0) =C j Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 74 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi VÍ DỤ 6.2 Bằng phương pháp lặp Jacobi, tìm nghiệm gần đúng của hệ phương trình với sai số theo công thức hậu nghiệm nhỏ hơn 10−4, chọn chuẩn vô cùng 4x1+0.24x2−0.08x3 = 8 0.09x1+3x2−0.15x3 = 9 0.04x1−0.08x2+4x3 = 20 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 75 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi Giải. Ta thấy |0.24|+ |−0.08| < |4|; |0.09|+ |−0.15| < |3|; |0.04|+ |−0.08| < |4| nên ma trận hệ số A của hệ là ma trận đường chéo trội nghiệm ngặt. Do đó phương pháp lặp Jacobi luôn hội tụ. Đưa hệ về dạng X = T jX +C j Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 76 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi Cách 1: T j =D−1.(L+U )= 4 0 00 3 0 0 0 4  −1 .  0 −0.24 0.08−0.09 0 0.15 −0.04 0.08 0 = 0 −0.06 0.02−0.03 0 0.05 −0.01 0.02 0  Cách 2: x1 = 14(8−0.24x2+0.08x3)= 2−0.06x2+0.02x3 x2 = 13(9−0.09x1+0.15x3)= 3−0.03x1+0.05x3 x3 = 14(20−0.04x1+0.08x2)= 5−0.01x1+0.02x2 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 77 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi ⇔  x1x2 x3 =  23 5 +  0 −0.06 0.02−0.03 0 0.05 −0.01 0.02 0   x1x2 x3  Khi đó công thức lặp có dạng X (m) = T jX (m−1)+C j , m = 1,2, . . . Chọn X (0) =  23 5  tính X (1),X (2), . . . Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 78 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi Bấmmáy. X = (8−0.24B +0.08C )÷4 : Y = (9−0.09A+0.15C )÷3 : C = (20−0.04A+0.08B)÷4 : A = X :B = Y CALC B=3, C=5, A=2 Nhấn tiếp dấu ”=” cho tới nghiệm x(3)1 ,x (3) 2 ,x (3) 3 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 79 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi m x(m)1 x (m) 2 x (m) 3 0 2 3 5 1 1.92 3.19 5.04 2 1.9094 3.1944 5.0446 3 1.909228 3.194948 5.044794 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 80 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp lặp Jacobi ||T j ||∞ = max i=1,2,3 3∑ j=1 |ti j | = max{|0|+ |−0.06|+ |0.02|, |−0.03|+ |0|+ |0.05|, |−0.01|+|0.02|+|0|}=max{0.08;0.08;0.03}= 0.08 Đánh giá sai số ||X (3)−X (2)||∞ = max i=1,2,3 |x(3)i −x(2)i | =max{ |−1.72.10−4|, |5.48.10−4|, |1.94.10−4|}= 5.48×10−4 và ||X (3)−X ||∞ É ||T j ||∞ 1−||T j ||∞ .||X (3)−X (2)||∞ = 0.08 1−0.08×5.48×10 −4≈ 0.4765×10−4 < 10−4 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 81 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Bài tập trắc nghiệm BÀI TẬP 7.1 Cho hệ phương trình { 12x1−4x2 = 4 2x1+13x2 = 3 Với x(0) = [0.5,0.3]T , véctơ x(3) tính theo phương pháp Jacobi là 1 ( 0.384 0.176 ) 2 ( 0.386 0.174 ) 3 ( 0.388 0.172 ) 4 ( 0.390 0.170 ) 5 Các câu kia sai.Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 82 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Bài tập trắc nghiệm x1 = 1 12 (4+4x2) x2 = 1 13 (3−2x1) Bấmmáy. X = (4+4B)÷12 :B = (3−2A)÷13 : A = X CALC B=0.3, A=0.5 Nhấn tiếp dấu ”=” cho tới nghiệm x(3)1 ,x (3) 2 Kết quả: x(3) = ( 0.388 0.172 ) ⇒ Câu 3 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 83 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel NỘI DUNG PHƯƠNG PHÁP Đưa hệ (1) về dạng tương đương X = T X +C . Chọn véctơ xấp xỉ ban đầu X (0) (thường thì ta chọn X (0) =C .) Trong đó T =  t11 t12 . . . t1n t21 t22 . . . t2n . . . . . . . . . . . . tn1 tn2 . . . tnn  và C =  c1 c2 . . . cn  Từ hệ phương trình (2) ta được (D−L)X = UX +B⇒ X = (D−L)−1UX + (D−L)−1B. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 84 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel Đặt Tg = (D−L)−1U ,Cg = (D−L)−1B ta được công thức lặp Gauss-Seidel có dạng X (m) = TgX (m−1)+Cg , m = 1,2, . . . Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 85 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel DẠNG TƯỜNG MINH CỦA CÔNG THỨC LẶP GAUSS-SEIDEL x(m)1 = c1+ n∑ j=2 t1 jx m−1 j , x(m)2 = c2+ t21x(m)1 + n∑ j=3 t2 jx m−1 j , . . . . . . . . . . . . . . . x(m)i = ci + i−1∑ j=1 ti jx (m) j + n∑ j=i+1 ti jx m−1 j , . . . . . . . . . . . . . . . x(m)n = cn+ n−1∑ j=1 tn jx (m) j , (m = 1,2, . . .) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 86 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel Phương pháp Gauss-Seidel có thể xem là 1 biến dạng của phương pháp lặp Jacobi, nhưng khác phương pháp Jacobi ở chỗ: khi tính thành phần thứ i của véctơ lặp X (m) thì ta sử dụng ngay những thành phần x(m)1 ,x (m) 2 , . . . ,x (m) i−1 vừa tính được. Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 87 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel SỰ HỘI TỤ CỦA PHƯƠNG PHÁP GAUSS-SEIDEL Điều kiện hội tụ của phương pháp Gauss-Seidel hoàn toàn giống với phương pháp Jacobi. Công thức đánh giá sai số của nghiệm gần đúng ||X (m)−X || É ||Tg || 1−||Tg || .||X (m)−X (m−1)|| hoặc ||X (m)−X || É ||Tg || m 1−||Tg || .||X (1)−X (0)|| Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 88 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel VÍ DỤ 6.3 Bằng phương pháp lặp Gauss-Seidel, tìm nghiệm gần đúng của hệ phương trình với sai số theo công thức hậu nghiệm nhỏ hơn 10−4, chọn chuẩn vô cùng 4x1+0.24x2−0.08x3 = 8 0.09x1+3x2−0.15x3 = 9 0.04x1−0.08x2+4x3 = 20 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 89 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel Giải. Ta thấy |0.24|+ |−0.08| < |4|; |0.09|+ |−0.15| < |3|; |0.04|+ |−0.08| < |4| nên ma trận hệ số A của hệ là ma trận đường chéo trội nghiệm ngặt. Đưa hệ về dạng X (m) = TgX (m−1)+Cg , m = 1,2, . . .  x(m)1 = 14(8−0.24x(m−1)2 +0.08x(m−1)3 ) x(m)2 = 13(9−0.09x(m)1 +0.15x(m−1)3 ) x(m)3 = 14(20−0.04x(m)1 +0.08x(m)2 ) Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 90 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel Tg = (D−L)−1.U = 4 0 00.09 3 0 0.04 0.08 4  −1 .  0 −0.24 0.080 0 0.15 0 0 0 =  0 −0.06 0.020 1.8.10−3 0.0494 0 5.64.10−4 −1.188.10−3  Khi đó công thức lặp có dạng X (m) = TgX (m−1)+Cg , m = 1,2, . . . Chọn X (0) =  23 5  tính X (1),X (2), . . . Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 91 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel x(1)1 = c1+ t12x(0)2 + t13x(0)3 , x(1)2 = c2+ t21x(1)1 + t23x(0)3 , x(1)3 = c3+ t31x(1)1 + t32x(1)2 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 92 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel Bấmmáy. A = (8−0.24B +0.08C )÷4 : B = (9−0.09A+0.15C )÷3 : C = (20−0.04A+0.08B)÷4 CALC B=3, C=5. (không nhập A) Nhấn tiếp dấu ”=” cho tới nghiệm x(3)1 ,x (3) 2 ,x (3) 3 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 93 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel m x(m)1 x (m) 2 x (m) 3 0 2 3 5 1 1.92 3.1924 5.044648 2 1.9093489 3.194952 5.0448056 3 1.909199 3.1949643 5.0448073 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 94 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Phương pháp Gauss-Seidel Đánh giá sai số ||X (3)−X (2)||∞ = max i=1,2,3 |x(3)i −x(2)i | = max{|−1.499.10−4|, |0.123.10−4|, |0.017.10−4|}= 1.499×10−4 ||Tg ||∞ =max{|0|+|−0.06|+|0.02|, |0|+|1.8.10−3|+|0.0494|, |0|+ |5.64.10−4|+ |−1.188.10−3|}= max{0.08,0.0512,1.744.10−3}= 0.08. ||X (3)−X ||∞ É ||Tg || 1−||Tg || · ||X (3)−X (2)||∞ = 0.08 1−0.08×1.499×10 −4≈ 0.1303×10−4 < 10−4 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 95 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Bài tập trắc nghiệm BÀI TẬP 9.1 Cho hệ phương trình { 15x1−6x2 = 5 −5x1+8x2 = 5 Với x(0) = [0.3,0.2]T , véctơ x(3) tính theo phương pháp Gauss-Seidel là 1 ( 0.753 1.099 ) 2 ( 0.755 1.097 ) 3 ( 0.757 1.095 ) 4 ( 0.759 1.093 ) . 5 Các câu kia sai.Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 96 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Bài tập trắc nghiệm A = (5+6B)÷15 :B = (5+5A)÷8 CALC B=0.2 (không nhập A) Nhấn tiếp dấu ”=” cho tới nghiệm x(3)1 ,x (3) 2 Vậy x(3) = ( 0.755 1.096875 ) .⇒ Câu 2 Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 97 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt Những phương pháp lặp Bài tập trắc nghiệm CÁMƠNCÁC EMĐÃCHÚ Ý LẮNGNGHE Nguyễn Thị Cẩm Vân (BK TPHCM) HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Ngày 12 tháng 2 năm 2018 98 / 98 CuuDuongThanCong.com https://fb.com/tailieudientucntt