29-相似矩阵和若尔当形

一、接着上一节说正定矩阵

 所谓正定,就是$x^TAx > 0$($except \space for \space x = 0$)成立,我们通常也可以通过特征值,主元,行列式来判断

 虽然我们知道了什么是正定矩阵,如何判断正定矩阵,那么正定矩阵是从何而来的呢?主要来自:最小二乘法

 实际上,大量的物理问题需要用长方形矩阵来描述,我们知道最小二乘法的关键是矩阵:$A^TA$,我们希望证明这是正定矩阵

 如果我们知道矩阵$A, B$都是正定的,那么矩阵$A+B$是否是正定的呢?我们只需要证明:$x^T(A+B)x > 0$,我们知道:

$x^TAx > 0$

$x^TBx > 0$

两式相加,即可证明$x^T(A+B)x > 0$,因此矩阵$A+B$正定

 下面我们来证明最小二乘法的关键矩阵是否正定:假设矩阵$A_{m*n}$(不是对称矩阵,所以不是正定的)但我们知道:$A^TA$是方阵,并且对称,我们只需要证明下面的式子恒大于0:

$x^T(A^TA)x > 0$

证明:$x^T(A^TA)x = (Ax)^T(Ax) = |Ax|^2 >= 0$,矩阵乘以向量得到向量,向量的平方相当于求向量的模长,只要$x$不等于0,那么上面的式子恒大于0

二、相似矩阵

 假设矩阵为:$A_{n*n}, B_{n*n}$,不要再把他当做对称矩阵,就是普通的方阵,$A, B$相似的意思是:存在某个可逆矩阵$M^{-1}$,使得下式成立

$B = M^{-1}AM$,可能你会觉得该式子来得比较突然,这样组合有什么意义

 其实上面的式子,你应该不陌生,还记得下面的公式吧,矩阵对角化

$A = S\Lambda S^{-1}, \Lambda = S^{-1}AS$

这说明$\Lambda$和$A$相似,$S$就是上面的$M$,但是如果$M$不取$S$,用其他矩阵的话,$M^{-1}AM$的结果就不是对角化矩阵$\Lambda$,而是新的相似于$A$的矩阵

其实改变$M$,只要可逆,就会得到很多与$A$相似的矩阵,他们和$\Lambda$一样与$A$相似,这些与$A$相似的矩阵应该归为一类,只是$\Lambda$是这些相似矩阵中最为简洁的一个:因为其为对角化矩阵

 我们为什么会对相似矩阵感兴趣呢?因为他们有共同的性质-他们的特征值都相同,我们来举例

$A = \left[\begin{array}{ll}{2} & {1} \\ {1} & {2}\end{array}\right]$,其特征值为$3, 1$

其对角化矩阵为:$\Lambda = \left[\begin{array}{ll}{3} & {0} \\ {0} & {1}\end{array}\right]$,其特征值为$3, 1$

我们再随便取一个$M = \left[\begin{array}{ll}{1} & {4} \\ {0} & {1}\end{array}\right]$

$M$的逆矩阵为:$M^{-1} = \left[\begin{array}{ll}{1} & {-4} \\ {0} & {1}\end{array}\right]$

$B = M^{-1}AM = \left[\begin{array}{ll}{1} & {-4} \\ {0} & {1}\end{array}\right]\left[\begin{array}{ll}{2} & {1} \\ {1} & {2}\end{array}\right]\left[\begin{array}{ll}{1} & {4} \\ {0} & {1}\end{array}\right]=\left[\begin{array}{ll}{-2} & {-15} \\ {1} & {6}\end{array}\right]$,其特征值为$3, 1$

再取其他$M$,会得到新的矩阵,和$A$相似,这些相似矩阵都有相同的特征值$3, 1$

 那么为什么相似矩阵会有相同的特征值呢?我们来证明一下

  首先,$Ax = \lambda x$,我们知道:$B = M^{-1}AM$,我们把$MM^{-1}$添加到$Ax$中间:

$AMM^{-1}x = \lambda x$

等式两侧同乘$M^{-1}$:$M^{-1}AMM^{-1}x = M^{-1} \lambda x$,$B$出现啦

$BM^{-1}x = \lambda M^{-1}x$,我们把$M^{-1}x$当作新的向量,并命名为$z$,则

