常见的特殊矩阵及分解(一)-程序员宅基地

技术标签: 算法推导  矩阵  线性代数  

参考:矩阵分析、维基百科、线代启示录、百度百科

在SLAM的学习过程中,经常会遇到矩阵的求解问题,这就需要用到矩阵分解,而使用一种分解方法很多时候都需要矩阵满足某种特性,故打算分三次博客分别介绍下这些特殊矩阵、分解方法及在Eigen库中的使用.
下面简单给出常见特殊矩阵的性质(必要条件)与判别方法(充分条件).


1、相似矩阵

两个 n × n n×n n×n 矩阵 A A A B B B 为相似矩阵,当且仅当存在一个 n × n n×n n×n 的可逆矩阵 P P P,使得:
P − 1 A P = B P^{-1}AP=B P1AP=B
P P P 被称为矩阵 A A A B B B 之间的相似变换矩阵.

1.1 特性

相似变换下的不变性质:
1)两者的相等
2)两者的行列式值相等
3)两者的迹数相等
4)两者拥有同样的特征值,尽管相应的特征向量一般不同

1.2 相似对角化

A ∈ C n × n A \in C^{n \times n} ACn×n,则 A A A 可以对角化的充要条件是 A A A n n n 个线性无关的特征向量. 且 Λ \Lambda Λ 的对角线元素都是 A A A 的特征值(包括重数).

P − 1 A P = Λ P^{-1}AP=\Lambda P1AP=Λ

证明:
设有可逆矩阵 P P P 使得
P − 1 A P = Λ , A P = P Λ P^{-1}AP=\Lambda , AP=P\Lambda P1AP=Λ,AP=PΛ

P = ( P 1 , P 2 , . . . , P n ) , Λ = d i a g ( λ 1 , λ 2 , . . . , λ n ) P=(P_1, P_2,..., P_n) , \Lambda=diag(\lambda_1, \lambda_2, ..., \lambda_n) P=(P1,P2,...,Pn),Λ=diag(λ1,λ2,...,λn)

( A P 1 , A P 2 , . . . , A P n ) = ( λ 1 P 1 , λ 2 P 2 , . . . , λ n P n ) (AP_1, AP_2,..., AP_n) =(\lambda_1P_1, \lambda_2P_2, ..., \lambda_nP_n) (AP1,AP2,...,APn)=(λ1P1,λ2P2,...,λnPn)
可见 P i P_i Pi A A A 的特征向量,而 P P P 又是可逆的,所以 A A A n n n 个线性无关的特征向量

对角化条件
一个 n × n n×n n×n的矩阵 M M M 是可对角化的,当且仅当 M M M 满足下列条件之一:

  • M M M实对称矩阵(参见对角分解
  • M M M n n n 个线性无关的特征向量。或者说, M M M 有一个由特征向量组成的基。(称作极大无关条件)
  • M M M 的所有特征值的几何重数(即相应特征子空间的维数)等于相应的代数重数(即特征多项式中 ( x − λ ) (x-\lambda ) (xλ)项的次数)。或者说, M M M 的所有几何重数之和等于 n n n。(称作重数相等条件)
  • M M M 的极小多项式经标准分解后,每一项都是一次项,且重数都是1。(称作互异单根条件)
1.3 约当标准型

有上可知,并不是每个矩阵都可以进行相似对角化,如果不能对角化,那么通过相似变换,矩阵能够化成最简单的形式就是约当标准型.
P − 1 A P = J P^{-1}AP=J P1AP=J
约当标准型是由若干个约当块构成的分块对角矩阵 J i J_i Ji 组成
J = [ J 1 ⋱ J p ] J = \left[ \begin{array} { l l } { J _ { 1 } } & { } \\ { } & { \ddots } & { } \\ { } & { } & { J _ { p } } \end{array} \right] J= J1Jp
每一个小块矩阵都具备一种很简单的形状:
J i = [ λ i 1 λ i ⋱ ⋱ 1 λ i ] J _ { i } = \left[ \begin{array} { c c c c } { \lambda _ { i } } & { 1 } & { } & { } \\ { } & { \lambda _ { i } } & { \ddots } & { } \\ { } & { } & { \ddots } & { 1 } \\ { } & { } & { } & { \lambda _ { i } } \end{array} \right] Ji= λi1λi1λi


2、正规矩阵

满足下列条件的矩阵都属于正规矩阵
A H A = A A H A^HA=AA^H AHA=AAH
注意,正规矩阵不一定可逆,对称矩阵就不一定可逆(见下).
正规矩阵包括但不限于下列矩阵:
在这里插入图片描述

如,正规矩阵 A = [ 1 1 − 1 1 ] A=\begin{bmatrix} 1 & 1 \\ -1 & 1 \\ \end{bmatrix} A=[1111] 就不属于上面几种特殊矩阵.

性质
1)正规矩阵一定可以对角化

