主題
Search

高斯消元法

下載 Mathematica 筆記本下載 Wolfram 筆記本

高斯消元法是求解形如 矩陣方程 的方法

Ax=b.
(1)

要執行高斯消元法，從方程組開始

[a_(11) a_(12) ... a_(1k); a_(21) a_(22) ... a_(2k); | | ... |; a_(k1) a_(k2) ... a_(kk)][x_1; x_2; |; x_k]=[b_1; b_2; |; b_k],
(2)

組成“增廣矩陣方程”

[a_(11) a_(12) ... a_(1k); a_(21) a_(22) ... a_(2k); | | ... |; a_(k1) a_(k2) ... a_(kk)|b_1; b_2; |; b_k][x_1; x_2; |; x_k].
(3)

這裡，變數 x 中的列向量被保留用於標記矩陣行。現在，執行初等行變換以將增廣矩陣變為上三角形式

[a_(11)^' a_(12)^' ... a_(1k)^'; 0 a_(22)^' ... a_(2k)^'; | | ... |; 0 0 ... a_(kk)^'|b_1^'; b_2^'; |; b_k^'].
(4)

求解第 k 行的方程得到 x_k，然後代回第 (k-1) 行的方程以獲得 x_(k-1) 的解，等等，根據公式

x_i=1/(a_(ii)^')(b_i^'-sum_(j=i+1)^ka_(ij)^'x_j).
(5)

在 Wolfram 語言中，RowReduce執行高斯消元法的一個版本，方程 mx=b 透過以下方式求解：

  GaussianElimination[m_?MatrixQ, v_?VectorQ] :=
    Last /@ RowReduce[Flatten /@ Transpose[{m, v}]]

矩陣的 LU 分解 經常用作高斯消元過程的一部分，用於求解矩陣方程。

經過高斯消元法的矩陣被稱為處於階梯形。

例如，考慮矩陣方程

[9 3 4; 4 3 4; 1 1 1][x_1; x_2; x_3]=[7; 8; 3].
(6)

以增廣形式，這變為

[9 3 4; 4 3 4; 1 1 1|7; 8; 3][x_1; x_2; x_3].
(7)

交換第一行和第三行（不交換右手列向量中的元素）得到

[1 1 1; 4 3 4; 9 3 4|3; 8; 7][x_1; x_2; x_3].
(8)

從第三行減去第一行的 9 倍得到

[1 1 1; 4 3 4; 0 -6 -5| 3; 8; -20][x_1; x_2; x_3].
(9)

從第二行減去第一行的 4 倍得到

[1 1 1; 0 -1 0; 0 -6 -5| 3; -4; -20][x_1; x_2; x_3].
(10)

最後，將第二行的 -6 倍加到第三行得到

[1 1 1; 0 -1 0; 0 0 -5| 3; -4; 4][x_1; x_2; x_3].
(11)

恢復變換後的矩陣方程得到

[1 1 1; 0 -1 0; 0 0 -5][x_1; x_2; x_3]=[ 3; -4; 4],
(12)

可以立即求解得到 x_3=-4/5，反向代入得到 x_2=4 （在本例中實際上是顯然的），然後再次反向代入找到 x_1=-1/5.

參見

增廣矩陣, 凝結法, 初等行和列運算, 階梯形, 高斯-約旦消元法, LU 分解, 矩陣方程, 平方根法

使用 探索

參考文獻

Bareiss, E. H. "Multistep Integer-Preserving Gaussian Elimination." Argonne National Laboratory Report ANL-7213, May 1966.Bareiss, E. H. "Sylvester's Identity and Multistep Integer-Preserving Gaussian Elimination." Math. Comput. 22, 565-578, 1968.Garbow, B. S. "Integer-Preserving Gaussian Elimination." Program P-158 (3600F), Applied Mathematics Division, Argonne National Laboratory, Nov. 21, 1966.Gentle, J. E. "Gaussian Elimination." §3.1 in Numerical Linear Algebra for Applications in Statistics. Berlin: Springer-Verlag, pp. 87-91, 1998.

在 中引用

高斯消元法

引用為

Weisstein, Eric W. “高斯消元法。” 來自 Web 資源。 https://mathworld.tw/GaussianElimination.html

主題分類