相似对角化与合同对角化

news/2024/12/19 18:58:41/文章来源:https://www.cnblogs.com/YipChipqwq/p/18617791

本博客仅做个人笔记。

相似矩阵

定义

在一组基 \(e = (e_1, e_2, \ldots, e_n)\) 下,向量 \(s = (s_1, s_2, \ldots, s_n)^T\) 具有坐标 \(x = (x_1, x_2, \ldots, x_n)^T\),此时 \(s = ex\)

考虑一个变换矩阵 \(A_{n \times n}\),使得 \(s'\)\(s\) 在基 \(e\) 下关于 \(A\) 的映射,即 \(s' = As\),我们也有 \(s = ex, s' = ex'\),我们可以把 \(A_{n \times n}\) 直接作用于坐标上,即 \(s' = eAx = ex'\)

现在我们有一组新基 \(e' = (e_1', e_2', \ldots, e_n')\),我们希望求得在 \(e'\)\(s\) 的坐标 \(y = (y_1, y_2, \ldots, y_n)^T\),这很容易,因为 \(s = ex = e'y\),我们解一个非齐次线性方程组即可得出 \(y\)

我们知道 \(s'\) 在基 \(e\) 下从 \(s\) 映射而来的变换矩阵 \(A\),如何求出在基 \(e'\) 下从 \(s\) 映射而来的变换矩阵 \(B\) 呢?

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

