1. 程式人生 > >矩陣的逆、偽逆、左右逆,最小二乘,投影矩陣

矩陣的逆、偽逆、左右逆,最小二乘,投影矩陣

主要內容:

  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)。於是

二、矩陣的左逆與最小二乘

其實,最小二乘就是一個超定方程組的求解問題,根據上述的瞭解,超定方程組的求解方法之一就是通過求偽逆的形式,具體來說就是求左逆。即:

最小二乘也可以從幾何的角度來考慮,那就是下面要說的投影矩陣。

三、左右逆與投影矩陣

左逆中, ,如果將左逆寫在A右邊將得不到單位矩陣了,那麼 是什麼?是在A矩陣列空間(A矩陣各列張成的子空間)投影的投影矩陣,它會盡量靠近單位矩陣,一個投影矩陣很想成為單位矩陣,但不可能做到。

右逆中, ,如果將右逆寫在A左邊也不是單位矩陣了,那是什麼?是在A矩陣行空間(A矩陣各行張成的子空間)投影的投影矩陣。

四、參考文章