-
2021-11-16 10:16:20
线性代数之 实对称矩阵,正交对角化,二次型与正定矩阵
前言
终于快到矩阵分解了。在矩阵分解前,最后一个内容是实对称矩阵,二次型和正定矩阵。这三个概念与矩阵分解相关。
实对称矩阵
对于矩阵 A ∈ R n × n A\in R^{n\times n} A∈Rn×n,如果 A T = A A^T=A AT=A,则称 A A A为实对称矩阵。
实对称矩阵不同特征值的特征向量正交,n重特征值有n个线性无关的特征向量。因此实对称矩阵必然能够对角化。
实对称矩阵是 n × n n\times n n×n矩阵能够正交对角化的充分必要条件。
正交对角化
如果存在一个正交矩阵 Q Q Q,使得方阵 A = Q Λ Q − 1 = Q Λ Q T A=Q\Lambda Q^-1=Q\Lambda Q^T A=QΛQ−1=QΛQT能够对角化,称为正交对角化。
能够正交对角化的矩阵都是对称矩阵。
证明:
A = Q Λ Q T A T = Q Λ T Q T = Q Λ Q T = A A=Q\Lambda Q^T \\ A^T=Q\Lambda^T Q^T=Q\Lambda Q^T=A A=QΛQTAT=QΛTQT=QΛQT=A二次型
A A A是实对称矩阵,将一个变量满足 f ( x ) = x T A x f(x)=x^TAx f(x)=xTAx函数称为二次型。
对于 f ( x ) = x T A x f(x)=x^TAx f(x)=xTAx,替换变量 x = P y , f ( x ) = f ( P y ) = y T P T A P y x=Py,f(x)=f(Py)=y^TP^TAPy x=Py,f(x)=f(Py)=yTPTAPy,而 A A A是实对称矩阵,因此存在正交矩阵 Q , f ( Q y ) = y T Λ y Q,f(Qy)=y^T\Lambda y Q,f(Qy)=yTΛy,使得二次型化为标准型。
正定矩阵
广义的正定矩阵:对于矩阵 A ∈ R n × n A\in R^{n\times n} A∈Rn×n,函数 f ( x ) = x T A x > 0 f(x)=x^TAx>0 f(x)=xTAx>0对任意非零向量 x ∈ R n x\in R^n x∈Rn都成立,则称 A A A为正定矩阵。如果 f ( x ) = x T A x ≥ 0 f(x)=x^TAx\ge0 f(x)=xTAx≥0,则称 A A A为半正定矩阵。
狭义的正定矩阵:对于对称矩阵 A ∈ R n × n A\in R^{n\times n} A∈Rn×n,函数 f ( x ) = x T A x > 0 f(x)=x^TAx>0 f(x)=xTAx>0对任意非零向量 x ∈ R n x\in R^n x∈Rn都成立,则称 A A A为正定矩阵。如果 f ( x ) = x T A x ≥ 0 f(x)=x^TAx\ge0 f(x)=xTAx≥0,则称 A A A为半正定矩阵。也把这种正定矩阵称为对称正定矩阵。
实对称矩阵的正定判断条件
如果实对称矩阵的特征值都大于0,则是对称正定矩阵;如果特征值都非负,则是对称半正定矩阵。
证明:
f ( x ) = x T A x Q Q T = E x = Q y f ( x ) = f ( Q y ) = y T Q T A Q y = y T Λ y = ∑ i = 1 n λ i y i 2 i f λ i > 0 , f ( x ) > 0 i f λ i ≥ 0 , f ( x ) ≥ 0 f(x)=x^TAx \\ QQ^T=E \\ x=Qy\\ f(x)=f(Qy)=y^TQ^TAQy=y^T\Lambda y=\sum_{i=1}^n \lambda_i y_i^2 \\ if \quad \lambda_i>0,f(x)>0 \\ if \quad \lambda_i\ge0,f(x)\ge 0\\ f(x)=xTAxQQT=Ex=Qyf(x)=f(Qy)=yTQTAQy=yTΛy=i=1∑nλiyi2ifλi>0,f(x)>0ifλi≥0,f(x)≥0一个常见的半正定矩阵
对于任意矩阵 A ∈ R m × n A\in R^{m\times n} A∈Rm×n,矩阵 A T A A^TA ATA是半正定矩阵。
证明:
x T A T A x = ( A x ) T A x = ( A x ) ⋅ ( A x ) = ∣ ∣ A x ∣ ∣ 2 2 ≥ 0 x^TA^TAx=(Ax)^TAx=(Ax)\cdot(Ax)=||Ax||^2_2\ge0 xTATAx=(Ax)TAx=(Ax)⋅(Ax)=∣∣Ax∣∣22≥0如果 A ∈ R m × n A\in R^{m\times n} A∈Rm×n列满秩,则 A T A A^TA ATA是正定矩阵。
证明:
A x = 0 → A T A x = 0 A T A x = 0 → x T A T A x = 0 → ∣ ∣ A x ∣ ∣ 2 2 = 0 → A x = 0 ∴ N U L ( A ) = N U L ( A T A ) ∵ d i m N U L ( A ) + r a n k ( A ) = n ∴ r a n k ( A T A ) = r a n k ( A ) r a n k ( A ) = n , d i m N U L ( A ) = 0 ∴ ∀ x ≠ 0 , A x ≠ 0 ∴ ∣ ∣ A x ∣ ∣ 2 2 > 0 i . e . x T A T A x > 0 Ax=0 \to A^TAx=0 \\ A^TAx=0 \to x^TA^TAx=0 \to ||Ax||^2_2=0\to Ax=0 \\ \therefore NUL(A)=NUL(A^TA) \\ \because dimNUL(A)+rank(A)=n \\ \therefore rank(A^TA)=rank(A) \\ \quad \\ rank(A)=n,dimNUL(A)=0 \\ \therefore \forall x\ne 0, A x\ne0 \\ \therefore ||Ax||_2^2>0 \\ \quad \\ i.e. \quad x^TA^TAx>0 Ax=0→ATAx=0ATAx=0→xTATAx=0→∣∣Ax∣∣22=0→Ax=0∴NUL(A)=NUL(ATA)∵dimNUL(A)+rank(A)=n∴rank(ATA)=rank(A)rank(A)=n,dimNUL(A)=0∴∀x=0,Ax=0∴∣∣Ax∣∣22>0i.e.xTATAx>0后记
线性代数的矩阵性质部分大概就记录完了。下一篇就进入到了矩阵计算的内容——矩阵分解。
更多相关内容 -
两个实对称阵同时正交对角化的MATLAB程序 (2013年)
2021-05-11 13:24:12使用Matlab 语言设计出实现两个实对称阵同时正交相似对角化的计算机程序。 -
将实对称矩阵正交对角化的流程
2020-11-15 16:45:21 -
两个实对称阵同时正交对角化的MATLAB程序.pdf
2021-07-10 10:22:23两个实对称阵同时正交对角化的MATLAB程序.pdf -
对称矩阵的正交对角化和它的奇异分解
2020-12-18 20:24:05对称矩阵的正交对角化就是他的奇异分解- 对称矩阵的正交对角化就是他的奇异分解
- 当然不是啦!
- 因为奇异值分解 Σ \Sigma Σ是个对角线都是正的啊!
- 而对称矩阵的正交对角化显然不是的!
- 如果
A
A
A是对称的,特征值为
λ
1
,
.
.
.
,
λ
n
\lambda_1,...,\lambda_n
λ1,...,λn
- 且 ∣ λ 1 ∣ ≥ ∣ λ 2 ∣ ≥ . . . ≥ ∣ λ r ∣ > ∣ λ r + 1 ∣ = . . . = ∣ λ n ∣ = 0 |\lambda_1|\ge |\lambda_2|\ge ...\ge|\lambda_r|>|\lambda_{r+1}|=...= |\lambda_n|=0 ∣λ1∣≥∣λ2∣≥...≥∣λr∣>∣λr+1∣=...=∣λn∣=0
- 那 A = W Σ W T A=W\Sigma W^T A=WΣWT, W = ( α 1 , . . . , α n ) W=(\alpha_1,...,\alpha_n) W=(α1,...,αn)
- A 2 = W Σ 2 W T A^2=W\Sigma^2W^T A2=WΣ2WT
- 这说明
A
2
A^2
A2也可对角化,特征值为
λ
1
2
,
.
.
.
,
λ
n
2
\lambda_1^2,...,\lambda_n^2
λ12,...,λn2
- 特征向量为 W = ( α 1 , . . . , α n ) W=(\alpha_1,...,\alpha_n) W=(α1,...,αn)
-
设 A = U Σ V T A=U\Sigma V^T A=UΣVT
-
Σ
\Sigma
Σ的对角线是
A
2
A^2
A2的特征值开根号
- 分别为 ∣ λ 1 ∣ , ∣ λ 2 ∣ , . . . , ∣ λ n ∣ |\lambda_1|,|\lambda_2|,...,|\lambda_n| ∣λ1∣,∣λ2∣,...,∣λn∣
-
Σ
\Sigma
Σ的对角线是
A
2
A^2
A2的特征值开根号
-
V V V的每一列是 A 2 A^2 A2的特征向量
- 也就是 ( α 1 , . . . , α n ) (\alpha_1,...,\alpha_n) (α1,...,αn)
-
U U U的每一列就是 ( A α 1 , . . . , A α n ) (A\alpha_1,...,A\alpha_n) (Aα1,...,Aαn)
-
所以
- 上面写的不对。
- 要记得把 A α i A\alpha_i Aαi单位化一下!
-
实对称矩阵必可正交对角化证明
2018-08-05 13:36:26n阶矩阵A可正交对角化的充分条件是A是实对称矩阵,即若A是实对称矩阵则A必可正交对角化。 首先,有以下定理: 若的特征值为,且,则存在正交矩阵Q,使A相似于如下三角矩阵: 证明如下(数学归纳法): 设n*n阶...n阶矩阵A可正交对角化的充分条件是A是实对称矩阵,即若A是实对称矩阵则A必可正交对角化。
首先,有以下定理:
若
的特征值为
,且
,则存在正交矩阵Q,使A相似于如下三角矩阵:
-
计算机视觉中的数学方法——7正交对角化——7.1 内积空间与正交矩阵
2019-10-17 09:25:137.1 内积空间与正交矩阵7.1.1 内积空间7.1.2 正交矩阵正交矩阵的几何意义2阶正交矩阵3阶正交矩阵 7.1.1 内积空间 定义一个函数<.,.>: 7.1.2 正交矩阵 单位正交:正交矩阵每个向量的模长都为1。 正交矩阵的... -
线性代数的问题:是否存在这样的矩阵,它满足正交对角化的条件,但它不是实对称矩阵呢?
2020-11-03 14:38:38想到了一点东西: 命题:实对称矩阵 ⇒\Rightarrow⇒ 正交对角化 逆否命题: 若一个矩阵不能正交对角化,那它一定不是对称矩阵。 然后想了想,是否存在 可以正交对角化,但它不是实对称的矩阵的呢? 首先,如果可... -
非正交联合对角化盲源分离算法 (2008年)
2021-05-09 21:42:58提出一种表征联合对角化近似程度的代价函数,给出优化该代价函数的非正交联合对角化算法。该代价函数对经典的最小二乘代价函数进行了改进,使原算法中关于混迭矩阵的四次函数转化为3 组待定参数的二次函数。因此,可... -
实对称阵的正交相似对角化
2022-04-10 12:54:35一、正交向量组与正交矩阵 正交向量组的定义,是一组非零向量,且两两正交,那么这组向量,则成为正交向量组。 两个向量正交的意思是,两个向量的内积为0,什么是两个向量的内积,就是向量内对应元素的积的和。 ... -
实对称矩阵必可正交相似对角化
2020-05-22 11:35:02k重特征值对应的k个特征向量可以施密特正交化得到新的k个特征向量,而依然对应同一个特征值;不同特征值的对应特征向量一定正交。而特征向量乘以常数依然对应同一个特征值,从而可以使n个特征向量是单位向量且两两... -
第四节 实对称矩阵的对角化
2021-04-20 14:06:19《第四节 实对称矩阵的对角化》由会员分享,可在线阅读,更多相关《第四节 实对称矩阵的对角化(9页珍藏版)》请在人人文库网上搜索。1、一、实对称矩阵的性质,1,2,1,2,1,2,1,2,6.4.1,A,引理,设,是实对称矩阵,的两个... -
对称矩阵的SVD分解和正交对角化的区别
2020-12-18 23:42:42对称矩阵的SVD分解和正交对角化的区别 -
非正交联合对角化的罚函数粒子群算法
2021-01-14 17:04:21针对多个矩阵近似联合对角化盲分离问题, 提出一种新的非正交近似联合对角化算法. 首先采用罚函数法 将联合对角化的非线性约束优化模型转化为无约束优化模型; 其次将粒子群优化算法引入无约束优化模型中实现目... -
实对称矩阵一定要用正交矩阵来对角化吗?
2020-11-12 22:12:50不一定,可以直接用一般矩阵的方法求其对角阵,即可以不用正交单位化,直接用【p逆Ap=A的对角阵】来做,书上一上来就说用正交阵来对角化就是淡村为了体现这个方法而已,但是还是有好处的,比如正交单位化后,要求p逆... -
MATLAB学习笔记:方阵的相似对角化
2021-04-18 08:14:510.1695 D = 3.0000 0 0 0 0 0 5.0000 0 0 0 0 0 5.0000 0 0 0 0 0 1.0000 0 0 0 0 0 1.0000 >> r=rank(V) r = 5A有5个线性无关的特征向量,所以矩阵A可以相似对角化,取P=V,则 P^(-1)*A*P=D >> B=[-2 1 -2;... -
EigenliftDiagonalization:Mathematica代码以证明Eignelift对角化
2021-04-13 23:21:43本征提升对角化 Mathematica代码以证明Eignelift对角化 -
实对称矩阵对角化含Matlab代码
2010-05-19 13:26:48实对称矩阵,对角化,有两篇论文,内含Matlab实现代码,在文章里的,可以直接写下来使用。测试过,还可以。 -
线性代数笔记:矩阵的对角化、SVD分解及应用
2019-07-31 20:45:38矩阵的对角化、SVD分解及应用矩阵的对角化、SVD分解及应用矩阵运算的总结矩阵对角化SVD分解张量(tensor)转置(transpose) 矩阵的对角化、SVD分解及应用 许多数学对象可以通过将它们分解成多个组成部分或者找到它们的... -
实对称阵正交相似对角化方法及其优化 (2013年)
2021-06-14 04:59:00对代数教材中实对称阵正交相似对角化的经典方法从三个方面进行了优化:线性方程组法省去施密特正交化步骤,特别当实对称阵有n-1重特征值时更为简化;积矩阵法改解线性方程组为求矩阵的极大无关列;对称变换法则颠覆... -
实对称阵可对角化的几种证明
2019-10-03 23:49:35实对称阵是一类常见的矩阵, 它与实二次型和实内积空间上的自伴随算子有着密切的联系. 任一实对称阵$A$ 均正交相似于对角阵, 即存在正交阵 $P$, 使得 $P'AP=\mathrm{diag\,}\{... 然而, 实对称阵可对角化这一性质可以... -
mathematica 矩阵对角化
2020-08-12 16:38:09mathematica 矩阵对角化 eigensystem给出的Vector默认是列向量,使用vectors做矩阵乘法的时候,就是利用向量的内积来做乘法 eg. A.B(AB都是矩阵),AB=(A1A2A3)(B1B2B3),其中A1A2A3是A的行向量,B1B2B3是B的列向量 -
某些正交矩阵的对角化 (2011年)
2021-05-18 01:50:57A是一个实正交矩阵,给出了n维复向量组线性相关和线性无关的定义,证明了二阶正交矩阵的对角化,以及某些三阶正交矩阵对角化. -
基于高斯-牛顿迭代的非正交联合对角化算法 (2014年)
2021-05-28 15:39:42针对传统的分离算法因迭代次数过多而不能满足通信信号分离时对信号实时处理的要求,将最佳权矩阵引入到联合对角化准则中,提出了一种改进的基于“高斯-牛顿”迭代法的非正交联合对角化算法(WEDGE),提高了算法的分离... -
基于非正交联合对角化的DOA估计新算法
2021-04-26 13:31:37基于非正交联合对角化的DOA估计新算法 -
考研数学基础 之线性代数通法——Chapter6:合同对角化与二次型
2021-05-19 22:10:24合同对角化与二次型 2022考研数学基础 主讲老师: 刘金峰 武忠祥 对称矩阵的对角化 考研范围内只考察实对称矩阵 以下内容所表述的对象均为实对称矩阵 性质 对称矩阵的特征值为实数 不同特征值对应的特征向量正交... -
实对称矩阵可对角化的证明
2020-04-01 19:38:20对于实对称矩阵 AAA,有 J=P−1APJ=P^{-1}APJ=P−1AP, JJJ为 AAA的若当标准型。 而 JT=PTAP−T=PTPP−1APP−1(P−1)T=(PTP)J(PTP)−1J^{T}=P^{T}AP^{-T}=P^TPP^{-1}APP^{-1}(P^{-1})^T=(P^TP)J(P^TP)... 即 AAA可对角化 -
线性代数笔记5.3实对称矩阵的对角化
2021-12-11 17:01:535.3实对称矩阵的对角化 正交矩阵 概念 A是一个n阶方阵 有A转置A == E,那么这个矩阵就是正交矩阵 性质 若A是正交矩阵,|A| == 1或 -1 若A是正交矩阵,A的逆 == A的转置 且A逆和A转置均为正交 A,B正交,AB也正交 若A... -
基于Lapack接口的矩阵双对角化:代码通过调用Lapack例程计算矩阵的双对角矩阵分解。-matlab开发
2021-06-01 09:22:36因此,此提交帮助人们使用矩阵的双对角化,A = Q * B * P',其中 A 是非平方的复数或实数,B 是双对角的实数,Q 和 P 是酉或正交的。 最初的 Lapack 接口例程来自 Tim Toolan 在“文件交换”中的工作。 此代码可... -
实对称矩阵 二次型 合同 相似对角化
2019-11-13 23:13:294.一定可以用正交矩阵相似对角化(满足的矩阵为正交阵),步骤如下 (1)求A的特征值λ1、λ2、λ3 (2)求特征向量α1、α2、α3 (3)改造特征向量 a. 如λi≠λj 只需要单位化 ...