伪逆矩阵的介绍&定义

矩阵的逆&伪逆

矩阵的伪逆推导

定义：令A是任意mXn矩阵，称矩阵G是A的广义逆矩阵，若G满足下述条件（Moore-penrose条件）：
（1）GAG = G；
（2）AGA = A；
（3）AG为hermitian矩阵，即（AG）^H=AG；
（4）GA为hermitian矩阵，即（GA）^H=GA；

伪逆矩阵的求法

直接求解：

求导，令导数为0，结果如下: InvA=(ATA)−1ATInvA=(A^TA)^{-1}A^TInvA=(ATA)−1AT

%直接求伪逆
InvA = inv(A’A)A’;

SVD求解

SVD分解求伪逆
- 原理和公式：
  - SVD分解得到的矩阵:U和V是正交阵，S是对角阵
  - 正交阵的逆=转置
  - 对角阵的逆=非零元素求倒
- Step1: 求解A的SVD分解
  - [U,S,V] = svd(A); % A = USV’
- Step2: 将S中的非零元素求倒
  - T=S;
  - T(find(S~=0)) = 1./S(find(S~=0));
- Step3: 求invA
  - svdInvA = V * T’ * U’;

QR求解

QR分解求伪逆,适用于稀疏矩阵
- 原理：
  - QR分解得到的矩阵:Q是正交阵，R是非奇异上三角阵
  - 正交阵的逆=转置
  - 上（下）三角矩阵的逆也仍然是上（下）三角矩阵。
  - 不必用高斯消去法，向前替换法解方程。
  - 但是具体的我不知道怎么用程序来写，这里仍旧用了matlab的函数。
- 公式：
  - [Q,R] = qr(A);
  - InvR = inv(R’*R)*R’;
  - qrInvA =InvR*Q’;

伪逆矩阵的应用

(1)信号的检测干扰消除

a=floor(10*rand(4,3))
a =
7 4 2
7 6 6
1 7 6
4 7 1

b=inv(a’*a)*a’
b =
0.1018 0.0650 -0.0959 -0.0180
-0.0263 -0.0767 0.0578 0.1658
-0.0296 0.1149 0.0903 -0.1719

b*a
ans =
1.0000 -0.0000 0.0000
-0.0000 1.0000 -0.0000
-0.0000 -0.0000 1.0000

伪逆矩阵在维基百科中的详细介绍

参考链接

伪逆矩阵
伪逆矩阵（广义逆矩阵）

伪逆矩阵（pseudo-inverse）相关推荐

逆矩阵、伪逆矩阵：数据的压缩和复原
逆矩阵.伪逆矩阵.数据的压缩和复原:这一块知识虽然很简单,但在光学各种实验情况下经常用到,特此总结. 矩阵的乘法 (观测矩阵):C = A* B C:M*1 A:M*N B:N*1 A可以理解为一个转 ...
pinv--求矩阵的伪逆矩阵
pinv--求矩阵的伪逆矩阵 [功能简介]用于求矩阵的伪逆矩阵. [语法格式] 1．B=pinv(A) 函数返回矩阵A的伪逆矩阵.如果矩阵A是可逆(非奇异)的,那么pinv(A)与inv(A)的结果是 ...
线性代数笔记：逆矩阵及伪逆矩阵，最小二乘估计，最小范数估计
逆矩阵及伪逆矩阵,最小二乘估计,最小范数估计逆矩阵的概念我们为什么需要逆矩阵? 伪逆矩阵和最小二乘估计正则化求伪逆矩阵逆矩阵的概念矩阵A的逆矩阵(matrix inversion)记作A−1 ...
【转录调控网络】典型的基因转录调控网络推导方法——伪逆矩阵模型
基因转录调控网络推导方法--伪逆矩阵模型在基因调控网络推导中,使用到的基因芯片数据通常具有样本个数少(通常小于10)而基因数目大(通常大于1000)的局限性,也就是说,实验样本个数远远小于基因 ...
MATLAB与线性代数--逆矩阵与伪逆矩阵
求逆矩阵与伪逆矩阵矩阵A的逆矩阵用表示,并且满足下面的关系: 看下面的矩阵方程: 如果A的逆矩阵存在,那么解可以写成: 在MATLAB输入下面的命令就可以计算矩阵A的逆矩阵但是逆矩阵并不一定存在 ...
关于数学库Eigen两个常用函数实现（伪逆矩阵和协方差矩阵)
关于数学库Eigen两个常用函数实现(伪逆矩阵和协方差矩阵) cov pinv 实现代码
一分钟读懂：矩阵的特征值分解、奇异值分解和伪逆矩阵
通过把矩阵运算分解成多个矩阵的乘法,可以简化矩阵运算,也可发现对应线性变换的一些内在规律和特性.根据不同的目的,有不同的分解策略.本文我们讨论最常用的特征值分解和奇异值分解. 1. 矩阵的乘方运算定 ...
奇异矩阵，非奇异矩阵，伪逆矩阵
奇异矩阵就是Singular Matrix 的中文翻译. Singular 就是唯一的,可以想成是单身狗,所以他没有对象逆矩阵. Non-singular的非奇异矩阵就是Couple 有逆矩阵. 奇 ...
奇异矩阵和非奇异矩阵，伪逆矩阵
奇异矩阵: 奇异矩阵是线性代数的概念,就是对应的行列式等于0的矩阵. 奇异矩阵的判断方法:首先,看这个矩阵是不是方阵(即行数和列数相等的矩阵.若行数和列数不相等,那就谈不上奇异矩阵和非奇异矩阵). 然 ...
c语言矩阵求伪逆算法pinv,pinv--求矩阵的伪逆矩阵
pinv--求矩阵的伪逆矩阵 [功能简介]用于求矩阵的伪逆矩阵. [语法格式] 1．B=pinv(A) 函数返回矩阵A的伪逆矩阵.如果矩阵A是可逆(非奇异)的,那么pinv(A)与inv(A)的结果是 ...

伪逆矩阵（pseudo-inverse）