2.1 实对称矩阵(symmetric matrix

A = A T A=A^T A=AT

特性
1)实对称矩阵都可以正交对角化
2)参见,2.4 厄米特矩阵
3)如果 X X X 是对称矩阵,那么 A X A T AXA^{T} AXAT 也是对称矩阵
4)一个矩阵同时为对称矩阵及斜对称矩阵当且仅当所有元素都是零
5)两个对称矩阵的积是对称矩阵,当且仅当两者的乘法可交换,两个实对称矩阵乘法可交换当且仅当两者的特征空间相同
6)对角矩阵都是对称矩阵
7)对于任何方形矩阵 X X X X + X T X+X^T X+XT 是对称矩阵

2.2 反实对称矩阵/ 斜对称矩阵(skew-symmetric matrix

A T = − A A^T=-A AT=A
行列式
det ⁡ ( A ) = det ⁡ ( A T ) = det ⁡ ( − A ) = ( − 1 ) n det ⁡ ( A ) \operatorname { det } ( A ) = \operatorname { det } \left( A ^ { T } \right) = \operatorname { det } ( - A ) = ( - 1 ) ^ { n } \operatorname { det } ( A ) det(A)=det(AT)=det(A)=(1)ndet(A)

  • n n n 是奇数,行列式等于0(雅可比定理),从这里看出,李代数 s o ( 3 ) \mathfrak { s o } ( 3 ) so(3)中的 [ ω ] × [ \omega ] _ { \times } [ω]× 的行列式就是0
  • n n n 是偶数,行列式可以写成部分元素的多项式的平方
  • 任意实斜对称矩阵的行列式非负数

特性
1)斜对称矩阵的主对角线元素必是0,所以其迹数为零
2)斜对称矩阵自身相乘的积是对称矩阵
3)任意矩阵 A A A A T − A A^T-A ATA是斜对称矩阵
4)若 A A A 是斜对称矩阵, x x x 是向量, x T A x x^T A x xTAx= 0,即二次型为0

2.3 正交矩阵(orthogonal matrix

正交矩阵是一个方块矩阵 Q Q Q,其元素为实数,而且行与列皆为正交的单位向量,使得该矩阵的转置矩阵为其逆矩阵:
Q T = Q − 1 ⇔ Q T Q = Q Q T = I Q ^ { T } = Q ^ { - 1 } \Leftrightarrow Q ^ { T } Q = Q Q ^ { T } = I QT=Q1QTQ=QQT=I
正交矩阵的行列式值必定为+1-1,因为:
1 = det ⁡ ( I ) = det ⁡ ( Q T Q ) = det ⁡ ( Q T ) det ⁡ ( Q ) = ( det ⁡ ( Q ) ) 2 ⇒ ∣ Q ∣ = ± 1 1 = \operatorname { det } ( I ) = \operatorname { det } \left( Q ^ { T } Q \right) = \operatorname { det } \left( Q ^ { T } \right) \operatorname { det } ( Q ) = ( \operatorname { det } ( Q ) ) ^ { 2 } \Rightarrow | Q | = \pm 1 1=det(I)=det(QTQ)=det(QT)det(Q)=(det(Q))2Q=±1

特性
1)作为一个线性映射(变换矩阵),正交矩阵保持距离不变(即,二范数的酉不变性),所以它是一个保距映射,具体例子为旋转与镜射
2)行列式值为+1的正交矩阵,称为特殊正交矩阵,即我们常说的旋转矩阵
3)行列式值为-1的正交矩阵,称为瑕旋转矩阵,瑕旋转就是旋转加上镜射。
4)所有 n × n n \times n n×n 正交矩阵形成一个群 O ( n ) O(n) O(n),称为正交群。亦即,正交矩阵与正交矩阵的乘积也是一个正交矩阵。
5)所有特殊正交矩阵形成一个子群 S O ( n ) SO(n) SO(n),称为特殊正交群。亦即,旋转矩阵与旋转矩阵的乘积也是一个旋转矩阵。

