LA@方阵相似@相似矩阵的性质@正交相似对角化
文章目录
方阵相似
引言
- 对角阵是矩阵中最简单的一类矩阵
- 对角阵相关的乘法运算是很高效的
- 相似方阵是和对角阵相关的概念
相似矩阵定义
- 设 A , B \bold{{A},\bold{B}} A,B是 n n n阶方阵,如果存在 n n n阶可逆方阵 P \bold{P} P,使得 P − 1 A P = B \bold{P^{-1}{AP}={B}} P−1AP=B,则称方阵 A , B \bold{A},\bold{B} A,B相似,记为 A ∼ B \bold{A}\sim{\bold{B}} A∼B
相似变换
- 对 A \bold{A} A进行运算 P − 1 A P \bold{P^{-1}AP} P−1AP称为对 A \bold{A} A进行相似变换
相似变换矩阵
- 矩阵 P \bold{P} P称为相似变换 P − 1 A P \bold{P^{-1}AP} P−1AP的相似变换矩阵
相似矩阵的矩阵多项式和特征值相同
-
若 n n n阶矩阵 A , B \bold{A,B} A,B相似,则 A , B \bold{A,B} A,B的特征多项式相同,从而 A , B \bold{A,B} A,B的特征值相同
-
证明:
- 由 A ∼ B \bold{A\sim{B}} A∼B,有 P \bold{P} P满足 P − 1 A P = B \bold{P^{-1}AP=B} P−1AP=B
- 所以 f B ( λ ) f_{\bold{B}}(\lambda) fB(λ)= ∣ B − λ E ∣ \bold{|B-\lambda{E}|} ∣B−λE∣= ∣ P − 1 A P − λ E ∣ |\bold{P^{-1}AP-\lambda{E}}| ∣P−1AP−λE∣
- 由于 P λ E P − 1 \bold{P\lambda{E}P^{-1}} PλEP−1= λ P E P − 1 \lambda\bold{PEP^{-1}} λPEP−1= λ P P − 1 \lambda\bold{PP^{-1}} λPP−1= λ E \lambda\bold{E} λE,因此,可以将 λ E \bold{\lambda{E}} λE变形为 P ( λ E ) P − 1 \bold{P(\lambda{E})P^{-1}} P(λE)P−1或 P − 1 ( λ E ) P \bold{P^{-1}(\lambda{E})P} P−1(λE)P
- f B ( λ ) = ∣ P − 1 A P − P − 1 ( λ E ) P ∣ f_{\bold{B}}(\lambda)=|\bold{P^{-1}AP-\bold{P^{-1}(\lambda{E})P}}| fB(λ)=∣P−1AP−P−1(λE)P∣= ∣ P − 1 ( A − λ E ) P ∣ |\bold{P^{-1}(A-\lambda{E})P}| ∣P−1(A−λE)P∣= ∣ P − 1 ∣ ∣ A − λ E ∣ ∣ P ∣ \bold{|P^{-1}||A-\lambda{E}||P|} ∣P−1∣∣A−λE∣∣P∣= ∣ P ∣ − 1 ∣ A − λ E ∣ ∣ P ∣ \bold{|P|^{-1}|A-\lambda{E}||P|} ∣P∣−1∣A−λE∣∣P∣= ∣ A − λ E ∣ \bold{|A-\lambda{E}|} ∣A−λE∣
- 显然 f A ( λ ) = f B ( λ ) f_{\bold{A}}(\lambda)=f_{\bold{B}}(\lambda) fA(λ)=fB(λ)
-
但是,特征值相同的方阵未必相似
推论:与对角阵相似的矩阵性质定理
-
与对角阵相似的矩阵的特征值就是对角阵对角元素
-
若 n n n阶矩阵 A ∼ Λ ( λ 1 , ⋯ , λ n ) \bold{A}\sim{\Lambda(\lambda_1,\cdots,\lambda_n)} A∼Λ(λ1,⋯,λn),则 λ 1 , ⋯ , λ n \lambda_1,\cdots,\lambda_n λ1,⋯,λn是 A A A的 n n n个特征值
-
证明:对角阵的特征值是对角元素,由本节定理可知, A \bold{A} A的特征值与 Λ \Lambda Λ相同,所以推论成立
相似矩阵性质
- A ∼ A \bold{A}\sim{\bold{A}} A∼A
-
A
∼
B
⇒
B
∼
A
\bold{A}\sim{\bold{B}}\Rightarrow{\bold{B}\sim{\bold{A}}}
A∼B⇒B∼A
- P − 1 A P = B , A = P B P − 1 \bold{P}^{-1}\bold{A}\bold P=\bold{B},\bold{A}=\bold P\bold{B}\bold P^{-1} P−1AP=B,A=PBP−1
-
A
∼
B
,
B
∼
C
⇒
A
∼
C
\bold{A}\sim{\bold{B}},\bold{B}\sim{\bold C}\Rightarrow{\bold{A}\sim{\bold C}}
A∼B,B∼C⇒A∼C
- P − 1 A P = B , Q − 1 B Q = C \bold{P^{-1}{A}P=\bold{B},Q^{-1}\bold{B}Q=C} P−1AP=B,Q−1BQ=C
- Q C Q − 1 = B = P − 1 A P \bold{QCQ^{-1}=\bold{B}=P^{-1}\bold{A}P} QCQ−1=B=P−1AP
- P Q C Q − 1 P − 1 = A \bold{PQCQ^{-1}P^{-1}=\bold{A}} PQCQ−1P−1=A
- ( P Q ) − 1 = Q − 1 P − 1 \bold{(PQ)^{-1}=Q^{-1}P^{-1}} (PQ)−1=Q−1P−1
- 因此 C ∼ A \bold{C\sim{{A}}} C∼A
- 单位矩阵只和自身相似
- 设方阵 A \bold{A} A和单位阵 E \bold{E} E相似
- P − 1 A P = E \bold{P^{-1}{A}P=E} P−1AP=E
- A = P E P − 1 = E \bold{A=PEP^{-1}=E} A=PEP−1=E
- 因此和单位阵 E \bold{E} E相似的矩阵是 E \bold{E} E本身
导出性质
-
设 A , B \bold{A,{B}} A,B相似 ( A = P − 1 B P ) (\bold{A=P^{-1}BP}) (A=P−1BP),则有以下性质
-
∣ A ∣ = ∣ B ∣ |\bold{A}|=|\bold{B}| ∣A∣=∣B∣
- ∣ A ∣ = ∣ P − 1 B P ∣ = ∣ P − 1 ∣ ∣ B ∣ ∣ P ∣ = ∣ P ∣ − 1 ∣ P ∣ ∣ B ∣ = ∣ B ∣ |\bold{A}|=|\bold P^{-1}\bold{B}\bold P|=|\bold P^{-1}||\bold{B}||\bold P|=|\bold P|^{-1}|\bold P||\bold{B}|=\bold{|B|} ∣A∣=∣P−1BP∣=∣P−1∣∣B∣∣P∣=∣P∣−1∣P∣∣B∣=∣B∣
- 或者,因为相似方阵有相同的特征值,且特征值之积等于方阵的行列式,所以相似方阵的行列式相同
-
t r ( A ) = t r ( B ) tr(\bold{A})=tr(\bold{B}) tr(A)=tr(B)
- A , B \bold{A},\bold{B} A,B具有相同的特征值
- t r ( A ) = ∑ i = 1 n a i i = ∑ i = 1 n λ i tr(\bold{A})=\sum\limits_{i=1}^{n}a_{ii}=\sum\limits_{i=1}^{n}\lambda_{i} tr(A)=i=1∑naii=i=1∑nλi
- t r ( B ) = ∑ i = 1 n b i i = ∑ i = 1 n λ i tr(\bold{B})=\sum\limits_{i=1}^{n}b_{ii}=\sum\limits_{i=1}^{n}\lambda_{i} tr(B)=i=1∑nbii=i=1∑nλi
- ∴ t r ( A ) = t r ( B ) \therefore tr(\bold{A})=tr(\bold{B}) ∴tr(A)=tr(B)
-
r ( A ) = r ( B ) r(\bold{A})=r(\bold{B}) r(A)=r(B)
-
A
=
P
−
1
B
P
\bold{A}=\bold P^{-1}\bold{B}\bold P
A=P−1BP
- P , P − 1 \bold{P,P^{-1}} P,P−1都是可逆矩阵,它们都可以表示为一系列的初等矩阵的乘积
- 因此, A \bold{A} A相当于有 B \bold{B} B经过初等变换得到的等价矩阵,它们的秩相等(初等变换不改变秩)
-
A
=
P
−
1
B
P
\bold{A}=\bold P^{-1}\bold{B}\bold P
A=P−1BP
-
A T ∼ B T \bold{A}^T\sim{\bold{B}^T} AT∼BT
- P − 1 A P = B \bold{P^{-1}{A}P={B}} P−1AP=B
-
(
P
−
1
A
P
)
T
=
B
T
\bold{(P^{-1}{A}P)^T={B}^T}
(P−1AP)T=BT
- P T A T ( P − 1 ) T = B T \bold{P^T\bold{A}^T(P^{-1})^T=\bold{B}^T} PTAT(P−1)T=BT
- P T A T ( P T ) − 1 = B T \bold{P^T\bold{A}^T(P^{T})^{-1}=\bold{B}^T} PTAT(PT)−1=BT
- 可见 A T ∼ B T \bold{A}^T\sim{\bold{B}^T} AT∼BT
-
A m ∼ B m \bold{A}^m\sim{\bold{B}^m} Am∼Bm
-
B m \bold{B}^m Bm= ( P − 1 A P ) m (\bold P^{-1}\bold{A}\bold P)^m (P−1AP)m
- = ( P − 1 A P ) ( P − 1 A P ) ⋯ ( P − 1 A P ) (\bold P^{-1}\bold{A}\bold P)(\bold P^{-1}\bold{A}\bold P)\cdots{(\bold P^{-1}\bold{A}\bold P)} (P−1AP)(P−1AP)⋯(P−1AP)
- = P − 1 A ( P P − 1 ) A ( P P − 1 ) ⋯ ( P P − 1 ) A P =\bold P^{-1}\bold{A}(\bold {PP}^{-1})\bold{A}(\bold P\bold P^{-1})\cdots{(\bold P\bold P^{-1})\bold{A}\bold P} =P−1A(PP−1)A(PP−1)⋯(PP−1)AP
- = P − 1 A m P =\bold P^{-1}\bold{A}^m\bold P =P−1AmP
-
-
设 A = P − 1 B P \bold{A=P^{-1}BP} A=P−1BP
<1>
,若 A − 1 \bold{A}^{-1} A−1存在,则 B − 1 \bold{B}^{-1} B−1存在,且 A − 1 ∼ B − 1 \bold{A}^{-1}\sim{\bold{B}^{-1}} A−1∼B−1, A ∗ ∼ B ∗ \bold{A}^*\sim{\bold{B}^*} A∗∼B∗- B \bold{B} B可逆:
- 方法1:
- A ∼ B ⇒ ∣ A ∣ = ∣ B ∣ = k \bold{A}\sim{\bold{B}}\Rightarrow{|\bold{A}|=|\bold{B}|}=k A∼B⇒∣A∣=∣B∣=k
- A − 1 \bold{A}^{-1} A−1存在, ∣ A ∣ ≠ 0 |\bold{A}|\neq{0} ∣A∣=0,则 ∣ B ∣ = ∣ A ∣ ≠ 0 |\bold{B}|=|\bold{A}|\neq{0} ∣B∣=∣A∣=0
- 方法2:
- 由于 A \bold{A} A可逆,则 P − 1 A P = B \bold P^{-1}\bold{A}\bold P=\bold{B} P−1AP=B表明, B \bold{B} B是可逆矩阵经过一系列初等变换得到的,所以 B \bold{B} B也可逆
- 对
<1>
两边取逆: A − 1 = P B − 1 P − 1 \bold{A}^{-1}=\bold P\bold{B}^{-1}\bold P^{-1} A−1=PB−1P−1<1.1>
,因此 A − 1 ∼ B − 1 \bold{A}^{-1}\sim{\bold{B}^{-1}} A−1∼B−1 - 因为 ∣ A ∣ = ∣ B ∣ = k \bold{|A|=|B|}=k ∣A∣=∣B∣=k, k ≠ 0 k\neq{0} k=0
- A − 1 = 1 ∣ A ∣ A ∗ = k − 1 A ∗ \bold{A}^{-1}=\frac{1}{|\bold{A}|}\bold{A}^*=k^{-1}\bold{A}^* A−1=∣A∣1A∗=k−1A∗
- B − 1 = 1 ∣ B ∣ B ∗ = k − 1 B ∗ \bold{B}^{-1}=\frac{1}{|\bold{B}|}\bold{B}^{*}=k^{-1}\bold{B}^* B−1=∣B∣1B∗=k−1B∗
- 将
<2>,<3>
代入<1.1>
,得 k − 1 A ∗ = P ( k − 1 B ∗ ) P − 1 k^{-1}\bold{A}^*=\bold{P}(k^{-1}\bold{B}^*)\bold{P}^{-1} k−1A∗=P(k−1B∗)P−1 - A ∗ = P − 1 B ∗ P \bold{A}^*=\bold P^{-1}\bold{B}^*\bold P A∗=P−1B∗P,即 A ∗ ∼ B ∗ \bold{A}^*\sim{\bold{B}^*} A∗∼B∗
相似矩阵的乘方性质
设
A
,
B
\bold{A,B}
A,B相似,即
A
=
P
B
P
−
1
\bold{A=PB}\bold{P}^{-1}
A=PBP−1<0>
,则:
A
k
\bold{A}^k
Ak=
P
B
k
P
−
1
\bold{P}\bold{B}^k{\bold{P}^{-1}}
PBkP−1<1>
- 证:对
<0>
两边同时取 k k k次幂运算: A k \bold{A}^k Ak= ( P B P − 1 ) ( P B P − 1 ) ⋯ ( P B P − 1 ) (\bold{P}\bold{B}{\bold{P}^{-1}})(\bold{P}\bold{B}{\bold{P}^{-1}})\cdots(\bold{P}\bold{B}{\bold{P}^{-1}}) (PBP−1)(PBP−1)⋯(PBP−1)- = P B ( P − 1 P ) B ( P − 1 P ) B ⋯ B ( P − 1 P ) B P − 1 \bold{P}\bold{B}(\bold{P}^{-1}\bold{P})\bold{B}{(\bold{P}^{-1}\bold{P})}\bold{B}\cdots\bold{B}(\bold{P}^{-1}\bold{P})\bold{B}{\bold{P}^{-1}} PB(P−1P)B(P−1P)B⋯B(P−1P)BP−1
- = P B k P − 1 \bold{P}\bold{B}^k{\bold{P}^{-1}} PBkP−1
相似矩阵和矩阵多项式
- 设矩阵多项式
f
(
A
)
=
∑
i
=
0
m
a
i
A
i
f(\bold{A})=\sum\limits_{i=0}^{m}a_i\bold{A}^i
f(A)=i=0∑maiAi
<2>
,将<1>
代入<2>
有: f ( A ) f(\bold{A}) f(A)= ∑ i = 0 m a i A i \sum\limits_{i=0}^{m}a_i\bold{A}^i i=0∑maiAi= ∑ i = 0 m a i ( P B i P − 1 ) ) \sum\limits_{i=0}^{m}a_i(\bold P\bold{B}^i{\bold P^{-1}})) i=0∑mai(PBiP−1))= ∑ i = 0 m P ( a i B i ) P − 1 \sum\limits_{i=0}^{m}\bold P(a_i\bold{B}^i){\bold P^{-1}} i=0∑mP(aiBi)P−1,根据矩阵乘法的分配律, f ( A ) f(\bold{A}) f(A)= P ( ∑ i = 0 m a i B i ) P − 1 \bold{P}(\sum_{i=0}^{m}a_i\bold{B}^{i})\bold{P}^{-1} P(∑i=0maiBi)P−1= P f ( B ) P − 1 \bold Pf(\bold{{B}})\bold P^{-1} Pf(B)P−1
相似对角阵
-
当 A \bold{A} A相似于某个对角阵 Λ \bold{\Lambda} Λ,则:
- A k \bold{A}^k Ak= P Λ k P − 1 \bold{P}\bold{\Lambda}^k{\bold{P}^{-1}} PΛkP−1
- f ( A ) f(\bold{A}) f(A)= P f ( Λ ) P − 1 \bold Pf(\bold{{\Lambda}})\bold P^{-1} Pf(Λ)P−1
-
由此可见,若矩阵 A \bold{A} A能够表示成 A = P Λ P − 1 \bold{A}=\bold P\Lambda{\bold P^{-1}} A=PΛP−1(相似对角化问题),矩阵 A \bold{A} A的多项式问题就能够被转换为对角阵的多项式
对角阵多项式的展开
-
f ( Λ ) = ∑ i = 0 m a i Λ i f(\bold\Lambda)=\sum_{i=0}^{m}a_{i}\bold\Lambda^{i} f(Λ)=∑i=0maiΛi= diag ( f ( λ 1 ) , f ( λ 2 ) , ⋯ , f ( λ n ) ) \text{diag}(f(\lambda_1),f(\lambda_2),\cdots,f(\lambda_n)) diag(f(λ1),f(λ2),⋯,f(λn)),
-
推导:
-
对角阵乘方运算性质:若 Λ = d i a g ( λ 1 , λ 2 , ⋯ , λ n ) \bold\Lambda=\mathrm{diag}(\lambda_1,\lambda_{2},\cdots,\lambda_{n}) Λ=diag(λ1,λ2,⋯,λn)为对角阵,则 Λ k \bold\Lambda^k Λk= d i a g ( λ 1 k , λ 2 k , ⋯ , λ n k ) \mathrm{diag}(\lambda_1^k,\lambda_{2}^k,\cdots,\lambda_{n}^k) diag(λ1k,λ2k,⋯,λnk)
-
f ( Λ ) = a 0 ( 1 1 ⋱ 1 ) + a 1 ( λ 1 λ 2 ⋱ λ n ) + ⋯ + a n ( λ 1 n λ 2 n ⋱ λ n n ) = ( ∑ i = 0 m a i λ 1 i ∑ i = 0 m a i λ 2 i ⋱ ∑ i = 0 m a i λ n i ) = ( f ( λ 1 ) f ( λ 2 ) ⋱ f ( λ n ) ) \begin{aligned} f(\Lambda) =&\small{a_0\begin{pmatrix} {{1}} & {} & {} & {} \cr {} & {{ 1}} & {} & {} \cr {} & {} & \ddots & {} \cr {} & {} & {} & {{1}} \cr \end{pmatrix} +a_1\begin{pmatrix} {{\lambda _1}} & {} & {} & {} \cr {} & {{\lambda _2}} & {} & {} \cr {} & {} & \ddots & {} \cr {} & {} & {} & {{\lambda _n}} \cr \end{pmatrix} +\cdots +a_n\begin{pmatrix} {{\lambda _1^n}} & {} & {} & {} \cr {} & {{\lambda _2^n}} & {} & {} \cr {} & {} & \ddots & {} \cr {} & {} & {} & {{\lambda _n^n}} \cr \end{pmatrix}} \\ =&\begin{pmatrix} \sum_{i=0}^{m}a_{i}\lambda_1^{i} & {} & {} & {} \cr {} & \sum_{i=0}^{m}a_{i}\lambda_2^{i} & {} & {} \cr {} & {} & \ddots & {} \cr {} & {} & {} & \sum_{i=0}^{m}a_{i}\lambda_n^{i} \cr \end{pmatrix} =\begin{pmatrix} {f({\lambda _1}}) & {} & {} & {} \cr {} & f({{\lambda _2}}) & {} & {} \cr {} & {} & \ddots & {} \cr {} & {} & {} & f({{\lambda _n}}) \cr \end{pmatrix} \end{aligned} f(Λ)==a0 11⋱1 +a1 λ1λ2⋱λn +⋯+an λ1nλ2n⋱λnn ∑i=0maiλ1i∑i=0maiλ2i⋱∑i=0maiλni = f(λ1)f(λ2)⋱f(λn)
-
这个展开式告诉我们,对角阵(矩阵)的多项式可以归结为数(标量)的多项式的计算
-
小结
- 上述结论说明,相似阵之间有很多共同点
- 特别是,当 A \bold{A} A有一个与之相似的对角阵时,许多关于 A \bold{A} A的计算就可以被简化,例如矩阵多项式的计算,这个问题归结为方阵相似对角化
强矩阵相似
对角化相似
- 若方阵 Q − 1 A Q = Λ \bold{Q^{-1}AQ=\Lambda} Q−1AQ=Λ, Λ \bold{\Lambda} Λ是对角阵,则称 A \bold{A} A对角化相似于 B \bold B B
- 这是一种重要的相似关系,"对角化"指的是原矩阵 A \bold{A} A被相似对角化为一个对角阵,这使得许多计算变得简单,特别是矩阵多项式的相关计算种
正交相似
- 若方阵 Q − 1 A Q = B \bold{Q^{-1}AQ=B} Q−1AQ=B, Q − 1 = Q T \bold{Q^{-1}=Q^{T}} Q−1=QT,则称 A \bold{A} A正交相似于 B \bold B B
- 这里的正交相似的正交体现在相似变换阵是一个正交阵上,而不是说将原矩阵 A \bold{A} A转换为一个正交阵
正交相似对角化
- 正交相似对角化:若方阵 Q − 1 A Q = Λ \bold{Q^{-1}AQ=\Lambda} Q−1AQ=Λ, Q − 1 = Q T \bold{Q^{-1}=Q^{T}} Q−1=QT,则 A \bold{A} A正交相似于对角阵 Λ \bold{\Lambda} Λ,称 A \bold{A} A正交相似对角化
- 在对称阵的对角化相关定理中,利用本概念,可以描述定理:对称阵一定可以正交相似对角化
实对称方阵A正交对角化方法
- 求出实对称阵A的全部特征值
λ
1
,
⋯
,
λ
m
\lambda_1,\cdots,\lambda_m
λ1,⋯,λm(对称阵一定可以正交对角化)
- 如果特征值 λ i \lambda_i λi是单根,则从 f ( λ i ) = 0 , 即 ( λ i E − A ) x = 0 f(\lambda_i)=0,即(\lambda_iE-A)x=0 f(λi)=0,即(λiE−A)x=0对应的求出一个特征向量 α i \alpha_i αi
- 如果特征值
λ
i
\lambda_i
λi是
n
i
n_i
ni重根,
∑
i
=
1
m
n
i
=
n
\sum_{i=1}^{m}n_i=n
∑i=1mni=n
- 则从 f ( λ i ) = 0 f(\lambda_i)=0 f(λi)=0求出 n i n_i ni个线性无关特征向量 A i = α i 1 , ⋯ , α n i A_i=\alpha_{i_1},\cdots,\alpha_{n_i} Ai=αi1,⋯,αni, i = 1 , ⋯ , m i=1,\cdots,m i=1,⋯,m
- 对 A i A_i Ai执行Schmidt正交化得到 B i B_i Bi, i = 1 , ⋯ , m i=1,\cdots,m i=1,⋯,m
- 对 B i B_i Bi在执行Normalization单位化得到 C i C_i Ci, i = 1 , ⋯ , m i=1,\cdots,m i=1,⋯,m
- 将得到的所有标准正交特征向量组 C i C_i Ci中的向量依此排列起来得到正交矩阵 C \bold C C= ( c 11 , ⋯ , c 1 n 1 , ⋯ , c m 1 , ⋯ , c m , n m ) (\bold{c}_{11},\cdots,\bold{c}_{1n_1},\cdots,\bold{c}_{m1},\cdots,\bold{c}_{m,n_m}) (c11,⋯,c1n1,⋯,cm1,⋯,cm,nm), Λ = diag ( λ 1 , ⋯ , λ n ) \Lambda=\text{diag}(\lambda_1,\cdots,\lambda_n) Λ=diag(λ1,⋯,λn)
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了