\[e' = eP \]

\[x = Py \]

容易求出 \(P = e^{-1}e'\),进而得出 \(y = Px\),带入 \(s' = e'By = eAx\) 有:

\[eAx = e'P^{-1}APy \]

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

\[B = P^{-1}AP \]

因此 \(A\)\(B\) 描述了 \(n\) 维线性空间上同一个线性变换,这种变换是客观存在的,但可以使用不同的基从不同的矩阵进行描述,在不同基下描述同一个线性变换的矩阵称作相似矩阵,即 \(A \sim B\),上述表达式 \(B = P^{-1}AP\) 就是相似矩阵的表达式。

已知某线性变换在基 \(C = (c_1, c_2, \ldots, c_n)\) 下的矩阵是 \(A\),求该线性变换在基 \(D\) 下的矩阵 \(B\)

我们首先求出 \(C\)\(D\) 的基变换矩阵 \(P\)

根据 \(D = CP \to P = C^{-1}D\),使用公式 \(B = P^{-1}AP\) 即可得到相似矩阵 \(B\)

性质

  • 反身性:\(A \sim A\)

  • 对称性:若 \(A \sim B\),则 \(B \sim A\)

  • 传递性:若 \(A \sim B\)\(B \sim C\),则 \(A \sim C\)

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

  • \(A \sim B\),则两者具有相同的特征多项式。

    • 证明:
      \(B = P^{-1}AP \to \lvert \lambda E - B \rvert = \lvert \lambda E - P^{-1}AP \rvert\)
      \(\begin{aligned} \lvert \lambda E - P^{-1}AP \rvert &= \lvert P^{-1} (\lambda E)P - P^{-1}AP \rvert \\ &= \lvert P^{-1} (\lambda E - A)P \rvert\\ &= \lvert P^{-1} \rvert \lvert \lambda E - A \rvert \lvert P \rvert \\ &= \lvert \lambda E - A \rvert \end{aligned}\)

    这进而推出了更多的性质,此时 \(A\)\(B\) 的特征多项式各个系数均相等,对比对应两者 \(\lambda^{n - 1}\) 的系数和常数项,可知:

    \[tr(A) = tr(B) \]

    \[\lvert A \rvert = \lvert B \rvert \]

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

    \[\lambda_A = \lambda_B \]

  • \(A \sim B\),如果一方可逆,另一方一定可逆

  • \(A \sim B\),则 \(A^{*} \sim B^{*}\)

    • 证明:
      \(B^{*} = (P^{-1} A P)^{*} = P^{*} A^{*} (P^{-1})^{*}\)
      \(P^{*} = P^{-1} \lvert P \rvert\)
      \(\to (P^{-1} \lvert P \rvert) A^{*} (P \lvert P^{-1} \rvert) = P^{-1} A^{*} P\)

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

  • \(A \sim B\),则 \(A^T \sim B^T\)

    • 证明:
      \(\begin{aligned} B = P^{-1}AP \to B^T &= (P^{-1}AP)^T \\ &= P^T A^T (P^{-1})T \\ &= P^T A^T (P^T)^{-1} \\ \end{aligned}\)

    该矩阵的基变换矩阵为 \(P^T\),两者转置矩阵也相似。

  • \(A \sim B\),则 \(kA \sim kB\)\(A^m \sim B^m\)\(p(A) \sim p(B)\)

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

相似对角化

定义

\(n\) 阶方阵 \(A\) 经过 \(n\) 阶可逆矩阵 \(X\) 的相似变换成为对角矩阵
\(U = X^{-1}AX\) 这一过程称作相似对角化

因为对角矩阵的优秀性质,他的线性运算相当于各个坐标关于基的线性运算,仅仅具有倍数关系,我们令 \(U = \operatorname{diag}(\lambda_1, \lambda_2, \ldots, \lambda_n), X = (x_1, x_2, \ldots, x_n)\)\(x_i\) 为列向量。

可以看出

\[AX = XU \to Ax_i = x_i\lambda_i \]

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

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

合同对角化

定义

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

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

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

二次型 基向量组 坐标
\(f = x^T A x\) \(C = (c_1, c_2, \ldots, c_n)\) \(x\)
\(f = y^T B y\) \(D = (d_1, d_2, \ldots, d_n)\) \(y\)

设基向量组 \(C\)\(D\) 的基变换矩阵为 \(P\),则 \(P\) 可逆且 \(x = Py\)

\[f = x^T A x = (Py)^T A (Py) = y^T (P^T A P) y \]

我们称 \(B = P^T A P\)合同变换

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/855526.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

记录一下:小华半导体HC32F448建立MDK工程

1.先到官网上下载文件 a>下载驱动库:HC32F448_DDL_Rev1.1.0.zip 驱动库中是包括了例程的。 b>下载样例:HC32F448_Template_Rev1.0.1.zip 可以直接复制官方的样例,就不用自己创建工程了。 c>下载芯片支持包:HC32F448_IDE_Rev1.0.1.zip 下载后双击安装即可。否则KEI…

4大应用场景揭秘:AI视频监控在养老院中的智能化管理与安全保障

随着人口老龄化的加剧,养老院的管理面临着越来越多的挑战。传统的人工巡查方式不仅难以做到全天候监控,而且存在响应迟缓、效率低下等问题。为了解决这些问题,思通数科推出的AI视频监控系统,利用人工智能技术提供了一种高效、智能化的解决方案。尤其在养老院的老人体征监控…

考勤管理系统

实现对用户信息的添加,查看,删除 实现对部门信息的添加,查看 实现用户登录功能,修改密码功能。 页面与之前的大相径庭。 还是一样的框架总体感觉换汤不换药,目前就是熟练度问题

浏览器怎么渲染数据的

突然发现自己对于css的样式规则一直都没有完全明白,今天写文好好整明白 浏览器渲染原理 1. 浏览器向服务器请求html文件 2. html文件返回浏览器 3. 浏览器解读html文件<!DOCTYPE html>//告诉浏览器,这是html5文件// html树<html lang="en">//语言<h…

在宝塔里添加反向代理

安装OnlyOffice时,需要配置反向代理; 完成以上操作之后,访问时有可能出现504报错,也可能就可以正常访问了。如果报错可以修改上图反向代理中的目标URL的内容将域名改成127.0.0.1。我当时改成http://127.0.0.1:9633就可以了本文来自博客园,作者:飞龙在生,转载请注明原文链…

2024 Clion安装使用教程(附激活以及常见问题处理)

第一步:下载Clion安装包 访问Clion官网,下载Clion第二步: 安装 Clion下载完成后,进行安装,next,安装完成点击xx 关掉程序! 第三步: 下载补丁 Clion补丁文件 点击获取补丁下载成功后,打开标注的文件文件夹 , 进入到文件夹 /jetbra 注意: 这个文件夹单独copy一份,所属文…

语言沟通中的设计实现

系统设计的目的是更好的支持需求 我们常说,只要业务能将你的需求描述清楚,能自圆其说,我们就有办法实现。 这其实是系统设计的最理想的状态, 如果业务没想清楚,那么在系统实现中,一定会把问题暴露出来。很多时候,问题的暴露源于没有考虑周全亦或都没有考虑这种场景! 比…

Python 解密 Navicat导出的数据库连接,Navicat数据库连接导入DBeaver。

最近公司收到Navicat律师告知书,让停止使用Navicat,用了那么久的数据库连接工具,不得不换其他的。 最终选择了开源的DBeaver。 安装完DBeaver后,把Navicat导出的connections.ncx文件直接导入DBeaver。直接访问提示连接失败,因为connections.ncx文件里的密码都是加密的。 如…

营销系统缺失投放概念

营销系统现阶段分为营销权益工具、玩法、招商提报三大块。其中权益工具主要负责创建券、促销、换购类优惠权益。玩法负责用户与平台之间的连接,通过任务式交互获得权益。招商提报主要是面向商家参与活动的连接。比如要做一个活动,需要从招商发布活动信息,提报统一创建营销权…

一款使用NET+MQTT+Arduino开发的智能浇花工具

最近闲来无事,对硬件控制产生了兴趣。看到家里的盆栽,我突然萌生了制作一个自动浇水工具的想法。通过在淘宝搜索并查找相关资料,我了解了需要的硬件和通信协议。接下来,我们先看看需要做哪些准备工作(如安装 Arduino、.NET、EMQX 工具等,请自行搜索并完成安装)。 准备工…

营销平台过去与展望

营销是什么 营销是商家给用户权益,让用户少花钱买到商品,以此达到商家的营销目的。 以此逻辑为基础,展开几个关键角色用户 平台 商家他们各自在营销侧的诉求是用户需要有获得优惠的途径 门店需要有发放优惠的能力 商家需要有优惠活动的功能营销工具权益建设 以这些诉求为切入…