实对称矩阵一定可以对角化
UTF8gbsn
实对称矩阵一定可以对角化.
最近看共轭梯度下降的时候看到有人的推导里面用到了这个命题. 虽然以前学过,
但是学得很渣, 所以没有自己想过这个命题怎么样成立的.
现在将这些证明过程梳理一下.
实对称矩阵含有n个实根
首先我们来证明一个命题, 实对称矩阵含有n个实根,
注意,n个实根并不一定都是不同的, 可能含有重根.
比如(r−1)2=0(r-1)^2=0(r−1)2=0就含有两个重根r=1r=1r=1.在计算根数目的时候这个方程的解算两个.
首先, 任意的矩阵A\mathbf{A}A,它的特征多项式
∣A−λI∣=0|\mathbf{A}-\lambda\mathbf{I}|=0∣A−λI∣=0
是一个nnn次多项式(这是很显然的).
由于nnn次多项式必定有nnn个根(在复数域上). 这个命题暂不证明,
直接使用. 我写过另外一篇文章简要的证明了一下这个定理.有了上一步的结论,
我们现在只需要证明每一个根λi\lambda_iλi是实根就可以了.
这个证明过程很简单. 假设λi\lambda_iλi是任意根之一,
并且αi\mathbf{\alpha}_iαi(当然也是在复数域),
那么根据特征值和特征向量的定义.我们可以得
Aαi=λiαi\mathbf{A}\mathbf{\alpha}_i=\lambda_i\mathbf{\alpha}_iAαi=λiαi 取共轭得
Aα‾i=λ‾iα‾i\mathbf{A}\mathbf{\overline{\alpha}}_i=\overline{\lambda}_i\mathbf{\overline{\alpha}}_iAαi=λiαi
再进行转置得, 注意AT=AA^T=AAT=A, 对称矩阵.
α‾iTA=λ‾iα‾iT\mathbf{\overline{\alpha}}_i^T\mathbf{A}=\overline{\lambda}_i\mathbf{\overline{\alpha}}_i^TαiTA=λiαiT
右边乘αi\mathbf{\alpha}_iαi得
α‾iTAαi=λ‾iα‾iTαi\mathbf{\overline{\alpha}}_i^T\mathbf{A}\mathbf{\alpha}_i=\overline{\lambda}_i\mathbf{\overline{\alpha}}_i^T\mathbf{\alpha}_iαiTAαi=λiαiTαi再看Aαi=λiαi\mathbf{A}\mathbf{\alpha}_i=\lambda_i\mathbf{\alpha}_iAαi=λiαi,
对它左边乘α‾iT\overline{\mathbf{\alpha}}_i^{T}αiT可得
α‾iTAαi=λiα‾iTαi\mathbf{\overline{\alpha}}_i^T\mathbf{A}\mathbf{\alpha}_i=\lambda_i\mathbf{\overline{\alpha}}_i^T\mathbf{\alpha}_iαiTAαi=λiαiTαi上面两个式子相减得
0=(λ‾i−λi)α‾iTαi0=(\mathbf{\overline{\lambda}_i-\lambda_i})\mathbf{\overline{\alpha}}_i^T\mathbf{\alpha}_i0=(λi−λi)αiTαi
因为,α‾iTαi\mathbf{\overline{\alpha}}_i^T\mathbf{\alpha}_iαiTαi是非0向量.所以我们可得λ‾i−λi=0\overline{\lambda}_i-\lambda_i=0λi−λi=0.
也就是说λi\lambda_iλi是实数.
又因为λi\lambda_iλi是任意的特征值,所以A\mathbf{A}A,
的所有特征值都是实数.
实对称矩阵属于不同特征值的特征向量正交
接下来我们再来证明一个命题,实对称矩阵属于不同特征值的特征向量正交.我们先假设两个不同的特征值位λi,λj\lambda_i,\lambda_jλi,λj,
他们对应的特征向量为αi,αj\mathbf{\alpha}_i, \mathbf{\alpha}_jαi,αj. 假如,
我们定义(αi,αj)(\mathbf{\alpha}_i, \mathbf{\alpha}_j)(αi,αj)表示点积.
那么我们可以按照下面的推导.
λi(αi,αj)=(λiαi,αj)=(Aαi,αj)=αjTAαi=αiTAαj\lambda_i(\mathbf{\alpha}_i, \mathbf{\alpha}_j)=(\lambda_i\mathbf{\alpha}_i, \mathbf{\alpha}_j)=(\mathbf{A\alpha}_i, \mathbf{\alpha}_j)=\mathbf{\alpha_j^TA\alpha_i}=\mathbf{\alpha_i^TA\alpha_j}λi(αi,αj)=(λiαi,αj)=(Aαi,αj)=αjTAαi=αiTAαj
λj(αi,αj)=(αi,λjαj)=(αi,Aαj)=αiTAαj\lambda_j(\mathbf{\alpha}_i, \mathbf{\alpha}_j)=(\mathbf{\alpha}_i, \lambda_j\mathbf{\alpha}_j)=(\mathbf{\alpha}_i, \mathbf{A\alpha}_j)=\mathbf{\alpha_i^TA\alpha_j}λj(αi,αj)=(αi,λjαj)=(αi,Aαj)=αiTAαj
下相减得(λi−λj)(αi,αj)=0(\lambda_i-\lambda_j)(\mathbf{\alpha}_i, \mathbf{\alpha}_j)=0(λi−λj)(αi,αj)=0,
又因为λi≠λj⇒αiT⋅αj=0\lambda_i\neq \lambda_j \Rightarrow \mathbf{\alpha}_i^T\cdot \mathbf{\alpha}_j=0λi=λj⇒αiT⋅αj=0,
也就是正交成立.
实对称矩阵可对角化
这里我们使用归纳法来证明.
首先假设n=1n=1n=1. A=a11\mathbf{A}=a_{11}A=a11. 这个不证自明.
假设n=k−1n=k-1n=k−1, 命题撤成立.
现在n=kn=kn=k, 我们假设其中一个特征值位λ1\lambda_1λ1,
那么我可以利用第一个特征值对应的特征向量构造一组RnR^nRn的正交基.
T=(η1,η2,⋯,ηn)T=(\eta_1, \eta_2,\cdots, \eta_n)T=(η1,η2,⋯,ηn). 那么我们可以得
T−1AT=(T−1λ1η1,T−1Aη2,⋯,T−1Aηn)T^{-1}AT=(T^{-1}\lambda_1\eta_1, T^{-1}A\eta_2, \cdots, T^{-1}A\eta_n)T−1AT=(T−1λ1η1,T−1Aη2,⋯,T−1Aηn)
又因为T−1T=IT^{-1}T=IT−1T=I, 那么,
我们可以得T−1η1=ε1T^{-1}\eta_1=\mathbf{\varepsilon}_1T−1η1=ε1. 那么可得
T−1AT=(λ1α0B)T^{-1}AT=\left( \begin{array}{cc} \lambda_1 & \mathbf{\alpha} \\ \mathbf{0} & \mathbf{B} \end{array} \right)T−1AT=(λ10αB)由于AAA, 是一个实对称矩阵,
那么T−1ATT^{-1}ATT−1AT也是一个实对称矩阵.进而α=0\mathbf{\alpha}=\mathbf{0}α=0.
由此可见BBB也是一个(k−1)×(k−1)(k-1)\times (k-1)(k−1)×(k−1)的实对称矩阵. 按照假设,
它是可以对角化的.现在假设.
T2−1BT2=diag{λ2,λ3,⋯,λn}T_2^{-1}BT_{2}=diag\{\lambda_2, \lambda_3, \cdots, \lambda_n\}T2−1BT2=diag{λ2,λ3,⋯,λn}
并设 Tf=T(100T2)T_f=T\left( \begin{array}{cc} 1 & 0 \\ 0 & T_2 \end{array} \right)Tf=T(100T2) 那么 Tf−1ATf=(100T2)−1T−1AT(100T2)=(100T2−1)(λ100B)(100T2)T^{-1}_fAT_{f}=\left( \begin{array}{cc} 1 & 0 \\ 0 & T_2 \end{array} \right)^{-1}T^{-1}AT\left( \begin{array}{cc} 1 & 0 \\ 0 & T_2 \end{array} \right)=\left( \begin{array}{cc} 1 & 0 \\ 0 & T_2^{-1} \end{array} \right)\left( \begin{array}{cc} \lambda_1 & 0 \\ 0 & B \end{array} \right)\left( \begin{array}{cc} 1 & 0 \\ 0 & T_2 \end{array} \right)Tf−1ATf=(100T2)−1T−1AT(100T2)=(100T2−1)(λ100B)(100T2) =(λ100T2−1BT2)=diag{λ1,λ2,⋯,λn}=\left( \begin{array}{cc} \lambda_1& 0 \\ 0 & T_2^{-1}BT_{2} \end{array} \right)=diag\{\lambda_1,\lambda_2, \cdots, \lambda_n\}=(λ100T2−1BT2)=diag{λ1,λ2,⋯,λn}
至此, 原式得证. 也就是说实对称矩阵一定可以对角化.
反过来也说明实对称矩阵的特征多项式的代叔重述等于几何重数.
实对称矩阵一定可以对角化相关推荐
- 线性代数的问题:是否存在这样的矩阵,它满足正交对角化的条件,但它不是实对称矩阵呢?
对称矩阵的对角化问题 定理 :对称矩阵的特征值是实数. 定理:设A是n阶对称阵,则必有正交阵P,使得P−1AP=PTAP=ΛP^{-1}AP=P^{T}AP= \LambdaP−1AP=PTAP=Λ, ...
- 【机器学习】【线性代数 for PCA】矩阵与对角阵相似、 一般矩阵的相似对角化、实对称矩阵的相似对角化
Note:PCA主成分分析用到实对称阵的相似对角化,用个文章复习一下相关概念和计算过程. 1.对角矩阵 如果一个矩阵满足如下条件,则它就是一个对角阵: (1)是一个方阵 (2)只有对角线元素是非零元素 ...
- 实对称矩阵的性质_浅谈矩阵的相似对角化(一)
森屿瑾年:浅谈线性变换和矩阵之间的关系zhuanlan.zhihu.com 通过前面的讨论,我们引出了线性变换在不同基下的矩阵之间的关系,知道了线性变换在不同基下的矩阵是相似的,进而我们可以通过选取 ...
- 参考答案:05 实对称矩阵与二次型
本篇图文为<线性代数及其应用>这本教材对应习题册的参考答案. 本章是特征值与特征向量知识的延续,根据谱定理可知实对称矩阵可以正交对角化,对角阵为其特征值,正交矩阵为其两两正交的单位特征向量 ...
- 线性代数:05 实对称矩阵与二次型
本讲义是自己上课所用幻灯片,里面没有详细的推导过程(笔者板书推导)只以大纲的方式来展示课上的内容,以方便大家下来复习. 本章是特征值与特征向量知识的延续,根据谱定理可知实对称矩阵可以正交对角化,对角阵 ...
- 矩阵的基础知识回顾:矩阵乘法,矩阵的逆,伴随矩阵,矩阵的转置,行列式,相似矩阵,实对称矩阵
Agenda 1. 矩阵matrix 1.1 矩阵运算matrix operations 1.1.1 矩阵乘法matrix multiplication 1.1.1.1 简化矩阵乘法(facilita ...
- 证明:对于实对称矩阵,不同特征值对应的特征向量相互正交
前言 不同特征值对应的特征向量相互正交,是实对称矩阵的一个重要属性,而且从这个属性出发可以证明实对称矩阵的另一个属性:实对称矩阵必可相似对角化.对于一个 n 维矩阵,其可相似对角化的充分且必要条件是- ...
- 线性代数学习笔记——第七十三讲——实对称矩阵的特征值与特征向量
1. 实对称矩阵的特征值都是实数 2. 实对称矩阵不同特征值的实特征向量相互正交
- 实对称矩阵为正定矩阵的一个充分必要条件
本文是为了在学习凸优化的时候遇到的一个问题展开讨论的.目的是能够明白凸优化的理论基础,或者尽可能的明白它的理论基础. 1,对称矩阵的特征值是实数. 证明如下:(我是用latex编辑的,这里不能显示公式 ...
- 对称矩阵与实对称矩阵
对称矩阵(Symmetric Matrices):是指以主对角线为对称轴,各元素对应相等的矩阵. [1] 在线性代数中,对称矩阵是一个方形矩阵,其转置矩阵和自身相等. 实对称矩阵:如果有n阶矩阵A,其 ...
最新文章
- php数独游戏开发,使用vue如何开发数独游戏
- Android 获取keystore SHA1方法
- 团队作业-Beta冲刺(周三)
- C语言编程蛟龙转身,[蓝桥杯][2014年第五届真题]排列序数-题解(C语言代码)
- Internet Explorer 无法显示网页
- Office Word2019您正试图运行的函数包含有宏或需要宏语言
- PyQt4开发环境搭建指导
- myecplise 完全开发
- Maven最佳实践:管理依赖
- 最大似然估计程序c语言,极大似然估计(示例代码)
- 跟父亲一样伟大的程序员,请照顾好自己!
- 概率论中的一些基础知识——条件概率 先验概率 后验概率 似然 概率分布函数 概率密度函数
- python词云生成与设计实现_python简单词云生成
- 复制英文和中文PDF文字到word并去除回车符
- 千峰python线下培训
- 利率风险结构和期限结构
- python正则表达式aabb式成语_python——正则表达式
- 算法的力量 李开复__转自李开复网易博客
- Window VMware MAC OS 安装
- 2021年安全生产模拟考试(建筑特种作业操作证-建筑架子工模拟考试题库)安考星