眩しさだけは、忘れなかった。|

相似、合同、等价

一、相似矩阵

1. 特征值与特征向量

(1)定义

n 阶矩阵 A 满足 Aα=λα(λ0),则 λA 的特征值,αA 的属于 λ 的特征向量,|λEA|=0A 的特征多项式。

【注】特征向量不能是零向量!

(2)特征值的性质

(2.1)设 An 阶矩阵,特征值为 λ1,λ2,...,λn,则:

  • λ1+λ2+...+λn=tr(A)=a11+a22+...+ann
  • λ1λ2···λn=|A|
  • A 有特征值 λ,则 λ 的重数 knr(λEA)
  • r(A)1,则 A 的特征值为 0,0,...,0,tr(A)(有 n10
  • A 为三角矩阵或对角矩阵,则 A 的特征值为主对角线上的元素
  • α0,则矩阵 αβT 的特征值为 0,0,...,0,βTα,其中特征值 βTα 对应的特征向量为 α

(2.2)设 An实对称矩阵,则:

  • A 的元素均为实数,且 AT=A
  • A 的特征值必为实数
  • A=αβT+βαT 为实对称矩阵

(3)特征向量的性质

(3.1)设 An 阶矩阵,则:

  • A 的不同特征值对应的特征向量线性无关
  • A 的不同特征值对应的特征向量之线性组合不是 A 的特征向量
  • Ak 重特征值 λ,则属于 λ 的线性无关的特征向量个数 s=nr(λEA)k
  • Ak 重特征值 λ,则属于 λ 的线性无关的特征向量之线性组合仍为 A 的特征向量

(3.2)设 An实对称矩阵,则:

  • A不同特征值的特征向量相互正交
  • Ak 重特征值 λ,则属于 λ 的线性无关的特征向量个数 s=nr(λEA)=k

(4)常用结论

矩阵 A kA An A+kE f(A) A1 A P1AP AT
特征值 λ kλ λn λ+k f(λ) 1λ |A|λ λ λ
特征向量 α α α α α α α P1α 不一定是 α

【注 1】关于 Af(A) 的几个要点:

  • f(A)=0,则 A 的每个特征值 λ 都满足 f(λ)=0
  • f(λ)=0 求得解 λ1,λ2,...,λt,则 A 的特征值可能有 λ1,λ2,...,λt 的其中几个(或一个都没有!),但不能确定 A 的特征值一定都有 λ1,λ2,...,λt

【注 2】一个可以快速计算矩阵 A 的特征值的技巧:将 A 拆分成数量矩阵 kE 加一个秩为 1 的矩阵 B,于是 A=B+kE,矩阵 B 的特征值容易写出,自然也就得到矩阵 A 的特征值了。

2. 相似关系

(1)相似的定义

A,Bn 阶矩阵,若存在可逆矩阵 P,使得 P1AP=B(即 AP=PB),则称 AB 相似,记作 AB

(2)相似的性质

(2.1)对称性和传递性:

  • 对称性ABBA
  • 传递性AB,BCAC

(2.2)若 ABP1AP=B 时,有:

  • ABATBT
  • ABAnBn
  • ABA1B1,并且 P1A1P=B1(当 A 可逆时)
  • ABAB,并且 P1AP=B(当 A 可逆时)
  • ABf(A)f(B),并且 P1f(A)P=f(B)
  • ηA 的特征向量 P1ηB 的特征向量
  • AB|A|=|B|
  • AB|λEA|=|λEB|
  • ABr(A)=r(B)
  • ABtr(A)=tr(B)
  • ABAB 有相同的特征值

【注 1】由 AB 可得出以上结论,但反过来,这些条件却并不能得到 AB
【注 2】以上条件只要有一个不满足,即可判断 A 不相似于 B
【注 3】若要判断 AB,则可尝试求出 AB 的对角矩阵,若它们都相似于同一个对角矩阵(即 AΛ,BΛ),则根据相似的传递性,可得 AB

(2.3)设 A,B 均为 n实对称矩阵,则:

  • A 必相似于实对角矩阵,即 AΛ
  • ABAB 有相同的特征值及重数 |λEA|=|λEB|

【注】|λEA|=|λEB| AB 有相同的特征值及重数AΛ,BΛAB

  • 存在正交矩阵 Q,使得 QTAQ=Q1AQ=Λ

3. 相似对角化

(1)相似对角化的定义

n 阶对角矩阵 Λ=diag(λ1,λ2,...,λn),其中 λiA 的特征值, 若存在可逆矩阵 P,使得 P1AP=Λ(即 AP=PΛ),则称 A相似对角化,简称为可对角化,记作 AΛ

(2)可对角化的判别

(设 An 阶矩阵)

  • An 个不同的特征值 AΛ
  • AΛAn 个线性无关的特征向量
  • AΛAk 重特征值 λk 个线性无关的特征向量 k=nr(λEA)

【注】当 k=1 时,k=nr(λEA) 一定成立,因此只需对 k2 的特征值进行判断。

  • A 满足 (AaE)(AbE)=0(ab)AΛ,且特征值满足 (λa)(λb)=0
  • 实对称矩阵必可相似对角化,且正交于对角矩阵

(3)相似对角化的步骤

(3.1)一般矩阵 A 的相似对角化的步骤:

  • |λEA|=0 求出 A 的特征值 λ1,λ2,...,λn
  • 对每个 λi,由 (λiEA)x=0 求出 A 的一组特征向量 α1,α2,...,αn
  • P=(α1,α2,...,αn),当 P 可逆时,有 P1AP=Λ=diag(λ1,λ2,...,λn)

【注】有些题目没有给出具体的矩阵 A,只给定一些已知的特征值和特征向量,要求反求出矩阵 A。有两种解法:

  • 传统的解法:求出对角矩阵 Λ=diag(λ1,λ2,λ3) 和可逆矩阵 P=(α1,α2,α3),然后由 P1AP=Λ 得到 A=PΛP1。这种解法需要对 P 求逆(需使用初等行变换求出),然后进行两次矩阵乘法。
  • 较快的解法:求解矩阵 A 的过程实际上是在求解一个矩阵方程。因为 P1AP=ΛAP=PΛA(α1,α2,α3)=(λ1α1,λ2α2,λ3α3),取转置即得矩阵方程 (α1,α2,α3)TAT=(λ1α1,λ2α2,λ3α3)T,于是求解 AT 只需进行初等行变换:(α1T,α2T,α3T|λ1α1T,λ2α2T,λ3α3T)(E|AT)。显然该法比传统解法要更快!

(3.2)实对称矩阵 A 的相似对角化的步骤:

  • |λEA|=0 求出 A 的特征值 λ1,λ2,...,λn
  • 对每个 λi,由 (λiEA)x=0 求出 A 的一组特征向量 α1,α2,...,αn
  • 对单重特征值的特征向量进行单位化;对多重特征值对应的特征向量进行施密特正交化单位化
  • 令正交矩阵 Q=(η1,η2,...,ηn),有 QTAQ=Q1AQ=Λ=diag(λ1,λ2,...,λn)

【注 1】有些题目没有给出具体的实对称矩阵 A,只给出其中一个或两个特征向量 α,若需要算出其他的特征向量,应使用“实对称矩阵的不同特征值的特征向量相互正交”这一性质来求解。

【注 2】有些题目没有给出具体的实对称矩阵 A,只给定一些已知的特征值和特征向量,要求反求出实对称矩阵 A,一种比较快的解法是使用谱分解定理,在本人专栏有涉及。

(3.3)施密特正交化

{β1=α1β2=α2(α2,β1)(β1,β1)β1β3=α3(α3,β1)(β1,β1)β1(α3,β2)(β2,β2)β2

【注】施密特正交化的推导过程以及其延伸出的一些解题思路,可见本人专栏的另一篇文章。

(3.4)单位化

{η1=β1/||β1||η2=β2/||β2||η3=β3/||β3||

||βi|| 为向量的长度)

令正交矩阵 Q=(η1,η2,η3),有 QTAQ=Q1AQ=Λ=diag(λ1,λ2,λ3)

二、合同矩阵

1. 二次型

(1)二次型的定义

  • 二次型f(x1,x2,...,xn)=xTAx,其中 A 是实对称矩阵。
  • 标准二次型:若交叉项的系数为 0,则得到标准二次型A 是实对角矩阵。
  • 规范二次型:若每一项去掉系数,只保留正负,则得到规范二次型A 是实规范对角矩阵,即 [EpEqO],其中 p正惯性指数(正平方项个数),q负惯性指数(负平方项个数)。

【注 1】标准二次型是不唯一的,规范二次型是唯一的。

【注 2】部分二次型所对应的矩阵不是实对称矩阵,则需要将其改成实对称矩阵:B=12(A+AT)。相关例题:f(x1,x2,x3)=(x1+2x2+3x3)(x12x2+3x3)

  • 可逆线性变量替换:对二次型 f(x1,x2,...,xn)=xTAx 引进新变量 y1,y2,...,yn 用来表示 x1,x2,...,xn

{x1=c11y1+c12y2+...+c1nynx2=c21y1+c22y2+...+c2nyn..................x3=cn1y1+cn2y2+...+cnnyn

将其中的系数矩阵记为 C,若 C 为可逆矩阵,则称为可逆线性变量替换,上式又可写成:x=Cy,所以二次型可化为:f(x1,x2,...,xn)=yTCTACy,可逆线性变量替换后的二次型为 g(y1,y2,...,yn)=CTAC

【注】坐标变换必须可逆,若不可逆则变换后的结果不是二次型 f 的标准型。相关例题:f(x1,x2,x3)=(x1+x2)2+(x2x3)2+(x3+x1)2

(2)惯性定理

标准二次型 f=xTAx 中,A实对称矩阵p 为正惯性指数(正平方项个数),q 为负惯性指数(负平方项个数),则 r(A)=p+q

【注】必须是在实对称矩阵的条件下!

【一类特殊的二次型】已知二次型

f(x1,x2,x3)=(a1x1+a2x2+a3x3)2+(b1x1+b2x2+b3x3)2+(c1x1+c2x2+c3x3)2

α=(a1,a2,a3)T,β=(b1,b2,b3)T,γ=(c1,c2,c3)T,则二次型 f(x1,x2,x3) 对应矩阵为 A=(α,β,γ)(α,β,γ)T,正惯性指数 p=r(α,β,γ),负惯性指数 q=0

(3)最大和最小值

n 元二次型 f=xTAx,其中实对称矩阵 A 的特征值 λ1,λ2,...,λn 中最大值为 λmax,最小值为 λmin,且 xTAx=M>0,则有:

MλminxTAxMλmax

(4)二次型的标准化(合同对角化)

(4.1)正交变换法

  • |λEA|=0 求出二次型矩阵 A 的特征值 λ1,λ2,...,λn
  • 对每个 λi,由 (λiEA)x=0 求出 A 的一组特征向量 α1,α2,...,αn
  • 对单重特征值的特征向量进行单位化;对多重特征值对应的特征向量进行施密特正交化单位化
  • 令正交矩阵 Q=(η1,η2,...,ηn),有可逆线性变量替换 x=Qy,把原二次型化为标准二次型,A 的特征值 λ1,λ2,...,λn 对应标准二次型中每一项的系数

(4.2)拉格朗日配方法

  • 若二次型中有平方项 xi2 和交叉项 xixj,则把含有 xi 的项集中起来进行配方
  • 若二次型中仅有交叉项 xixj,则进行以下换元,此时将产生出平方项,按第一种方法进行配方:

{xi=yi+yjxj=yiyjxk=yk(k=1,2,...,n)(ki,j)

  • 也可使用公式 ab=(a+b)2(ab)24 产生出平方项

【例】用配方法将二次型 f(x1,x2,x3)=x1x2+x2x3+x1x3 化为标准型。

【解】根据交叉项 x1x2 可进行以下换元(当然也可以挑选其他交叉项进行换元):

{x1=y1+y2x2=y1y2x3=y3

所以 f(x1,x2,x3)=y12+y22+2y1y3,配方得 f(x1,x2,x3)=(y1+y3)2y22y32

2. 合同关系

(1)合同的定义

n 阶矩阵 A,B,若存在可逆实矩阵 C,使得 B=CTAC,则称 AB 合同,记为 AB

【注】在矩阵合同的定义中,并没有要求合同的矩阵一定是实对称矩阵。

(2)合同的性质

(2.1)一般矩阵的性质(设 A,B 为一般矩阵):

  • 两个二次型(分别对应实对称矩阵 A,B)可用可逆线性变量替换互相转化 AB
  • AB 正、负惯性指数相同,即 pA=pB,qA=qB
  • AB 正、负特征值个数相同
  • ABr(A)=r(B)

【注】ABpA=pB,qA=qBr(A)=r(B),但 r(A)=r(B)AB,只能推出 AB 等价(若 A,B 同型)。

  • ABAB

【注】这是在一般矩阵下的得出的结论:

  • 合同不一定相似:很容易理解,合同只能推出矩阵 A,B 所对应的对角矩阵元素的正负个数相等,但无法推出对角矩阵元素均相等。
  • 相似不一定合同:由 AB 可得 P1AP=B,但无法保证 P1=PT

(2.2)实对称矩阵的性质(设 A,B 为实对称矩阵):

  • 实对称矩阵必能合同对角化,即 CTAC=Λ
  • A 为实对称矩阵,则:ABB 为实对称矩阵

【证明】ABCTAC=B(CTAC)T=(B)TCTATC=BT=B,说明 B 也为实对称矩阵。

  • ABAB,但 ABAB

【注】这是在实对称矩阵下的得出的结论:

  • 相似是合同的特例:实对称矩阵必与对角矩阵相似,可得 ABABΛ,所以 A,B 有相同的特征值,即 A,B 有相同的正、负惯性指数,由惯性定理知 AB
  • 合同不一定相似:很容易理解,合同只能推出矩阵 A,B 所对应的对角矩阵元素的正负个数相等,但无法推出对角矩阵元素均相等。

3. 正定矩阵

(1)正定的定义

设二次型 f(x1,x2,...,xn)=xTAx,其中 A 是实对称矩阵。若对任意 x0,都有 f(x1,x2,...,xn)=xTAx>0,则称二次型 f 正定,称 A正定矩阵

【注 1】判定矩阵 A 正定时,需要检验 A 是否为实对称矩阵。
【注 2】若二次型 f 正定,则仅当 x=0 时,f(x1,x2,...,xn)=xTAx=0

(2)正定的性质

  • A 正定 A 为实对称矩阵
  • A 正定 A 的特征值全部大于 0
  • A 正定 A 的顺序主子式全大于 0
  • A 正定 P,使得 A=PTPAE
  • A 正定 aii>0
  • A 正定 |A|>0
  • A 正定 Ak,A1,A 均正定
  • [AOOB] 正定 A,B 均正定
  • 对于实矩阵 ATA
    • ATA 的负惯性指数为 0
    • r(A)=n,则 ATA 正定

【证明】(1)首先证明矩阵 ATA 为实对称矩阵。因为 (ATA)T=AT(AT)T=ATA,所以 ATA 为实对称矩阵。

(2)现在用特征值证明其正定。设 λ 是矩阵 ATA 的特征值,所对应的特征向量为 α,则有:ATAα=λα,等式两边同乘 αT 得:αTATAα=λαTα,化为内积形式即:(Aα,Aα)=λ(α,α),显然 λ0,矩阵 ATA 负惯性指数为 0

(3)当 r(A)=n 时,表示 Ax=0 仅有非零解,所以 Aα0(Aα,Aα)>0,显然 λ>0,矩阵 ATA 正定。

三、等价关系

1. 等价的定义

若矩阵 A 经过有限次初等变换变成矩阵 B,则称矩阵 A 与矩阵 B 等价,记为 AB

2. 等价的判定

A,B 是同型矩阵,则:

  • ABA 经过初等变换得到 B
  • ABPAQ=B,其中 P,Q 可逆
  • ABr(A)=r(B)

本文作者:漫舞八月(Mount256)

本文链接:https://www.cnblogs.com/Mount256/p/17633705.html

版权声明:本作品采用CC 4.0 BY-SA许可协议进行许可。

posted @   漫舞八月(Mount256)  阅读(347)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· DeepSeek 开源周回顾「GitHub 热点速览」
· 记一次.NET内存居高不下排查解决与启示
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· .NET10 - 预览版1新功能体验(一)
历史上的今天:
2022-08-16 【组成原理-数据】浮点数的编码与运算
2021-08-16 STM32学习笔记(9)——(I2C续)读写EEPROM
点击右上角即可分享
微信分享提示
评论
收藏
关注
推荐
深色
回顶
展开
  1. 1 Main Menu Theme Syd Matters
  2. 2 Luminous Memory (Acyanxi Remix) Acyanxi
  3. 3 夏影 麻枝准
  4. 4 潮騒の香り 水月陵
  5. 5 stand still 井口裕香 (いぐち ゆか)
  6. 6 流星雨 麻枝准
  7. 7 Summer Fantasy 傅许
  8. 8 失う 米白
  9. 9 epilogue 霜月はるか
  10. 10 夏に君を待ちながら 小原好美
  11. 11 桜のような恋でした 鹿乃 (かの)
  12. 12 風は微かに、熱を残し… 水月陵
  13. 13 夏凪ぎ 麻枝准/やなぎなぎ
  14. 14 空に光る 戸越まごめ
  15. 15 木漏れ日 riya
  16. 16 Songbirds Homecomings (ホームカミングス)
  17. 17 宝物になった日 麻枝准/やなぎなぎ
  18. 18 夏影~あの飛行機雲を超えた、その先へ~ 雪桜草 (雪樱草)
  19. 19 快晴 Orangestar (蜜柑星P),初音未来 (初音ミク)
  20. 20 永遠 霜月はるか
  21. 21 Sion 天門
  22. 22 遙かな年月-piano- 麻枝准
  23. 23 夏恋慕 kobasolo/春茶
  24. 24 夏凪ぎ-piano ver.- MANYO/麻枝准
  25. 25 Goodbye Seven Seas -piano ver.- MANYO/麻枝准
  26. 26 Light Years 麻枝准/やなぎなぎ
  27. 27 優しさの記憶 鹿乃 (かの)
夏影 - 麻枝准
00:00 / 00:00
An audio error has occurred, player will skip forward in 2 seconds.