Quy tắc Cramer: Phương pháp trực tiếp để giải hệ tuyến tính

Quy tắc Cramer cung cấp một công thức trực tiếp để giải hệ phương trình tuyến tính bằng định thức. Mặc dù các phương pháp tính toán hiện đại thường hiệu quả hơn quy tắc Cramer khi giải các hệ lớn, việc hiểu kỹ thuật này giúp chúng ta nắm bắt mối quan hệ giữa hệ tuyến tính, định thức và lý thuyết ma trận trong đại số tuyến tính. 

Bài viết này giới thiệu quy tắc Cramer, công thức, nền tảng toán học và diễn giải hình học của nó. Chúng tôi cũng minh họa cách áp dụng qua nhiều ví dụ, các giới hạn, lựa chọn thay thế, và khi nào nó hiệu quả nhất trong bối cảnh khoa học dữ liệu.

Quy tắc Cramer là gì?

Quy tắc Cramer phát biểu rằng mỗi biến trong hệ có thể được tìm bằng cách chia hai định thức: định thức của một ma trận hệ số đã được thay đổi cho định thức của ma trận hệ số gốc.

Ma trận đã thay đổi này được tạo bằng cách thay thế một cột hệ số bằng các hằng số từ vế phải của các phương trình.

Để hiểu các điều kiện của định lý, hãy xét một hệ phương trình tuyến tính được biểu diễn dưới dạng ma trận là Ax = b, trong đó:

• A là ma trận hệ số
• x là vectơ ẩn số
• b là vectơ hằng số

Quy tắc Cramer chỉ áp dụng khi A là ma trận vuông (số phương trình bằng số ẩn) và có định thức khác 0, đảm bảo tồn tại nghiệm duy nhất.

Nền tảng toán học và công thức

Với hệ tổng quát n×n, quy tắc Cramer phát biểu rằng thành phần thứ i của vectơ nghiệm là:

Công thức quy tắc Cramer. (Hình ảnh và tất cả hình ảnh bên dưới bởi Tác giả)

Trong đó:

• x_i là biến thứ i mà chúng ta cần tìm
• det(A_i) là định thức của ma trận A khi thay cột thứ i bởi vectơ b
• det(A) là định thức của ma trận hệ số gốc

Để hiểu rõ hơn, hãy xét một hệ 2×2:

Ma trận hệ số A và vectơ hằng số b là:

Theo quy tắc Cramer:

Trong đó A₁ là A với cột thứ nhất được thay bằng b:

Và A₂ là A với cột thứ hai được thay bằng b:

Định thức của ma trận 2×2 được tính như sau:

Khi đã hiểu khuôn mẫu nghiệm ở trên, chúng ta có thể áp dụng trực tiếp cho bất kỳ hệ tương thích nào mà không cần giải hệ đồng thời. Chúng ta cũng sẽ mở rộng quy tắc cho hệ 3x3 trong một ví dụ sau.

Diễn giải hình học của quy tắc Cramer

Ngoài việc là công thức trực tiếp để giải hệ tuyến tính, quy tắc Cramer còn có một diễn giải hình học phong phú.

Diện tích hình bình hành trong 2D

Hãy coi các cột của ma trận A là các vectơ trong mặt phẳng. Các vectơ này tạo thành một hình bình hành, và định thức của A biểu thị diện tích có dấu của hình bình hành đó.

Khi chúng ta thay một cột bằng vectơ b, ta tạo ra một hình bình hành mới. Định thức của ma trận mới này, det⁡(A_i), biểu thị diện tích của hình đã thay đổi, trong đó một trong các vectơ phương ban đầu được thay bằng b.

Do đó:

Mỗi biến x_i cho biết b nằm bao nhiêu theo hướng của vectơ cột thứ i của A, được chuẩn hóa theo diện tích tổng thể do các cột của A tạo ra.

Thể tích hình hộp trong 3D

Lập luận tương tự mở rộng sang 3D. Ở đây, định thức của A biểu thị thể tích của hình hộp song song được tạo bởi ba vectơ cột của A.

• Thay một cột bằng b cho ta một hình hộp mới.
• Tỷ số det⁡(A_i)/det⁡(A) cho biết mức đóng góp thể tích tương đối của cột đó vào hình tổng thể.

Trong cả 2D và 3D, quy tắc Cramer đo lường mức độ b “lấp đầy” thể tích do các cột của A xác định.

Hướng dẫn từng bước áp dụng quy tắc Cramer

Hãy thực hiện hai ví dụ để minh họa cách quy tắc Cramer áp dụng cho các hệ có kích thước khác nhau.

Giải hệ 2x2

Xét hệ gồm hai phương trình tuyến tính sau:

