隐藏页面特效

相似对角化与合同对角化的关系

本博客仅做个人经验性总结,非学习笔记,是在学习完两者后做的一个区分。

1|0相似矩阵


1|1定义


在一组基 C=(c1,c2,,cn) 下,向量 s=(s1,s2,,sn)T 具有坐标 x=(x1,x2,,xn)T,此时 s=Cx

考虑一个变换矩阵 An×n,使得 ss 在基 C 下关于 A 的映射,即 s=As,我们也有 s=Cx,s=Cx,我们可以把 An×n 直接作用于坐标上,即 s=CAx=Cx

现在我们有一组新基 D=(d1,d2,,dn),我们希望求得在基 Ds 的坐标 y=(y1,y2,,yn)T,这很容易,因为 s=Cx=Dy,我们解一个非齐次线性方程组即可得出 y

我们知道 s 在基 C 下从 s 映射而来的变换矩阵 A,如何求出在基 D 下从 s 映射而来的变换矩阵 B 呢?

P 为基变换矩阵,根据基变换公式:

D=CP

x=Py

容易求出 P=C1D,进而得出 y=Px,带入 s=DBy=CAx 有:

CAx=DP1APy

所以我们可以得出重要结论:

B=P1AP

因此 AB 描述了 n 维线性空间上同一个线性变换,这种变换是客观存在的,但可以使用不同的基从不同的矩阵进行描述,在不同基下描述同一个线性变换的矩阵称作相似矩阵,即 AB,上述表达式 B=P1AP 就是相似矩阵的表达式。

已知某线性变换在基 C=(c1,c2,,cn) 下的矩阵是 A,求该线性变换在基 D 下的矩阵 B

我们首先求出 CD 的基变换矩阵 P

根据 D=CPP=C1D,使用公式 B=P1AP 即可得到相似矩阵 B

1|2性质


  • 反身性:AA

  • 对称性:若 AB,则 BA

  • 传递性:若 ABBC,则 AC

  • 单位阵相似唯一性:单位阵 E 只与自身相似

  • AB,则两者具有相同的特征多项式。

    • 证明:
      B=P1AP|λEB|=|λEP1AP|
      |λEP1AP|=|P1(λE)PP1AP|=|P1(λEA)P|=|P1||λEA||P|=|λEA|

      这进而推出了更多的性质,此时 AB 的特征多项式各个系数均相等,对比对应两者 λn1 的系数和常数项,可知:

      tr(A)=tr(B)

      |A|=|B|

      更进一步,两者特征方程相同,特征根相同, 有:

      λA=λB

  • AB,如果一方可逆,另一方一定可逆

  • AB,则 AB

    • 证明:
      B=(P1AP)=PA(P1)
      P=P1|P|
      (P1|P|)A(P|P1|)=P1AP

      容易发现伴随矩阵的基变换矩阵也是 P,两者伴随矩阵相似。

  • AB,则 ATBT

    • 证明:
      B=P1APBT=(P1AP)T=PTAT(P1)T=PTAT(PT)1

      该矩阵的基变换矩阵为 PT,两者转置矩阵也相似。

  • AB,则 kAkBAmBmp(A)p(B)

    • 证明较为显然,故略过
  • AB,则 r(B)=r(A)
    因为基变换矩阵相当于行和列的线性运算,不改变其空间大小,所以他们的秩是相同的。

我们列出一个表格来总结一些性质:

矩阵 A AT A1 A p(A) 行向量 特征根和特征向量
基变换矩阵 P (PT)1 P1 P P - - -
性质 tr(A)=tr(B) |A|=|B| λA=λB,XA=XB

2|0相似对角化


2|1定义


n 阶方阵 A 经过 n 阶可逆矩阵 X 的相似变换成为对角矩阵 U=X1AX 这一过程称作相似对角化

因为对角矩阵的优秀性质,他的线性运算相当于各个坐标关于基的线性运算,仅仅具有倍数关系,我们令 U=diag(λ1,λ2,,λn),X=(x1,x2,,xn)xi 为列向量。

可以看出

AX=XUAxi=xiλi

如果 A 通过可逆矩阵 X 的相似变换成为了对角阵 U,那么 U 的各个元素就是的 A 的特征值,X 的各列就是 A 的特征向量。