在矩阵分解中应用
一些重要的矩阵分解(Golub & Van Loan, 1996)涉及到了正交矩阵,包括:
1)QR分解: M = Q R M = QR M=QR, Q Q Q 正交, R R R 上三角。
2)奇异值分解: M = U Σ V T M = UΣV^T M=UΣVT, U U U V V V 正交, Σ Σ Σ 非负对角。
3)谱分解: S = Q Λ Q T S = QΛQ^T S=QΛQT, S S S 对称, Q Q Q 正交, Λ Λ Λ 对角。
4)极分解: M = Q S M = QS M=QS, Q Q Q 正交, S S S 对称半正定

2.4 埃尔米特矩阵/厄米特矩阵(Hermitian matrix/self-adjoint matrix

A = A H A=A^H A=AH
其中 A H A^H AH 表示 A A A 的共轭转置(conjugate transpose),埃尔米特矩阵中每一个第 i i i 行第 j j j 列的元素都与第 j j j 行第 i i i 列的元素的复共轭,埃尔米特矩阵主对角线上的元素都是实数的

特别注意,线性代数还有另一种伴随矩阵。共轭转置 A ∗ = A ‾ T A^{\ast}=\overline{A}^T A=AT也称为 A A A 的伴随(adjoint),所以厄米特矩阵又叫自伴随矩阵.

特性
1)埃尔米特矩阵特征值是实数
2)若 A A A 是Hermitian矩阵,对于任意向量 x ∈ C n \mathbf{x}\in\mathbb{C}^n xCn x ∗ A x \mathbf{x}^{\ast}A\mathbf{x} xAx 是实数
3)Hermitian矩阵对应相异特征值的特征向量互为正交
4)Hermitian矩阵的特征值的代数重数等于几何重数

2.5 反厄米特矩阵/反埃尔米特矩阵

同上

2.6 酉矩阵

3、正定矩阵

广义定义
M M M n n n 阶方阵,如果对任何非零向量z,都有 z T M z > 0 z^TMz> 0 zTMz>0,其中 z T z^T zT 表示 z z z 的转置,就称 M M M 为正定矩阵

狭义定义
正定矩阵是埃尔米特矩阵的一种. 维基百科使用的是这种定义方式,默认正定矩阵属于埃尔米特矩阵.
A A A 为一个 n × n n\times n n×n实对称矩阵。若每一 n n n 维非零实向量 x \mathbf{x} x 皆使得
x T A x > 0 \mathbf{x}^TA \mathbf{x}>0 xTAx>0
我们称 A A A 为正定(positive definite);若将上述条件放松为
x T A x ≥ 0 \mathbf{x}^TA \mathbf{x} \geq 0 xTAx0
A A A 称为半正定(positive semidefinite)。改变正定和半正定的不等式方向就有 A A A 是负定或半负定的概念,也可以说 A A A 是正定或半正定。如果 x T A x \mathbf{x}^TA\mathbf{x} xTAx 可能是正值也可能是负值,则称A是未定的(indefinite)
上述实正定矩阵的定义同样适用于复矩阵,但转置 ( ⋅ ) T (\cdot)^T ()T 要改成共轭转置 ( ⋅ ) ∗ (\cdot)^\ast (),因此实对称矩阵 A T = A A^T=A AT=A 替换为Hermitian矩阵 A ∗ = A A^\ast=A A=A.

任何一个实方阵 A A A 必可表示为一个实对称矩阵与一个反对称矩阵之和, A = B + C A = B + C A=B+C,其中,
B = 1 2 ( A + A T ) C = 1 2 ( A − A T ) \begin{aligned} B & = \frac { 1 } { 2 } \left( A + A ^ { T } \right) \\ C & = \frac { 1 } { 2 } \left( A - A ^ { T } \right) \end{aligned} BC=21(A+AT)=21(AAT)