Bước 1: Xác định ma trận hệ số và vectơ hằng

Ta có thể viết hệ dưới dạng ma trận hệ số A và vectơ hằng b như sau:

Bước 2: Tính định thức của ma trận hệ số

Định thức của A sẽ là:

det(A) = (3)(-1) - (2)(2)
det(A) = -3 - 4
det(A) = -7Được hỗ trợ bởi

Vì det(A) = -7 ≠ 0, hệ có nghiệm duy nhất và chúng ta có thể áp dụng quy tắc Cramer.

Bước 3: Tạo các ma trận đã thay đổi và tính định thức của chúng

Đối với x, thay cột thứ nhất của A bằng b:

Bây giờ hãy tính định thức của A₁:

det(A₁) = (13)(-1) - (2)(4)
det(A₁) = -13 - 8
det(A₁) = -21Được hỗ trợ bởi

Đối với y, thay cột thứ hai của A bằng b:

Bây giờ hãy tính định thức của A₂:

det(A₂) = (3)(4) - (13)(2)
det(A₂) = 12 - 26
det(A₂) = -14Được hỗ trợ bởi

Bước 4: Áp dụng công thức quy tắc Cramer

Chúng ta có thể tính các giá trị theo công thức:

x = det(A₁) / det(A) = -21 / -7 = 3
y = det(A₂) / det(A) = -14 / -7 = 2Được hỗ trợ bởi

Cuối cùng, ta có thể kiểm chứng kết quả bằng cách thay các giá trị vào phương trình:

• Phương trình thứ nhất: 3(3) + 2(2) = 9 + 4 = 13 
• Phương trình thứ hai: 2(3) — 1(2) = 6–2 = 4 

Như vậy, chúng ta đã giải hệ tuyến tính 2x2 bằng quy tắc Cramer.

Giải hệ 3x3

Bây giờ hãy mở rộng sang hệ 3x3 phức tạp hơn:

Chúng ta sẽ lặp lại các bước như trên cho ví dụ này.

Bước 1: Xác định ma trận hệ số và vectơ hằng

Từ hệ tuyến tính, ta có thể viết ma trận hệ số A và vectơ hằng b như sau:

Bước 2: Tính định thức của ma trận hệ số

Dùng khai triển theo hàng đầu tiên:

det(A) =       2|2  1| - 1 |1  1| + 1 |1  2|
        		 |1  2|     |1  2|       |1  1|Được hỗ trợ bởi

Tính từng định thức 2×2:

|2  1| = (2)(2) - (1)(1) = 4 - 1 = 3
|1  2|

|1  1| = (1)(2) - (1)(1) = 2 - 1 = 1
|1  2|

|1  2| = (1)(1) - (2)(1) = 1 - 2 = -1
|1  1|Được hỗ trợ bởi

Vì vậy, định thức của A là:

det(A) = 2(3) - 1(1) + 1(-1)
det(A) = 6 - 1 - 1
det(A) = 4Được hỗ trợ bởi

Vì det(A) = 4 ≠ 0, hệ có nghiệm duy nhất và chúng ta có thể áp dụng quy tắc Cramer.

Bước 3: Tạo các ma trận đã thay đổi và tính định thức của chúng

Đối với x, thay cột thứ nhất của A bằng b:

Ta có thể khai triển theo hàng đầu tiên:

det(A₁) = 7|2  1| - 1|8  1| + 1|8  2|
   	        |1  2|     |9  2|      |9  1|Được hỗ trợ bởi

Tính định thức của A₁:

det(A₁) = 7(3) - 1(16 - 9) + 1(8 - 18)
det(A₁) = 21 - 1(7) + 1(-10)
det(A₁) = 21 - 7 - 10
det(A₁) = 4Được hỗ trợ bởi

Đối với y, thay cột thứ hai của A bằng b:

Ta có thể tính định thức của A₂ như sau:

det(A₂) = 2|8  1| - 7|1  1| + 1|1  8|
                     |9  2|     |1  2|    |1  9|

det(A₂) = 2(16 - 9) - 7(2 - 1) + 1(9 - 8)
det(A₂) = 2(7) - 7(1) + 1(1)
det(A₂) = 14 - 7 + 1
det(A₂) = 8Được hỗ trợ bởi

Đối với z, thay cột thứ ba của A bằng b:

Ta có thể tính định thức của A₃ như sau:

det(A₃) = 2|2  8| - 1|1  8| + 7|1  2|
                     |1  9|      |1  9|       |1  1|

det(A₃) = 2(18 - 8) - 1(9 - 8) + 7(1 - 2)
det(A₃) = 2(10) - 1(1) + 7(-1)
det(A₃) = 20 - 1 - 7
det(A₃) = 12Được hỗ trợ bởi

