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

这里我们可以预报一下:

  1. 对称矩阵只有实数特征值
  2. 对称矩阵特征向量可以选择正交单位向量 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)的更多相关文章

  1. 【线性代数】6-6:相似矩阵(Similar Matrices)

    title: [线性代数]6-6:相似矩阵(Similar Matrices) categories: Mathematic Linear Algebra keywords: Similar Matr ...

  2. 【线性代数】6-5:正定矩阵(Positive Definite Matrices)

    title: [线性代数]6-5:正定矩阵(Positive Definite Matrices) categories: Mathematic Linear Algebra keywords: Po ...

  3. 线性代数导论 | Linear Algebra 课程

    搞统计的线性代数和概率论必须精通,最好要能锻炼出直觉,再学机器学习才会事半功倍. 线性代数只推荐Prof. Gilbert Strang的MIT课程,有视频,有教材,有习题,有考试,一套学下来基本就入 ...

  4. Properties in Algebra

    附录-Properties in Algebra 部分证明转载自标注,仅作个人整理查阅用. 范数 (norm) \(^{[1]}\) 要更好的理解范数,就要从函数.几何与矩阵的角度去理解,我尽量讲的通 ...

  5. Linear Algebra Lecture5 note

    Section 2.7     PA=LU and Section 3.1   Vector Spaces and Subspaces   Transpose(转置) example: 特殊情况,对称 ...

  6. Moore-Penrose Matrix Inverse 摩尔-彭若斯广义逆 埃尔米特矩阵 Hermitian matrix

    http://mathworld.wolfram.com/Moore-PenroseMatrixInverse.html 显然,埃尔米特矩阵主对角线上的元素都是实数的,其特征值也是实数.对于只包含实数 ...

  7. A.Kaw矩阵代数初步学习笔记 4. Unary Matrix Operations

    “矩阵代数初步”(Introduction to MATRIX ALGEBRA)课程由Prof. A.K.Kaw(University of South Florida)设计并讲授. PDF格式学习笔 ...

  8. We Recommend a Singular Value Decomposition

    We Recommend a Singular Value Decomposition Introduction The topic of this article, the singular val ...

  9. <转>机器学习笔记之奇异值分解的几何解释与简单应用

    看到的一篇比较好的关于SVD几何解释与简单应用的文章,其实是有中文译本的,但是翻译的太烂,还不如直接看英文原文的.课本上学的往往是知其然不知其所以然,希望这篇文能为所有初学svd的童鞋提供些直观的认识 ...

随机推荐

  1. Solr介绍 入门练习

    1.1 Solr是什么 Solr是一个基于全文检索的企业级应用服务器. 全文检索:可以输入一段文字,通过分词检索数据!!(复习) 应用服务器:它是单独的服务. 1.2 Solr能做什么 它就是用于做全 ...

  2. java异常那些事

    异常的基本定义: 异常情形是指阻止当前方法或者作用域继续执行的问题.在这里一定要明确一点:异常代码某种程度的错误,尽管Java有异常处理机制,但是我们不能以“正常”的眼光来看待异常,异常处理机制的原因 ...

  3. 动手篇:简单的注册界面与防SQL注入(续)

    注册前先判断用户名是否已经存在,通过if (SqlHelper.Exists(strSql))查询判断是否已经存在.没存在则进行一个数据插入数据库操作 string strSql1 = string. ...

  4. css 小图标 & iconfont 字体图标

    前言:这是笔者学习之后自己的理解与整理.如果有错误或者疑问的地方,请大家指正,我会持续更新! iconfont 字体图标 我们的需求中,很多时候会看到一些小的图形,或者叫图标,比如天猫网站中:     ...

  5. hdu 1114需要装满的完全背包 重点是背包初始化的问题

    .,. 最近在看背包九讲 所以就刷了一下背包的题目 这道题目是一个典型的完全背包问题 而且要求满包 在这里 我就简单整理一下背包初始化问题吧 对于没有要求满包的问题 也就是背包可以不取满的问题 在背包 ...

  6. (三)自定义Realm

    一.Realm概念 Realm:域,Shiro从从Realm获取安全数据(如用户.角色.权限),就是说SecurityManager要验证用户身份,那么它需要从Realm获取相应的用户进行比较以确定用 ...

  7. ASP.NET Core 入门(1)(搭建环境CentOS)

    一.CentOS 7 安装 下载CentOS http://isoredirect.centos.org/centos/7/isos/x86_64/  选择其中下载即可. 下载完成后打开vmware准 ...

  8. Asp.Net Core Linux环境下 找不到配置文件、静态文件的问题

    直接发布会找不到配置文件,和静态文件.需要先cd到项目文件夹,然后在发布.

  9. FPGA上外挂DDR2&DDR3&MIG IP的使用记录

    前言 当需要大容量数据存储及处理的时候,FPGA内部自带的存储资源是远远不够的,所以问题来了,怎么使用外带的DDR3? 首要问题在于DDR3是什么?有没有协议?当然只是需要用Xilinx MIG IP ...

  10. iOS圆弧渐变进度条的实现

    由于项目需要一个环形渐变进度条显示课程,这方便网上的确有很多相关资料但是,都是比较零散的而且,大多数只是放一堆代码就算完了.这里我想详细写一篇我自己实现这个进度条的过程. 实现一个圆弧进度条主要分为三 ...