高等工程数学 —— 第三章(2)奇异值分解和A的加号逆

文章目录

  • 高等工程数学 —— 第三章(2)奇异值分解和A的加号逆
    • 奇异值分解
    • 广义逆矩阵
    • A + A^{+} A+的直接计算方法
      • 奇异值分解计算 A + A^{+} A+
      • 满秩分解计算 A + A^{+} A+
    • A + A^{+} A+的迭代计算方法
    • A + A^{+} A+的基本性质
    • 广义逆矩阵的应用

奇异值分解

首先来看什么是奇异值


也别管什么原理了,直接看方法和例题。盘它!

奇异值分解步骤:

  • 这里就是先求 A H A A^{H}A AHA的特征值,然后求其特征向量并将每一个特征向量进行单位化得 V V V
  • 然后看有几个非零特征向量就分出来几列当 V 1 V_1 V1​
  • 求出 U 1 U_1 U1​后将其补全成方阵,因为是酉矩阵所以补的列向量要满足都两两单位正交
  • 最后求 A A A时对 Σ \Sigma Σ的补全满足矩阵乘法就行,让补全后的 Σ \Sigma Σ矩阵的行数和 U U U矩阵的列数一样;补全后的 Σ \Sigma Σ矩阵的列数和 U U U矩阵的行数一样

看道例题好喽!

首先求 A H A A^{H}A AHA的特征值和特征向量来求 V V V

因为非零特征值有两个,所以取 V V V的前两列作为 V 1 V_1 V1​来参与 U 1 U_1 U1​的计算。


这里 U 1 U_1 U1​已经是方阵了,所以不用补全了。对 Σ \Sigma Σ的补全满足矩阵乘法即可。

广义逆矩阵

  • 这里其实就是当 A − 1 A^{-1} A−1不存在时,我们做的一个逆矩阵的推广

  • 这里的四个公式我们把 X X X当作 A − A^{-} A−代入发现都成立。可见推广后的逆矩阵 A + A^{+} A+满足这四个条件用来某种程度上代替 A − A^{-} A−

A + A^{+} A+的直接计算方法

奇异值分解计算 A + A^{+} A+

  • 注意奇异值分解得到的 U U U和 V V V可以是不一样的,但是相乘后得到 A + A^{+} A+的值是一样的。

例:

  • 这里看 U U U补全的最后一列可知,我们补全的原则就是要满足两两正交的单位向量。

  • 这里 V V V的列向量都是标准化后的单位向量

  • 最后解出来的 A + A^{+} A+值是唯一的

满秩分解计算 A + A^{+} A+

对于行满秩矩阵和列满秩矩阵而言:

  • 对于行满秩矩阵而言,满秩分解后 F = E F = E F=E,代入得上式。
  • 对于列满秩矩阵而言,满秩分解后 G = E G = E G=E,代入得上式。
  • 也不用刻意记,满秩分解完顺水推舟的化成这样了

例题:

  • 注意这里要先求出括号里的 ( G G T ) − 1 (GG^T)^{-1} (GGT)−1和 ( F T F ) − 1 (F^TF)^{-1} (FTF)−1的值,然后将四个矩阵从左到右依次相乘。
  • 看似满秩分解法求 A + A^{+} A+比奇异值分解求 A + A^{+} A+要简单一点,但是要求两个逆矩阵然后矩阵相乘到最后一大堆分数,恶心的很,呜呜呜呜。
  • 注意求逆的时候只做初等行变换,笔者做题的时候傻了,有行有列的最后一直算不对,服了。

A + A^{+} A+的迭代计算方法

这部分我只能说不考,哈哈哈哈。这迭代应该是那种写计算机程序用的,而且这玩意书上一道题的没有。所以我们就简单列一下公式算了,嘻嘻。

Greville递推法


A + A^{+} A+的基本性质

A + A^{+} A+存在且唯一

同时满足以下性质:


广义逆矩阵的应用

例:

先用满秩分解求得 A + A^{+} A+

