欢迎访问 生活随笔!

生活随笔

当前位置: 首页 >

20211130 正定矩阵的几个不等式

发布时间:2025/3/15 45 豆豆
生活随笔 收集整理的这篇文章主要介绍了 20211130 正定矩阵的几个不等式 小编觉得挺不错的,现在分享给大家,帮大家做个参考.

1

Let A∈Rn×nA\in\mathbb{R}^{n\times n}ARn×n be positive definite and B∈Rn×nB\in\mathbb{R}^{n\times n}BRn×n be symmetric. Then, for any x∈Rnx \in \mathbb{R}^{n}xRn, λ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(A1B)xAxxBxλmax(A1B)xAx

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}=IPP=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=P1diag{λi(A)}P=P1diag{λi(A)}PP1diag{λi(A)}P
Let M=P−1diag{λi(A)}PM=P^{-1}\text{diag}\{\sqrt{\lambda_i(A)}\}PM=P1diag{λ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}M1 is positive definite (obviously, also symmetric), which, together with that BBB is symmetric, leads to M−1BM−1M^{-1}BM^{-1}M1BM1 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)IM1BM1λmax(M1BM1)I=λmax(M1M1B)I=λmax(A1B)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(A1B)MM=λmax(A1B)A

Followng the same line, we have λmin⁡(A−1B)A≤B\lambda_{\min}(A^{-1}B)A \leq Bλmin(A1B)AB

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(A1B)xAxxBxλmax(A1B)xAx holds.

2

AAAm×nm\times nm×n维实矩阵,则有 λ(A⊤A)⩾0\lambda (A^{\top}A) \geqslant 0 λ(AA)0

3

A1A_1A1A2A_2A2m×nm\times nm×n维实矩阵,A3A_3A3m×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} A1A2+A2A1A1A3A1+A2A31A2

证明:
因为A3A_3A3m×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}A1A3A1+A2A31A2A1A2+A2A1

4

AAABBB 是正定矩阵,ABABAB 的特征值均正;如果满足如下三个等价条件中的一个:① ABABAB 是对称的;② BABABA是对称的;③AB=BAAB=BAAB=BA,那么 ABABAB 是正定矩阵。

证明:
AAABBB 是正定矩阵,所以存在正定矩阵 PPPQQQ ,使得A=P⊤PB=Q⊤QA = P^{\top}P\quad\quad B=Q^{\top}QA=PPB=QQ,同时有PPPQQQ的逆存在,那么Q(AB)Q−1=QP⊤PQ⊤=(PQ⊤)⊤PQ⊤Q(AB)Q^{-1} = QP^{\top}PQ^{\top}=(PQ^{\top})^{\top}PQ^{\top}Q(AB)Q1=QPPQ=(PQ)PQ因为 PQ⊤x=0PQ^{\top}x=0PQx=0只有一个零解,所以 (PQ⊤)⊤PQ⊤(PQ^{\top})^{\top}PQ^{\top}(PQ)PQ 是正定矩阵,相似于 ABABAB,所以 ABABAB 的特征值均正。很明显,如果 ABABAB 为对称阵 ,那么 ABABAB 是正定矩阵。

接下来考虑三个条件的等价:

①➡②: 如果 ABABAB 为对称阵 ,又因为 AAABBB 是正定矩阵,那么 ATBT=AB=BTAT=BAA^TB^T=AB=B^TA^T=BAATBT=AB=BTAT=BA

②➡③: 如果 BABABA 为对称阵 ,又因为 AAABBB 是正定矩阵,那么 BA=ATBT=ABBA=A^TB^T=ABBA=ATBT=AB

③ ➡①:如果 AB=BAAB=BAAB=BA,又因为 AAABBB 是正定矩阵,所以可知 BA=BTAT=(AB)T=ABBA=B^{T}A^{T}=(AB)^T=ABBA=BTAT=(AB)T=AB,说明 ABABAB 是对称的。

5

给定任意 x∈Rnx\in\mathbb{R}^nxRnσ>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}maxiATA
假设 AAA 是对称阵,那么存在正交阵 PPP 使得 A=PTXPA=P^T XPA=PTXP,X是以AAA的特征值为对角线的对角矩阵。ATA=PTXXPA^TA=P^T XXPATA=PTXXP ,那么 ATAA^TAATA 的最大特征值就是 AAA 的最大特征值的平方。

总结

以上是生活随笔为你收集整理的20211130 正定矩阵的几个不等式的全部内容,希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错,欢迎将生活随笔推荐给好友。