1. 程式人生 > >矩陣求逆常見演算法

矩陣求逆常見演算法

前言

    不知道從哪天開始,看到矩陣就頭疼,特別是矩陣的運算更是蛋疼,都不好意思說自己是數學專業的,哈哈。這兩天在搞opencv影象處理,又涉及到這一塊,無語之,乾脆收集整理下,以饗同痛苦者。

一、逆矩陣的概念

利用矩陣的乘法和矩陣相等的含義,可以把線性方程組寫成矩陣形式。對於線性方程組

A   X   B

則方程組可寫成AX=B.

方程AX=B是線性方程組的矩陣表達形式,稱為矩陣方程。其中A稱為方程組的係數矩陣X稱為未知矩陣B稱為常數項矩陣

這樣,解線性方程組的問題就變成求矩陣方程中未知矩陣X的問題。類似於一元一次方程ax=ba≠0)的解可以寫成x=a-1b,矩陣方程AX

=B的解是否也可以表示為X=A-1B的形式?如果可以,則X可求出,但A-1的含義和存在的條件是什麼呢?下面來討論這些問題。

定義11  對於n階方陣A,如果存在n階方陣C,使得AC=CA=EEn階單位矩陣),則把方陣C稱為A逆矩陣(簡稱逆陣)記作A-1,即C=A-1

例如     

因為  AC==

CA==

所以C是的A逆矩陣,即C=A-1

由定義可知,AC=CA=ECA的逆矩陣,也可以稱AC的逆矩陣,即A=C-1。因此,AC稱為互逆矩陣。

可以證明,逆矩陣有如下性質:

(1)若A是可逆的,則逆矩陣唯一。

(2)若A可逆,則(A-1)-1=A.

(3)若AB為同階方陣且均可逆,則AB

可逆,且(AB)-1=B-1A-1

(4)若A可逆,則detA≠0。反之,若detA≠0,則A是可逆的。

  (1)如果BC都是A的逆矩陣,則

C=CE=C(AB)=(CA)=EB=B

即逆矩陣唯一。

其它證明略。

二、逆矩陣的求法

1、用伴隨矩陣求逆矩陣

定義12  設矩陣

A=

所對應的行列式detA中元素aij的代數餘子式矩陣

稱為A伴隨矩陣,記為A*

顯然,AA=

仍是一個n階方陣,其中第i行第j列的元素為

由行列式按一行(列)展開式可知

=

所以       AA==detAE    (1)

同理  AA=detAE=AA

定理3  n階方陣A可逆的充分必要條件是A為非奇異矩陣,而且

A-1=A=

  必要性:

如果A可逆,則A-1存在使AA-1=E,兩邊取行列式det(AA-1)= detE,即detAdetA-1=1,因而detA≠0,即A為非奇異矩陣。

充分性:

A為非奇異矩陣,所以detA≠0,由(1)式可知A(A)= (A)A=E

所以A是可逆矩陣。

A-1=A

例1  求矩陣A=的逆矩陣。

  因為detA=,所以A是可逆的。又因為

    

    

    

所以A-1=A=

=

2、用初等變換求逆矩陣

用初等變換求一個可逆矩陣A的逆矩陣,其具體方法為:把方陣A和同階的單位矩陣E,寫成一個長方矩陣,對該矩陣的行實施初等變換,當虛線左邊的A變成單位矩陣E時,虛線右邊的E變成了A-1

從而可求A-1

例2  用初等變換求

的逆矩陣。

  因為 =

所以 A-1=

例3 解線性方程組

  方程組可寫成

=

A=  X=  B=  則AX=B

由例2知A可逆,且A-1=

所以X=A-1B,即=A-1B==

於是,方程組的解是

習題 12--6

1、用伴隨矩陣求下列矩陣的逆矩陣:

(1)   (2)   (3)

2、用初等變換求逆矩陣:

(1)   (2)   (3)

(4)

3、解矩陣方程

(1)X

(2)

(3)

(4)

4、解線性方程組

(1)

(2)