对角占优矩阵是指一矩陣的每一橫行,對角線上元素的大小大於或等於同一橫行其他元素大小的和,一矩陣A為对角占优矩阵若
:|a_{ii}| \geq \sum_{j\neq i} |a_{ij}| \quad\text{for all } i, \,
其中aij為第i行第j列的元素。
上述的定義中用到大於等於,其條件較鬆,因此有時會稱為弱对角占优矩阵,若上述的定義用大於代替大於等於,則稱為強对角占优矩阵。对角优势矩阵可以指弱对角占优矩阵,也可以指強对角占优矩阵,視上下文而定。
變體
第一段的定義是考慮同一橫行其他元素大小的和,有時也稱為行对角优势矩阵,若是考慮同一直列其他元素大小的和,則稱為列对角优势矩阵。
若一矩阵是弱对角优势矩陣,但至少一橫行(或一直列)符合強对角优势的條件,則此矩陣稱為不可约对角优势矩阵。
例子
矩陣
:\mathbf A = \begin{bmatrix}
3 & -2 & 1\\
1 & -3 & 2\\
-1 & 2 & 4\end{bmatrix}
可得
:|a_{11}| \ge |a_{12}| + |a_{13}| 因為 |3| \ge |-2| + |1|
:|a_{22}| \ge |a_{21}| + |a_{23}| 因為 |-3| \ge |1| + |2|
:|a_{33}| \ge |a_{31}| + |a_{32}| 因為 |4| \ge |-1| + |2|.
因為任一對角線元素大小都大於等於同一行其他元素的和,因此A為对角优势矩阵。
矩陣
:\mathbf B = \begin{bmatrix}
-2 & 2 & 1\\
1 & 3 & 2\\
1 & -2 & 0\end{bmatrix}
但是
:|b_{11}| 因為 |-2|
:|b_{22}| \ge |b_{21}| + |b_{23}| 因為 |3| \ge |1| + |2|
:|b_{33}| 因為 |0| .
因為|b_{11}|和|b_{33}|都小於同一列其他元素大小的和,因此B不是对角优势矩阵。
矩陣
:\mathbf C = \begin{bmatrix}
-4 & 2 & 1\\
1 & 6 & 2\\
1 & -2 & 5\end{bmatrix}
可得
:|c_{11}| \ge |c_{12}| + |c_{13}| 因為 |-4| > |2| + |1|
:|c_{22}| \ge |c_{21}| + |c_{23}| 因為 |6| > |1| + |2|
:|c_{33}| \ge |c_{31}| + |c_{32}| 因為 |5| > |1| + |-2|.
因為任一對角線元素大小都大於同一行其他元素的和,因此C為強对角优势矩阵。
應用及性質
強对角优势矩阵(或不可约对角优势矩阵)是非奇異方陣,此結果即為Levy–Desplanques定理,針對強对角优势矩阵的結果,可以用證明。
若埃爾米特对角优势矩阵 A ,其對角線為非負值,即為正定矩陣。
若不考慮對稱性的條件,上述的矩陣不一定會是半正定矩陣。(例如, \begin{bmatrix}-\sqrt 5&2&1\end{bmatrix}\begin{bmatrix}
1&1&0\\
1&1&0\\
1&0&1\end{bmatrix}\begin{bmatrix}-\sqrt 5\\2\\1\end{bmatrix}=10-5\sqrt 5),但其特徵值的實部為非負數(參見对角优势矩阵的結果,可以用。)
類似的,若埃爾米特強对角优势矩阵的對角線元素為正,此矩陣為正定矩陣,此矩陣等於某個對角線元素為非負值實數的埃爾米特強对角优势矩阵A加上xI,其中x為正的實數(也是正定矩陣)。
若高斯消去法(LU分解)的矩陣為強对角优势矩阵,不需要進行尋找主元的過程。
若一線性聯立方程的矩陣為強对角优势矩阵或不可约对角优势矩阵,利用雅可比法及高斯-賽德爾迭代的計算結果會收斂。
許多從有限元素法中產生的矩陣都是对角优势矩阵。
參考資料
- Gene H. Golub & Charles F. Van Loan. Matrix Computations, 1996. ISBN 0-8018-5414-8
- Roger A. Horn & Charles R. Johnson. Matrix Analysis, Cambridge University Press, 1985. ISBN 0-521-38632-2 (paperback).
外部連結
- [http://planetmath.org/?op=getobj&from=objects&id=4512 PlanetMath: Diagonal dominance definition]
- [http://planetmath.org/?op=getobj&from=objects&id=7483 PlanetMath: Properties of diagonally dominant matrices]
- [https://mathworld.wolfram.com/DiagonallyDominantMatrix.html Mathworld]
评论 (0)