matlab求伪逆,逆、伪逆、左右逆、最小二乘、投影矩阵
主要内容:
矩阵的逆、伪逆、左右逆
矩阵的左逆与最小二乘
左右逆与投影矩阵
一、矩阵的逆、伪逆、左右逆
1、矩阵的逆
定义:
设A是数域上的一个n阶方阵,若在相同数域上存在另一个n阶矩阵B,使得: AB=BA=I。 则我们称B是A的逆矩阵,而A则被称为可逆矩阵。
可逆条件:
A是可逆矩阵的充分必要条件是,即可逆矩阵就是非奇异矩阵。(当
时,A称为奇异矩阵)
性质:
矩阵A可逆的充要条件是A的行列式不等于0。
可逆矩阵一定是方阵。
如果矩阵A是可逆的,A的逆矩阵是唯一的。
可逆矩阵也被称为非奇异矩阵、满秩矩阵。
两个可逆矩阵的乘积依然可逆。
可逆矩阵的转置矩阵也可逆。
矩阵可逆当且仅当它是满秩矩阵。
求逆方法:
伴随矩阵法、初等变换法
2、矩阵的伪逆和左右逆
伪逆矩阵:
伪逆矩阵是逆矩阵的广义形式。由于奇异矩阵或非方阵的矩阵不存在逆矩阵,但在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)花费更少的时间。
伪逆矩阵求法:
A 为m*n矩阵,r代表矩阵的秩:
若矩阵A是方阵,且|A|!=0,则存在AA-1=E;
若A不是方阵,或者|A|=0,那么只能求A的伪逆,所谓伪逆是通过SVD计算出来的;
pinv(A)表示A是伪逆:
如果A列满秩,列向量线性无关,r=n,Ax=b为超定方程组,存在0个或1个解,那么
,因为
,因此也称为左逆;
如果A行满秩,行向量线性无关,Ax=b为欠定方程组,存在0个或无穷个解,那么
,因为
,因此也称为右逆;
如果秩亏损,那么只好先做奇异值分解
,U,V是正交阵,D是对角阵;然后取对角阵S,如果D(i,i)=0,那么S(i,i)=0,如果D(i,i)<>0,那么S(i,i)=1/D(i,i)。于是
;
二、矩阵的左逆与最小二乘
其实,最小二乘就是一个超定方程组的求解问题,根据上述的了解,超定方程组的求解方法之一就是通过求伪逆的形式,具体来说就是求左逆。即:
最小二乘也可以从几何的角度来考虑,那就是下面要说的投影矩阵。
三、左右逆与投影矩阵
左逆中,
,如果将左逆写在A右边将得不到单位矩阵了,那么
是什么?是在A矩阵列空间(A矩阵各列张成的子空间)投影的投影矩阵,它会尽量靠近单位矩阵,一个投影矩阵很想成为单位矩阵,但不可能做到。
右逆中,
,如果将右逆写在A左边也不是单位矩阵了,那
是什么?是在A矩阵行空间(A矩阵各行张成的子空间)投影的投影矩阵。
四、参考文章
matlab求伪逆,逆、伪逆、左右逆、最小二乘、投影矩阵相关推荐
- matlab求矩阵违逆,两种求矩阵伪逆的方法
伪逆矩阵是逆矩阵的广义形式.由于奇异矩阵或非方阵的矩阵不存在逆矩阵,但在matlab里可以用函数pinv(A)求其伪逆矩阵.基本语法为X=pinv(A),X=pinv(A,tol),其中tol为误差, ...
- OpenCV求逆(伪逆)矩阵函数
转自 double invert(InputArray src, OutputArraydst, int flags=DECOMP_LU); 功能:用以求取一个矩阵的逆或者伪逆. src: 输入,浮点 ...
- matlab:inv,pinv逆与伪逆
对于方阵A,如果为非奇异方阵,则存在逆矩阵inv(A) 对于奇异矩阵或者非方阵,并不存在逆矩阵,但可以使用pinv(A)求其伪逆 inv: inv(A)*B 实际上可以写成A\B B*inv(A) 实 ...
- 矩阵的逆、伪逆、左右逆,最小二乘,投影矩阵
主要内容: 矩阵的逆.伪逆.左右逆 矩阵的左逆与最小二乘 左右逆与投影矩阵 一.矩阵的逆.伪逆.左右逆 1.矩阵的逆 定义: 设A是数域上的一个n阶方阵,若在相同数域上存在另一个n阶矩阵B,使得: A ...
- (数学概念)矩阵的逆、伪逆、左右逆,最小二乘,投影矩阵
主要内容: 矩阵的逆.伪逆.左右逆 矩阵的左逆与最小二乘 左右逆与投影矩阵 一.矩阵的逆.伪逆.左右逆 1.矩阵的逆 定义: 设A是数域上的一个n阶方阵,若在相同数域上存在另一个n阶矩阵B,使得: A ...
- 奇异矩阵,非奇异矩阵,矩阵的逆和伪逆
奇异矩阵/非奇异矩阵 首先看矩阵是不是方阵,只有是方阵了,才有这两个概念. 方阵A的行列式等于零,记为|A|=0,A是奇异矩阵 方阵A的行列式不等于零,记为|A|~=0,A是非奇异矩阵 一些性质: 如 ...
- 线性代数学习笔记10-4:左右逆、伪逆/M-P广义逆(从四个子空间和SVD角度理解)
下面讨论m×nm\times nm×n的秩为rrr的矩阵 对于不同情况,讨论逆矩阵 两侧逆矩阵 2-sided inverse 这也是一般所说的"逆矩阵"的含义 方阵A\bolds ...
- 线性代数系列(八)--基变换、左逆、右逆和伪逆
正文 这里关于基变换和伪逆做的都是简单的介绍,关于他们的更深入的理论介绍和更深入的应用介绍还需参考其他资料,然后补充. 基变换 基变换是图像压缩.信号压缩等应用的理论基础,通俗来讲就是对于给定的数据矩 ...
- 如何用Matlab求矩阵的秩、乘积、逆、行列式的值、转置
https://jingyan.baidu.com/article/a65957f495b3ab24e67f9bc2.html 如何用Matlab求矩阵的秩.乘积.逆.行列式的值.转置_Tracy_L ...
- matlab求逆序数
matlab求逆序数 逆序数概念: 因为没时间详细介绍逆序数概念,上传图片仅作参考. 逆序数matlab代码: clc clear %author:猪猪侠 %date:2018-7-18 x=inpu ...
最新文章
- 一些量化(quantization)技巧
- 学习3dmax(四)
- AI换脸技术再创新高度,DeepMind发布的VQ-VAE二代算法有多厉害?
- 目标检测 最新源码大放送
- Javascript中的typeof() 与undefined 与undefined
- 图片在线转换base64
- Linux中环境变量文件及配置
- Apple 如何知道你使用了私有API
- STN(spatial transformer network)论文与源码理解
- 企业入职测试题目及答案_[企业入职性格测试题以及答案信息] 新员工入职性格测试题...
- 华为畅享max什么时候鸿蒙,华为鸿蒙系统正式推送 具体怎么样及支持哪些手机带你前先看...
- hdu 5755 Gambler Bo 高斯消元
- 【新手上路常见问答】关于物联网传输协议MQTT
- Kubernetes v1.10.4 安装记录
- 综合布线 、 子网划分
- 二元一次在线计算机,(二元一次)方程组求解在线计算器_三贝计算网_23bei.com
- 网上不需要实名的物联卡能用吗?有什么安全隐患?
- CentOS安装Eclipse,Eclipse启动时报错
- 短学期微机接口课程设计
- 范德堡大学排名计算机,【usnews世界大学排名2019】2019范德堡大学排名(USNews排名)...
热门文章
- 猜字游戏的各种玩法(附详细源码+演示结果),你都玩过吗?
- openjudge 1.12.2 短信计费
- 送男孩子情人节礼物送什么?情人节必备礼物合集
- NovelAi + Webui + Stable-diffusion本地配置
- php西游记源码,基于Python第三方插件实现西游记章节标注汉语拼音的方法
- Microsoft.Office.Interop.Word 分页操作(文本和表格分开)
- Python转UTC世界标准时间(包含T和Z) 成为北京时间
- RPC - 麻雀虽小,五脏俱全
- 无偿加班竟然被开除,程序员坐不住了
- python制作简单网页_怎么用python简单的制作一个网页