当前位置:   article > 正文

线性代数笔记2--矩阵消元

线性代数笔记2--矩阵消元

0. 简介

矩阵消元

1. 消元过程

实例方程组
{ x + 2 y + z = 2 3 x + 8 y + z = 12 4 y + z = 2

{x+2y+z=23x+8y+z=124y+z=2
x+2y+z=23x+8y+z=124y+z=2
矩阵化
A = [ 1 2 1 3 8 1 0 4 1 ] X = [ x y z ] A=
[121381041]
\\ X=
[xyz]
A= 130284111 X= xyz

B = [ 2 12 2 ] B=
[2122]
B= 2122

消元

[ 1 2 1 3 8 1 0 4 1 ] ⟶ ( 2 , 1 ) [ 1 2 1 0 2 − 2 0 4 1 ] ⟶ ( 3 , 2 ) [ 1 2 1 0 2 − 2 0 0 5 ]

[121381041]
\stackrel{(2,1)}{\longrightarrow}
[121022041]
\stackrel{(3,2)}\longrightarrow
[121022005]
130284111 (2,1) 100224121 (3,2) 100220125
回代
[ 2 12 2 ] ⟶ r o w 2 − 3 r o w 1 [ 2 6 2 ] ⟶ r o w 3 − 2 r o w 2 [ 2 6 − 10 ]
[2122]
\stackrel{row_2-3row_1}\longrightarrow
[262]
\stackrel{row_3-2row_2}\longrightarrow
[2610]
2122 row23row1 262 row32row2 2610

求解
[ 1 2 1 0 2 − 2 0 0 5 ] [ x y z ] = [ 2 6 − 10 ]
[121022005]
[xyz]
=
[2610]
100220125 xyz = 2610

结果
[ x y z ] = [ 2 1 − 2 ]
[xyz]
=
[212]
xyz = 212

2. 消元矩阵

将上述消元的过程变为矩阵相乘的形式。

向量式思考
矩阵乘列向量
[ . . . . . . . . . ] [   x y z ] = [ x ∗ c o l 1 + y ∗ c o l 2 + z ∗ c o l 3 ]

[.........]
[ xyz]
=
[xcol1+ycol2+zcol3]
.........  xyz =[xcol1+ycol2+zcol3]
行向量乘矩阵
[   x   y   z ] [ . . . . . . . . . ] = [ x ∗ r o w 1 + y ∗ r o w 2 + z ∗ r o w 3 ]
[ x y z]
[.........]
=
[xrow1+yrow2+zrow3]
[ x y z] ......... = xrow1+yrow2+zrow3

一个矩阵左边乘一个单位矩阵并不改变其值
A = [ 1 2 1 3 8 1 0 4 1 ] = [ 1 0 0 0 1 0 0 0 1 ] [ 1 2 1 3 8 1 0 4 1 ] A=
[121381041]
=
[100010001]
[121381041]
A= 130284111 = 100010001 130284111

而做行的加减则可以
A = [ 1 2 1 3 8 1 0 4 1 ] A ′ = [ 1 0 0 − 3 1 0 0 0 1 ] [ 1 2 1 3 8 1 0 4 1 ] = [ 1 2 1 0 2 − 2 0 4 1 ] A=
[121381041]
\\ A'=
[100310001]
[121381041]
=\\
[121022041]
A= 130284111 A= 130010001 130284111 = 100224121

实际上这个过程就是,我们在之前的消元过程中的第二行减去三倍第一行的过程。我们继续下去将这个矩阵对角化。
A ′ ′ = [ 1 0 0 0 1 0 0 − 2 1 ] A ′ = [ 1 0 0 0 1 0 0 − 2 1 ] [ 1 2 1 0 2 − 2 0 4 1 ] = [ 1 2 1 0 2 − 2 0 0 5 ] A''=
[100010021]
A'\\=
[100010021]
[121022041]
\\=
[121022005]
A′′= 100012001 A= 100012001 100224121 = 100220125

我们令最后的上三角矩阵为
U = [ 1 2 1 0 2 − 2 0 0 5 ] U=

[121022005]
U= 100220125
两个变换矩阵为
E 21 = [ 1 0 0 − 3 1 0 0 0 1 ] E 32 = [ 1 0 0 0 1 0 0 − 2 1 ] E_{21}=
[100310001]
\\ E_{32}=
[100010021]
E21= 130010001 E32= 100012001

E 32 ( E 21 A ) = U E_{32}(E_{21}A)=U E32(E21A)=U
而矩阵乘法满足结合律证明
E 32 E 21 A = E 32 ( E 21 A ) E_{32}E_{21}A=E_{32}(E_{21}A) E32E21A=E32(E21A)
所以最终消元的过程变成了寻找矩阵E的过程
E = E 32 E 21 E=E_{32}E_{21} E=E32E21
这一过程。

3. 置换矩阵

在上述的消元矩阵中,我们并没有进行列的交换。那么如何进行交换呢?

我们知道在原矩阵基础左边乘单位矩阵,矩阵不会发生变化。
A = [ 1 2 3 4 ] =   [ 1 0 0 1 ] [ 1 2 3 4 ] A=

[1234]
=\
[1001]
[1234]
A=[1324]= [1001][1324]

如何交换两行呢,将单位矩阵变形
A ′ = [ 0 1 1 0 ] [ 1 2 3 4 ] = [ 3 4 1 2 ] A'=

[0110]
[1234]
=
[3412]
A=[0110][1324]=[3142]
推广到多行
A = [ 1 2 3 4 5 6 7 8 9 ] A=
[123456789]
A= 147258369

  • 行变换
    交换第一行和第三行
    A ′ = [ 0 0 1 0 1 0 1 0 0 ] [ 1 2 3 4 5 6 7 8 9 ] A'=
    [001010100]
    [123456789]
    A= 001010100 147258369

    交换第一行和第二行
    A ′ ′ = [ 0 1 0 1 0 0 0 0 1 ] [ 1 2 3 4 5 6 7 8 9 ] A''=
    [010100001]
    [123456789]
    A′′= 010100001 147258369

所以交换任意两行,只需将单位矩阵中对应行 1 1 1的位置进行交换。

  • 列变换

在矩阵左边乘是对原矩阵行变换,而在矩阵右边则是列变换
交换矩阵两列
A = [ 1 2 3 4 ] A ′ = [ 1 2 3 4 ] [ 0 1 1 0 ] = [ 2 1 4 3 ] A=

[1234]
\\ A'=
[1234]
[0110]
=
[2143]
A=[1324]A=[1324][0110]=[2413]

交换多列也是一样的效果
交换第 1 1 1 2 2 2
A = [ 1 2 3 4 5 6 7 8 9 ] A ′ = [ 1 2 3 4 5 6 7 8 9 ] [ 0 1 0 1 0 0 0 0 1 ] = [ 2 1 3 5 4 6 8 7 9 ] A=

[123456789]
\\ A'=
[123456789]
[010100001]
=
[213546879]
A= 147258369 A= 147258369 010100001 = 258147369

所以交换任意两列,只需将单位矩阵中对应行 1 1 1的位置进行交换。
与行交换的不同地方在于,矩阵乘的在右边了。

4. 矩阵的逆

A = [ 1 0 0 − 3 1 0 0 0 1 ] A − 1 = [ 1 0 0 3 1 0 0 0 1 ] A − 1 A = [ 1 0 0 3 1 0 0 0 1 ] [ 1 0 0 − 3 1 0 0 0 1 ] = [ 1 0 0 0 1 0 0 0 1 ] A=

[100310001]
\\ A^{-1}=
[100310001]
\\ A^{-1}A=
[100310001]
[100310001]
=
[100010001]
A= 130010001 A1= 130010001 A1A= 130010001 130010001 = 100010001

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/AllinToyou/article/detail/89330
推荐阅读
相关标签
  

闽ICP备14008679号