$Bz = \lambda z$,所以$\lambda$也是$B$的特征值,同时我们知道:$M^{-1}$乘以$A$的特征向量是$B$的特征向量

注意:相似矩阵的特征值相同,但特征向量不同,比如一个是$x$,一个是$z$,即$M^{-1}x$,特征向量不会相同的,如果特征值和特征向量都相同,那么不叫相似矩阵了,他们就是相同的矩阵

三、特殊情况

 上面二中所讲的例子,正好是两个特征值不相等(一个是3,一个是1)的情况,也就是原矩阵$n$个特征向量线性无关,可以对角化。

 但是我们不能排除如果特征值相等(也就是特征向量存在相关性),特征值相等的矩阵求相似矩阵,还可以分为两种情况,我们还是以$2*2$矩阵为例:

 1)相似矩阵只有原矩阵一个(孤零零的一个人),也可以说不存在相似矩阵,因为只有自己一个,如:

$A = \left[\begin{array}{ll}{4} & {0} \\ {0} & {4}\end{array}\right]$,其特征值为$4, 4$

我们发现:$B = M^{-1}\left[\begin{array}{ll}{4} & {0} \\ {0} & {4}\end{array}\right]M=M^{-1}4IM=4I=A$

无论$M$取何矩阵,最后与$A$相似的矩阵只有原矩阵自己一个

 2)另外一种特征值相同,不可以对角化,如:

$\left[\begin{array}{ll}{4} & {1\space or \space something \space other} \\ {0} & {4}\end{array}\right]$

  我们发现上面的矩阵特征值相同,均为$4$,该家族中最好的矩阵是右上角那个值是1,,而这种形式被称为若尔当标准型(Jordan Form),这是该家族中最简洁的,最接近对角阵的一个。也就是说虽然原矩阵不可以对角化,但是可以找到最接近对角化矩阵的一个

  我们来看看这个家族(相似矩阵)的几个成员:

$\left[\begin{array}{ll}{4} & {1} \\ {0} & {4}\end{array}\right],\left[\begin{array}{ll}{5} & {1} \\ {-1} & {3}\end{array}\right],\left[\begin{array}{ll}{4} & {0} \\ {17} & {4}\end{array}\right],......$

是不是第一个最接近对角化矩阵呢

这些矩阵都是相似的,只要找到合适的$M$,都可以证明一个矩阵相似于另一个

但要注意:一般矩阵很难化成若尔当标准型,因为若尔当标准型依赖于特征值严格相等

 3)下面还有一点没看懂,有空再补充

  请参考

原文地址:https://www.cnblogs.com/always-fight/p/12339616.html

时间: 2024-10-13 00:56:34

29-相似矩阵和若尔当形的相关文章

线性代数导论35——线性代数全总结(麻省理工公开课:线性代数)

