20211130 正定矩阵的几个不等式
1
Let A∈Rn×nA\in\mathbb{R}^{n\times n}A∈Rn×n be positive definite and B∈Rn×nB\in\mathbb{R}^{n\times n}B∈Rn×n be symmetric. Then, for any x∈Rnx \in \mathbb{R}^{n}x∈Rn, λmin?(A?1B)x?Ax≤x?Bx≤λmax?(A?1B)x?Ax\lambda_{\min }\left(A^{-1} B\right) x^{\top} A x\leq x^{\top} B x \leq \lambda_{\max }\left(A^{-1} B\right) x^{\top} A xλmin?(A?1B)x?Ax≤x?Bx≤λmax?(A?1B)x?Ax
Proof:
Because AAA is positive definite, it is symmetric, which shows that there exists some orthogonal matrices PPP (P?P=PP?=IP^{\top}P=PP^{\top}=IP?P=PP?=I) such that
A=P?1diag{λi(A)}P=P?1diag{λi(A)}PP?1diag{λi(A)}PA=P^{-1}\text{diag}\{\lambda_i(A)\}P=P^{-1}\text{diag}\{\sqrt{\lambda_i(A)}\}PP^{-1}\text{diag}\{\sqrt{\lambda_i(A)}\}PA=P?1diag{λi?(A)}P=P?1diag{λi?(A)?}PP?1diag{λi?(A)?}P
Let M=P?1diag{λi(A)}PM=P^{-1}\text{diag}\{\sqrt{\lambda_i(A)}\}PM=P?1diag{λi?(A)?}P, we have
A=MMA=MMA=MM where MMM is positive definite.
Because MMM is positive definite (obviously, also symmetric), it is shown that M?1M^{-1}M?1 is positive definite (obviously, also symmetric), which, together with that BBB is symmetric, leads to M?1BM?1M^{-1}BM^{-1}M?1BM?1 is symmetric. Therefore, we have
M?1BM?1≤λmax?(M?1BM?1)I=λmax?(M?1M?1B)I=λmax?(A?1B)IM^{-1}BM^{-1}\leq \lambda_{\max}(M^{-1}BM^{-1})I=\lambda_{\max}(M^{-1}M^{-1}B)I=\lambda_{\max}(A^{-1}B)IM?1BM?1≤λmax?(M?1BM?1)I=λmax?(M?1M?1B)I=λmax?(A?1B)I leading to B≤λmax?(A?1B)MM=λmax?(A?1B)AB\leq \lambda_{\max}(A^{-1}B) MM=\lambda_{\max}(A^{-1}B)AB≤λmax?(A?1B)MM=λmax?(A?1B)A
Followng the same line, we have λmin?(A?1B)A≤B\lambda_{\min}(A^{-1}B)A \leq Bλmin?(A?1B)A≤B
As a result, λmin?(A?1B)x?Ax≤x?Bx≤λmax?(A?1B)x?Ax\lambda_{\min }\left(A^{-1} B\right) x^{\top} A x\leq x^{\top} B x \leq \lambda_{\max }\left(A^{-1} B\right) x^{\top} A xλmin?(A?1B)x?Ax≤x?Bx≤λmax?(A?1B)x?Ax holds.
2
令AAA是m×nm\times nm×n維實矩陣,則有 λ(A?A)?0\lambda (A^{\top}A) \geqslant 0 λ(A?A)?0
3
令A1A_1A1?和A2A_2A2?是m×nm\times nm×n維實矩陣,A3A_3A3?是m×mm \times mm×m維的正定陣,則有 A1?A2+A2?A1≤A1?A3A1+A2?A3?1A2A_{1}^{\top} A_{2}+A_{2}^{\top} A_{1} \leq A_{1}^{\top} A_{3} A_{1}+ A_{2}^{\top} A_{3}^{-1} A_{2} A1??A2?+A2??A1?≤A1??A3?A1?+A2??A3?1?A2?
證明:
因為A3A_3A3?是m×mm \times mm×m維的正定陣,存在正交陣PPP使得
A3=P?ΛPA_3=P^{\top}\Lambda PA3?=P?ΛP,其中Λ\LambdaΛ為對角陣,且為正實數,那么存在Λ=Λ~Λ~\Lambda=\tilde\Lambda\tilde\LambdaΛ=Λ~Λ~另外,(Λ~PA1?Λ~?1PA2)?(Λ~PA1?Λ~?1PA2)≥0(\tilde\Lambda P A_1-\tilde\Lambda^{-1}PA_2)^{\top}(\tilde\Lambda P A_1-\tilde\Lambda^{-1}PA_2)\geq0(Λ~PA1??Λ~?1PA2?)?(Λ~PA1??Λ~?1PA2?)≥0展開后可得
A1?A3A1+A2?A3?1A2≥A1?A2+A2?A1A_{1}^{\top} A_{3} A_{1}+A_{2}^{\top} A_{3}^{-1} A_{2}\geq A_{1}^{\top} A_{2}+A_{2}^{\top} A_{1}A1??A3?A1?+A2??A3?1?A2?≥A1??A2?+A2??A1?
4
AAA 和 BBB 是正定矩陣,ABABAB 的特征值均正;如果滿足如下三個等價條件中的一個:① ABABAB 是對稱的;② BABABA是對稱的;③AB=BAAB=BAAB=BA,那么 ABABAB 是正定矩陣。
證明:
AAA 和 BBB 是正定矩陣,所以存在正定矩陣 PPP 和 QQQ ,使得A=P?PB=Q?QA = P^{\top}P\quad\quad B=Q^{\top}QA=P?PB=Q?Q,同時有PPP 和 QQQ的逆存在,那么Q(AB)Q?1=QP?PQ?=(PQ?)?PQ?Q(AB)Q^{-1} = QP^{\top}PQ^{\top}=(PQ^{\top})^{\top}PQ^{\top}Q(AB)Q?1=QP?PQ?=(PQ?)?PQ?因為 PQ?x=0PQ^{\top}x=0PQ?x=0只有一個零解,所以 (PQ?)?PQ?(PQ^{\top})^{\top}PQ^{\top}(PQ?)?PQ? 是正定矩陣,相似于 ABABAB,所以 ABABAB 的特征值均正。很明顯,如果 ABABAB 為對稱陣 ,那么 ABABAB 是正定矩陣。
接下來考慮三個條件的等價:
①?②: 如果 ABABAB 為對稱陣 ,又因為 AAA 和 BBB 是正定矩陣,那么 ATBT=AB=BTAT=BAA^TB^T=AB=B^TA^T=BAATBT=AB=BTAT=BA;
②?③: 如果 BABABA 為對稱陣 ,又因為 AAA 和 BBB 是正定矩陣,那么 BA=ATBT=ABBA=A^TB^T=ABBA=ATBT=AB;
③ ?①:如果 AB=BAAB=BAAB=BA,又因為 AAA 和 BBB 是正定矩陣,所以可知 BA=BTAT=(AB)T=ABBA=B^{T}A^{T}=(AB)^T=ABBA=BTAT=(AB)T=AB,說明 ABABAB 是對稱的。
5
給定任意 x∈Rnx\in\mathbb{R}^nx∈Rn 和 σ>0\sigma>0σ>0 ,那么 σI+xxT\sigma I+ x x^TσI+xxT 是正定的,其逆是 1σ(I?xxTσ+xTx)\frac{1}{\sigma}(I-\frac{xx^T}{\sigma+x^Tx})σ1?(I?σ+xTxxxT?)。
證明:
1σ(σI+xxT)(I?xxTσ+xTx)=I\frac{1}{\sigma}(\sigma I+ x x^T)(I-\frac{xx^T}{\sigma+x^Tx})=Iσ1?(σI+xxT)(I?σ+xTxxxT?)=I
6
對稱矩陣的譜半徑等于其譜范數。
證明:
譜半徑是 max?i∣λ(A)∣\max_i{|\lambda(A)|}maxi?∣λ(A)∣,譜范數是 max?iATA\sqrt{\max_i A^TA}maxi?ATA?。
假設 AAA 是對稱陣,那么存在正交陣 PPP 使得 A=PTXPA=P^T XPA=PTXP,X是以AAA的特征值為對角線的對角矩陣。ATA=PTXXPA^TA=P^T XXPATA=PTXXP ,那么 ATAA^TAATA 的最大特征值就是 AAA 的最大特征值的平方。
總結
以上是生活随笔為你收集整理的20211130 正定矩阵的几个不等式的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 20211126 为什么转动惯量矩阵是正
- 下一篇: 20211201 二范数的导数小于等于导