Ma trận xác định dương là gì? Ý nghĩa và ứng dụng chuyên sâu

Ma trận xác định dương (positive definite matrix) là một ma trận đối xứng thực mà dạng toàn phương tương ứng luôn nhận giá trị thực dương nghiêm ngặt đối với mọi vectơ khác không, đóng vai trò nền tảng trong tối ưu hóa lồi và học máy.

Định nghĩa Ma trận Xác định Dương

Ma trận vuông $A \in \mathbb{R}^{n\times n}$ gọi là xác định dương nếu với mọi vectơ cột không không $x \in \mathbb{R}^n$ ta luôn có

 $x^T A x > 0$

Yêu cầu này ngụ ý rằng biểu thức song tuyến tính do $A$ sinh ra đóng vai trò tích lũy năng lượng dương trên mọi hướng của không gian vectơ. Kết quả là hàm số $f(x)=x^T A x$ chỉ có giá trị lớn hơn không, trừ trường hợp $x=0$ .

Định nghĩa này đòi hỏi hai điều kiện: $A$ phải đối xứng (tức $A=A^T$ ) và mọi cặp phần tử trong phép tính $x^T A x$ đều góp giá trị dương. Hầu hết các ma trận trong ứng dụng khoa học tính toán, tối ưu hóa và thống kê đều được giả thiết xác định dương để bảo đảm tính ổn định và khả năng giải hệ.

Tính chất Cơ bản

Một ma trận $A$ xác định dương có các tính chất then chốt sau đây:

Đối xứng: $A = A^T$ . Mọi ma trận xác định dương buộc phải đối xứng để đảm bảo $x^T A x$ luôn thực.
Giá trị riêng dương: Nếu $A v = \lambda v$ với vectơ riêng $v\neq0$ , thì $\lambda>0$ .
Định thức và các principal minors: $\det(A)>0$ và mọi định thức con chính đều dương.
Không có vectơ không dương tính: Không tồn tại $x\neq0$ sao cho $Ax=0$ hoặc $x^T A x \le 0$ .

Nhờ những tính chất này, $A$ cho phép định nghĩa căn bậc hai $A^{1/2}$ và nghịch đảo $A^{-1}$ cũng là xác định dương, đồng thời các phép toán như phân tích Cholesky và phân tích phổ luôn khả thi và ổn định về mặt số.

Điều kiện Sylvester

Theo tiêu chí Sylvester, một ma trận đối xứng $A$ xác định dương nếu và chỉ nếu tất cả các định thức con chính từ cấp 1 đến cấp $n$ đều dương. Cụ thể, ký hiệu $\Delta_k$ là định thức của ma trận con $A_{1:k,1:k}$ , ta có:


