王品超上一道关于矩阵不等式的证明及其推广
1.设A为2阶正定矩阵,0<∣∣→α∣∣≤1,则有(αTAα)(αTA‐1α)αTα≤(λ1+λ2)24λ1λ2,其中λ1,λ2是A特征值
证明:由条件知,存在正交矩阵T使TTAT=(λ100λ2),λi>0,i=1,2
令TTα=(y1y2),则(αTAα)(αTA‐1α)αTα=(λ1y21+λ2y22)(λ‐11y21+λ‐12y22)y21+y22=y41+y42+(λ2λ1+λ1λ2)y21y22y21+y22=(y21+y22)2−2y21y22+(λ2λ1+λ1λ2)y21y22y21+y22=y21+y22+(√λ2λ1−√λ1λ2)2y21y22y21+y22=(y21+y22)⎡⎢
⎢⎣1+(√λ2λ1−√λ1λ2)2y21y22(y21+y22)2⎤⎥
⎥⎦≤1+(√λ2λ1−√λ1λ2)2y21y22(y21+y22)2
因为y21+y22≥2|y1y2|,故12≥|y1y2|y21+y22,即14≥y21y22(y21+y22)2,⇒(αTAα)(αTA‐1α)αTα≤1+14(√λ2λ1−√λ1λ2)2=(λ1+λ2)24λ1λ2
更进一步,可以得到更强的不等式
2.设A为n(n≥2)阶正定矩阵,0<∣∣→α∣∣≤1,则有(αTAα)(αTA‐1α)αTα≤(n−1)(λ1+λ2)22nλ1λ2,其中λ1,λ2分别是A最小及最大特征值
证明:由条件知,存在正交矩阵T使TTAT=⎛⎜
⎜⎝λ10⋱0λm⎞⎟
⎟⎠,λi>0,i=1,2,⋯,n,不妨设λ1≤λ2≤⋯≤λn。
令TTα=⎛⎜
⎜⎝y1⋮yn⎞⎟
⎟⎠,有(αTAα)(αTA‐1α)αTα=(αTT)(TTAT)(TTα)(αTT)(TTA‐1T)(TTα)(αTT)(TTα)=(λ1y21+⋯+λny2n)(λ‐11y21+⋯+λ‐1ny2n)y21+⋯+y2n
令△=y21+⋯+y2n,则原式=n∑i=1y4i+∑1≤i<j≤n(λiλj+λjλi)y2iy2j△=△2+∑1≤i<j≤n(√λiλj−√λjλi)2y2iy2j△
由(√λ1λn−√λnλ1)2−(√λiλj−√λjλi)2=λ21+λ2nλ1λn−λ2i+λ2jλiλj=(λnλj−λiλ1)(λiλn−λ1λj)λ1λnλiλj
因为λnλj−λiλ1≥0,λiλn−λ1λj≥0,则(√λ1λn−√λnλ1)2≥(√λiλj−√λjλi)2
⇒原式≤△+(∑1≤i<j≤ny2iy2j)(√λ1λn−√λnλ1)2△=△⎡⎢
⎢
⎢
⎢
⎢
⎢⎣1+(∑1≤i<j≤ny2iy2j)(√λ1λn−√λnλ1)2△2⎤⎥
⎥
⎥
⎥
⎥
⎥⎦
由△2=(y21+⋯+y2n)2=y41+⋯+y4n+2∑1≤i<j≤ny2iy2j
=(y21√n−1−y22√n−1)2+⋯+(y21√n−1−y2n√n−1)2+⋯+(y2n−1√n−1−y2n√n−1)2+2n−1∑1≤i<j≤ny2iy2j+2∑1≤i<j≤ny2iy2j≥2nn−1∑1≤i<j≤ny2iy2j
⇒n−12n≥1△2∑1≤i<j≤ny2iy2j,且因为△≤1
⇒原式≤1+n−12n(√λ1λn−√λnλ1)2=(n−1)λ21+2λ1λn+(n−1)λ2n2nλ1λn≤(n−1)λ21+2(n−1)λ1λn+(n−1)λ2n2nλ1λn=(n−1)(λ1+λn)22nλ1λn
若当—谢瓦莱(Jordan-Chevalley)分解定理证明
引理1:已知g1(x),g2(x),⋯,gk(x)两两互素,r1(x),r2(x),⋯,rk(x)为k个非零多项式,并且∂(ri(x))<∂(gi(x))(i=1,2,⋯,k),求存在一个多项式f(x),被gi(x)除余式是ri(x)
证明:令Gi(x)=g1(x)⋯gi−1(x)gi+1(x)⋯gk(x),因为g1(x),g2(x),⋯,gk(x)两两互素
⇒(gi(x),Gi(x))=1⇒存在ui(x)gi(x)+vi(x)Gi(x)=1
⇒ui(x)gi(x)ri(x)+vi(x)Gi(x)ri(x)=ri(x)⇒vi(x)Gi(x)ri(x)=ri(x)−ui(x)gi(x)ri(x),所以g1(x),⋯,gi−1(x),gi+1(x),⋯,gk(x)都整除vi(x)Gi(x)ri(x),而gi(x)除v(x)Gi(x)ri(x)余ri(x),取
f(x)=k∑i=1vi(x)Gi(x)ri(x)=v1(x)G1(x)r1(x)+⋯+vk(x)Gk(x)rk(x)
⇒gi(x)除f(x)余式ri(x)
引理2:已知A,B都是n级方阵,AB=BA,且A,B均可对角化,证明A,B可以同时对角化。
证明:由于A可对角化。所以存在可逆矩阵P使得P−1AP=diag{λ1E1,λ2E2,⋯,λsEs},其中λ1,λ2,⋯,λs是A所有互异的特征值,E1,E2,⋯,Es分别是r1,r2,⋯,rs级单位矩阵,且r1+r2+⋯+rs=n,由AB=BA得P−1BPP−1AP=P−1APP−1BP
⇒P−1BP=diag{B1,B2,⋯,Bs},B1,B2,⋯,Bs分别是r1,r2,⋯,rs级方阵
由B可对角化,所以P−1BP可以对角化 ⇒B1,B2,⋯,Bs可以对角化
⇒对任意的i=1,2,⋯,s,存在可逆矩阵Qi使得Q−1iBiQi为对角矩阵。
取Q=diag{Q1,Q2,⋯,Qs},则Q−1P−1BPQ=diag{Q−11B1Q1,Q−12B2Q2,⋯,Q−1sBsQs}为对角矩阵
⇒取T=PQ,有T−1BT,T−1BT同时为对角矩阵。
(来自复旦大学高代教材的证明过程)
若当—谢瓦莱(Jordan-Chevalley)分解定理:设A是n阶复矩阵,则A可分解为A=B+C,其中B,C适合下面的条件:
(1)B是一个可对角化的矩阵
(2)C是一个幂零矩阵
(3)BC=CB
(4)B,C均可表示为A的多项式
不仅如此,上述满足条件(1)——(3)分解是唯一的。
证明:先对A的若当标准形J证明结论,
设A的全体不同特征值为λ1,λ2,⋯,λs,且J=⎛⎜
⎜
⎜
⎜
⎜⎝J1J2⋱Js⎞⎟
⎟
⎟
⎟
⎟⎠,其中Ji是属于特征值λi的根子空间对应的块,其阶设为mi.显然,对每个i均有Ji=Mi+Ni,其中Mi=λiI是对角阵,Ni是一个幂零矩阵。又MiNi=NiMi。
令M=⎛⎜
⎜
⎜
⎜
⎜⎝M1M2⋱MS⎞⎟
⎟
⎟
⎟
⎟⎠,N=⎛⎜
⎜
⎜
⎜
⎜⎝N1N2⋱NS⎞⎟
⎟
⎟
⎟
⎟⎠,
则J=M+N,MN=NM,M是对角阵,N是幂零矩阵。
因为(Ji−λiI)mi=0,所以Ji适合多项式(λ−λi)mi.而λi互不相同,因此多项式(λ−λ1)m1,(λ−λ2)m2,⋯,(λ−λs)ms两两互素。由引理1,存在多项式g(λ)满足条件:
g(λ)=hi(λ)(λ−λi)mi+λi,对所有的i=1,2,⋯,s成立(这里hi(λ)也是多项式)。
代入Ji得到g(Ji)=hi(Ji)(Ji−λiI)mi+λiI=λiI=Mi
于是g(J)=⎛⎜
⎜
⎜
⎜
⎜⎝g(J1)g(J2)⋱g(Js)⎞⎟
⎟
⎟
⎟
⎟⎠=⎛⎜
⎜
⎜
⎜
⎜⎝M1M2⋱Ms⎞⎟
⎟
⎟
⎟
⎟⎠=M.
又因为N=J−M=J−g(J),所以N也是J的多项式。
现考虑一般情况:设P−1AP=J,则A=PJP−1=P(M+N)P−1
令B=PMP−1,C=PNP−1,则B是可对角化矩阵,而C是幂零矩阵,又g(A)=g(PJP−1)=Pg(J)P−1=PMP−1=B
又易证明BC=CB(∵MN=NM⇒PMP−1PNP−1=PNP−1PMP−1),C=A−g(A)
(注意方法)最后证明唯一性,假设A有另一满足条件的分解A=B1+C1,则B1−B=C1−C。由B1C1=C1B1不难验证AB1=B1A,AC1=C1A(注:在A=B1+C1右乘B1,得AB1=B1B1+C1B1=B1B1+B1C1=B1A,同理AC1=C1A).
因为B=g(A),故BB1=B1B。同理CC1=C1C。设Cr=0,Ct1=0,用二项式定理知道(C1−C)r+t=0。于是
(B1−B)r+t=(C1−C)r+t=0.
因为BB1=B1B,他们都是可对角化的矩阵,由引理2知道它们可以同时对角化,即存在可逆阵Q,使得Q−1BQ和Q−1B1Q都是对角阵。注意到(Q−1BQ−Q−1B1Q)r+t=[Q−1(B−B1)Q]r+t=Q−1(B−B1)r+tQ=0,
两个对角阵之差依旧是对角阵,这个差的幂要等于零矩阵,这两个矩阵必相等,由此即得
B=B1,于是C=C1.证毕□
一道似乎不好下手的高代题
此题是我发在论坛里面的一个问题,在论坛好友jiangjun7116的指导下,现在把过程完整的总结分享给大家。
设A为n阶实对称正定矩阵,β1,β2,⋯,βn∈Rn是n个关于A共轭的非零列向量,即βi≠0(i=1,2,⋯,n),且βTiAβj=0(i≠j;i,j=1,2,⋯,n),证明:A=n∑i=1AβiβTiAβTiAβi和A−1=n∑i=1βiβTiβTiAβi
证明:根据题目注意到n∑i=1AβiβTiAβTiAβiβj=Aβj,n∑i=1βiβTiβTiAβiAβj=βj
(因为i≠j时βTiAβj=0,所以累加里面只有i=j那一项不等于0)
下证β1,β2,⋯,βn是一组基,即β1,β2,⋯,βn是线性无关的
假设存在一组实数k1,k2,⋯,kn使得k1β1+k2β2+⋯+knβn=0,式子两边左乘βTiA,
可以得出kiβTiAβi=0,然而A是实对称的正定矩阵且βi≠0,可以得出ki=0(i=1,2,⋯,n)
所以β1,β2,⋯,βn是一组线性无关的向量组
⇒(n∑i=1AβiβTiAβTiAβi−A)X=0,(n∑i=1βiβTiβTiAβiA−E)X=0有n个线性无关的解
⇒n∑i=1AβiβTiAβTiAβi−A=0,n∑i=1βiβTiβTiAβiA−E=0
即n∑i=1AβiβTiAβTiAβi=A,n∑i=1βiβTiβTiAβiA=E⇒n∑i=1βiβTiβTiAβi=A−1
“许多问题解决之后都会一阵暗爽,这可能就是我喜欢数学的原因吧!”
周不通问题:
一个方阵总相似于它自己的转置。那么,方阵总合同于自己的转置吗?
即:对于任意方阵A,证明存在可逆方阵P,使得 P'AP=A'。
容易证明,对于A=E+B,B反对称,是对的。猜想一般情形也对。
各位看看。
证明:两个实反对称矩阵彼此合同充要条件是秩相等。
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 开发者必知的日志记录最佳实践
· SQL Server 2025 AI相关能力初探
· Linux系列:如何用 C#调用 C方法造成内存泄露
· AI与.NET技术实操系列(二):开始使用ML.NET
· 记一次.NET内存居高不下排查解决与启示
· Manus重磅发布:全球首款通用AI代理技术深度解析与实战指南
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· 没有Manus邀请码?试试免邀请码的MGX或者开源的OpenManus吧
· 园子的第一款AI主题卫衣上架——"HELLO! HOW CAN I ASSIST YOU TODAY
· 【自荐】一款简洁、开源的在线白板工具 Drawnix
2017-11-21 LaTeX算法排版