伪逆矩阵是逆矩阵的广义形式。由于奇异矩阵或非方阵的矩阵不存在逆矩阵,但在matlab里可以用函数pinv(A)求其伪逆矩阵。基本语法为X=pinv(A),X=pinv(A,tol),其中tol为误差,pinv为pseudo-inverse的缩写:max(size(A))*norm(A)*eps。函数返回一个与A的转置矩阵A' 同型的矩阵X,并且满足:AXA=A,XAX=X.此时,称矩阵X为矩阵A的伪逆,也称为广义逆矩阵。pinv(A)具有inv(A)的部分特性,但不与inv(A)完全等同。 如果A为非奇异方阵,pinv(A)=inv(A),但却会耗费大量的计算时间,相比较而言,inv(A)花费更少的时间。

机器视觉中,finite camera的投影矩阵P为一个3*4,秩为3 的矩阵,矩阵P的伪逆如下:

矩阵P的伪逆为:

如果A列满秩,那么pinv(A)=(A'*A)^{-1}*A'。

如果A行满秩,那么pinv(A)=pinv(A')'。
如果秩亏损,那么只好先做奇异值分解A=UDV',U,V是正交阵,D是对角阵。
然后取对角阵S,如果D(i,i)=0,那么S(i,i)=0,如果D(i,i)<0或者D(i,i)>0,那么S(i,i)=1/D(i,i)。于是pinv(A)=VSU'。
-----------------------------------------------------------------------------------------------------------

原文:https://blog.csdn.net/u014260892/article/details/38581175

两种求矩阵伪逆的方法相关推荐

  1. matlab求矩阵违逆,两种求矩阵伪逆的方法

    伪逆矩阵是逆矩阵的广义形式.由于奇异矩阵或非方阵的矩阵不存在逆矩阵,但在matlab里可以用函数pinv(A)求其伪逆矩阵.基本语法为X=pinv(A),X=pinv(A,tol),其中tol为误差, ...

  2. 【MATLAB】几种求矩阵伪逆的方法

    几种求矩阵伪逆的方法 A*x=y当A不可逆如何解x 求伪逆有五种方法,左右同时乘以A','chol()分解,qr() 分解,svd()分解和pinv()求伪逆,下面来比较一下哪个最优秀. 如果A是个3 ...

  3. 两种求集合全部子集的方法

    如果我们有一个求集合的所有子集(包括集合自身)的需求,即有一个集合s,包括两个元素 <a,b>,则其所有的子集为<a,ab,b>. 不难求得,子集个数sn与原集合元素个数n之间 ...

  4. 模为2的逆元是什么_两种求模m逆元的方法

    在a|b(a能整除b)的前提下,计算(b/a)mod m的时候转化为 计算(b*x)mod m ; 这时的x就是a的逆元(a模m的逆元): 此时x满足  (a*x mod m == 1): 这个x的求 ...

  5. 求矩阵伪逆的matlab方法,手把手教学

    1.定义矩阵a 2.求b,c b为a的线性无关的两列 c为a的行最简的线性无关的两行 3.求b的左逆和c的右逆 左逆和右逆的具体见一下例题: 4.求出结果即可 或者一步到位 有时候会出现小数,转变表达 ...

  6. 求矩阵的逆的三种方法

    我们知道求矩阵的逆具有非常重要的意义,本文分享给大家如何针对3阶以内的方阵,求出逆矩阵的3种手算方法:待定系数法.伴随矩阵法.初等变换法(只介绍初等行变换) 待定系数法求逆矩阵 1 首先,我们来看如何 ...

  7. Matlab求矩阵的逆(三种方法)

    Matlab求矩阵的逆(三种方法) 说明:若所求矩阵为非奇异矩阵(可逆矩阵),则可以精确求得其逆矩阵:若所求矩阵为奇异矩阵,则所求出的逆矩阵是近似的(不精确). 下面以矩阵A为例. inv()方法 A ...

  8. Matlab两种求相位的方法

    在matlab中有两种求相位方法. 1.使用phase函数. 比如a = 1+j*2 phase(a),则可以就a的相位.这个函数取实部为x轴,虚部为y轴. 2.使用函数atan2. atan2(y, ...

  9. 机器人运动学---雅克比矩阵伪逆

    Problem 使用雅克比矩阵可以通过关节速度,得到EE的速度和角速度.如果雅克比矩阵是方阵且满稚可逆,这时,可以通过EE的线速度和角速度求出关节的速度. 以上是比较理想的情况. 但是,实际情况,往往 ...

  10. C语言求矩阵的逆(高斯法)

    初等变换法是常用的矩阵求逆方法之一 相对于伴随法,初等行变换法有着较低的时间复杂度,可以进行相对高维的矩阵运算,但同时也会损失一点点精度. 伴随法可参考之前的博客:C语言求矩阵的逆(伴随法) 目录 数 ...

最新文章

  1. 计算机二级公共基础知识总结百度云,计算机二级公共基础知识总结详细版本[精]...
  2. python模板是什么意思_什么是最快的Python模板系统?
  3. POJ 3734 Blocks
  4. wpf listview mysql_Kivy:使用MySQL的Kivy页面的Listview实现
  5. 说说我为什么看好Spring Cloud Alibaba
  6. MFC六大核心机制之五、六:消息映射和命令传递
  7. Python中os库的使用
  8. python123输出hello world_Python基础:输入与输出
  9. 页面错误!请稍后再试_微信内嵌H5页面授权和分享
  10. DriverMessageBean配置详解
  11. swt中关于Text.setSelection()的记录
  12. python静态检查工具哪个好用_介绍几款 Python 类型检查工具
  13. python爬取股票历史数据_爬取股票历史数据并绘制K线图
  14. 【信息系统项目管理师】2018下半年系统集成项目管理工程师案例分析
  15. (usb转串口)CH340驱动安装失败的解决方案
  16. SBG Ellipse2最强替代型号推荐?AHRS INS/GNSS
  17. C++修正儒略日转UTC时间
  18. 景元利老师--沪师经纪
  19. 嵌入式开发要学习哪些知识
  20. IBM笔记本选购注意事项(笔记本是一个大物品,所以会三思而后行的)

热门文章

  1. 正则表达式汇总--小鱼儿
  2. 哪家的云游戏服务器好?如何选择云游戏服务器?
  3. SCHMERSAL AZM 161 Z ST1-AS PT西溪,且留下
  4. windows+clion Process finished with exit code -1073741515 (0xC0000135)
  5. 2023厦门大学计算机考研信息汇总
  6. 前端 环形统计_用Canvas+js 实现简单圆形百分比统计图表
  7. 英特尔神经计算棒二代(NCS2)安装、配置及使用
  8. 全流程东方时尚C1考试经历
  9. 杭州地铁四期尘埃落定和我预测差异,看地铁空白地,遥想地铁五期规划
  10. IT公民:293个公司人压力和心理调查