码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • 矩阵分析与应用-18-Moore-Penrose逆矩阵02


    Moore-Penrose 逆矩阵的计算

    方程求解法

    第一步:分别求解矩阵方程得到X^H,Y.

    AA^HX^H=A

    A^HAY=A^H

    第二步:计算广义逆矩阵A^\dagger=XAY

    若矩阵A 为Hermitian矩阵,上述方法可以简化,因为以上两个矩阵方程等价为一个矩阵方程

    A^2X^H=A

    若A^H=A,Moore-Penrose逆矩阵可以计算为

    A^\dagger=XAX^H

    算法1

    步骤1计算矩阵B=AA^H。

    步骤2求解矩阵方程B^2X^H=B得到矩阵X^H。

    步骤3计算B的Moore-Penrose逆矩阵B^\dagger=(AA^H)^\dagger=XBX^H。

    步骤4计算矩阵A的Moore-Penrose逆矩阵A^\dagger=A^H(AA^H)^\dagger=A^HB^\dagger。

    算法2

    步骤1计算矩阵B=AA^H。

    步骤2求解矩阵方程B^2X^H=B得到矩阵X^H。

    步骤3计算的Moore-Penrose逆矩阵B^\dagger=(A^HA)^\dagger=XBX^H。

    步骤4计算矩阵A的Moore-Penrose逆矩阵A^\dagger=(A^HA)^\dagger A^H=B^\dagger A^H。

    KL分解法

    若A=KL是矩阵Amxn的满秩分解,则

    G=L^H(K^HAL^H)^{-1}K^H

    递推法

    算法3

    初始值A_1^\dagger=a_1^\dagger=(a_1^Ha_1)^{-1}a_1^H。

    递推令k= 2,3,…- ,n,进行以下计算:

    d_k=A^\dagger_{k-1}a_k

    b_k=\begin{cases} (1+d_k^Hd_k)^{-1}d_k^HA_{k-1}^\dagger & d_k^Hd_k\neq -1 \\ (a_k-A_{k-1}d_k)^\dagger &d_k^Hd_k= -1 \end{cases}

    A_k^\dagger=\begin{bmatrix} A^\dagger_{k-1}d_kb_k\\ b_k \end{bmatrix}

    迹方法

    步骤1计算矩阵B=AA^H。

    步骤2令C_1=I。

    步骤3计算

    C_{i+1}=\frac{1}{i}tr(C_iB)I-C_iB,i=1,2,...,r-1

    步骤4计算

    A^\dagger=\frac{r}{tr(C_iB)}C_iA^T

    注:C_{i+1}B=O,tr(C_iB)\neq 0.

  • 相关阅读:
    【面试题】——JavaIO篇(23题)
    Blender做一个小凳子学习笔记
    匈牙利算法讲解
    Codeforces Round #816 (Div. 2)
    【Linux】基本指令(三)
    分享几个常问的vue3面试题!
    【JavaGuide学习笔记】Day.1
    Mybatis(整合spring)
    友盟+|如何通过阿里云Flink+Paimon实现流式湖仓落地方案
    鸿蒙原生应用开发-DevEco Studio本地模拟器的使用
  • 原文地址:https://blog.csdn.net/Xiao__fly/article/details/126280172
  • 最新文章
  • 攻防演习之三天拿下官网站群
    数据安全治理学习——前期安全规划和安全管理体系建设
    企业安全 | 企业内一次钓鱼演练准备过程
    内网渗透测试 | Kerberos协议及其部分攻击手法
    0day的产生 | 不懂代码的"代码审计"
    安装scrcpy-client模块av模块异常,环境问题解决方案
    leetcode hot100【LeetCode 279. 完全平方数】java实现
    OpenWrt下安装Mosquitto
    AnatoMask论文汇总
    【AI日记】24.11.01 LangChain、openai api和github copilot
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1
正则表达式工具 cron表达式工具 密码生成工具

京公网安备 11010502049817号