课程介绍 "线性代数",同微积分一样,是高等数学中两大入门课程之一,不仅是一门非常好的数学课程,也是一门非常好的工具学科,在很多领域都有广泛的用途.本课程讲述了矩阵理论及线性代数的基本知识,侧重于那些与其他学科相关的内容,包括方程组.向量空间.行列式.特征值.相似矩阵及正定矩阵. [第1集] 方程组的几何解释    [第2集] 矩阵消元    [第3集] 乘法和逆矩阵    [第4集] A的LU分解    [第5集] 转置-置换-向量空间R    [第6集] 列空间和零空间    [第

BadVPN详解之--组网原理剖析

这可能是第一篇关于BadVPN原理的中文介绍,所以我恪守职责,希望能将BadVPN的原理阐述清楚,而不仅仅是为了欺世盗名取天下之先. 序:格尔上市 这一段插在这里有点突兀,与BadVPN无关,只是我的一点实感,忽略即可.       明天小小要去厦门演出,今天是什么日子?是小小准备且出发的日子.但是同样在今天,也就是2017年4月21日(写完本文时估计已经到了22日...),上海格尔软件股份有限公司成功上市了.也许你不知道这个公司,确实这公司的业务是很小众,不像BAT布局那般宏大,但这并不意味着

暗黑符文之语1.10

公式要求的物品类型 确定你使用了正确的物品类型. 衣服 = 任何衣服部位都有效, 但是不包括例如头盔等其他部位, 在英文版中这个非常容易有误解. 头盔 = 任何头部物品都有效, 包括一般头盔, 野蛮人专用头盔, 德鲁依专用皮盔以及头饰. 盾牌 = 任何盾牌类型都有效, 包括盾牌, 死灵法师萎缩头颅, 以及圣骑士盾牌. 棍棒 = 木棒, 狼牙棒, 以及扩展/精华类型. 锤类 = 巨战铁锤, 大木棍, 卓越巨棍, 以及扩展/精华类型. 它和钉头锤是不同的. 钉头锤 = 钉头锤, 流星锤, 连枷, 以

Power Map

推荐64位版本Office,但会遇到以下问题 建议大多数用户使用 32 位 Office 我们建议使用 32 版本的 Office,因为它与大多数其他应用程序更加兼容,尤其是第三方加载项.这也是默认安装 32 位版本的 Office 2013 的原因所在,即使在 64 位 Windows 操作系统上也是如此.在这些系统上,支持将 32 位 Office 客户端作为 Windows-32-on-Windows-64 (WOW64) 安装.WOW64 是 x86 模拟器,它使基于 Windows 的

若尔当标准形 弗罗贝尼乌斯标准形

["两个相似的矩阵可以看做是同一个线性变换的'两面',即在两个不同的基下的表现":故特征向量不一定相同外,其他都相同, 其中"其他"指矩阵的其他描述量,如行列式的值.秩.迹数.特征值.特征多项式.初等因子:根据逆否命题的成立关系,可出 2个矩阵为相似矩阵的必要条件.] https://zh.wikipedia.org/wiki/相似矩陣 两个相似的矩阵有许多相同的性质: 两者的秩相等. 两者的行列式值相等. 两者的迹数相等. 两者拥有同样的特征值,尽管相应的特征向量

一文搞懂HMM(隐马尔可夫模型)

本文转自于:http://www.cnblogs.com/skyme/p/4651331.html 隐马尔可夫模型(Hidden Markov Model,HMM)是统计模型,它用来描述一个含有隐含未知参数的马尔可夫过程.其难点是从可观察的参数中确定该过程的隐含参数.然后利用这些参数来作进一步的分析,例如模式识别. 是在被建模的系统被认为是一个马尔可夫过程与未观测到的(隐藏的)的状态的统计马尔可夫模型. 下面用一个简单的例子来阐述: 假设我手里有三个不同的骰子.第一个骰子是我们平常见的骰子(称这

[官方软件] Easy Sysprep v4.3.29.602 【系统封装部署利器】(2016.01.22)--skyfree大神

[官方软件] Easy Sysprep v4.3.29.602 [系统封装部署利器](2016.01.22) Skyfree 发表于 2016-1-22 13:55:55 https://www.itsk.com/forum.php?mod=viewthread&tid=362766&highlight=Easy%2BSysprep [官方软件] Easy Sysprep v4.3.29.602 [系统封装部署利器](2016.01.22) [Easy Sysprep]概述:Easy Sy

隐马尔可夫模型的前向算法(java实现),今天奉上

隐马尔可夫模型的前向算法(手动实现),今天奉上,由于研究生期间,实现的时候没有多加注释,这里为了让更好的人进入自然语言处理领域,特此,将前向算法奉上,具体公式可参考52nlp的HMN系列博客. 参考了大部分网站公式和借鉴.在此表示感谢. 后向算法和维特比算法,后续更新. HMM类: 1 package jxutcm.edu.cn.hmm.model;  2   3 import jxutcm.edu.cn.hmm.bean.HMMHelper;  4   5 /**  6  * 实现了 HMM(

2月29日(编程之美2015资格赛)

时间限制:2000ms 单点时限:1000ms 内存限制:256MB 描述 给定两个日期,计算这两个日期之间有多少个2月29日(包括起始日期). 只有闰年有2月29日,满足以下一个条件的年份为闰年: 1. 年份能被4整除但不能被100整除 2. 年份能被400整除 输入 第一行为一个整数T,表示数据组数. 之后每组数据包含两行.每一行格式为"month day, year",表示一个日期.month为{"January", "February",