这样,对角化问题转化为了求特征值和特征向量问题。但是 AU 未必相似,因为 X 未必可逆,也就是说,相似特征化的前提条件是 X 可逆

3|0合同变换


3|1定义


对二次型 f=xTAx,其中 A 为实对称矩阵(因为实对称矩阵是唯一的)。

二次型的秩定义为 r(f)=r(A)

对两组二次型做基变换,如下图所示:

二次型 基向量组 坐标
f=xTAx C=(c1,c2,,cn) x
f=yTBy D=(d1,d2,,dn) y

设基向量组 CD 的基变换矩阵为 P,则 P 可逆且 x=Py

f=xTAx=(Py)TA(Py)=yT(PTAP)y

我们称 B=PTAP合同变换

不同基向量组下表示同一个二次型的实对称矩阵是一组合同矩阵,其中对应的合同变换矩阵是对应的基变换矩阵,它一定可逆

对于 n 维线性空间下同一个二次型不同基向量下,对应的实对称矩阵具有合同关系

4|0实对称矩阵


4|1定义


如果方阵 A 满足 A=AT,则称 A 为对称矩阵。

如果一个对称矩阵内的所有元素都是实数,那么称其为实对称矩阵,也就是 A=AT=¯A

4|2性质


复数运算可以扩展到矩阵:

¯z1+z2=¯z1+¯z2¯λA=¯λ¯A

¯z1z2=¯z1¯z2¯AB=¯A¯B

  • 实对称矩阵的特征值都是实数
    实对称矩阵的特征向量都是实向量

    • 证明:

      根据 A=AT=¯AAx=λx(x0),我们推导下列式子,因为点乘运算等价于转置乘原行向量:

    λ(¯x·x)=λ(¯xTx)=¯xT(λx)=¯xT(Ax)=¯xTATx=(A¯x)Tx=(¯A¯x)Tx=(¯Ax)Tx=(¯λx)Tx=(¯λ¯x)Tx=¯λ(¯xTx)=¯λ(¯x·x)

    这说明 ¯λ=λ,即 λR

  • 实对称矩阵的特征向量两两正交

    • 证明:
      考虑 λ1x1=Ax1λ2x2=Ax2(λ1λ2),有:

      λ1(x1·x2)=λ1(xT1x2)=(λ1x1)Tx2=(Ax1)Tx2=xT1ATx2=xT1(Ax2)=xT1(λ2x2)=λ2xT1x2=λ2(x1·x2)

      (λ1λ2)(x1·x2)=0x1·x2=0

      因此实对称矩阵的特征向量两两正交,也说明了实对称矩阵的几何重数等于代数重数,实对称矩阵一定可以相似对角化,这个性质比线性无关更强,下面我们会用到。

5|0合同对角化


选择合适的可逆矩阵 P,使 A 经过合同变换后变为对角阵 U=PTAP,这样就可以使某个一般形态的二次型化为只含有平方项的特殊二次型,这个过程叫合同对角化

如下列案例:

f=x21+3x22+2x1x2=(x1,x2)[1113][x1x2]=xTAx

=(x21+2x1x2+x22)+2x22=(x1+x2)2+2x22

此时取 y1=x1+x2y2=x2 可以得到:

f=(y1,y2)[1002][y1y2]=yTUy

使用配方法可以将任意一个二次型化为只含平方项的形式,这种形式叫做二次型的标准型,即任意二次型都可化为标准型。

利用基变换矩阵:

[y1y2]=[1101][x1x2]

[x1x2]=[1101][y1y2]

根据基变换公式 x=Py,利用合同变换可以得到 U=PTAP

如果我们令 f=1,容易发现,上述两个方程均代表了一个椭圆方程:

f=x21+3x22+2x1x2=1

image

f=y21+2y22=1

image

我们成功将一个斜椭圆转换为了正椭圆,因此可以通过合同变换使得圆锥曲线的长轴和短轴与坐标轴共线,从而获得标准方程以简化研究。

但是我们发现刚才的变化有一个显著缺陷,也就是在转换之后,椭圆的形状改变了,这破坏了他原有的性质。

image

