【线性代数】6-4:对称矩阵(Symmetric Matrices)
title: 【线性代数】6-4:对称矩阵(Symmetric Matrices)
categories:
- Mathematic
- Linear Algebra
keywords: - Eigenvalues
- Eigenvectors
- Symmetric Matrices
- Projection Matrices
- Spectral Theorem
- Principal Axis Theorem
toc: true
date: 2017-11-22 15:18:03
Abstract: 本篇继续线性代数的高潮部分,对称矩阵,以及对称矩阵的一些性质
Keywords: Eigenvalues,Eigenvectors,Symmetric Matrices,Projection Matrices,Spectral Theorem,Principal Axis Theorem
开篇废话
这几篇在难度上确实要比前面的内容大很多,所以看书理解和总结都变得不那么流畅了,但是慢慢看下来收获还是有很大的,而且我发现不管学的多认真,还是会有遗漏,所以我觉得之前的想法就是一次性把什么什么学会是不可能的,只能学到自己觉得达到自己能发现的最大限度,等到应用之时还是要回来查阅,这样往往会有进一步的更大发现,
Symmetric Matrices
对称矩阵我们在最早的知识里面就学过 AT=AA^T=AAT=A 的矩阵叫做对称矩阵,我们也学过投影矩阵,但是当时我们并没有强调过一点就是投影矩阵都是对称的,这个性质今天在这里会有很大的用途。
我们继续说投影矩阵,所谓投影矩阵,就是在和向量 c⃗\vec{c}c 相乘的时候,投影到矩阵A的列空间内,那么其中,投影 ppp 和 原向量 c⃗\vec{c}c 的差 e⃗=c⃗−p⃗\vec{e} =\vec{c}-\vec{p}e=c−p 与子空间正交。
举个例子,在三维空间内,A的列空间是一个二维平面那么,A对应的投影矩阵P能够把任何方向的向量投影到平面上,那么如果向量本身属于平面那么 Px=xPx=xPx=x 显然是不用质疑的(我们之前在投影那篇文章中也讲过) 但是,同志们,看看这个有木有很面熟啊,这个明显就是投影矩阵 PPP 的特征值和特征向量么?没错,PPP 有一个平面的特征向量,可以随便选!能选多少个呢、当然是无数个,但是问题又来了,这无数多个并不是独立的,因为一共就二维,选出来三个线性独立的向量都是不可能的,所以这个平面能选出两个线性独立的特征向量,并且对应的特征值都是1,这里有人可能疑惑为啥要选两个,因为我们6-2的时候说过只有特诊向量足够的情况下才能对角化,投影矩阵明显是个3x3的矩阵,那么特征向量也应该有三个呀!我们的子空间是二维的,所以理论上应该有两个特征向量在上面,剩下一维存在一个,那么这一个也能很好找,e⃗\vec{e}e 就是 也就是和子空间正交的向量都行 Px=0xPx=0xPx=0x 表明 x⃗\vec{x}x 和子空间正交,那么这是个特征值为0的特征向量,这样我们又进一步规范一下,选择三个特征向量相互正交,这个也是可以做到的,也就是对于矩阵P我们找到了三个相互正交的特征向量,并且长度缩放到单位长度。
以上三维投影到二维平面可以通过几何来解释,但为了能让大家从线性空间来理解,就没用几何方法,大家可以自己脑补。
得出结论,对称矩阵 PT=PP^T=PPT=P 的特征向量相互正交并且为单位向量。
对称矩阵"It is no exaggeration to say that these are the most important matrices the world will ever see – in the theory of linear algebra and alos in the applications" 翻译成中文:“对称矩阵是史上最牛B的矩阵,无论在理论还是应用”
这个我们目前还无法考证,还没做过应用呢?不是么,但是我知道PCA中确实用了对称矩阵,SVD等一些列相关技术。
一个矩阵能被如此称赞,不外乎几点原因,首先是其本身拥有较好的性质,其次这个矩阵在自然生活中经常出现,就像正态分布,那么难的公式,却能准确的描述自然届的现象。最后就是如果表现形式简单,那么这个就是非常有用的东西啦。
下面我们开始探索对称矩阵的性质。
如果一个对称矩阵满足:
suppose:AT=AA=SΛS−1then:AT=(S−1)TΛTST
suppose:\\
A^T=A\\
A=S\Lambda S^{-1}\\
then:\\
A^T=(S^{-1})^T\Lambda^T S^T
suppose:AT=AA=SΛS−1then:AT=(S−1)TΛTST
这种情况下就有下面这种***可能***了,也就是对应的 S=(S−1)TS=(S^{-1})^TS=(S−1)T 注意我们这里说的是可能,并不排除不可能的情况,原文书上用的也是possibly,也就是说我们目前假设:
ST=S−1STS=I
S^T=S^{-1}\\
S^TS=I
ST=S−1STS=I
这里我们可以预报一下:
- 对称矩阵只有实数特征值
- 对称矩阵特征向量可以选择正交单位向量 orthonormal
对于 STS=IS^TS=ISTS=I 面熟么?还有印象么?我们认识啊,正交矩阵
QTQ=IQ^TQ=IQTQ=I 矩阵Q中每列之间相互正交,也就是我们对于对称矩阵可以写成:
A=QΛQ−1=QΛQTwith:Q−1=QT
A=Q\Lambda Q^{-1}=Q\Lambda Q^T\\
with:\\
Q^{-1}=Q^T
A=QΛQ−1=QΛQTwith:Q−1=QT
这个就是著名的普定理 “Spectral Theorem”:
Every symmetric matrix has the factorization A=QΛQTA=Q\Lambda Q^TA=QΛQT with real eigenvalues in Λ\LambdaΛ and orthonormal eigenvectors in S=QS=QS=Q
对于所有对称矩阵都能分解成 A=QΛQTA=Q\Lambda Q^TA=QΛQT 的形式并且在 Λ\LambdaΛ 中的所有特征值都是实数,其对应的特征向量是正交单位矩阵,即 S=QS=QS=Q
本文为节选,完整内容地址:https://www.face2ai.com/Math-Linear-Algebra-Chapter-6-4转载请标明出处
【线性代数】6-4:对称矩阵(Symmetric Matrices)的更多相关文章
- 【线性代数】6-6:相似矩阵(Similar Matrices)
title: [线性代数]6-6:相似矩阵(Similar Matrices) categories: Mathematic Linear Algebra keywords: Similar Matr ...
- 【线性代数】6-5:正定矩阵(Positive Definite Matrices)
title: [线性代数]6-5:正定矩阵(Positive Definite Matrices) categories: Mathematic Linear Algebra keywords: Po ...
- 线性代数导论 | Linear Algebra 课程
搞统计的线性代数和概率论必须精通,最好要能锻炼出直觉,再学机器学习才会事半功倍. 线性代数只推荐Prof. Gilbert Strang的MIT课程,有视频,有教材,有习题,有考试,一套学下来基本就入 ...
- Properties in Algebra
附录-Properties in Algebra 部分证明转载自标注,仅作个人整理查阅用. 范数 (norm) \(^{[1]}\) 要更好的理解范数,就要从函数.几何与矩阵的角度去理解,我尽量讲的通 ...
- Linear Algebra Lecture5 note
Section 2.7 PA=LU and Section 3.1 Vector Spaces and Subspaces Transpose(转置) example: 特殊情况,对称 ...
- Moore-Penrose Matrix Inverse 摩尔-彭若斯广义逆 埃尔米特矩阵 Hermitian matrix
http://mathworld.wolfram.com/Moore-PenroseMatrixInverse.html 显然,埃尔米特矩阵主对角线上的元素都是实数的,其特征值也是实数.对于只包含实数 ...
- A.Kaw矩阵代数初步学习笔记 4. Unary Matrix Operations
“矩阵代数初步”(Introduction to MATRIX ALGEBRA)课程由Prof. A.K.Kaw(University of South Florida)设计并讲授. PDF格式学习笔 ...
- We Recommend a Singular Value Decomposition
We Recommend a Singular Value Decomposition Introduction The topic of this article, the singular val ...
- <转>机器学习笔记之奇异值分解的几何解释与简单应用
看到的一篇比较好的关于SVD几何解释与简单应用的文章,其实是有中文译本的,但是翻译的太烂,还不如直接看英文原文的.课本上学的往往是知其然不知其所以然,希望这篇文能为所有初学svd的童鞋提供些直观的认识 ...
随机推荐
- Arm-Linux 移植 QT5.9 带 tslib,QT-creator配置
平台 : Ubuntu 16.04 QT :5.9.8 tslib : 1.4arm-gcc ...
- (二)手动配置第一个HelloWorld程序
上例的HelloWorld是由Android sutudio 自动生成的,现在我们手动来配置. 1. 重新创建工程 2. 创建空的Activity 生成的MainActivity.java 文件: p ...
- Windows群集之NLB【转】
本文转自:http://www.talkwithtrend.com/Article/31746 网络负载平衡群集(Network Load balancing) 在Internet快速发展的今天,为了 ...
- B树Java代码实现以及测试
B树定义 B 树又叫平衡多路查找树.一棵m阶的B 树 (m叉树)的特性如下: 根节点至少有两个孩子 每个非根节点至少有M/2(上取整)个孩子,至多有M个孩子. 每个非根节点至少有M/2-1(上取整)个 ...
- nginx搭建反向代理服务器详解
一.反向代理:Web服务器的“经纪人” 1.1 反向代理初印象 反向代理(Reverse Proxy)方式是指以代理服务器来接受internet上的连接请求,然后将请求转发给内部网络上的服务器,并将从 ...
- SSL/TLS 受诫礼(BAR-MITZVAH)攻击漏洞(CVE-2015-2808)
最近发现SSL/TLS漏洞已经修改过,但是绿盟扫描器还可以扫描出来,网上看了很多文章,但是能用的比较少,今天刚好有空,就自己写一下. 方法一: 控制面板--->系统和安全--->管理工具- ...
- Django drf:幂等性
一.什么叫做幂等性 用户对于同一操作发起的一次请求或者多次请求的结果是一致的,不会因为多次点击而产生了副作用.举个最简单的例子,那就是支付,用户购买商品使用约支付,支付扣款成功,但是返回结果的时候网络 ...
- inotify+rsync文件实时同步
原文转自http://dl528888.blog.51cto.com/2382721/771533/ 之前做了“ssh信任与scp自动传输脚本”的技术文档,此方案是作为公司里备份的方法,但在实际的运行 ...
- SQL性能优化思路
1. 尽可能把数据的存储和计算放入Memory而不是Disk,且减少IO操作,比如运用Redis等缓存技术 2. 对数据表进行精心设计,特别是大数据表,对常用数据字段进行适当的冗余,尽可能避免分表导致 ...
- sql 发生死锁
SELECT request_session_id spid , OBJECT_NAME(resource_associated_entity_id) tableName FROM sys.dm_tr ...
