Math Labs

Đại số

Ma trận và định thức

Các bảng số hình chữ nhật và đại lượng vô hướng đo mức độ chúng làm co giãn thể tích.

Trực giácTrực giác: ma trận là một cỗ máy biến dạng không gian

Ma trận là một bảng số hình chữ nhật, nhưng công dụng lớn nhất của nó là một quy tắc biến một vectơ này thành một vectơ khác theo một công thức tuyến tính cố định: kéo dãn ở đây, xoay ở đó, có khi lật ngược. Định thức là một số duy nhất gắn với một ma trận vuông, trả lời một câu hỏi rất cụ thể: ma trận này làm diện tích (trong mặt phẳng) hay thể tích (trong không gian) co giãn theo hệ số nào, và trong quá trình đó nó có lật ngược hướng (từ tay trái thành tay phải) hay không? Định thức bằng 00 nghĩa là ma trận làm không gian bị dẹt hẳn lại, mất đi hẳn một chiều — chỉ một sự kiện đó lại chi phối gần như mọi thứ khác về ma trận, kể cả việc nó có thể "hoàn tác" được hay không.

Ma trận 2×2 biến ô vuông đơn vị thành hình bình hành; định thức hiển thị là hệ số diện tích có dấu.
Cách ma trận 2×22\times 2 A=(a11a12a21a22)A = \begin{pmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{pmatrix} biến dạng mặt phẳng: ô vuông đơn vị biến thành hình bình hành có diện tích đại số det⁡A=a11a22−a12a21\det A = a_{11}a_{22} - a_{12}a_{21}.

Phổ thôngPhép nhân ma trận và công thức định thức

Định nghĩa: Phép nhân ma trận

Nếu AA có cỡ m×nm\times n và BB có cỡ n×pn\times p, tích C=ABC=AB là ma trận cỡ m×pm\times p có các phần tử Cij=∑k=1nAikBkjC_{ij} = \sum_{k=1}^{n} A_{ik}B_{kj}: phần tử (i,j)(i,j) của CC là tích vô hướng của dòng ii của AA với cột jj của BB.

Cij=∑k=1nAikBkjC_{ij} = \sum_{k=1}^{n} A_{ik}B_{kj}

Trong đó AikA_{ik} là phần tử ở dòng ii, cột kk của AA, và BkjB_{kj} là phần tử ở dòng kk, cột jj của BB; tổng chạy theo chỉ số trong chung kk, đó là lý do vì sao số cột của AA phải bằng số dòng của BB để tích được định nghĩa.

Định nghĩa: Định thức (công thức Leibniz)

Với ma trận vuông n×nn\times n là AA, định thức được định nghĩa bởi det⁡(A)=∑σ∈Snsgn⁡(σ)∏i=1nai,σ(i)\det(A)=\sum_{\sigma\in S_n}\operatorname{sgn}(\sigma)\prod_{i=1}^n a_{i,\sigma(i)}, lấy tổng theo mọi hoán vị σ\sigma của {1,…,n}\{1,\dots,n\}, với sgn⁡(σ)=+1\operatorname{sgn}(\sigma)=+1 hoặc −1-1 tùy theo σ\sigma là hoán vị chẵn hay lẻ. Với n=2,3n=2,3, công thức này khai triển thành các công thức khai triển theo phần bù đại số quen thuộc dưới đây.

det⁡(abcd)=ad−bc\det\begin{pmatrix}a&b\\ c&d\end{pmatrix}=ad-bc

Với ma trận 2×22\times 2, det⁡(abcd)=ad−bc\det\begin{pmatrix}a&b\\ c&d\end{pmatrix}=ad-bc chính là diện tích có dấu của hình bình hành sinh bởi các vectơ dòng (a,b)(a,b) và (c,d)(c,d): hai hoán vị của {1,2}\{1,2\} là hoán vị đồng nhất (góp +ad+ad) và hoán vị đổi chỗ (góp −bc-bc).

det⁡(a11a12a13a21a22a23a31a32a33)=a11(a22a33−a23a32)−a12(a21a33−a23a31)+a13(a21a32−a22a31)\det\begin{pmatrix}a_{11}&a_{12}&a_{13}\\ a_{21}&a_{22}&a_{23}\\ a_{31}&a_{32}&a_{33}\end{pmatrix}=a_{11}(a_{22}a_{33}-a_{23}a_{32})-a_{12}(a_{21}a_{33}-a_{23}a_{31})+a_{13}(a_{21}a_{32}-a_{22}a_{31})

Với ma trận 3×33\times 3, det⁡(a11a12a13a21a22a23a31a32a33)=a11(a22a33−a23a32)−a12(a21a33−a23a31)+a13(a21a32−a22a31)\det\begin{pmatrix}a_{11}&a_{12}&a_{13}\\ a_{21}&a_{22}&a_{23}\\ a_{31}&a_{32}&a_{33}\end{pmatrix}=a_{11}(a_{22}a_{33}-a_{23}a_{32})-a_{12}(a_{21}a_{33}-a_{23}a_{31})+a_{13}(a_{21}a_{32}-a_{22}a_{31}) là khai triển theo phần bù đại số dọc theo dòng đầu: khai triển theo dòng 11, nhân mỗi phần tử a1ja_{1j} với định thức 2×22\times 2 của ma trận còn lại sau khi xóa dòng 11 và cột jj, và đổi dấu xen kẽ +,−,++,-,+ theo j=1,2,3j=1,2,3.

Các phép biến đổi dòng cơ bản ảnh hưởng thế nào đến định thức
Phép biến đổi dòngẢnh hưởng lên det⁡\det
Đổi chỗ hai dòngnhân det⁡\det với −1-1
Nhân một dòng với kknhân det⁡\det với kk
Cộng kk lần một dòng vào dòng khácdet⁡\det không đổi
Hai dòng bằng nhau hoặc tỉ lệdet⁡=0\det=0

Đại họcĐịnh lý: tính nhân của định thức và điều kiện khả nghịch

Với hai ma trận vuông n×nn\times n bất kỳ A,BA,B, ta có det⁡(AB)=det⁡(A)det⁡(B)\det(AB)=\det(A)\det(B).

Vì sao đúng?

Vì định thức đo mức mà một ánh xạ tuyến tính làm co giãn thể tích, việc áp dụng ánh xạ BB trước rồi ánh xạ AA sau lẽ ra phải làm thể tích co giãn theo hệ số của BB rồi theo hệ số của AA, tức là theo tích hai hệ số đó — tính nhân chính là phát biểu đại số cho việc hợp hai ánh xạ tuyến tính thì hợp luôn cả hai hệ số co giãn thể tích của chúng.

Chứng minh

Trước hết, nhắc lại ảnh hưởng của ba phép biến đổi dòng cơ bản lên định thức (được thiết lập trực tiếp từ tổng Leibniz det⁡(A)=∑σ∈Snsgn⁡(σ)∏i=1nai,σ(i)\det(A)=\sum_{\sigma\in S_n}\operatorname{sgn}(\sigma)\prod_{i=1}^n a_{i,\sigma(i)}, vì mỗi số hạng của tổng đó là tuyến tính theo từng dòng riêng biệt): đổi chỗ hai dòng nhân det⁡\det với −1-1, nhân một dòng với kk nhân det⁡\det với kk, và cộng một bội của một dòng vào dòng khác thì det⁡\det không đổi.

Mỗi phép biến đổi dòng cơ bản trên AA giống như nhân AA ở bên trái với một ma trận cơ bản EE tương ứng (thu được bằng cách áp dụng đúng phép biến đổi đó lên ma trận đơn vị II). So sánh với đoạn trên, det⁡(E)\det(E) chính bằng hệ số co giãn của phép biến đổi đó (−1-1, kk, hoặc 11 tương ứng), nên với mọi ma trận cơ bản EE, det⁡(EA)=det⁡(E)det⁡(A)\det(EA)=\det(E)\det(A).

Trường hợp AA khả nghịch: phép khử Gauss đưa mọi ma trận khả nghịch về ma trận đơn vị bằng một dãy hữu hạn các phép biến đổi dòng cơ bản, nghĩa là Em⋯E1A=IE_m\cdots E_1A=I với các ma trận cơ bản E1,…,EmE_1,\dots,E_m, nên A=E1−1⋯Em−1A=E_1^{-1}\cdots E_m^{-1}, tức chính nó là tích các ma trận cơ bản (nghịch đảo của một ma trận cơ bản cũng là ma trận cơ bản cùng loại). Áp dụng liên tiếp đẳng thức ở đoạn trước, det⁡(A)=det⁡(E1−1)⋯det⁡(Em−1)\det(A)=\det(E_1^{-1})\cdots\det(E_m^{-1}), và áp dụng lại cho AB=E1−1⋯Em−1BAB=E_1^{-1}\cdots E_m^{-1}B ta được det⁡(AB)=det⁡(E1−1)⋯det⁡(Em−1)det⁡(B)=det⁡(A)det⁡(B)\det(AB)=\det(E_1^{-1})\cdots\det(E_m^{-1})\det(B)=\det(A)\det(B).

Trường hợp AA suy biến: khi đó det⁡(A)=0\det(A)=0 (một ma trận suy biến không thể đưa hoàn toàn về ma trận đơn vị, và các quy tắc biến đổi dòng trên cho thấy mọi dạng bậc thang đạt được vẫn còn một dòng bằng không, buộc det⁡=0\det=0 theo khai triển phần bù đại số theo dòng đó). Suy biến cũng có nghĩa rank⁡(A)<n\operatorname{rank}(A)<n, và vì rank⁡(AB)≤rank⁡(A)<n\operatorname{rank}(AB)\leq\operatorname{rank}(A)<n với mọi BB, tích ABAB cũng suy biến, nên det⁡(AB)=0\det(AB)=0. Do đó det⁡(AB)=0=0⋅det⁡(B)=det⁡(A)det⁡(B)\det(AB)=0=0\cdot\det(B)=\det(A)\det(B), và đẳng thức vẫn đúng trong trường hợp này, hoàn tất chứng minh cho mọi AA.

Một ma trận vuông AA khả nghịch nếu và chỉ nếu det⁡(A)≠0\det(A)\neq 0; khi đó, A−1=1det⁡(A)adj⁡(A)A^{-1}=\dfrac{1}{\det(A)}\operatorname{adj}(A), trong đó adj⁡(A)\operatorname{adj}(A) là ma trận phụ hợp (chuyển vị của ma trận các phần bù đại số).

Vì sao đúng?

Một ánh xạ tuyến tính có thể "hoàn tác" được đúng khi nó không làm không gian sụp xuống một chiều thấp hơn — vì định thức đã đo chính xác câu hỏi "nó có làm thể tích sụp về không", nên tính khả nghịch và det⁡≠0\det\neq 0 thực ra là cùng một phát biểu nhìn từ hai góc: một góc đại số (có tồn tại ma trận nghịch đảo hay không), một góc hình học (hệ số co giãn thể tích có khác không hay không).

Chứng minh

(⇒\Rightarrow) Giả sử AA khả nghịch, nên tồn tại ma trận A−1A^{-1} với AA−1=IAA^{-1}=I. Áp dụng định lý tính nhân đã chứng minh trên, det⁡(A)det⁡(A−1)=det⁡(AA−1)=det⁡(I)=1\det(A)\det(A^{-1})=\det(AA^{-1})=\det(I)=1. Vì tích của hai số det⁡(A)\det(A) và det⁡(A−1)\det(A^{-1}) bằng 11, không số nào trong chúng có thể bằng 00 — nên det⁡(A)≠0\det(A)\neq 0.

(⇐\Leftarrow), thực hiện cụ thể với 2×22\times 2: giả sử A=(abcd)A=\begin{pmatrix}a&b\\ c&d\end{pmatrix} có det⁡(A)=ad−bc≠0\det(A)=ad-bc\neq 0. Đặt adj⁡(A)=(d−b−ca)\operatorname{adj}(A)=\begin{pmatrix}d&-b\\ -c&a\end{pmatrix}. Nhân trực tiếp, A⋅adj⁡(A)=(abcd)(d−b−ca)=(ad−bc−ab+abcd−cd−cb+ad)=(ad−bc00ad−bc)=det⁡(A) IA\cdot\operatorname{adj}(A)=\begin{pmatrix}a&b\\ c&d\end{pmatrix}\begin{pmatrix}d&-b\\ -c&a\end{pmatrix}=\begin{pmatrix}ad-bc&-ab+ab\\ cd-cd&-cb+ad\end{pmatrix}=\begin{pmatrix}ad-bc&0\\ 0&ad-bc\end{pmatrix}=\det(A)\,I, đúng là đẳng thức A⋅adj⁡(A)=det⁡(A) IA\cdot\operatorname{adj}(A)=\det(A)\,I. Cùng phép tính theo thứ tự ngược lại cũng cho adj⁡(A)⋅A=det⁡(A) I\operatorname{adj}(A)\cdot A=\det(A)\,I.

Vì det⁡(A)≠0\det(A)\neq 0, chia cả hai vế của A⋅adj⁡(A)=det⁡(A)IA\cdot\operatorname{adj}(A)=\det(A)I cho số vô hướng det⁡(A)\det(A) ta được A⋅adj⁡(A)det⁡(A)=IA\cdot\dfrac{\operatorname{adj}(A)}{\det(A)}=I, và tương tự theo thứ tự ngược lại adj⁡(A)det⁡(A)⋅A=I\dfrac{\operatorname{adj}(A)}{\det(A)}\cdot A=I. Điều này cho một nghịch đảo hai phía tường minh A−1=1det⁡(A)adj⁡(A)A^{-1}=\dfrac{1}{\det(A)}\operatorname{adj}(A), nên AA khả nghịch.

Với nn tổng quát (minh họa ở n=3n=3), lập luận tương tự vẫn đúng vì khai triển Laplace theo phần bù đại số cho ∑jaijCij=det⁡(A)\sum_j a_{ij}C_{ij}=\det(A) (khai triển theo dòng ii) còn ∑jaijCkj=0\sum_j a_{ij}C_{kj}=0 với i≠ki\neq k (tổng này là khai triển phần bù đại số của một ma trận có dòng ii được chép vào dòng kk, tức có hai dòng bằng nhau nên định thức bằng 00 theo bảng trên). Ghép các đẳng thức theo từng dòng này lại thành dạng ma trận cho đúng A⋅adj⁡(A)=det⁡(A)IA\cdot\operatorname{adj}(A)=\det(A)I với mọi cỡ nn, nên bất cứ khi nào det⁡(A)≠0\det(A)\neq 0, lập luận chia như ở trường hợp 2×22\times 2 cho ra A−1=adj⁡(A)/det⁡(A)A^{-1}=\operatorname{adj}(A)/\det(A); ngược lại, nếu det⁡(A)=0\det(A)=0, phần phản đảo của đoạn đầu loại hẳn khả năng khả nghịch, vì det⁡(A)det⁡(A−1)=1\det(A)\det(A^{-1})=1 sẽ buộc det⁡(A)≠0\det(A)\neq 0. Điều này chứng minh điều kiện theo cả hai chiều.

Đại họcỨng dụng thực tiễn và Ví dụ minh họa

Định thức cho phép kỹ sư giải các hệ tuyến tính nhỏ bằng một công thức tường minh (quy tắc Cramer) thay vì khử dần, và cho phép các phần mềm đồ họa máy tính kiểm tra chỉ bằng một số xem một phép biến đổi có giữ nguyên hay lật hướng không và co giãn diện tích theo hệ số nào — cả hai đều là ứng dụng trực tiếp của việc det⁡≠0\det\neq 0 nghĩa là "khả nghịch" và ∣det⁡∣|\det| nghĩa là "hệ số co giãn."

Ví dụ: Kỹ thuật: giải mạch điện hai vòng bằng quy tắc Cramer

Áp dụng các định luật Kirchhoff cho một mạch điện hai vòng cho hệ tuyến tính 3x+y=113x+y=11 và x+2y=8x+2y=8 với hai dòng điện vòng x,yx,y (đơn vị ampere). Hãy giải xx và yy bằng định thức.

Lời giải

Viết hệ dưới dạng phương trình ma trận (3112)(xy)=(118)\begin{pmatrix}3&1\\1&2\end{pmatrix}\begin{pmatrix}x\\y\end{pmatrix}=\begin{pmatrix}11\\8\end{pmatrix}. Ma trận hệ số có det⁡=3×2−1×1=6−1=5≠0\det=3\times2-1\times1=6-1=5\neq 0, nên theo điều kiện khả nghịch, hệ có nghiệm duy nhất, và quy tắc Cramer áp dụng được.

Quy tắc Cramer cho x=det⁡(11182)det⁡(3112)x=\dfrac{\det\begin{pmatrix}11&1\\8&2\end{pmatrix}}{\det\begin{pmatrix}3&1\\1&2\end{pmatrix}}: thay cột thứ nhất của ma trận hệ số bằng vế phải, lấy định thức, rồi chia cho định thức ban đầu. Tử số là 11×2−1×8=22−8=1411\times2-1\times8=22-8=14, nên x=14/5=2.8x=14/5=2.8 A.

Tương tự y=det⁡(31118)5y=\dfrac{\det\begin{pmatrix}3&11\\1&8\end{pmatrix}}{5}: thay cột thứ hai thay vì cột thứ nhất. Tử số là 3×8−11×1=24−11=133\times8-11\times1=24-11=13, nên y=13/5=2.6y=13/5=2.6 A.

Kiểm tra lại: 3(2.8)+2.6=8.4+2.6=113(2.8)+2.6=8.4+2.6=11 ✓ và 2.8+2(2.6)=2.8+5.2=82.8+2(2.6)=2.8+5.2=8 ✓, xác nhận x=2.8x=2.8 A và y=2.6y=2.6 A giải đúng các phương trình mạch.

Ví dụ: Đồ họa máy tính: một phép biến đổi có lật hình không, và co giãn diện tích theo hệ số nào?

Một bộ máy đồ họa 2D áp dụng ma trận biến đổi M=(2143)M=\begin{pmatrix}2&1\\4&3\end{pmatrix} lên mọi điểm của một hình có diện tích 55 đơn vị vuông. Hãy xác định hình sau biến đổi có bị lật gương so với hình ban đầu không, và tính diện tích của hình sau biến đổi.

Lời giải

Tính định thức của MM: det⁡(M)=2×3−1×4=6−4=2\det(M)=2\times3-1\times4=6-4=2.

Dấu của det⁡(M)\det(M) cho biết về hướng: vì det⁡(M)=2>0\det(M)=2>0, phép biến đổi giữ nguyên hướng, nên hình không bị lật gương (định thức âm sẽ cho biết có lật, vì nghĩa là có một số lẻ phép phản xạ được 'nhúng' trong ma trận).

Độ lớn ∣det⁡(M)∣|\det(M)| cho biết hệ số co giãn diện tích: diện tích của mọi vùng đều được nhân với ∣det⁡(M)∣=2|\det(M)|=2 sau phép biến đổi, bất kể hình dạng của vùng, vì đó chính là ý nghĩa hình học của định thức đã nêu ở trên.

Vậy hình sau biến đổi có diện tích 5×∣det⁡(M)∣=5×2=105\times|\det(M)|=5\times2=10 đơn vị vuông, và giữ nguyên hướng (không bị lật) so với hình ban đầu.

Tính det⁡(2513)\det\begin{pmatrix}2&5\\1&3\end{pmatrix}?

Nếu det⁡(A)=3\det(A)=3 và det⁡(B)=4\det(B)=4 với hai ma trận vuông n×nn\times n, hỏi det⁡(AB)\det(AB) bằng bao nhiêu?

Ma trận nào dưới đây không có ma trận nghịch đảo?

Một ma trận biến đổi đồ họa có det⁡=−3\det=-3. Điều này cho biết gì?

Tài liệu tham khảo

  1. Eric W. Weisstein (MathWorld) (2024). Determinant
  2. Gilbert Strang (2016). Introduction to Linear Algebra