合同变换时,我们需要保证变换前后椭圆是全等状态,即椭圆上的每一个点都单纯发生旋转变换,我们需要从矩阵的角度去探讨这个问题。

设二次型 f=xTAx,可逆合同变换矩阵 Q 使得 U=QTAQf=yTUy 成立,需要 U 是对角矩阵且变换前后任意一点到原点的距离不变,可以得出 x=yxTx=yTy,带入 x=Qy 有:

(Qy)TQy=yTyyT(QTQ)y=yTy

这要求 QTQ=E,我们称满足这样的 Q正交矩阵

因此我们想要维持原有二次型的性质,我们要求他的基变换矩阵必须要是正交矩阵。

5|1正交矩阵的性质


因为 QTQ=E|QT|=|Q|,可以推出 |Q|2=1,故正交矩阵的行列式只能为 11

更优秀的,因为 Q1Q=E,所以 Q1=QT,正交矩阵的逆矩阵与转置矩阵相等。

Q 分块为 n 组列向量 Q=(q1,,qn),我们有:

QTQ=[qT1q1qT1qnqTnq1qTnqn]=[1001]

这意味正交矩阵任意一列(行)向量都是单位向量且两两正交,全体列(行)向量可以作为一组单位正交基向量组

二次型 f=xTAx 中,使用正交矩阵 Q 对实对称矩阵 A 经过合同变换化为对角阵 U,则有:

U=QTAQ=Q1AQ

由于实对称矩阵 A 到他的特征值矩阵 U合同相似的,它必有正交的特征向量组,我们很多时候不需要求出可逆的变换矩阵 Q,因为这个矩阵是客观存在的。

通过这种方式化成的标准型不仅保留原先二次型的类型,还使得二次型得以全等地变换,并保持各个方向的伸缩因子不变,更重要的是建立起了二次型合同对角化实对称矩阵的相似对角化的关系,这样就可以使用对实对称矩阵相似对角化的方法题来解决化二次型为标准型的问题。

5|2与曲线的关系


不妨考虑上面 x21+2x1x2+3x2=1 的性质,他的二次型可以表示为:

(x1,x2)[1113][x1x2]

我们可以求解它的特征值:

|1λ113λ|=λ2+4λ+2=0{λ1=22λ2=2+2

所以经过合同变换后的二次型变为:

(y1,y2)[λ100λ2][y1y2]

也就是 (22)y21+(2+2)y22=1,我们现在可以直观地求出这个椭圆的长短轴和面积,无疑简化了我们的运算。

更一般的,由于 ax21+bx1x2+cx22=1 可以转换为 λ1y21+λ2y22=1,由于两个二次型矩阵的是相似的,我们可以得到:

λ1+λ2=a+c

λ1λ2=14(b24ac)

Δ=b24ac,我们甚至可以直接看出二次曲线的类型:

  • a+c>0Δ<0 时,λ 均为正数,对应曲线为椭圆。

  • a+c<0Δ<0 时,对应曲线在 f=1 时无意义。

  • Δ>0 时,λ 一正一负,对应曲线为双曲线。

  • Δ=0,即 |A|=0 时,必有特征值为 0,所以原二次型可写作 (a+c)y21=1,为两条平行直线。

更平凡的,如果为三元二次型,我们也能得到类似的结论,我们列出一个表格,但不再继续说明。

f(x1,x2,x3)=xTAX=m(m>0)

特征值 曲面类型
三正 椭球面
三负 无意义
两正一负 单叶单曲面
两负一正 双叶双曲面
两正一零 椭圆柱面
两负一零 无意义
一正一负一零 双曲柱面
一非零,两零 两平面平行

__EOF__

本文作者Yip.Chip
本文链接https://www.cnblogs.com/YipChipqwq/p/18617791.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   YipChip  阅读(236)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 无需6万激活码!GitHub神秘组织3小时极速复刻Manus,手把手教你使用OpenManus搭建本
· C#/.NET/.NET Core优秀项目和框架2025年2月简报
· Manus爆火,是硬核还是营销?
· 终于写完轮子一部分:tcp代理 了,记录一下
· 【杭电多校比赛记录】2025“钉耙编程”中国大学生算法设计春季联赛(1)
点击右上角即可分享
微信分享提示