x T A x = x T ( B + C ) x = x T B x + x T C x \mathbf { x } ^ { T } A \mathbf { x } = \mathbf { x } ^ { T } ( B + C ) \mathbf { x } = \mathbf { x } ^ { T } B \mathbf { x } + \mathbf { x } ^ { T } C \mathbf { x } xTAx=xT(B+C)x=xTBx+xTCx
由反对称矩阵性质可得,
x T C x = ( x T C x ) T = x T C T x = − x T C x \mathbf { x } ^ { T } C \mathbf { x } = \left( \mathbf { x } ^ { T } C \mathbf { x } \right) ^ { T } = \mathbf { x } ^ { T } C ^ { T } \mathbf { x } = - \mathbf { x } ^ { T } C \mathbf { x } xTCx=(xTCx)T=xTCTx=xTCx
即得, x T C x = 0 \mathbf { x } ^ { T } C \mathbf { x }=0 xTCx=0,进一步得,
x T A x = x T B X \mathbf { x } ^ { T } A \mathbf { x } = \mathbf { x } ^ { T } B \mathbf { X } xTAx=xTBX
即二次型 x T A x \mathbf{x}^TA\mathbf{x} xTAx 可用对称部分表示。若 A A A 不为对称矩阵,则 A A A 为正定矩阵(即 x T A x > 0 \mathbf{x}^TA\mathbf{x}>0 xTAx>0,对于所有 x ≠ 0 ) \mathbf{x}\neq\mathbf{0}) x=0) 等价于 1 2 ( A + A T ) \frac{1}{2}(A+A^T) 21(A+AT) 为正定矩阵

性质
1)正定矩阵的每一个主子阵都是正定的
2)正定矩阵的特征值皆为正数,所以,利用矩阵特征值、行列式与迹数(trace)性质可得,正定矩阵必可逆
λ \lambda λ 为正定矩阵 A A A 的一个特征值, x \mathbf{x} x 为对应的特征向量, x ≠ 0 \mathbf{x}\neq\mathbf{0} x=0,则
x ∗ A x = x ∗ λ x = λ x ∗ x \mathbf { x } ^ { * } A \mathbf { x } = \mathbf { x } ^ { * } \lambda \mathbf { x } = \lambda \mathbf { x } ^ { * } \mathbf { x } xAx=xλx=λxx
所以, λ = ( x ∗ A x ) / x ∗ x \lambda=(\mathbf{x}^{\ast}A\mathbf{x})/\mathbf{x}^{\ast}\mathbf{x} λ=(xAx)/xx,分子与分母都是正数,故 λ > 0 \lambda>0 λ>0
3)正定矩阵的轴(pivot)都是正数
@TODO
4)正定矩阵 A A A 可表示为 A = B ∗ B A=B^{\ast}B A=BB B B B 是一个可逆矩阵(这里的证明也使用的正定矩阵的狭义定义)

判别正定阵
1)求出A的所有特征值。若A的特征值均为正数,则A是正定的;若A的特征值均为负数,则A为负定的。

正定或半正定矩阵与奇异值分解有密切的关系,对于任意 m × n m\times n m×n阶矩阵 A A A,交互乘积 A T A A^TA ATA A A T AA^T AAT特征值都不为负值(即,SVD分解中的奇异值),故 A T A A^TA ATA A A T AA^T AAT半正定矩阵
如果 A T A A^TA ATA 可逆,那么 A T A A^TA ATA 一定是正定矩阵,因为 A T A A^TA ATA 的行列不为0,所以SVD分解中的 D D D 矩阵可逆!

2)计算A的各阶顺序主子式。若A的各阶顺序主子式均大于零,则A是正定的;若A的各阶顺序主子式中,奇数阶主子式为负,偶数阶为正,则A为负定的。

3)若 n × n n\times n n×n 阶Hermitian矩阵 A A A 的轴皆为正数,则 A A A 是正定矩阵
4)若 n × n n\times n n×n 阶Hermitian矩阵 A A A 可表示为 A = B ∗ B A=B^{\ast}B A=BB B B B 是一个可逆矩阵,则 A A A 是正定矩阵

判别半正定
1)所有的主子式非负,顺序主子式非负并不能推出矩阵是半正定的


4、伴随矩阵

A A A余子矩阵是一个 n × n n×n n×n 的矩阵 C C C,使得其第 i i i 行第 j j j 列的元素是 A A A 关于第 i i i 行第 j j j 列的代数余子式.
矩阵 A A A 的伴随矩阵是 A A A 的余子矩阵的转置矩阵. 即
adj ⁡ ( A ) = C T \operatorname { adj } ( \mathbf { A } ) = \mathbf { C } ^ { T } adj(A)=CT

如果矩阵可逆,那么它的逆矩阵和它的伴随矩阵之间只差一个系数,但是伴随矩阵对不可逆的矩阵也有定义的!

