原文 | https://mp.weixin.qq.com/s/zdQttJfuubyztiVplScbwA

对称矩阵

  对称矩阵是最重要的矩阵之一,对于对称矩阵来说,A=AT。矩阵的特殊性也表现在特征值和特征向量上,比如马尔可夫矩阵的有一个值为1的特征值,对称矩阵的特征值又有哪些特性呢?

  本文的相关知识:

    正交向量和正交矩阵  (线性代数20——格拉姆-施密特正交化)

    投影矩阵 (线性代数18——投影矩阵和最小二乘)

    复数 (闲话复数(1))

谱定理

  对于实对称矩阵来说,它的特征值也为实数,并且能够挑选出完全正交的特征向量。

  单位矩阵是对称矩阵,它的特征值都是1,并且单位矩阵的每一个列向量都是特征向量,它们完全正交,因此单位矩阵肯定符合实对称矩阵特征值和特征向量的性质。

  P是投影矩阵也是单位矩阵,x是一个二维向量,如果x在平面的投影是x本身,即Px=x,那么平面内的所有向量都是P的特征向量。更一般的情况是,在重特征值的情况下,可能一个平面内的所有向量都能作为特征向量,因此我们说实对称矩阵“能够挑选出完全正交的特征向量”,下面是一个例子:

  A的特征值全部是λ=a,对于任何向量都有Axx,因此任何向量都是特征向量,但这些特征向量并不都是正交的,所以是从中选出一套正交向量。

  如果A有n个线性无关的特征向量,那么A可以对角化为A=S∧S-1,如果A是对称矩阵,那么A对角化后有更好的性质:

  QA的特征向量矩阵,同时也是正交矩阵,列向量是标准正交基。对于一个列向量标准正交的矩阵来说,矩阵的逆等于矩阵的转置,因此:

  上式是说,如果给定一个对称矩阵,那么这个矩阵就可以分解成正交矩阵乘以特征值矩阵再乘以正交矩阵的转置的形式,这种分解在数学上称为“谱定理”,将特征值的集合视为谱,力学上称为“主轴定理”。

  谱定理展示了对称矩阵的对称性,Q∧QT的转置还是原矩阵:

  是对角矩阵,它的转置还是

为什么是实特征值?

  矩阵的特征值可能为虚数,但实对称矩阵的特征值一定是实数,这又是什么原理?

  先解释一下共轭复数(conjugate complex number):两个实部相等,虚部互为相反数的复数互为共轭复数。当虚部不为零时,共轭复数就是实部相等,虚部相反;如果虚部为零,其共轭复数就是自身。复数z的共轭复数用z上面加一横表示。

  “轭”的本意是两头牛背上的架子,轭使两头牛同步行走。共轭是指按一定的规律相配的一对。

  如果实矩阵A有特征值λ和特征向量x,则Axx。如果λ是复数,则λ的共轭复数满足:

  如果等号两侧同时转置:

  λ是对角矩阵,其共轭仍然是对角矩阵,因此:

  两侧同时乘以x

  另一方面,将Axx两侧同时乘以x共轭的转置:

  现在假设A是对称矩阵,则①和②相等,即:

  根据共轭复数的定义,如果一个复数的共轭等于这个数本身,那么其虚部为0,即这个数是实数。

  现在的问题是如何证明

  对于虚数来说,i2=-1,(bi)2=-b2。对于复数来说,z=a+bi来说,它的模几何意义是复平面上一点(a,b)到原点的距离,模长的计算公式是:

  因此:

  对于复矩阵来说,若A是共轭对称复矩阵,即,那么上面的推导仍然成立,A的特征值也是实数。

朝向正交向量的投影矩阵

  对于一个实对称矩阵A=AT,可以写成下面的形式:

  根据投影矩阵的定义:向量b的在向量a上的投影是一个矩阵作用在b上得到的,这个矩阵就叫做投影矩阵(Projection Matrix),用大写的P表达:

  由于Q中的向量是正交向量,因此:

  所以qkqkT是某个向量在qk上的投影矩阵,因此每一个对称矩阵也是朝向正交向量的投影矩阵的线性组合。