\Delta_1 > 0,\ \Delta_2 > 0,\ \dots,\ \Delta_n = \det(A) > 0
\end{script>

Cấp $k Định thức con chính <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi mathvariant="normal">Δ</mi><mi>k</mi></msub></mrow><annotation encoding="application/x-tex">\Delta_k</annotation></semantics></math> Yêu cầu 1 <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi mathvariant="normal">Δ</mi><mn>1</mn></msub><mo>=</mo><msub><mi>a</mi><mn>11</mn></msub></mrow><annotation encoding="application/x-tex">\Delta_1 = a_{11}</annotation></semantics></math> > 0 2 <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi mathvariant="normal">Δ</mi><mn>2</mn></msub><mo>=</mo><mrow><mo fence="true">∣</mo><mtable rowspacing="0.16em" columnalign="center center" columnspacing="1em"><mtr><mtd><mstyle scriptlevel="0" displaystyle="false"><msub><mi>a</mi><mn>11</mn></msub></mstyle></mtd><mtd><mstyle scriptlevel="0" displaystyle="false"><msub><mi>a</mi><mn>12</mn></msub></mstyle></mtd></mtr><mtr><mtd><mstyle scriptlevel="0" displaystyle="false"><msub><mi>a</mi><mn>21</mn></msub></mstyle></mtd><mtd><mstyle scriptlevel="0" displaystyle="false"><msub><mi>a</mi><mn>22</mn></msub></mstyle></mtd></mtr></mtable><mo fence="true">∣</mo></mrow></mrow><annotation encoding="application/x-tex">\Delta_2 = \begin{vmatrix}a_{11} & a_{12}\\a_{21}&a_{22}\end{vmatrix}</annotation></semantics></math> > 0 … … … <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>n</mi></mrow><annotation encoding="application/x-tex">n</annotation></semantics></math> <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi mathvariant="normal">Δ</mi><mi>n</mi></msub><mo>=</mo><mi>det</mi><mo>⁡</mo><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">\Delta_n = \det(A)</annotation></semantics></math> > 0 Việc kiểm tra tuần tự từng <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi mathvariant="normal">Δ</mi><mi>k</mi></msub></mrow><annotation encoding="application/x-tex">\Delta_k</annotation></semantics></math> cho phép sớm phát hiện ma trận không thỏa mãn điều kiện, tiết kiệm chi phí tính toán so với phân tích toàn bộ phổ giá trị riêng. Thông tin này thường được áp dụng trong các thư viện số như LAPACK (LAPACK Users’ Guide). <h2>Phân tích Giá trị Riêng</h2> Vì <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi></mrow><annotation encoding="application/x-tex">A</annotation></semantics></math> đối xứng nên tồn tại phân tích trực giao <pre><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi><mo>=</mo><mi>Q</mi><mi mathvariant="normal">Λ</mi><msup><mi>Q</mi><mi>T</mi></msup><mo separator="true">,</mo><mspace width="1em"/><msup><mi>Q</mi><mi>T</mi></msup><mi>Q</mi><mo>=</mo><mi>I</mi><mo separator="true">,</mo><mtext> </mtext><mi mathvariant="normal">Λ</mi><mo>=</mo><mrow><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">g</mi></mrow><mo stretchy="false">(</mo><msub><mi>λ</mi><mn>1</mn></msub><mo separator="true">,</mo><mo>…</mo><mo separator="true">,</mo><msub><mi>λ</mi><mi>n</mi></msub><mo stretchy="false">)</mo><mi mathvariant="normal">.</mi></mrow><annotation encoding="application/x-tex">A = Q \Lambda Q^T,\quad Q^TQ=I,\ \Lambda=\mathrm{diag}(\lambda_1,\dots,\lambda_n).</annotation></semantics></math></pre> Một ma trận xác định dương tương đương với việc mọi giá trị riêng <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi>λ</mi><mi>i</mi></msub></mrow><annotation encoding="application/x-tex">\lambda_i</annotation></semantics></math> trong <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi mathvariant="normal">Λ</mi></mrow><annotation encoding="application/x-tex">\Lambda</annotation></semantics></math> đều thỏa mãn <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi>λ</mi><mi>i</mi></msub><mo>></mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\lambda_i>0</annotation></semantics></math>. Việc kiểm tra phổ giá trị riêng được thực hiện bằng các thuật toán QR hoặc divide-and-conquer với độ chính xác cao (MathWorld). Giá trị riêng cũng cho phép định nghĩa condition number <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>κ</mi><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><msub><mi>λ</mi><mi>max</mi><mo>⁡</mo></msub><msub><mi>λ</mi><mi>min</mi><mo>⁡</mo></msub></mfrac></mrow><annotation encoding="application/x-tex">\kappa(A)=\frac{\lambda_{\max}}{\lambda_{\min}}</annotation></semantics></math>, đại lượng quan trọng đánh giá độ nhạy của việc giải hệ <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi><mi>x</mi><mo>=</mo><mi>b</mi></mrow><annotation encoding="application/x-tex">Ax=b</annotation></semantics></math> trước sai số làm tròn. Ma trận có <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>κ</mi><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">\kappa(A)</annotation></semantics></math> thấp (gần 1) cho thấy tính ổn định số học cao, trong khi <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>κ</mi><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">\kappa(A)</annotation></semantics></math> rất lớn dễ dẫn đến sai số tích lũy. <h2>Phân tích Cholesky</h2> Phân tích Cholesky phân tách ma trận xác định dương <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi></mrow><annotation encoding="application/x-tex">A</annotation></semantics></math> thành tích của ma trận tam giác dưới <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>L</mi></mrow><annotation encoding="application/x-tex">L</annotation></semantics></math> và chuyển vị của nó: <pre><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi><mo>=</mo><mi>L</mi><mtext> </mtext><msup><mi>L</mi><mi>T</mi></msup><mo separator="true">,</mo><mspace width="1em"/><mi>L</mi><mo>∈</mo><msup><mi mathvariant="double-struck">R</mi><mrow><mi>n</mi><mo>×</mo><mi>n</mi></mrow></msup><mrow><mtext> tam gi</mtext><mover accent="true"><mtext>a</mtext><mo>ˊ</mo></mover><mtext>c dưới, </mtext></mrow><msub><mi>L</mi><mrow><mi>i</mi><mi>i</mi></mrow></msub><mo>></mo><mn>0.</mn></mrow><annotation encoding="application/x-tex">A = L\,L^T,\quad L\in\mathbb{R}^{n\times n}\text{ tam giác dưới, }L_{ii}>0.</annotation></semantics></math></pre> Giải thuật Cholesky lặp qua hàng và cột, tính thành phần <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi>L</mi><mrow><mi>i</mi><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">L_{ii}</annotation></semantics></math> và <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi>L</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">L_{ji}</annotation></semantics></math> theo công thức: <pre><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi>L</mi><mrow><mi>i</mi><mi>i</mi></mrow></msub><mo>=</mo><msqrt><mrow><msub><mi>A</mi><mrow><mi>i</mi><mi>i</mi></mrow></msub><mo>−</mo><msubsup><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msubsup><msubsup><mi>L</mi><mrow><mi>i</mi><mi>k</mi></mrow><mn>2</mn></msubsup></mrow></msqrt><mo separator="true">,</mo><mspace width="1em"/><msub><mi>L</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>L</mi><mrow><mi>i</mi><mi>i</mi></mrow></msub></mfrac><mo fence="true" stretchy="true" minsize="1.2em" maxsize="1.2em">(</mo><msub><mi>A</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub><mo>−</mo><msubsup><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msubsup><msub><mi>L</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><msub><mi>L</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo fence="true" stretchy="true" minsize="1.2em" maxsize="1.2em">)</mo><mi mathvariant="normal">.</mi></mrow><annotation encoding="application/x-tex"> L_{ii} = \sqrt{A_{ii} - \sum_{k=1}^{i-1} L_{ik}^2},\quad L_{ji} = \frac{1}{L_{ii}}\bigl(A_{ji} - \sum_{k=1}^{i-1}L_{jk}L_{ik}\bigr). </annotation></semantics></math></pre> Độ phức tạp tính toán của Cholesky xấp xỉ <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mstyle scriptlevel="0" displaystyle="false"><mfrac><mn>1</mn><mn>3</mn></mfrac></mstyle><msup><mi>n</mi><mn>3</mn></msup></mrow><annotation encoding="application/x-tex">\tfrac{1}{3}n^3</annotation></semantics></math> phép toán, nhanh hơn phân tích phổ giá trị riêng hay khử Gaussian. Phân tích Cholesky thường được sử dụng để giải hệ phương trình tuyến tính <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi><mi>x</mi><mo>=</mo><mi>b</mi></mrow><annotation encoding="application/x-tex">Ax=b</annotation></semantics></math> nhanh chóng và ổn định, đồng thời phát hiện ngay lập tức nếu <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi></mrow><annotation encoding="application/x-tex">A</annotation></semantics></math> không xác định dương khi biểu thức dưới dấu căn bị âm (LAPACK Users’ Guide). <table> <thead> <tr> <th>Bước</th> <th>Phép toán chính</th> <th>Độ phức tạp</th> </tr> </thead> <tbody> <tr> <td>Tính <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi>L</mi><mrow><mi>i</mi><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">L_{ii}</annotation></semantics></math></td> <td>Căn bậc hai của số thực dương</td> <td><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi mathvariant="script">O</mi><mo stretchy="false">(</mo><msup><mi>n</mi><mn>2</mn></msup><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">\mathcal{O}(n^2)</annotation></semantics></math> tổng thể</td> </tr> <tr> <td>Tính <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msub><mi>L</mi><mrow><mi>j</mi><mi>i</mi></mrow></msub></mrow><annotation encoding="application/x-tex">L_{ji}</annotation></semantics></math></td> <td>Thương của hiệu hai tổng tích</td> <td><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi mathvariant="script">O</mi><mo stretchy="false">(</mo><msup><mi>n</mi><mn>3</mn></msup><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">\mathcal{O}(n^3)</annotation></semantics></math> tổng thể</td> </tr> </tbody> </table> <h2>Ứng dụng trong Tối ưu hóa</h2> Trong tối ưu hóa lồi, hàm bậc hai <pre><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle scriptlevel="0" displaystyle="false"><mfrac><mn>1</mn><mn>2</mn></mfrac></mstyle><msup><mi>x</mi><mi>T</mi></msup><mi>A</mi><mi>x</mi><mo>−</mo><msup><mi>b</mi><mi>T</mi></msup><mi>x</mi><mo>+</mo><mi>c</mi></mrow><annotation encoding="application/x-tex">f(x)=\tfrac12x^TAx - b^Tx + c</annotation></semantics></math></pre> là lồi toàn cục nếu <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi></mrow><annotation encoding="application/x-tex">A</annotation></semantics></math> xác định dương, đảm bảo tồn tại nghiệm tối ưu duy nhất giải <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi mathvariant="normal">∇</mi><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mi>A</mi><mi>x</mi><mo>−</mo><mi>b</mi><mo>=</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\nabla f(x)=Ax-b=0</annotation></semantics></math>. Phương pháp gradient descent và Newton tận dụng đặc tính này để hội tụ nhanh chóng: <ul> <li>Gradient Descent: bước nhảy <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>α</mi></mrow><annotation encoding="application/x-tex">\alpha</annotation></semantics></math> chọn trong khoảng <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mo stretchy="false">(</mo><mn>0</mn><mo separator="true">,</mo><mtext> </mtext><mn>2</mn><mi mathvariant="normal">/</mi><msub><mi>λ</mi><mi>max</mi><mo>⁡</mo></msub><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">(0,\,2/\lambda_{\max})</annotation></semantics></math> để đảm bảo giảm hàm mục tiêu mỗi bước.</li> <li>Newton’s Method: sử dụng phân tích Cholesky để giải hệ <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi><mtext> </mtext><mi>p</mi><mo>=</mo><mo>−</mo><mi mathvariant="normal">∇</mi><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">A\,p = -\nabla f(x)</annotation></semantics></math>, hội tụ bậc hai gần nghiệm tối ưu (MIT OpenCourseWare).</li> </ul> Độ hội tụ tuyến tính của gradient descent tỉ lệ nghịch với condition number <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>κ</mi><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">\kappa(A)</annotation></semantics></math>, trong khi Newton ổn định hơn nhưng yêu cầu tính nghịch đảo hoặc phân tích Cholesky của <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi></mrow><annotation encoding="application/x-tex">A</annotation></semantics></math> tại mỗi bước, chi phí <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi mathvariant="script">O</mi><mo stretchy="false">(</mo><msup><mi>n</mi><mn>3</mn></msup><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">\mathcal{O}(n^3)</annotation></semantics></math>. <h2>Ứng dụng trong Thống kê</h2> Ma trận hiệp phương sai <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi mathvariant="normal">Σ</mi></mrow><annotation encoding="application/x-tex">\Sigma</annotation></semantics></math> của phân phối chuẩn đa biến <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>N</mi><mo stretchy="false">(</mo><mi>μ</mi><mo separator="true">,</mo><mi mathvariant="normal">Σ</mi><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">N(\mu,\Sigma)</annotation></semantics></math> phải xác định dương để mật độ xác suất có nghĩa: <pre><math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><msqrt><mrow><mo stretchy="false">(</mo><mn>2</mn><mi>π</mi><msup><mo stretchy="false">)</mo><mi>n</mi></msup><mi>det</mi><mo>⁡</mo><mi mathvariant="normal">Σ</mi></mrow></msqrt></mfrac><mi>exp</mi><mo>⁡</mo><mo fence="true" stretchy="true" minsize="1.2em" maxsize="1.2em">(</mo><mo>−</mo><mstyle scriptlevel="0" displaystyle="false"><mfrac><mn>1</mn><mn>2</mn></mfrac></mstyle><mo stretchy="false">(</mo><mi>x</mi><mo>−</mo><mi>μ</mi><msup><mo stretchy="false">)</mo><mi>T</mi></msup><msup><mi mathvariant="normal">Σ</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo stretchy="false">(</mo><mi>x</mi><mo>−</mo><mi>μ</mi><mo stretchy="false">)</mo><mo fence="true" stretchy="true" minsize="1.2em" maxsize="1.2em">)</mo><mi mathvariant="normal">.</mi></mrow><annotation encoding="application/x-tex"> f(x)=\frac{1}{\sqrt{(2\pi)^n\det\Sigma}} \exp\bigl(-\tfrac12(x-\mu)^T\Sigma^{-1}(x-\mu)\bigr). </annotation></semantics></math></pre> Xác định dương của <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi mathvariant="normal">Σ</mi></mrow><annotation encoding="application/x-tex">\Sigma</annotation></semantics></math> đảm bảo tồn tại <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msup><mi mathvariant="normal">Σ</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><annotation encoding="application/x-tex">\Sigma^{-1}</annotation></semantics></math> và <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>det</mi><mo>⁡</mo><mi mathvariant="normal">Σ</mi><mo>></mo><mn>0</mn></mrow><annotation encoding="application/x-tex">\det\Sigma>0</annotation></semantics></math>. Trong phân tích thành phần chính (PCA), phân tích giá trị riêng của <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi mathvariant="normal">Σ</mi></mrow><annotation encoding="application/x-tex">\Sigma</annotation></semantics></math> cho biết phương sai theo từng hướng, cho phép giảm chiều dữ liệu bằng cách chọn các thành phần chính với giá trị riêng lớn nhất (CMU Lecture). Trong hồi quy tuyến tính <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>y</mi><mo>=</mo><mi>X</mi><mi>β</mi><mo>+</mo><mi>ε</mi></mrow><annotation encoding="application/x-tex">y = X\beta + \varepsilon</annotation></semantics></math>, ma trận thiết kế <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msup><mi>X</mi><mi>T</mi></msup><mi>X</mi></mrow><annotation encoding="application/x-tex">X^TX</annotation></semantics></math> phải xác định dương để phép ước lượng bình phương nhỏ nhất <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mover accent="true"><mi>β</mi><mo>^</mo></mover><mo>=</mo><mo stretchy="false">(</mo><msup><mi>X</mi><mi>T</mi></msup><mi>X</mi><msup><mo stretchy="false">)</mo><mrow><mo>−</mo><mn>1</mn></mrow></msup><msup><mi>X</mi><mi>T</mi></msup><mi>y</mi></mrow><annotation encoding="application/x-tex">\hat\beta=(X^TX)^{-1}X^Ty</annotation></semantics></math> tồn tại và duy nhất. Đa cộng tuyến dẫn đến <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msup><mi>X</mi><mi>T</mi></msup><mi>X</mi></mrow><annotation encoding="application/x-tex">X^TX</annotation></semantics></math> suy biến hoặc gần suy biến, gây khối ill-conditioning. <h2>Kiểm tra Số và Ổn định Tính Toán</h2> Condition number <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>κ</mi><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><msub><mi>λ</mi><mi>max</mi><mo>⁡</mo></msub><msub><mi>λ</mi><mi>min</mi><mo>⁡</mo></msub></mfrac></mrow><annotation encoding="application/x-tex">\kappa(A)=\frac{\lambda_{\max}}{\lambda_{\min}}</annotation></semantics></math> đánh giá độ nhạy của lời giải hệ <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi><mi>x</mi><mo>=</mo><mi>b</mi></mrow><annotation encoding="application/x-tex">Ax=b</annotation></semantics></math> đối với sai số trong <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi></mrow><annotation encoding="application/x-tex">A</annotation></semantics></math> và <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>b</mi></mrow><annotation encoding="application/x-tex">b</annotation></semantics></math>. Ma trận có <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>κ</mi><mo stretchy="false">(</mo><mi>A</mi><mo stretchy="false">)</mo></mrow><annotation encoding="application/x-tex">\kappa(A)</annotation></semantics></math> cao dễ dẫn đến mất chính xác do sai số làm tròn. Sử dụng phân tích Cholesky thay cho khử Gauss có thể giảm sai số tích lũy nhờ hạn chế số phép chia và nhân (NIST Matrix Computation). Thuật toán kiểm tra Cholesky breakdown dừng khi gặp <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><msubsup><mi>L</mi><mrow><mi>i</mi><mi>i</mi></mrow><mn>2</mn></msubsup><mo>=</mo><msub><mi>A</mi><mrow><mi>i</mi><mi>i</mi></mrow></msub><mo>−</mo><msub><mo>∑</mo><mi>k</mi></msub><msubsup><mi>L</mi><mrow><mi>i</mi><mi>k</mi></mrow><mn>2</mn></msubsup><mo>≤</mo><mn>0</mn></mrow><annotation encoding="application/x-tex">L_{ii}^2 = A_{ii} - \sum_{k}L_{ik}^2 \le 0</annotation></semantics></math>, báo hiệu <math xmlns="http://www.w3.org/1998/Math/MathML"><semantics><mrow><mi>A</mi></mrow><annotation encoding="application/x-tex">A</annotation></semantics></math> không xác định dương. Đây là phương pháp nhanh và ổn định hơn so với tính toàn bộ phổ giá trị riêng.</pre></div></article><section class="sh-section topic-faq" data-v-f6ff6a5e><div class="sh-section-heading sh-section-heading--bar" data-v-f6ff6a5e><h2 data-v-f6ff6a5e>Câu hỏi thường gặp</h2></div><details class="faq-item" open data-v-f6ff6a5e><summary data-v-f6ff6a5e><h3 class="faq-question" data-v-f6ff6a5e>Tiêu chuẩn Sylvester kiểm tra tính xác định dương của ma trận như thế nào?</h3></summary>Một ma trận đối xứng thực là xác định dương khi và chỉ khi tất cả các định thức con chính góc (Leading Principal Minors) của nó đều mang giá trị dương nghiêm ngặt.</details><details class="faq-item" data-v-f6ff6a5e><summary data-v-f6ff6a5e><h3 class="faq-question" data-v-f6ff6a5e>Phân tích Cholesky áp dụng cho ma trận xác định dương có ý nghĩa gì?</h3></summary>Phân tích ma trận A thành tích L nhân L chuyển vị (với L là ma trận tam giác dưới có đường chéo dương), giúp giải hệ phương trình tuyến tính nhanh gấp đôi so với phép khử Gauss.</details><details class="faq-item" data-v-f6ff6a5e><summary data-v-f6ff6a5e><h3 class="faq-question" data-v-f6ff6a5e>Tại sao ma trận hiệp phương sai (Covariance Matrix) luôn là ma trận nửa xác định dương?</h3></summary>Vì phương sai của bất kỳ tổ hợp tuyến tính nào của các biến ngẫu nhiên luôn không âm, bảo đảm dạng toàn phương tương ứng luôn lớn hơn hoặc bằng không.</details></section><section class="sh-section topic-references" data-v-f6ff6a5e><div class="sh-section-heading sh-section-heading--bar" data-v-f6ff6a5e><h2 data-v-f6ff6a5e>Tài liệu tham khảo</h2></div><ol class="reference-list" data-v-f6ff6a5e><li data-v-f6ff6a5e>Martin (1971). Symmetric Decomposition of a Positive Definite Matrix. Linear Algebra.<a href="https://doi.org/10.1007/978-3-662-39778-7_1" target="_blank" rel="noopener nofollow" class="reference-link" data-v-f6ff6a5e> DOI: 10.1007/978-3-662-39778-7_1</a></li><li data-v-f6ff6a5e>Seaks (1972). SYMINV: An Algorithm for the Inversion of a Positive Definite Matrix by the Cholesky Decomposition. Econometrica.<a href="https://doi.org/10.2307/1912094" target="_blank" rel="noopener nofollow" class="reference-link" data-v-f6ff6a5e> DOI: 10.2307/1912094</a></li><li data-v-f6ff6a5e>Berenhaut (2005). Monotone convex sequences and Cholesky decomposition of symmetric Toeplitz matrices. Linear Algebra and its Applications.<a href="https://doi.org/10.1016/j.laa.2005.01.021" target="_blank" rel="noopener nofollow" class="reference-link" data-v-f6ff6a5e> DOI: 10.1016/j.laa.2005.01.021</a></li></ol></section></div><aside class="sh-aside" data-v-f6ff6a5e><section class="sidebar-box-container" data-v-f6ff6a5e data-v-6e6149cf><header class="sidebar-box-title" data-v-6e6149cf>Tra cứu thuật ngữ</header><div class="sidebar-box-content" data-v-6e6149cf><form class="topic-search-fallback" action="/topic/browse" method="get" role="search" data-v-f6ff6a5e><input type="search" name="terms" placeholder="Nhập thuật ngữ tiếng Việt hoặc tiếng Anh" aria-label="Tra cứu thuật ngữ" data-v-f6ff6a5e></form><a href="/topic/browse" class="topic-search-more" data-v-f6ff6a5e> Mở từ điển học thuật </a></div></section><section class="sidebar-box-container" data-v-f6ff6a5e data-v-6e6149cf><header class="sidebar-box-title" data-v-6e6149cf>Mục từ liên quan</header><div class="sidebar-box-content" data-v-6e6149cf><div class="topic-relate-list" data-v-f6ff6a5e data-v-cafed715><a href="/topic/l%E1%BA%ADp%20tr%C3%ACnh%20b%E1%BA%ADc%20hai" class="topic-relate-item" data-v-cafed715>Lập trình bậc haiquadratic programming</a><a href="/topic/thu%E1%BA%ADt%20to%C3%A1n%20%C4%91i%E1%BB%81u%20khi%E1%BB%83n" class="topic-relate-item" data-v-cafed715>thuật toán điều khiểncontrol algorithm</a><a href="/topic/b%C3%A0i%20to%C3%A1n%20b%E1%BA%ADc%20hai" class="topic-relate-item" data-v-cafed715>bài toán bậc haiquadratic programming problem</a><a href="/topic/kho%E1%BA%A3ng%20c%C3%A1ch%20bregman" class="topic-relate-item" data-v-cafed715>khoảng cách bregmanBregman divergence</a><a href="/topic/proton" class="topic-relate-item" data-v-cafed715>protonproton</a><a href="/topic/nhi%E1%BB%87t%20%C4%91%E1%BB%99%20n%C3%B3ng%20ch%E1%BA%A3y" class="topic-relate-item" data-v-cafed715>nhiệt độ nóng chảynhiet do nong chay</a><a href="/topic/h%E1%BA%B1ng%20s%E1%BB%91%20t%E1%BB%91c%20%C4%91%E1%BB%99" class="topic-relate-item" data-v-cafed715>hằng số tốc độrate constant</a><a href="/topic/%C4%91%E1%BB%99ng%20v%E1%BA%ADt%20linh%20tr%C6%B0%E1%BB%9Fng" class="topic-relate-item" data-v-cafed715>động vật linh trưởngprimates</a><a href="/topic/%C4%91%E1%BB%93ng%20v%E1%BB%8B" class="topic-relate-item" data-v-cafed715>đồng vịisotope</a><a href="/topic/chu%20tr%C3%ACnh%20sinh%20%C4%91%E1%BB%8Ba%20h%C3%B3a" class="topic-relate-item" data-v-cafed715>chu trình sinh địa hóaBiogeochemical cycle</a><a href="/topic/browse" class="topic-relate-more" data-v-cafed715> Khám phá từ điển học thuật </a></div></div></section></aside></div></div></div></div></div></main><footer class="site-footer" data-v-f975c723 data-v-4e17e457><svg class="footer-grid" aria-hidden="true" data-v-4e17e457><defs data-v-4e17e457><pattern id="footer-dot-grid" width="28" height="28" patternUnits="userSpaceOnUse" data-v-4e17e457><circle cx="1.5" cy="1.5" r="1.5" fill="#ffffff" opacity="0.08" data-v-4e17e457></circle></pattern></defs><rect width="100%" height="100%" fill="url(#footer-dot-grid)" data-v-4e17e457></rect></svg><div class="footer-glow" aria-hidden="true" data-v-4e17e457></div><div class="el-row is-justify-center footer-row" style="" type="flex" data-v-4e17e457><div class="el-col el-col-24 el-col-xs-24 el-col-sm-24 el-col-md-22 el-col-lg-20" style="" data-v-4e17e457><div class="footer-grid-layout" data-v-4e17e457><div class="footer-brand" data-v-4e17e457><a href="/" class="footer-logo" aria-label="Scholar Hub - Trang chủ" data-v-4e17e457><img src="/logo.svg" alt="Logo ScholarHub" width="40" height="40" data-v-4e17e457>Scholar Hub</a>Công cụ hỗ trợ trích dẫn và phân tích khoa học Việt Nam Scholar Hub hỗ trợ trích dẫn tài liệu tham khảo, xác minh tính chính xác của tài liệu tham khảo, tính toán chỉ số ảnh hưởng của các tạp chí tham gia Scholar Hub, cùng từ điển học thuật và phân tích xu hướng nghiên cứu. ScholarHub không đăng thông tin tổng hợp, không đăng lại nội dung từ các trang báo chí Việt Nam hoặc trang thông tin điện tử khác tại Việt Nam. </div><div class="footer-column" data-v-4e17e457><h3 class="footer-heading" data-v-4e17e457>Khám phá</h3><ul class="footer-links" data-v-4e17e457><li data-v-4e17e457><a href="/publication/browse" class="" data-v-4e17e457>Công bố khoa học</a></li><li data-v-4e17e457><a href="/publisher/browse" class="" data-v-4e17e457>Tạp chí khoa học</a></li><li data-v-4e17e457><a href="/topic/browse" class="" data-v-4e17e457>Từ điển học thuật</a></li><li data-v-4e17e457><a href="/user/citation" class="" data-v-4e17e457>Công cụ trích dẫn</a></li><li data-v-4e17e457><a href="/post/browse" class="" data-v-4e17e457>Bài viết, nội dung cập nhật</a></li></ul><h3 class="footer-heading footer-heading--spaced" data-v-4e17e457>Tham gia</h3><ul class="footer-links" data-v-4e17e457><li data-v-4e17e457><a href="https://docs.google.com/forms/d/e/1FAIpQLScPJI8r9B2cFAOO-00E0cAsU1cVnZ1_pk3YaXn5c21vgNrdfQ/viewform?usp=sf_link" target="_blank" rel="noopener" data-v-4e17e457> Đăng ký Tạp chí tham gia Scholar Hub </a></li><li data-v-4e17e457><a href="https://forms.gle/bWrbmfR549VNF1UP7" target="_blank" rel="noopener" data-v-4e17e457>Phản hồi ý kiến về Scholar Hub</a></li></ul></div><div class="footer-column" data-v-4e17e457><h3 class="footer-heading" data-v-4e17e457>Website liên kết</h3><ul class="footer-links" data-v-4e17e457><li data-v-4e17e457><a href="https://metis.vn" target="_blank" data-v-4e17e457>Hệ thống CSDL Khoa học & Công nghệ SciBase</a></li><li data-v-4e17e457><a href="https://kiemtratailieu.vn" target="_blank" data-v-4e17e457>Phần mềm kiểm tra trùng lặp Kiểm Tra Tài Liệu</a></li><li data-v-4e17e457><a href="https://vojs.vn" target="_blank" data-v-4e17e457>Phần mềm xuất bản tạp chí điện tử VOJS</a></li><li data-v-4e17e457><a href="https://conference.io.vn" target="_blank" data-v-4e17e457>Hệ thống hội thảo khoa học Việt Nam</a></li><li data-v-4e17e457><a href="https://letqa.vn" target="_blank" data-v-4e17e457>Nền tảng trắc nghiệm và đề thi đa lĩnh vực LetQA</a></li><li data-v-4e17e457><a href="https://unidoc.vn" target="_blank" data-v-4e17e457>Công cụ soạn thảo văn bản chuẩn thể thức UniDoc</a></li></ul></div><div class="footer-column footer-column--contact" data-v-4e17e457><h3 class="footer-heading" data-v-4e17e457>Liên hệ & hỗ trợ</h3><ul class="footer-contact" data-v-4e17e457><li data-v-4e17e457>Đơn vị chủ quản, phát triển và vận hành: Công ty Cổ phần Metis</li><li data-v-4e17e457>Số 8 B14, Khu đô thị Mỹ Đình, Phường Từ Liêm, Thành phố Hà Nội</li><li data-v-4e17e457>Hotline: <a href="tel:0566685688" data-v-4e17e457>0566.685.688</a></li><li data-v-4e17e457>Email: <a href="mailto:hotro@metis.vn" data-v-4e17e457>hotro@metis.vn</a></li><li data-v-4e17e457>Số giấy chứng nhận ĐKKD: 0109293202 cấp ngày 03/08/2020 tại Sở Kế hoạch và Đầu tư thành phố Hà Nội</li><li data-v-4e17e457>Người quản lý và chịu trách nhiệm nội dung: Nguyễn Ngọc Sơn</li><li data-v-4e17e457>Facebook: <a href="https://www.facebook.com/scholarhub.vn" target="_blank" rel="noopener" data-v-4e17e457>facebook.com/scholarhub.vn</a></li></ul></div></div><div class="footer-bottom" data-v-4e17e457>Copyright © 2023-2026 ScholarHub - Metis., Jsc. All rights reserved</div></div></div></footer></div></div><div id="teleports"></div><script>window.__NUXT__={};window.__NUXT__.config={public:{baseAPIUrl:"/api",baseInternalAPIUrl:"http://scholarhub-public-backend/api",baseFrontendUrl:"https://scholarhub.vn",googleAnalyticId:"G-TP4MKRFFZQ"},app:{baseURL:"/",buildId:"515db0c1-fb1e-4670-ad97-72850217475b",buildAssetsDir:"/_nuxt/",cdnURL:""}}$