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⊤A3A1+A2⊤A3−1A2
证明:
因为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⊤A3A1+A2⊤A3−1A2≥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
对称矩阵的谱半径等于其谱范数。
证明:
谱半径是 maxi∣λ(A)∣\max_i{|\lambda(A)|}maxi∣λ(A)∣,谱范数是 maxiATA\sqrt{\max_i A^TA}maxiATA。
假设 AAA 是对称阵,那么存在正交阵 PPP 使得 A=PTXPA=P^T XPA=PTXP,X是以AAA的特征值为对角线的对角矩阵。ATA=PTXXPA^TA=P^T XXPATA=PTXXP ,那么 ATAA^TAATA 的最大特征值就是 AAA 的最大特征值的平方。
总结
以上是生活随笔为你收集整理的20211130 正定矩阵的几个不等式的全部内容,希望文章能够帮你解决所遇到的问题。
- 上一篇: 20211126 为什么转动惯量矩阵是正
- 下一篇: 20211201 二范数的导数小于等于导