特征值的符号

  我们已经知道对称矩阵的特征值是实数,下一个问题是弄清它们的符号,这对微分方程来说意味着状态是否稳定。

  我们可以通过计算的方式求解特征值,然后回答特征值的符号问题,但对于一个大型矩阵来说,通过计算det(AI) = 0来求解特征值并不容易。值得庆幸的是,对于对称矩阵来说,主元和特征值存在着相关性:主元和特征值的个数一样,且正负主元的个数都和正负特征值的个数相同。

正定矩阵

  正定矩阵是一类特殊的实对称矩阵,如果一个矩阵M满足对于任何非零向量z,都有zTMz> 0,那么这个矩阵是正定矩阵。

  正定矩阵有很多重要的性质,其中一个是:正定矩阵的特征值和主元都是正数。

  来看一个正定矩阵:

  首先A是一个对称矩阵,现在来计算一下它的主元。可以通过化简行阶梯矩阵的形式求得主元,在经过变换后,矩阵表示的“数表”改变了,但是如果将矩阵看方程组,那么方程组的本质没有变,可以将初等变换看成方程组的消元过程。

  两个主元是5和11/5。另一种方式或许更为简单,原矩阵中转换成上三角矩阵的时候,一定能变成下面的形式:

  它的行列式是主对角线元素的乘积,行列式的值又和原矩阵的行列式相等,因此a=det(A)/5=11/5。

  类似地,可对角化的矩阵可也以化成对角元素都是其特征值的对角矩阵,而行列式的值不变,对角矩阵的行列式就是对角元素相乘,因此A的行列式也等于A的特征值的乘积。

  特征值:

与行列式的关系

  正定矩阵的主元和特征值都是正数,因此可以确定其行列式也是正数(行列式等于主元的乘积,也等于所有特征值的乘积),但行列式是正数的矩阵不一定是正定矩阵,比如下面这个矩阵,虽然行列式是正值,但并不是一个正定矩阵,甚至都不是对称矩阵:

  如果把行列式作为正定矩阵的判定依据,则对于n阶矩阵来说,需要矩阵左上角的所有k×k (1≤k≤n)子行列式均为正,才能判定矩阵是正定矩阵。

正定矩阵的性质

  正定矩阵都是可逆的。

  矩阵可逆的条件是行列式不等于0,行列式等于特征值的乘积,正定矩阵的性质又规定特征值是正数,因此正定矩阵的行列式一定大于0,是可逆矩阵。

  

  只有正定的投影矩阵才是单位矩阵。

  如果P是投影矩阵,那么P的特征值要么是0,要么是1。如果P是正定的,那么根据定义,它的特征值只能是1,特征值矩阵是单位矩阵,因此:

  

  如果D是一个对角元素都是正数的对角矩阵,那么D一定是个正定矩阵。

  对角矩阵肯定是对称的,对于任何非零向量x来说:

  满足正定矩阵的定义。

  

  A是正定矩阵,则A的逆矩阵也是正定矩阵。

  首先证明矩阵A的逆是对称矩阵。因为A是正定的,所以:

  再证明xTA-1x > 0

  A是正定矩阵,对于任意向量u来说,uTAu > 0,因此xTA-1x > 0,A-1也是正定矩阵。

  

  两个正定矩阵的和是正定矩阵。

  

  正实数与正定矩阵的乘积是正定矩阵。


  作者:我是8位的

  出处:https://mp.weixin.qq.com/s/zdQttJfuubyztiVplScbwA

  本文以学习、研究和分享为主,如需转载,请联系本人,标明作者和出处,非商业用途!

  扫描二维码关注作者公众号“我是8位的”

