逆矩阵、伪逆矩阵、数据的压缩和复原:这一块知识虽然很简单,但在光学各种实验情况下经常用到,特此总结。

矩阵的乘法 (观测矩阵):C = A* B

C:M*1

A:M*N

B:N*1

A可以理解为一个转化矩阵,或者说 观测矩阵。

矩阵B在观测矩阵A上,观察到的效果是矩阵C

这个思想在物理思想上非常重要,尤其是光学。

实际生活中,往往我们需要通过观测矩阵(A)和观察到的结果(C),来求出原矩阵(B)

也就是解非线性齐次方程,M个方程,N个未知数。

求逆矩阵当矩阵A必须是方阵的时候,A才是可逆矩阵,这时M=N,有唯一解,B = inv(A) * C,matlab中测试如下。

%测试逆矩阵
A = rand(10,10);
B = rand(10,1);
C = A*B;
BB = inv(A)*C;
B'
BB'结果:ans =0.0986    0.1420    0.1683    0.1962    0.3175    0.3164    0.2176    0.2510    0.8929    0.7032ans =0.0986    0.1420    0.1683    0.1962    0.3175    0.3164    0.2176    0.2510    0.8929    0.7032

求伪逆矩阵:伪逆矩阵是逆矩阵的广义形式。由于奇异矩阵或非方阵的矩阵不存在逆矩阵,但在matlab里可以用函数pinv(A)求其伪逆矩阵。理解其实就是求不存在逆矩阵的逆矩阵。

M>N的时候,也就是说我们探测到的数据量大于原数据量,此时方程数量大于未知数数量,可以直接通过求伪逆求出矩阵B,B = pinv(A) * C

%测试伪逆,M > N
A = rand(10,8);
B = rand(8,1);
C = A*B;
BB = pinv(A)*C;
B'
BB'结果:
ans =0.1375    0.3900    0.9274    0.9175    0.7136    0.6183    0.3433    0.9360ans =0.1375    0.3900    0.9274    0.9175    0.7136    0.6183    0.3433    0.9360

