矩阵消元
三行三元的消元步骤1是这样的:
用矩阵来表示就是找到 E21 E 21 矩阵,使得第2行第1列的元变为0:
Ax=b A x = b 变为了 E21Ax=E21b E 21 A x = E 21 b :
矩阵乘法
从上面的例子可以知道
E(Ax)=(EA)x
E
(
A
x
)
=
(
E
A
)
x
:
结合律
A(BC)=(AB)C A ( B C ) = ( A B ) C
交换律是不成立的:
AB≠BA A B ≠ B A
矩阵乘法
矩阵A乘以矩阵B的第j列,作为结果矩阵的第j列
AB=A[b1b2b3]=[Ab1Ab2Ab3] A B = A [ b 1 b 2 b 3 ] = [ A b 1 A b 2 A b 3 ]
使用置换矩阵(Permutation matrix) Pij P i j 交换行
P23 P 23 表示置换矩阵的一种,行交换矩阵(Row Exchange Matrix),可以用于交换第2和第3行
Row Exchange Matrix:
Pij P i j 是通过将单位矩阵的i行和j行交换得到,用该矩阵乘以其他矩阵,会交换其他矩阵的i行和j行
增广矩阵(Augmented Maxtrix)
可以向矩阵A增加一列b:
消元的过程变为:
其中第二行0,1,1,4表示 x2+x3=4 x 2 + x 3 = 4
矩阵乘法
E[Ab]
E
[
A
b
]
的两种计算方式:
Rows:E的每行acts on
[Ab]
[
A
b
]
,作为
[EAEb]
[
E
A
E
b
]
的行
Columns:E acts on
[Ab]
[
A
b
]
的每列,作为
[EAEb]
[
E
A
E
b
]
的列
每一步的Elimination操作都是影响一行的改变, E32E31E21A E 32 E 31 E 21 A 最终得到了三角矩阵