伴随矩阵的秩
当矩阵 A A A 可逆时,它的伴随矩阵也可逆,因此两者的秩一样,都是 n n n,并且
A − 1 = 1 det ⁡ A adj ⁡ A A ^ { - 1 } = \frac { 1 } { \operatorname { det } A } \operatorname { adj } A A1=detA1adjA

当矩阵 A A A 不可逆时, A A A 的伴随矩阵的秩通常并不与 A A A 相同,

  • A A A 的秩为 n − 1 n-1 n1 时,其伴随矩阵的秩为1
  • A A A 的秩小于 n − 1 n-1 n1 时,其伴随矩阵为零矩阵

性质
1) adj ⁡ ( I ) = I \operatorname {adj} ( \mathbf { I } ) = \mathbf { I } adj(I)=I
2) adj ⁡ ( A B ) = a d j ( B ) adj ⁡ ( A ) \operatorname { adj } ( \mathbf { A } \mathbf { B } ) = \mathbf { a } \mathrm { d } \mathbf { j } ( \mathbf { B } ) \operatorname { adj } ( \mathbf { A } ) adj(AB)=adj(B)adj(A),注意顺序
3) adj ⁡ ( A T ) = adj ⁡ ( A ) T \operatorname { adj } \left( \mathbf { A } ^ { T } \right) = \operatorname { adj } ( \mathbf { A } ) ^ { T } adj(AT)=adj(A)T
4)如果 A A A 可逆,那么 adj ⁡ ( A − 1 ) = adj ⁡ ( A ) − 1 = A det ⁡ A \operatorname { adj } \left( \mathbf { A } ^ { - 1 } \right) = \operatorname { adj } ( \mathbf { A } ) ^ { - 1 } = \frac { A } { \operatorname { det } A } adj(A1)=adj(A)1=detAA
5)如果 A A A对称矩阵,那么其伴随矩阵也是对称矩阵;如果 A A A反对称矩阵,那么当 n n n 为偶数时, A A A 的伴随矩阵也是反对称矩阵, n n n 为奇数时则是对称矩阵
6)如果 A A A(半)正定矩阵,那么其伴随矩阵也是(半)正定矩阵
7)如果矩阵 A A A B B B 相似,那么 a d j ( A ) \mathrm{adj}(\mathbf{A}) adj(A) a d j ( B ) \mathrm{adj}(\mathbf{B}) adj(B) 也相似
8)如果 n > 2 n>2 n>2,那么非零矩阵 A A A正交矩阵当且仅当 a d j ( A ) = ± A T \mathrm{adj}(\mathbf{A}) = \pm A^T adj(A)=±AT注意,可使用该方法判别正交矩阵


5、可交换矩阵

一般来说,矩阵的乘法是不可交换的,但当矩阵满足一定条件时,其乘法变成可交换,即
A ⋅ B = B ⋅ A A·B=B·A AB=BA
判定:
1)设A , B 至少有一个为零矩阵,则A , B 可交换;
2)设A , B 至少有一个为单位矩阵, 则A , B可交换;
3)设A , B 至少有一个为数量矩阵, 则A , B可交换;
4)设A , B 均为对角矩阵,则A , B 可交换;
5)设A , B 均为准对角矩阵(准对角矩阵是分块矩阵概念下的一种矩阵。即除去主对角线上分块矩阵不为零矩阵外,其余分块矩阵均为零矩阵),且对角线上的子块均可交换,则A , B 可交换
6)设 A ∗ A* A A A A 的伴随矩阵,则 A ∗ A* A A A A 可交换
7)设 A A A 可逆,则 A A A 与其逆矩阵可交换
注: A A A 的逆矩阵经过数乘变换所得到的矩阵也可以与 A A A 进行交换
8) A n ( n = 0 , 1... , n ∈ N ) A^n(n=0,1..., n \in N) An(n=0,1...,nN)可与 A m ( n = 0 , 1... , m ∈ N ) A^m(n=0,1..., m \in N) Am(n=0,1...,mN) 交换. 这一点由矩阵乘法的结合律证明


6、置换矩阵

置换矩阵是一种系数只由0和1组成的方块矩阵,置换矩阵的每一行和每一列都恰好有一个1,其余的系数都是0.
当一个矩阵乘上一个置换矩阵时,所得到的是原来矩阵的横行(置换矩阵在左)或纵列(置换矩阵在右)经过置换后得到的矩阵.
一个置换矩阵必然是正交矩阵,并且它的逆也是置换矩阵