N>M的时候,也就是说我们探测到的数据量大于原数据量,此时方程数量小于未知数数量,方程没有唯一解,不能通过B = pinv(A) * C 求出B(压缩感知情况另外考虑 :https://blog.csdn.net/tyfwin/article/details/88902091)

%测试伪逆 M<N
A = rand(8,10);
B = rand(10,1);
C = A*B;
BB = pinv(A)*C;
B'
BB'结果:
ans =0.7360    0.7947    0.5449    0.6862    0.8936    0.0548    0.3037    0.0462    0.1955    0.7202ans =0.7045    0.7584    0.6073    0.5426    1.0107    0.0920    0.2619    0.1060    0.2999    0.6264

逆矩阵、伪逆矩阵:数据的压缩和复原相关推荐

  1. 【转录调控网络】典型的基因转录调控网络推导方法——伪逆矩阵模型

    基因转录调控网络推导方法--伪逆矩阵模型 在基因调控网络推导中,使用到的 基因芯片数据通常具有 样本个数少(通常小于10)而 基因数目大(通常大于1000)的局限性,也就是说,实验样本个数远远小于基因 ...

  2. 最小二乘法与伪逆矩阵

    一.简介 最小二乘法是一种数学优化技术.通过最小化误差的平方和寻找数据的最佳函数匹配.利用最小二乘法可以简便地求得未知的数据,并是得这些求得的数据与实际数据之间误差的平方和最小. 二.最小二乘法拟合直 ...

  3. pinv--求矩阵的伪逆矩阵

    pinv--求矩阵的伪逆矩阵 [功能简介]用于求矩阵的伪逆矩阵. [语法格式] 1.B=pinv(A) 函数返回矩阵A的伪逆矩阵.如果矩阵A是可逆(非奇异)的,那么pinv(A)与inv(A)的结果是 ...

  4. 线性代数笔记:逆矩阵及伪逆矩阵,最小二乘估计,最小范数估计

    逆矩阵及伪逆矩阵,最小二乘估计,最小范数估计 逆矩阵的概念 我们为什么需要逆矩阵? 伪逆矩阵和最小二乘估计 正则化求伪逆矩阵 逆矩阵的概念 矩阵A的逆矩阵(matrix inversion)记作A−1 ...

  5. MATLAB与线性代数--逆矩阵与伪逆矩阵

    求逆矩阵与伪逆矩阵 矩阵A的逆矩阵用 表示,并且满足下面的关系: 看下面的矩阵方程: 如果A的逆矩阵存在,那么解可以写成: 在MATLAB输入下面的命令就可以计算矩阵A的逆矩阵 但是逆矩阵并不一定存在 ...

  6. 关于数学库Eigen两个常用函数实现(伪逆矩阵和协方差矩阵)

    关于数学库Eigen两个常用函数实现(伪逆矩阵和协方差矩阵) cov  pinv 实现代码

  7. 一分钟读懂:矩阵的特征值分解、奇异值分解和伪逆矩阵

    通过把矩阵运算分解成多个矩阵的乘法,可以简化矩阵运算,也可发现对应线性变换的一些内在规律和特性.根据不同的目的,有不同的分解策略.本文我们讨论最常用的特征值分解和奇异值分解. 1. 矩阵的乘方运算 定 ...

  8. 奇异矩阵,非奇异矩阵,伪逆矩阵

    奇异矩阵就是Singular Matrix 的中文翻译. Singular 就是唯一的,可以想成是单身狗,所以他没有对象 逆矩阵. Non-singular的非奇异矩阵就是Couple 有逆矩阵. 奇 ...

  9. 奇异矩阵和非奇异矩阵,伪逆矩阵

    奇异矩阵: 奇异矩阵是线性代数的概念,就是对应的行列式等于0的矩阵. 奇异矩阵的判断方法:首先,看这个矩阵是不是方阵(即行数和列数相等的矩阵.若行数和列数不相等,那就谈不上奇异矩阵和非奇异矩阵). 然 ...

最新文章

  1. Linux终端运行fasterrcnn,对yolo与fasterrcnn anchors的理解
  2. Linux常用命令——fsck
  3. 网络即服务(NaaS)是什么???
  4. Clinicast让癌症治疗不那么昂贵
  5. Gym - 100917F Find the Length-用最小路径树求最小环
  6. LaTeX对公式字体加粗
  7. 网上购物安全防范很重要
  8. java html写入到word文档_java 使用jacob将html页面写入word
  9. AWS表示:我们将用Lambda函数让移动应用程序再次变得伟大
  10. matlab怎么根据公式建立三维模型,紧急!求助!用matlab建立三维copula模型,正用极大似然估计方法求copula函数参数...
  11. NLP 语料库 大全
  12. Fedora14硬盘安装
  13. passenger 部署
  14. 按步搭建简单IoT微服务(2)
  15. 利用公有云cvm 公网IP实现内网穿透,远程登录内网服务器
  16. 科比最新超酷广告《门徒》
  17. fifa18怎么改服务器位置,fifa18 球员职业生涯怎么改位置 | 手游网游页游攻略大全...
  18. 脑的计算模型能带我们走多远
  19. winform程序使用clickonce方式发布之后点击安装没反应
  20. 神经网络——最易懂最清晰的一篇文章

热门文章

  1. FreeBSD5.3下安装Apache+PHP+MySQL+Tomcat
  2. 下载Hibernate中文参考手册2.1版
  3. memcached+keepalived+magent高群集
  4. 用VLC读取摄像头产生RTSP流,DSS主动取流转发(一)
  5. 设置sqlplus环境变量
  6. 对话框中加入标签页的5种方法
  7. Java:Linux上java -jar xxx.jar命令执行jar包时出现Error: Invalid or corrupt jarfile xxx.jar解决方案...
  8. 用vector实现一个变长数组
  9. C#实现(递归和非递归)快速排序和简单排序
  10. Jquery性能优化(转自蓝色理想)