Bước 4: Áp dụng công thức quy tắc Cramer

Vì chúng ta đã tính đủ các giá trị cần thiết, có thể tính trực tiếp các giá trị của x, y và z theo công thức:

x = det(A₁) / det(A) = 4 / 4 = 1
y = det(A₂) / det(A) = 8 / 4 = 2
z = det(A₃) / det(A) = 12 / 4 = 3Được hỗ trợ bởi

Cuối cùng, ta có thể kiểm chứng kết quả bằng cách thay các giá trị vào phương trình:

• Phương trình thứ nhất: 2(1) + 1(2) + 1(3) = 2 + 2 + 3 = 7 
• Phương trình thứ hai: 1(1) + 2(2) + 1(3) = 1 + 4 + 3 = 8 
• Phương trình thứ ba: 1(1) + 1(2) + 2(3) = 1 + 2 + 6 = 9 

Lưu ý quy trình vẫn nhất quán: chúng ta thay các cột bằng vectơ hằng và tính định thức. Trường hợp 3×3 đòi hỏi nhiều phép tính hơn, nhưng phương pháp cơ bản giống hệt trường hợp 2×2.

Giới hạn, lựa chọn thay thế và khi nào dùng quy tắc Cramer

Giống như mọi thuật toán, quy tắc Cramer cũng có một số giới hạn. Hãy cùng tìm hiểu:

Độ phức tạp tính toán

Với hệ n×n, quy tắc Cramer cần tính (n+1) định thức. Vì tính định thức bằng khai triển phần bù đại số có độ phức tạp O(n!), độ phức tạp tổng thể trở nên quá lớn đối với các hệ kích thước lớn. 

Ngay cả với các thuật toán định thức hiệu quả hơn, phương pháp này vẫn không thực tế vượt quá các hệ nhỏ.

Vấn đề ổn định số

Quy tắc Cramer có thể gặp mất ổn định số khi định thức gần bằng không. 

Các sai số làm tròn nhỏ khi tính det(A) có thể dẫn đến sai số lớn ở nghiệm khi chia cho một giá trị gần 0. Điều này khiến phương pháp không đáng tin cậy với các hệ điều kiện kém, vốn phổ biến trong dữ liệu thực tế.

Các lựa chọn thay thế quy tắc Cramer

Do hạn chế chính về độ phức tạp tính toán, chúng ta thường dùng một trong các phương án sau:

• Khử Gauss với chọn phần tử trụ một phần biến đổi hệ về dạng bậc thang hàng với độ phức tạp O(n³), giúp hiệu quả hơn rất nhiều so với quy tắc Cramer cho các hệ lớn hơn 3×3.
• Phân rã LU phân tích ma trận hệ số thành ma trận tam giác dưới và tam giác trên, cho phép giải hiệu quả nhiều hệ có cùng ma trận hệ số.
• Phân rã QR mang lại độ ổn định số tốt hơn cho các hệ điều kiện kém bằng cách phân rã ma trận thành thành phần trực giao và tam giác trên.

Khi nào dùng quy tắc Cramer

Bất chấp các giới hạn này, quy tắc Cramer vẫn hữu ích trong các bối cảnh cụ thể như:

• Phân tích và chứng minh lý thuyết trong đại số tuyến tính
• Tính toán ký hiệu khi cần nghiệm đại số chính xác
• Các hệ nhỏ (2×2 hoặc 3×3) nơi công thức tường minh dễ xử lý
• Bối cảnh giáo dục để hiểu mối liên hệ giữa định thức và hệ tuyến tính
• Suy ra nghiệm dạng đóng trong công việc phân tích

Tuy nhiên, đáng lưu ý rằng trong thực hành khoa học dữ liệu, quy tắc Cramer đóng vai trò như một công cụ khái niệm và lý thuyết để hiểu các hệ tuyến tính hơn là một phương pháp tính toán.

Kết luận

Quy tắc Cramer là một cách tiếp cận trực tiếp dựa trên công thức để giải hệ phương trình tuyến tính thông qua định thức. Chúng ta đã khám phá cách phương pháp này tính giá trị từng biến, thực hành qua các ví dụ 2×2 và 3×3, và hiểu diễn giải hình học của nó. Cuối cùng, chúng ta cũng tìm hiểu về các giới hạn, lựa chọn thay thế và khi nào nên dùng quy tắc Cramer.

Để đào sâu hiểu biết của bạn về các hệ tuyến tính và nghiệm của chúng, hãy xem khóa học Đại số tuyến tính cho Khoa học dữ liệu, nơi bạn sẽ học nhiều kỹ thuật đại số tuyến tính như quy tắc Cramer.