7、顺序主子式

顺序主子式是n阶方阵的n个行列式按顺序排列而成,第k个行列式是由该方阵的前k行和k列组成。对于n阶方阵 A A A,其共有n个顺序主子式
A A A i i i 阶顺序主子式为:
D i = ∣ a 11 a 12 … a 1 i a 21 a 22 … a 2 i ⋮ ⋮ : a i 1 a i 2 … a i i ∣ ( i = 1 , 2 , … , n ) D _ { i } = \left| \begin{array} { c c c c } { a _ { 11 } } & { a _ { 12 } } & { \dots } & { a _ { 1 i } } \\ { a _ { 21 } } & { a _ { 22 } } & { \dots } & { a _ { 2 i } } \\ { \vdots } & { \vdots } & { } & { : } \\ { a _ { i 1 } } & { a _ { i 2 } } & { \dots } & { a _ { i i } } \end{array} \right| ( i = 1,2 , \dots , n ) Di= a11a21ai1a12a22ai2a1ia2i:aii (i=1,2,,n)

应用
通过计算方阵A的所有顺序主子式,可以来判断一个实二次型是否正定或方阵A是否为正定矩阵,也可以判断方阵A是否可以进行唯一LU分解


8、分块矩阵

通过将大的矩阵通过分块的方式划分,并将每个分块看做另一个矩阵的元素,这样之后再参与运算,通常可以让计算变得清晰甚至得以大幅简化。例如,有的大矩阵可以通过分块变为对角矩阵或者是三角矩阵等特殊形式的矩阵.
分块矩阵中,位在同一行(列)的每一个子矩阵,都拥有相同的列数(行数).

8.1 矩阵的分块乘法

假设存在 C m × n = A m × p B p × n C^{m \times n} = A^{m \times p}B^{p \times n} Cm×n=Am×pBp×n,如果将 A m × p A^{m \times p} Am×p B p × n B^{p \times n} Bp×n 写成分块矩阵再相乘,那么,
1) A m × p A^{m \times p} Am×p的列分块形式必须与 B p × n B^{p \times n} Bp×n行分块形式相同,
2) A m × p A^{m \times p} Am×p的行分块形式随便,
3) B p × n B^{p \times n} Bp×n的列分块形式随便,举例如下:
[ A 1 3 × 4 A 2 3 × 5 A 3 6 × 4 A 4 6 × 5 ] [ B 1 4 × 1 B 2 4 × 8 B 3 5 × 1 B 4 5 × 8 ] \begin{bmatrix} A_1^{3 \times 4} & A_2^{3 \times 5}\\ A_3^{6 \times 4} & A_4^{6 \times 5} \end{bmatrix} \begin{bmatrix} B_1^{4 \times 1} & B_2^{4 \times 8}\\ B_3^{5 \times 1} & B_4^{5 \times 8} \end{bmatrix} [A13×4A36×4A23×5A46×5][B14×1B35×1B24×8B45×8]

8.2 分块矩阵的行列式

证明参考:分块矩阵的行列式

公式一:
分块对角矩阵的行列式:
∣ A 0 0 D ∣ = ∣ A 0 0 I ∣ ∣ I 0 0 D ∣ = ( det ⁡ A ) ( det ⁡ D ) \left| \begin{array} { c c } { A } & { 0 } \\ { 0 } & { D } \end{array} \right| = \left| \begin{array} { c c } { A } & { 0 } \\ { 0 } & { I } \end{array} \right| \left| \begin{array} { c c } { I } & { 0 } \\ { 0 } & { D } \end{array} \right| = ( \operatorname { det } A ) ( \operatorname { det } D ) A00D = A00I I00D =(detA)(detD)

公式二:
分块上(下)三角矩阵的行列式也等于主对角分块行列式之积,其中 A A A D D D 是方阵(但大小可以相异),即:
∣ A B 0 D ∣ = ( det ⁡ A ) ( det ⁡ D ) \left| \begin{array} { c c } { A } & { B } \\ { 0 } & { D } \end{array} \right| = ( \operatorname { det } A ) ( \operatorname { det } D ) A0BD =(detA)(detD)