高等工程数学 —— 第三章(2)奇异值分解和A的加号逆相关推荐

  1. 高等工程数学 —— 第四章 (1)线性方程组的直接解法与广义逆矩阵求解矛盾方程组

    高等工程数学 -- 第四章 (1)线性方程组的直接解法与广义逆矩阵求解矛盾方程组 文章目录 高等工程数学 -- 第四章 (1)线性方程组的直接解法与广义逆矩阵求解矛盾方程组 线性方程组的直接解法 Ga ...

  2. 高等工程数学 —— 第一章 (1)距离与范数

    前言 研一生活开始了,看了大家对我之前博客的鼓励让我知道写博客是一件多么有意义的事情.写这些让我遇见许多陌生的有缘人,有老骥伏枥的大叔.也有可爱温暖的学妹-- 这里将高等工程数学的笔记留给不爱吃香菜的 ...

  3. 高等工程数学 —— 第二章 (1) Smith标准型与Jordan标准型

    高等工程数学 -- 第二章 (1) Smith标准型与Jordan标准型 文章目录 高等工程数学 -- 第二章 (1) Smith标准型与Jordan标准型 不变因子与Smith标准型 行列式因子 行 ...

  4. 【高等工程数学】南理工研究生课程 突击笔记3 不变因子与Jordan标准型

    高等工程数学 突击笔记3 文章目录 高等工程数学 突击笔记3 一.标准型 λ矩阵 行列式因子D 不变因子d 初等因子 Jordan标准型 二.盖尔圆 特征值隔离 总结 第二章内容大致分成三个部分 标准 ...

  5. 【高等工程数学】南理工研究生课程 突击笔记4 幂迭代

    幂迭代法,和逆幂迭代法 文章目录 幂迭代法,和逆幂迭代法 写在前面 一.幂迭代法 二.逆幂迭代法 三.规范化迭代方式 四.A分解 例 总结 写在前面 承接笔记3,先补一个盖尔圆的题目 如果特征值是复数 ...

  6. 【高等工程数学】南理工研究生课程 突击笔记2 谱半径,条件数

    高等工程数学 突击笔记2 文章目录 高等工程数学 突击笔记2 写在前面 一.谱半径 什么是正规矩阵 矩阵级数敛散性判断 二.条件数 举个栗子 条件数求误差 总结 写在前面 提示:这里可以添加本文要记录 ...

  7. 华科计算机博士2201,2018年华中科技大学2201高等工程数学考博大纲

    资源大小:0.1-10.0 MB 资源类型: rar 发布时间: 2020/5/14 23:35:45 资源评分: ★★★ 资源简介: 2018年华中科技大学2201高等工程数学考博大纲 下载流程: ...

  8. 【高等工程数学】南理工研究生课程 突击笔记1 距离与范数1

    高等工程数学 突击笔记1 距离与范数1 文章目录 高等工程数学 突击笔记1 距离与范数1 写在开头 一.什么是范数 二.诱导范数 总结 写在开头 高等工程数学这个课真的很恶心,感谢B站UP ASH丶零 ...

  9. 工程伦理第三章学习笔记2020最新

    工程伦理第三章学习笔记2020最新 因为之前自己在网上找答案总是觉得费劲,一道一道的找,很慢,突然找到了前两章的答案,感觉有一种前人种树后人乘凉的感觉,于是自己在艰难找完第三章习题并全对的情况下,将题 ...

最新文章

  1. linux命令telnet
  2. shell条件测试操作 if分支 for循环 while 循环
  3. auuc 评估指标_广告中增益模型理解
  4. html坐标定位图解,HTML5地理定位实例
  5. jasmine fixture.debugElement.query的单步调试
  6. STL 容器中的元素必须满足的条件
  7. 码农回家过个年,真心不容易!
  8. C语言数组及相关函数
  9. SAP License:FI-CO集成
  10. Ubuntu 16.04安装 Nmap 6.46.1
  11. Oracle 中 start with 递归查询、case when 条件表达式、rowid 伪列去重
  12. 计算机io设备是cpu控制的吗,操作系统--IO设备
  13. 分享几个免费的根据IP 获取地址的链接
  14. 计算机网络知识点概括
  15. 7-14 电话聊天狂人 (25分) 【map】
  16. 揭秘世界福利的主要来源
  17. redis锁和分布式锁的实现
  18. 库存管理软件免费版软件哪个好用?
  19. mysql内存机制_浅析Oracle 的体系架构及内存分配机制
  20. HTML:初次见面,先做个自我介绍

热门文章

  1. 算法唯手熟尔4月(PAT剩余清单 or leetcode)---希望可以日更
  2. 由南非世界杯足球赛的用球“JABULANI”说起
  3. Leetcode 738. Monotone Increasing Digits
  4. python画素描画_Python素描画的两种程序解析
  5. MATLAB_GUI项目 自然场景下限速标志识别
  6. 惯性导航之四元数与欧垃角互转(五)
  7. mnt/hgfs中不显示共享文件夹解决方法
  8. java并发编程第一课 线程的创建、停止和状态变更
  9. Python 将控制台输出另存为日志文件
  10. 查看数据库日志文件及路径