线性代数笔记27——对称矩阵及正定性相关推荐

  1. matlab向量与x正方向的夹角_MIT—线性代数笔记25 对称矩阵和正定性

    第25讲 对称矩阵和正定性 Symmetric matrices and positive definiteness 网易公开课​open.163.com 进入第三单元学习,主题是正定矩阵及其应用.对 ...

  2. MIT线性代数笔记二十五讲 对称矩阵和正定性

      对称矩阵是最重要的矩阵之一.那么它的特征值和特征向量有什么特别之处吗? 文章目录 1. 对称矩阵 Symmetric matrices 2. 实特征值 Real eigenvalues 3. 正定 ...

  3. 线性代数之——对称矩阵及正定性

    当 AAA 是对称的时候,Ax=λxAx=\lambda xAx=λx 有什么特殊的呢? 1. 对称矩阵的分解 A=SΛS−1A = S\Lambda S^{-1}A=SΛS−1 AT=(S−1)TΛ ...

  4. CP.26对称矩阵及正定性

    本部分从对称矩阵的特征值.特征向量为出发点,并引出正定矩阵 1.对称矩阵 对称矩阵有如下两个性质: (1)A=AT(1)\textbf{A}\:=\:\textbf{A}^T(1)A=AT (2)有正 ...

  5. 26对称矩阵及正定性

    一.知识概要 本节从对称矩阵的特征值,特征向量入手,介绍对称矩阵在我们之前学习的一些内容上的特殊性质.并借此引出了正定矩阵. 二.对称矩阵 正如我们之前学习的很多特殊矩阵一样(如马尔科夫矩阵),对称矩 ...

  6. MIT公开课18.06 Gilbert Strang 线性代数 笔记3 - 正定矩阵及其应用

    文章目录 第26讲:对称矩阵及正定性 对称矩阵 性质 描述 证明 谱定理,对称矩阵的分解 对称矩阵特征值的符号 正定性 定义 例子 第27讲:复数矩阵和快速傅里叶变换 复数矩阵运算 计算复向量的模 计 ...

  7. 矩阵论(五):矩阵的正定性

    矩阵论专栏:专栏(文章按照顺序排序) 矩阵的正定性是矩阵理论的基础,在机器学习的很多算法中都有它的身影.正定矩阵是很基础的工具,当涉及到诸如协方差矩阵.海森矩阵.不等式时就可能会用到正定性.另外,控制 ...

  8. 线性代数——矩阵正定性及二次型的矩阵表示

    一. 判断矩阵正定性的方法(4种) 1.矩阵所有特征值为正 即λi>0 2.矩阵的所有主元为正数 3.矩阵的顺序主子式均为正数 4.矩阵表示的二次型为正 二.二次型矩阵形式及代数形式的转化 二次 ...

  9. [笔记][总结] MIT线性代数 Gilbert Strang 对称矩阵

    作者水平有限,欢迎大家提出文中错误 正定性与对称矩阵 对称矩阵 对称矩阵的对角化 正定性 正定矩阵判据 正定矩阵的性质 正定性与最小二乘法 二次型 对称矩阵的LU分解与二次型的配方 连续多元函数在某点 ...

最新文章

  1. JavaScript_上
  2. ContentProvider再探——Document Provider
  3. springboot配置连接rds_java – 在AWS-EC2上运行的SpringBoot应用程序无法连接到MySQL AWS-RDS数据库...
  4. 前端学习(509):水平和垂直居中第二种方式
  5. linux创建了目录复制不存在,Linux:复制并创建目标目录,如果它不存在
  6. React与ES6(一)开篇介绍
  7. 助力南京打造创新名城 第三届未来网络发展大会将召开
  8. 游戏开发之类和对象的基本概念(C++基础)
  9. 帆软决策报表全页面刷新或局部组件定时刷新
  10. spring的核心模块
  11. android打印机字体大小,打印自定义文档  |  Android 开发者  |  Android Developers
  12. 利用定义式求圆锥曲线方程
  13. pip使用清华镜像源安装库
  14. c语言编译defined,#if defined(__GNUC__)的意思是不是如果使用的是GCC编译器?
  15. drawableTop,button图片显示在文字上面,文字不显示的问题
  16. [spring boot] 我的图床我做主,自己动手写一个shareX图床
  17. 教程:BootToGecko系统虚机与实机开发环境搭建及FirefoxOS设备使用和开发网站推荐;
  18. [视觉概述] 机器视觉应用方向、项目流程及学习思路总结
  19. 问题1:修改日志文件(redo log)的位置的方法
  20. typescript完成日期转换

热门文章

  1. 目前网罗的各种有关Kinect的有趣应用
  2. 短视频是屌丝逆袭的一个绝好的机会
  3. MySQL高集群方案
  4. HDU 2212 DFS (伪·DFS)
  5. *迟来的爱*——《Foursquare》应用源码学习(一) 下载、编译、运行
  6. Hibernate中发生Session is closed 的另一种可能
  7. 幼儿园计算机课件制作评比表,幼儿园课件制作基本原则与评价标准_表格评价.doc...
  8. 如何选择U盘加密软件?
  9. 教育培训网站源码,大班课/小班课线上直播都能用
  10. Altium Designer 鼠标移动到元件高亮显示连接网络