公式三:
设A和D是方阵。若A是可逆的,则
∣ A B C D ∣ = ( det ⁡ A ) ( det ⁡ ( D − C A − 1 B ) ) \left| \begin{array} { c c } { A } & { B } \\ { C } & { D } \end{array} \right| = ( \operatorname { det } A ) \left( \operatorname { det } \left( D - C A ^ { - 1 } B \right) \right) ACBD =(detA)(det(DCA1B))
若D可逆:
∣ A B C D ∣ = ( det ⁡ D ) ( det ⁡ ( A − B D − 1 C ) ) \left| \begin{array} { c c } { A } & { B } \\ { C } & { D } \end{array} \right| = ( \operatorname { det } D ) \left( \operatorname { det } \left( A - B D ^ { - 1 } C \right) \right) ACBD =(detD)(det(ABD1C))

公式四:
设A, B, C, D是 n × n n\times n n×n阶矩阵。若A, B, C, D其中之一是零矩阵
∣ A B C D ∣ = det ⁡ ( A D − B C ) \left| \begin{array} { c c } { A } & { B } \\ { C } & { D } \end{array} \right| = \operatorname { det } ( A D - B C ) ACBD =det(ADBC)

还有公式五和公式六,貌似不常用,先占个空,用到再补^~^


下一篇:常见的特殊矩阵及分解(二)


@leatherwang
二零一八年十二月八日

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/hzwwpgmwy/article/details/84848042

智能推荐

详解音视频直播平台软件开发中的低延时-程序员宅基地

文章浏览阅读427次。详解音视频直播平台软件开发中的低延时音视频实时通讯的应用场景已经随处可见,从“吃鸡”的语音对讲、直播连麦、直播答题组队开黑,再到银行视频开户等。对于开发者来讲,除了关注如何能快速实现不同应用场景重点额音视频通讯,另一个更需要关注的可能就是“低延时”。但是,到底实时音视频传输延时应该如何“低”,才能满足你的应用场景呢?延时的产生与优化在聊低延时之前,我们先要讲清延时是如何产生的。由于音视频的传输路径一样,我们可以通过一张图来说明延时的产生:在音视频传输过程中,在不同阶段都...

网络精英赛模拟练习(5)_在企业网络安全应急响应过程中,对当前网络安全情况进行评估时,可以使用()等方法。-程序员宅基地

文章浏览阅读2.3k次,点赞3次,收藏18次。1、作为全方位的、整体的信息安全防范体系是分层次的,以下关于企业信息系统层次划分的描述,错误的是_________(单选题,1分)A.越接近内部的网络安全要求等级越低,越接近外部的网络安全要求等级越高B.业务专用网是企业为了特殊工作需要而建造的专用网络C.互联网区域用于日常的互联网业务,安全防护等级要求最低D.企业内网是企业的核心网络,拥有最高的安全防护等级回答正确(得分: 1分)正确答案A解析暂无解析2、信息安全风险评估是信息系统安全工程的重要组成部分,以下数据的操作与安_在企业网络安全应急响应过程中,对当前网络安全情况进行评估时,可以使用()等方法。

无敌破坏王2之大闹互联网观后感_拉尔夫网络观后感-程序员宅基地

文章浏览阅读5.4k次。因为同事要离职的缘故,而且我们平时又玩的很好,所以决定今天晚上请去看电影,电影院中毒液和无名之辈又都看过,而无敌破坏王的口碑还不错,所以就问是否想看这部电影,结果一拍即合,我和高工两个人陪他去看了这场动画电影。依然是美轮美奂的电影,拉尔夫这个无敌破坏王依然有着第一季中的个性,云泥洛普依然调皮可爱,生动活泼。第一季中由于拉尔夫总是不断的破坏城堡导致自己不能受到同伴们的喜爱而悲伤低落,想要像阿修一样..._拉尔夫网络观后感

c++ 学习中遇到的问题: error C2065: 'ifstream' : undeclared identifier_ifstream undeclared-程序员宅基地

文章浏览阅读2.2k次,点赞2次,收藏4次。出现这个错误有两种可能:1. 没有包含文件流类的头文件 #include2. 没有使用命名空间 using namespace std 或者 std::ifstream#include #include using namespace std;int main(){ifstream infile("C:\\Users\\Administrat_ifstream undeclared

Python中的eval函数_real signature unknown-程序员宅基地

文章浏览阅读1.1w次,点赞9次,收藏12次。代码版本:3.6.3 文档:3.6.6eval() eval()是Python内置函数,点进去是这样的def eval(*args, **kwargs): # real signature unknown """ Evaluate the given source in the context of globals and locals. Th..._real signature unknown

