高斯消除矩阵

矩阵运算

初始矩阵

( a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 ) = ( b 1 b 2 b 3 ) \begin{pmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \\ \end{pmatrix}= \begin{pmatrix} b_{1}\\ b_{2}\\ b_{3}\\ \end{pmatrix} a11a21a31a12a22a32a13a23a33=b1b2b3

运算流程

规定从第一个元素 a11 开始, 如果a11 !=0 ,那么第一行除以a11然后减去相应的倍数从每个其他行,第一行归零法第一个条目的所有行(条件a11>0), 其他项自身 减去矩阵第一行的值*当前行的第一个值

  • 操作步骤
    • a11 变成1
    • a21 变成0
    • a31 变成0
    • 第一行每一个除a11
    • 第二行 自身 - 第一个数*第一行值
    • 第三行 自身 - 第一个数*第一行值

( 1 a 12 a 11 a 13 a 11 0 a 22 − a 21 ∗ a 12 a 11 a 23 − a 21 ∗ a 13 a 11 0 a 32 − a 31 ∗ a 12 a 11 a 33 − a 31 ∗ a 13 a 11 ) = ( b 1 a 11 b 2 − a 21 ∗ b 1 a 11 b 3 − a 31 ∗ b 1 a 11 ) \begin{pmatrix} 1 & \cfrac{a_{12}}{a_{11}} & \cfrac{a_{13}}{a_{11}} \\ 0 & a_{22}-a_{21}* \cfrac{a_{12}}{a_{11}} & a_{23} -a_{21}* \cfrac{a_{13}}{a_{11}} \\ 0 & a_{32}-a_{31}* \cfrac{a_{12}}{a_{11}} & a_{33} -a_{31}* \cfrac{a_{13}}{a_{11}}\\ \end{pmatrix} = \begin{pmatrix} \cfrac{b_{1}}{a_{11}} \\ b_{2}-a_{21}* \cfrac{b_{1}}{a_{11}} \\ b_{3}- a_{31}*\cfrac{b_{1}}{a_{11}}\\ \end{pmatrix} 100a11a12a22a21a11a12a32a31a11a12a11a13a23a21a11a13a33a31a11a13=a11b1b2a21a11b1b3a31a11b1

实践

{ x + 2 y + z = 2 3 x + 8 y + z = 12 4 y + z = 2 \begin{cases}x+2y+z=2\\3x+8y+z=12\\4y+z=2\\\end{cases} x+2y+z=23x+8y+z=124y+z=2 ( 1 2 1 3 8 1 0 4 1 ) = ( 2 12 2 ) \begin{pmatrix} 1& 2 & 1 \\ 3& 8 & 1 \\ 0 & 4 & 1 \\ \end{pmatrix} = \begin{pmatrix} 2\\ 12\\ 2\\ \end{pmatrix} 130284111=2122

第一步

  • 计算最外圈
    ( 1 2 1 1 1 ( 3 ) 8 − 3 ∗ 2 1 1 − 3 ∗ 1 1 ( 0 ) 4 − 0 ∗ 2 1 1 − 0 ∗ 1 1 ) = ( 2 1 12 − 3 ∗ 2 1 2 − 0 ∗ 2 1 ) \begin{pmatrix} 1 & \cfrac{2}{1} & \cfrac{1}{1} \\ (3) & 8-3* \cfrac{2}{1} & 1 -3* \cfrac{1}{1} \\ (0) & 4-0* \cfrac{2}{1} & 1 -0* \cfrac{1}{1}\\ \end{pmatrix} = \begin{pmatrix} \cfrac{2}{1}\\ 12 - 3*\cfrac{2}{1}\\ 2-0*\cfrac{2}{1}\\ \end{pmatrix} 1(3)(0)12831240121113111011=12123122012
  • 计算结果
    ( 1 2 1 0 2 − 2 0 4 1 ) = ( 2 6 2 ) \begin{pmatrix} 1 & 2 & 1 \\ 0 & 2 & -2 \\ 0 & 4 & 1 \\ \end{pmatrix}= \begin{pmatrix} 2\\ 6 \\ 2 \\ \end{pmatrix} 100224121=262

第二步

  • 计算下一个环
    ( 2 − 2 4 1 ) = ( 6 2 ) \begin{pmatrix} 2& -2\\ 4&1 \end{pmatrix}= \begin{pmatrix} 6 \\ 2 \\ \end{pmatrix} (2421)=(62) ( ( 2 ) 2 − 2 ( 4 ) 1 − 4 ∗ 2 − 2 ) = ( 6 2 2 − 4 ∗ 6 2 ) \begin{pmatrix} (2) & \cfrac{2}{-2} \\ (4) & 1 -4* \cfrac{2}{-2} \\ \end{pmatrix}= \begin{pmatrix} \cfrac{6}{2} \\ 2 -4*\cfrac{6}{2} \\ \end{pmatrix} (2)(4)221422=262426
  • 计算结果 ( 0 − 1 0 5 ) = ( 3 − 10 ) \begin{pmatrix} 0& -1 \\ 0 & 5 \\ \end{pmatrix}= \begin{pmatrix} 3\\ -10 \end{pmatrix} (0015)=(310)

第三步

( 1 2 1 0 2 − 2 0 0 5 ) = ( 2 6 − 10 ) \begin{pmatrix} 1& 2 & 1 \\ 0& 2 & -2 \\ 0& 0 & 5 \\ \end{pmatrix}= \begin{pmatrix} 2 \\ 6\\ -10\\ \end{pmatrix} 100220125=2610 { x + 2 y + z = 2 2 y − 2 z = 6 5 z = − 10 \begin{cases}x+2y+z=2\\ 2y-2z=6\\5z=-10\\\end{cases} x+2y+z=22y2z=65z=10 { x = 2 y = 1 z = − 2 \begin{cases}x=2\\y=1\\z=-2\\\end{cases} x=2y=1z=2

你可能感兴趣的:(map)