馬格尼託哥爾斯克

[拼音]:guangyi nijuzhen

[英文]:generalized inverse matrix

逆矩陣概念的推廣。若

A

為非奇異矩陣,則線性方程組

A

尣=b的解為尣=

A

_1b,其中

A

的逆矩陣

A

_1滿足

A

A

_1=

A

_1

A

=I(I為單位矩陣)。若

A

是奇異陣或長方陣,

A

尣=b可能無解或有很多解。若有解,則解為尣=Xb+(I-X

A

)у,其中у是維數與

A

的列數相同的任意向量,X是滿足AXA=

A

的任何一個矩陣,通常稱X為

A

的廣義逆矩陣,用

A

g、

A

_或

A

等符號表示,有時簡稱廣義逆。當

A

非異時,

A

_1也滿足

A

A

_1

A

=

A

,且

。故非異陣的廣義逆矩陣就是它的逆矩陣,說明廣義逆矩陣確是通常逆矩陣概念的推廣。

1955年R.彭羅斯證明了對每個m×n階矩陣

A

,都存在惟一的n×m階矩陣 X,它滿足:

(1)AXA=

A

(2)XAX=X;

(3)(

A

X)*=

A

X;

(4)(X

A

)*=X

A

。通常稱X為

A

的穆爾-彭羅斯廣義逆矩陣,簡稱M-P逆,記作

A

+。當

A

非異時,

A

_1也滿足①~④,因此M-P逆也是通常逆矩陣的推廣。在矛盾線性方程組

A

尣=b的最小二乘解中,尣=

A

+b是範數最小的一個解。

A

是n階方陣,k為滿足

的最小正整數(rank為矩陣秩的符號),記作k=Ind(

A

),則存在惟一的n階方陣X,滿足:

(1) AkX

A

=

A

k;(2) X

A

X=X; (3) AX=XA。

通常稱X為

A

的德雷津廣義逆矩陣,簡稱D逆,記作Ad,

A

(d)或

A

D等。雖然它和線性代數方程組的解無關,但它線上性差分方程、線性微分方程、最優控制等方面都有應用。例如,設

A

B

是n階方陣,齊次差分方程

,如果存在一個數λ,使

存在,則它的一般解為

式中q為任意n維向量;

;

根據實際問題需要還定義了其他各種型別的廣義逆矩陣,如網路理論中用到的博特-達芬逆矩陣等。一般說來,它們都具有下列一些性質:當

A

非異時,廣義逆矩陣就是

A

_1;廣義逆矩陣必存在;廣義逆矩陣具有逆矩陣的某些性質(或適當修改後的性質),如(

A

_1)_1=

A

,(

A

_1)*=(

A

*)_1等等。

廣義逆的思想可追溯到1903年(e.)i.弗雷德霍姆的工作,他討論了關於積分運算元的一種廣義逆(他稱之為偽逆)。1904年,d.希爾伯特在廣義格林函式的討論中,含蓄地提出了微分運算元的廣義逆。而任意矩陣的廣義逆定義最早是由E.H.穆爾在1920年提出的,他以抽象的形式發表在美國數學會會刊上。當時人們對此似乎很少注意。這一概念在以後30年中沒有多大發展。曾遠榮在1933年,F.J.默裡和j.馮·諾伊曼在1936年對希爾伯特空間中線性運算元的廣義逆作過討論。20世紀50年代圍繞著某些廣義逆的最小二乘性質的討論重新引起了人們對這個課題的興趣。1951年瑞典人A.布耶爾哈梅爾重新發現了穆爾所定義的廣義逆,並注意到廣義逆與線性方程組的關係。T.N.E.格雷維爾、C.R.拉奧和其他人也作出了重要的貢獻。1955年,彭羅斯證明了存在惟一的X=

A

+滿足前述性質①~④,並以此作為 A+的定義。1956年,R.拉多證明了彭羅斯定義的廣義逆與穆爾定義的廣義逆是等價的,因此通稱

A

+為穆爾-彭羅斯廣義逆矩陣。

廣義逆的計算方法大致可分為三類:以滿秩分解和奇異值分解為基礎的直接法,迭代法和其他一些常用於低階矩陣的特殊方法。

A

+的計算為例。若

A

是一個秩為r的m×n階非零矩陣,記作

,有滿秩分解

A

=F·G,其中

,則

,即將廣義逆矩陣的計算化為通常逆矩陣的計算。常用LU分解和

Q

R分解等方法實現滿秩分解,然後求出

A

+。

A

有奇異值分解

A

=UDV*,其中U、V為m階和n階酉矩陣,

是m×n階矩陣,

是r階對角陣,對角元

A

的r個非零奇異值(

A

A

*的非零特徵值的平方根),則

A

+=VD+U*,其中

是n×m階矩陣。也可用豪斯霍爾德變換先將 A化為上雙對角陣J0=

P

*

A

Q

,然後再對J0使用

Q

R演算法化為矩陣D=G*J0h,於是

A

=(

P

G)D(

Q

h)*,故

A

+1=(

Q

h)D+(

P

G)*。

設λ1是

A

A

*的最大非零特徵值,若0<α<2/λ1,則計算

A

+的一個迭代法是x0=α

A

*,xn+1=(2I-

A

xn),當n→∞時,xn收斂於

A

+。

格雷維爾逐次遞推法也是計算

A

+的常用方法。設

A

的第k列為αk(k=1,2,…,n),

A

1=α1,

A

k=(

A

k-1,αk)(k=2,3,…,n),則

式中

1955年以後,出現了大量的關於廣義逆矩陣的理論、應用和計算方法的文獻。70年代還出版了一些專著和會議錄,指出廣義逆矩陣在控制論、系統辨識、規劃論、網路理論、測量、統計和計量經濟學等方面的應用。

參考書目

S.L.Campbell and C.D.Meyer,Jr.,Generalized Inverses of Linear TransforMations,Pitman,London, 1979.

M.Z.Nashed, ed.,Generalized Inverses and Applications,Academic Press,New York,1976.