码迷,mamicode.com
首页 > 其他好文 > 详细

(数学概念)矩阵的逆、伪逆、左右逆,最小二乘,投影矩阵

时间:2015-12-28 15:40:52      阅读:286      评论:0      收藏:0      [点我收藏+]

标签:

主要内容:

  1. 矩阵的逆、伪逆、左右逆
  2. 矩阵的左逆与最小二乘
  3. 左右逆与投影矩阵

一、矩阵的逆、伪逆、左右逆

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)。于是技术分享

二、矩阵的左逆与最小二乘

关于最小二乘可以参考:最小二乘的几何意义及投影矩阵http://www.cnblogs.com/AndyJee/p/5053354.html

其实,最小二乘就是一个超定方程组的求解问题,根据上述的了解,超定方程组的求解方法之一就是通过求伪逆的形式,具体来说就是求左逆。即:

技术分享

最小二乘也可以从几何的角度来考虑,那就是下面要说的投影矩阵。

三、左右逆与投影矩阵

左逆中, 技术分享,如果将左逆写在A右边将得不到单位矩阵了,那么 技术分享是什么?是在A矩阵列空间(A矩阵各列张成的子空间)投影的投影矩阵,它会尽量靠近单位矩阵,一个投影矩阵很想成为单位矩阵,但不可能做到。

右逆中, 技术分享,如果将右逆写在A左边也不是单位矩阵了,那技术分享是什么?是在A矩阵行空间(A矩阵各行张成的子空间)投影的投影矩阵。

四、参考文章

http://baike.baidu.com/link?url=whnNGl6wlBJ7bIzn-ldxZ3KfXj03WlXxuJvLw2VPLcCjLvFtSU_7csPUyNQ57cMzk9zz-y6sG_7hrt88NHcg2a

http://baike.baidu.com/link?url=9BBn2Hc2IgUjr2bwr8CGOFvNRfSWZB3AW6_p5DjTxY74OtZJJYvXIMQPmQ3zDpDsX36HLkEbeskvVczEruqHFa

http://shijuanfeng.blogbus.com/logs/206966888.html

http://www.blogbus.com/shijuanfeng-logs/238839798.html

http://blog.sina.com.cn/s/blog_438e26440102vsm8.html

(数学概念)矩阵的逆、伪逆、左右逆,最小二乘,投影矩阵

标签:

原文地址:http://www.cnblogs.com/AndyJee/p/5082508.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!