ĐẠI HỌC QUỐC GIA TP HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA BÁO CÁO BÀI TẬP LỚN MÔN ĐẠI SỐ TUYẾN TÍNH ĐỀ TÀI 2 CƠ SỞ LÝ THUYẾT VÀ ỨNG DỤNG PHÂN TÍCH A = LU GVHD Nguyễn Xuân Mỹ Lớp DT01 Nhóm số 2 TP HỒ CHÍ MINH,[.]
ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA BÁO CÁO BÀI TẬP LỚN MÔN ĐẠI SỐ TUYẾN TÍNH ĐỀ TÀI CƠ SỞ LÝ THUYẾT VÀ ỨNG DỤNG PHÂN TÍCH A = LU GVHD: Nguyễn Xuân Mỹ Lớp: DT01 Nhóm số: TP HỒ CHÍ MINH, tháng 05 năm 2023 Tra ng I ĐẠI HỌC QUỐC GIA TP.HCM TRƯỜNG ĐẠI HỌC BÁCH KHOA BÁO CÁO BÀI TẬP LỚN MƠN ĐẠI SỐ TUYẾN TÍNH ĐỀ TÀI CƠ SỞ LÝ THUYẾT VÀ ỨNG DỤNG PHÂN TÍCH A = LU GVHD Nguyễn Xuân Mỹ Lớp: DT01 Nhóm số: Danh sách thành viên MSSV Nguyễn Tiến Đạt 2210704 Lê Minh Giang 2210821 Nguyễn Hương Giang 2210824 Từ Ngọc Hân 2113317 Phạm Ngọc Bảo Hằng 2210928 Nguyễn Minh Hảo 2013092 Tô Phước Hào 1913229 Mai Huy Hoàng 1913432 Phạm Trương Tung Hoành 2211127 10 2211344 Nguyễn Xuân Hùng Trang MỤC LỤC Trang MỤC LỤC _ VỀ ĐỀ TÀI CHƯƠNG 1: CƠ SỞ LÝ THUYẾT I.1 CÁC KHÁI NIỆM VỀ MA TRẬN I.1.1 Ma trận _ I.1.2 Các loại ma trận I.1.3 Phép nhân hai ma trận I.1.4 Hệ phương trình tuyến tính I.2 PHÂN TÍCH A = LU I.2.1 Giới thiệu _ I.2.2 Phân tích A = LU _ I.3 PHÂN TÍCH PA = LU _ CHƯƠNG 2: CHƯƠNG TRÌNH PHÂN TÍCH A = LU II.1 THUẬT TOÁN CƠ BẢN II.2 CHẠY CHƯƠNG TRÌNH II.3 NHẬN XÉT CHƯƠNG TRÌNH _ CHƯƠNG 3: ỨNG DỤNG CỦA PHÂN TÍCH A = LU _ 10 III.1 GIẢI PHƯƠNG TRÌNH TUYẾN TÍNH _ 10 III.2 TÌM MA TRẬN NGHỊCH ĐẢO _ 11 III.3 TÍNH ĐỊNH THỨC CỦA MA TRẬN _ 12 TÀI LIỆU THAM KHẢO _ 13 Trang VỀ ĐỀ TÀI ĐỀ 1/ Nêu sở lý thuyết phân tích A = LU phân tích PA = PLU 2/ Viết chương trình dùng để phân tích A = LU 3/ Tìm ứng dụng phân tích A = LU Trang CHƯƠNG 1: CƠ SỞ LÝ THUYẾT I.1 CÁC KHÁI NIỆM VỀ MA TRẬN I.1.1 Ma trận Trong toán học, ma trận mảng chữ nhật – số, ký hiệu, biểu thức, xếp theo hàng cột – mà ma trận tuân theo quy tắc định trước Từng ô ma trận gọi phần tử hoặ c mụ c Ví dụ ma trận có hàng cột: � −3 � I.1.2 Các loại ma trận Ma trận tam giác ma trận đường chéo Nếu phần tử A bên đường chéo 0, A gọi ma trận tam giác Tương tự, phần tử A bên đường chéo 0, A gọi ma trận tam giác Nếu phần tử nằm bên ngồi đường chéo 0, A gọi ma trận đường chéo Ma trận đơn vị Ma trận đơn vị In có số chiều n ma trận n × n phần tử đường chéo tất những1 ph0ần 0tử khác 0, ví dụ �010� 0 Ma trận bậc thang Ma trận A ma trận bậc thang thoả mãn hai điều kiện: + Các hàng (nếu có) ma trận + Phần tử sở hàng phía bên phải so với phần tử sở hàng −3 Biến đổi ma trận bậc thang thơng qua phương pháp khử Gauss Ví dụ 0 −1 0 0 � + � Ma trận nghịch đảo Ma trận vuông A gọi khả nghịch hay không suy biến tồn ma trận B cho AB = BA = In + Nếu B tồn tại, gọi ma trận nghịch đảo A, ký −1 hiệu A + Một số tính chất ma trận khả nghịch: −1 −1 (A ) =A (AB) −1 −1 −1 =B A T −1 (A ) −1 T = (A ) Trang I.1.3 Phép nhân hai ma trận Cho hai ma trận Am×n Bn×p ma trận A có số cột số dịng ma trận B Tích ma trận A ma trận B ma trận cấp m×p, kí hiệu AB Ví dụ phép tính ma trận 3×2 với ma trận 2×3: a× g + b× j ��×� a× h + b× k � = �c × g + d × j a× i + b × l c × h + d × k c × i + d × l� I.1.4 Hệ phương g+f×j h+f×k i+f×l trình tuyến tính Một hệ phương trình tuyến tính tập hợp phương trình tuyến tính có x1 − 2x2 − x3 + 5x4 = biến số, ví dụ: −x1 + 3x2 + 4x3 − 3x4 � −x1 + 4x2 + 7x3 2x1 − 5x2 − 5x3 − x4 = −1 = −1 + 8x4 = Giải hệ phương trình cách đưa dạng AX = b sử dụng phương pháp khử Gauss Nếu biến số hệ phương trình tuyến tính số thực số phức, có ba trường hợp xảy ra: hệ vơ nghiệm, hệ có nghiệm, hệ có vơ số nghiệm I.2 PHÂN TÍCH A = LU I.2.1 Giới thiệu Một ứng dụng phép biến đổi sơ cấp dùng để phân tích ma trận A = LU Phân tích có nhiều ứng dụng: dùng để tính định thức, giải phương trình tuyến tính, tìm ma trận nghịch đảo Trong đó: + A ma trận vng cấp n + L ma trận tam giác (L lower lower triangle) + U ma trận tam giác (U upper upper triangle) Sử dụng phép biến đổi sơ cấp đưa A ma trận phía U Để tính ma trận L cần dùng phép biến đổi sơ cấp ngược lại với biến đổi biến đổi I thành L I.2.2 Phân tích A = LU Tìm ma trận U Phép phân tích ma trận đơn giản, ta thực phép biến đổi dòng để đưa A thành ma trận bậc thang Lúc đó, ma trận bậc thang ma A= trận tam giác U Lấy ma trận có kích thước 3×3 để làm ví dụ: �3 −2 4� Trang Đầu tiên, ta cần tính tốn hệ số để nhân dịng ( =3 A21 thực loại trừ Có thể tính 31 l =A31 = −2 phép biến r2 A11 r2 �3 r33 = − r r r phép r A11 l 31 l để l − = + Đến đây, ta có ký hiệu chung Tiếp tục l21 r1 r) Vậy ọ i hệ số đổi dòng thứ hai r2 21 = lij biến đổi dòng ) trừ cho dòng ( r = G hệ số nhân cho dòng j trừ cho dòng i − 3r −2 4�� �� −9 −2 bước biến đổi , ta có l32 = −2� 12 −9 =− −9 +4� r2 −9 12 −9 1 11 � −9 � 0 −2 � � 25 −2� 25 Vậy kết thu ma trận tam giác U: −2 11 U= � � �0 � � −2� 12 � �0 + 2r r3 r3 � Tiếp theo cần phải xác định ma trận tam giác L Tìm ma trận L 5 −1 −1 Ta có phân tích A = LU nên tìm L thơng qua công thức L = AU L=�3 4� −2 L= �0 −9 −2 � 0 0 Trên thực tế không cần thực phép tính trên, ma tr 0� �l21 31 l � 32 l = ận L có t hể tìm sau: 25 �3 0� −4 Những phần tử dòng L giúp đảo � −2 21 Ví dụ với dịng ma trận U có từ việc trừ lại cá c phép t oán dịng mà ta làm 1 dịng ma trận L ta khơi phục lại, lấy dòng U cộng thêm vào dòng để lấy l ại dòng ban đầu: l r (l 21 0) � −9 −2 � 25 = + 0 � l21(1 2) (0 −9 −2)=(3 4) Trang l21 L biểu diễn sau: L = � ⋮ lm1 lm2 0 1 � � ⋯ 1 −9 0 �� −4 A=LU= ⋱ Vậy kết phân tích ⋯ ⋯ −2 � 25 � I.3 PHÂN TÍCH PA = LU Đôi cần phải thay đổi dịng với thực loại trừ Trong trường hợp ta khơng thể sử dụng cách phân tích A = LU nói LU kết phép phân tích ma trận A đổi dịng (tạm gọi A′) Lúc ta phải sử dụng⇔ thêm ma trận P để thực phép hốn đổi dịng A Lúc này: A′ = LU PA = LU Trang CHƯƠNG 2: CHƯƠNG TRÌNH PHÂN TÍCH A = LU − II.1 THUẬT TỐN CƠ BẢN Cho ma trận m×m A = (an,n), đặt A (0) :=A lặp với n = 1, 2, , m Khử phần tử bên đường chéo cột thứ n A (n–1) cách cộng vào dòng thứ i ma trận với dòng thứ n nhân thêm hệ số l i,n, với i = n + 1, n + 2, …, m Sau m bước, ta khử tất phần tử bên đường chéo chính, nhận tam giác A(m–1) A(m–1) ma trận − ma trận tam giác U Có thể sử dụng ma −1 trận nghịch đảo để tìm L đơn giản công thức L = AU A=�−3 � II.2 CHẠY CHƯƠNG TRÌNH Ví dụ 1: Cho A ma trận cấp 2×3, Phân tích A = LU Hình Chương trình báo lỗi Do ma trận A ma trận cấp 2×3 nên chương trình khơng thể chạy Để chương trình phân tích chạy cần nhập vào ma trận vng cấp n Ví dụ 2: Cho A ma trận vuông cấp 3, A=�2 −3� 1 −1 Kết phân tích Hình Cửa sổ Command Window phân tích ma trận A = LU Trang 2354 A=�3 6� 4593 Ví dụ 3: Cho A ma trận vng cấp 4, 1121 Phân tích A = LU Hình Cửa sổ Command Window phân tích ma trận A = LU II.3 NHẬN XÉT CHƯƠNG TRÌNH Chương trình phân tích A = LU Matlab giúp phân tích nhanh ma trận A thành ma trận tam giác U ma trận tam giác L Từ ứng dụng vào tốn khác Thuận tiện cho việc làm tập tính tốn Tuy nhiên để viết chương trình cần có kiến thức lập trình code Matlab Vì cần có thời gian tìm hiểu tìm tịi thuật tốn Nên sử dụng phân tích cho tập có khối lượng tính tốn nhiều Ngồi chương trình cịn chưa giải số trường hợp đặc biệt cần hoàn thiện Trang CHƯƠNG 3: ỨNG DỤNG CỦA PHÂN TÍCH A = LU III.1 GIẢI HỆ PHƯƠNG TRÌNH TUYẾN TÍNH Từ phương trình tuyến tính, ta có Ax = b (với A b cho trước) Từ A ta phân tích thành A = LU Khi nghiệm phương trình tính tốn theo bước sau: Giải Ly = b => Tìm y Giải Ux = y => Tìm x −x1 + 1x2 − 1x3 = −5 Ví dụ: Giải hệ phương trình 17x1 Từ hệ phương trình −1 A=�17 10 + 10x 19x1 + 6x −1 ta bi ến đổi t hành: −5 = −19 b = �−1919� x = �x2 � 0� 19 = 19 0 −1 � −1 Phân tích ma trận A = LU: x3 −17 19� A=LU=�−19 0 �� −−2817�27 0 27 Hình Kết phân tích ma trận A Matlab y Giải Ly = b: Ly = b ⇔ � 0 −17 � =�19� �y3� y 27 y1 −19 −5 y2 ⇔ � � 19 −19 −5 −66 � � = � −608 � −1 Giải Ux = y: Ux = y ⇔ � −1 x1 −5 y 27 −17 ⇔ x1 � x2 x3 x −66 �=� �� −28 Vậy nghiệm hệ cho x �27 �= � �9 � −151 −608 456 x1 = −151 , x2 = 2707 , x3 = 4567 Trang 10 III.2 TÌM MA TRẬN NGHỊCH ĐẢO Khi giải hệ phương trình, thường b xem vector có chiều dài số dịng A Nếu thay vector b, ta có ma trận B, với B ma trận kích thước n×p, ta phải tìm ma trận X (cũng có kích thước n×p): AX = LUX = B Có thể sử dụng phương pháp Ly = b, Ux = y để giải cho cột ma trận X Với giả sử B ma trận đơn vị với kích thước n X nghịch đảo A A=�2 Ví dụ cho ma trận Gọi ma trận A, B, X: 6� 1 Tìm ma trận nghịch đảo 11 B=�0 0� A= �2 � 2 x21 X=� Phân 1 2 0 2 0� �0 2� 0 x x12 x13� x22 x23 32 33 31 tí ch A = LU: A=LU⇔�2 6�=�2 1 i Giải Ly = B ột ma trận B , với i c 1 Ly = B = �0� − ⇔�2 0�� 1 −1 ⇔ y11 � �y21� = y31 Giải Ux = y UX = y ⇔ � − 2 x11 −−32 �0 2� �x21� = 0 ⇔ x113 �x21� = x31 �0 x31 x22 x23 Thực tương tự với B2 B3, tìm −4 12 x32 � �= �3 �1 � 1 13 x � 33 � X − �−22� 2� ⇒ Vậy −4 ma trận nghịch đảo A Trang 11 III.3 TÍNH ĐỊNH THỨC CỦA MA TRẬN Các ma trận L U dùng để tính định thức ma trận A hiệu det(A) = det(L) det(U) định thức ma trận tam giác đơn giản tích phần tử đường chéo Đặc biệt, L ma trận tam giác đơn vị thì: �3 Ví dụ: Tính định thức ma trận A = Ta có A = LU = �3 1 0 1 0 1 −1 �� 0 0 0 � 1 � −3 Hình Kết phân tích ma trận A Matlab − ⇒ det(A) = det(L) × det(U) = (1×1×1×1) × (1×1×2×( 3)) = Vậy det(A) = − − Trang 12 TÀI LIỆU THAM KHẢO [1] Đặng Văn Vinh, Giáo trình Đại Số Tuyến Tính, Nhà Xuất Bản Đại Học Quốc Gia, 2020 [2] https://vi.wikipedia.org/wiki/Phân_tích_LU [3] https://tailieuvnu.com/wp-content/uploads/2020/11/19/Slide-2-Dai-so-Tuyen-TinhMa-Tran-nghich-dao-va-phan-tich-LU-Le-Xuan-Thanh-UET.pdf [4] https://rootonchair.blogspot.com/2019/06/phep-phan-tich-ma-tran-alu.html [5] Vũ Thị Hương Trang, Phân Tích Ma Trận Và Một Số Ứng Dụng, Đại Học Thái Nguyên, 2016 Trang 13