Ubuntu Samba Server_2019 server wsllaunchinteractive /usr/sbin/adduser-程序员宅基地

文章浏览阅读979次。## Sample configuration file for the Samba suite for Debian GNU/Linux.### This is the main Samba configuration file. You should read the# smb.conf(5) manual page in order to understand the option_2019 server wsllaunchinteractive /usr/sbin/adduser --quiet --gecos failed wi

随便推点

iptables防火墙_iptables -t nat-程序员宅基地

文章浏览阅读1.1k次。IP信息包过滤系统,它实际上由两个组件 netfilter_和 iptables组成。主要工作在网络层,针对IP数据包,体现在对包内的IP地址、端口等信息的处理。iptables的作用是为包过滤机制的实现提供规则,通过各种不同的规则,告诉netfilter对来自某些源,前往某些目的或具有某些协议特征的数据包应该如何处理,为了更加方便的组织和管理防火墙规则。iptables采用了表和链的分层结构,所以它会对请求的数据包的包头数据进行分析,根据我们预先设定的规则进行匹配来决定是否可以进入主机。_iptables -t nat

sql级联删除_plsql 两表 级联删除-程序员宅基地

文章浏览阅读905次。sql级联删除功能:在删除主表时,自动删除副表(外键约束)相应内容删除包含主键值的行的操作,该值由其它表的现有行中的外键列引用。在级联删除中,还删除其外键值引用删除的主键值的所有行。如:create database tempgouse tempgocreate table UserInfo(UserId int id_plsql 两表 级联删除

c语言常用的运行速度优化方法,嵌入式C语言性能优化方法-程序员宅基地

文章浏览阅读1.2k次。嵌入式C语言性能优化方法嵌入式系统是指完成一种或几种特定功能的计算机系统,具有自动化程度高,响应速度快等优点,目前已广泛应用于消费电子,工业控制等领域.嵌入式系统受其使用的硬件以及运行环境的限制,非常注重代码的时间和空间效率,因此选择一种合适的开发语言十分重要.嵌入式C语言性能优化方法有哪些?下面是相关的知识,欢迎阅读。使用宏定义在C语言中,宏是产生内嵌代码的唯一方法。对于嵌入式系统而言,为了能达..._c语言编程之运行速度优化方法汇总

复活的 C4C Linux 发行版_linux版本复活-程序员宅基地

文章浏览阅读97次。导读 Computers4Christians 项目以定制发行版的形式进行了改革,该发行版为有基督教信仰的人提供了软件。当我刚开始在这里写作时,我介绍了一个 基督徒的Linux发行版,距离现在已经有 6 个年头了,让我们来速览一下这个项目在 6 年的时间里都有哪些变化吧。 Computers4Christians 项目以定制发行版的形式进行了改革,该发行版为有基督教信仰的人提供了软件。当我刚开始在这里写作时,我介绍了一个 基督徒的 Linux 发行版,距离现在已经有 6 个年头了,让.._linux版本复活

RDIFramework.NET ━ .NET快速信息化系统开发框架 V3.2->新增记录SQL执行过程-程序员宅基地

文章浏览阅读57次。有时我们需要记录整个系统运行的SQL以作分析,特别是在上线前这对我们做内部测试也非常有帮助,当然记录SQL的方法有很多,也可以使用三方的组件。3.2版本我们在框架底层新增了记录框架运行的所有SQl过程保存到用户指定的地方以便分析查看,只需要在配置文件把配置项”LogSQL”设置为True即可。框架会自动记录各常用数据库如:Oracle、SqlServer..._rdi指令 执行sql

如何查看主板的型号和名称-程序员宅基地

文章浏览阅读1.9k次。 电脑刚开机时第一屏,找到HOME键右上角的Pause Break键按下,电脑屏幕会暂停 ,找到屏幕最下边显示的一行就有主板的型号和名称。 下面我讲讲主板选购四大盲点: 对于电脑来说,主板上承载着电脑中最重要的元件,它的重要性是不言而喻的,同时,主板也是各大厂商竞争的舞台——放眼当今主板市场,各大厂商竞相登台,各种型号相继亮相。面对如此热闹的市场,你是否对主板的选购感觉到有一些说不出来的盲..._怎么查看主板型号和品牌

推荐文章

热门文章

相关标签