趙姣珍，許道云

(1.貴陽人文科技學(xué)院 大數(shù)據(jù)與信息工程學(xué)院，貴州 貴陽 550025；2.貴州大學(xué) 計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院，貴州 貴陽 550025)

矩陣的行列式值是矩陣的一個(gè)不變量，它是代數(shù)中一個(gè)重要的基礎(chǔ)概念，對(duì)矩陣各種性質(zhì)研究以及矩陣在其他領(lǐng)域中的應(yīng)用，幾乎都與行列式相關(guān)。本文考慮的矩陣均指n階方陣，給定一個(gè)n階矩陣，可用列向量表示為

通常，矩陣的行列式定義以如下公式給出[1]：

(1)

其中，π是集合{1,2,…,n}上的一個(gè)置換。置換π可以被分解為對(duì)換的乘積，π的奇偶性由它被分解為最少對(duì)換個(gè)數(shù)的奇偶性決定，此與序列(π(1),π(2),…,π(n))的逆序數(shù)(記為τ(π))的奇偶性一致。式(1)中的求和則是取遍{1,2,…,n}上的所有置換。

從式(1)可以找到二階和三階行列式的計(jì)算規(guī)律，并用于計(jì)算。但對(duì)于四階以上的n階行列式，利用式(1)作為計(jì)算公式是不現(xiàn)實(shí)的。實(shí)際計(jì)算時(shí)，主要是基于行列式的性質(zhì)不斷地降階，或化為特殊矩陣，或按性質(zhì)進(jìn)行計(jì)算。

一般，有關(guān)行列式的性質(zhì)是由式(1)推導(dǎo)，部分性質(zhì)推導(dǎo)較為復(fù)雜?；诠砘x的行列式，從公理形式出發(fā)進(jìn)行推導(dǎo)，在邏輯思路、推導(dǎo)過程、簡(jiǎn)潔性等方面都有其優(yōu)點(diǎn)。

文獻(xiàn)[2]中以矩陣的列向量(A1,A2,…,An)作為變量，引入行列式函數(shù)，以公理形式給出了行列式定義。

本文基于公理本身，討論了公理的等價(jià)性和獨(dú)立性，從公理化定義行列式出發(fā)，直接導(dǎo)出行列式常見的基本性質(zhì)和普通定義計(jì)算式(1)，并給出相關(guān)重要性質(zhì)的推導(dǎo)。

1 行列式的公理化定義及等價(jià)公理

定義[2]設(shè)矩陣A=(A1,A2,…,An)，考慮一個(gè)實(shí)數(shù)函數(shù)det(A1,A2,…,An)，如果函數(shù)滿足如下公理?xiàng)l件，則稱det(A1,A2,…,An)為矩陣A的行列式。

公理1對(duì)于任意固定的(1≤k≤n),以Ak為變量，其余列不變的情況下誘導(dǎo)出的函數(shù)detk(Ak)具有齊次線性性質(zhì)。即，對(duì)于任意常數(shù)a，b，detk(aAk+bBk)=adetk(Ak)+bdetk(Bk)。

公理2若存在相鄰兩列相等，其值為0。即，如果存在某個(gè)1≤k

公理3對(duì)于單位矩陣U=(U1,U2,…,Un)，det(U1,U2,…,Un)=1。其中Uk為第k個(gè)單位向量。

首先，可由公理?xiàng)l件直接推出如下基本性質(zhì)：

性質(zhì)1如果矩陣中有一列全為0，則行列式為0。

事實(shí)上，如果矩陣A的第k列全為0，由公理1中的齊次性，有detk(0)=detk(0·0)=0·detk(0)=0。

性質(zhì)2將矩陣中一列的c倍加到相鄰一列后，則行列式不變。

假設(shè)由矩陣中的第k+1列的c倍加到第k列，則新矩陣的第k列為Ak+cAk+1。由公理2，detk(Ak+1)=0。再由公理1,有detk(Ak+cAk+1)=detk(Ak)+c·detk(Ak+1)=detk(Ak)。

性質(zhì)3將矩陣中相鄰兩列互換后，行列式改變符號(hào)。

證明設(shè)A=(A1,…,Ak,Ak+1,…,An)，兩列互換后得到矩陣A′=(A1,…,Ak+1,Ak,…,An)。

將矩陣A中第k列加到第k+1列，得到B=(A1,…,Ak,Ak+Ak+1,…,An)；將矩陣B中第k+1列的(-1)倍加到第k列，得到C=(A1,…,(-1)Ak+1,Ak+Ak+1,…,An)；將矩陣C中第k列加到第k+1列，得到D=(A1,…,(-1)Ak+1,Ak,…,An)。

由性質(zhì)2以及公理1中的齊次性質(zhì)，有

det(A1,…,Ak,Ak+1,…,An)

=det(A1,…,Ak,Ak+Ak+1,…,An)

=det(A1,…,(-1)Ak+1,Ak+Ak+1,…,An)

=det(A1,…,(-1)Ak+1,Ak,…,An)

=(-1)·det(A1,A2,…,Ak+1,Ak,…,An)

性質(zhì)1～3完全由公理本身得到。很顯然：①性質(zhì)1和性質(zhì)2可以導(dǎo)出公理2；②性質(zhì)1和性質(zhì)3也可以導(dǎo)出公理2。換言之，分別以性質(zhì)1和性質(zhì)2取代公理2，以性質(zhì)1和性質(zhì)3取代公理2，可以得到行列式函數(shù)的另外兩個(gè)等價(jià)公理定義。

請(qǐng)注意：公理定義中公理1和公理3是本質(zhì)的。

2 行列式的其他性質(zhì)

下面的性質(zhì)表明：公理2、性質(zhì)2和性質(zhì)3中的“相鄰”條件可以去掉。

性質(zhì)4將矩陣中不同兩列互換后，行列式改變符號(hào)。

證明指定兩個(gè)不同列號(hào)i,j(i

因此，一共作了奇數(shù)次相鄰列互換。由性質(zhì)3，det(A)=(-1)det(C)。即：矩陣中不同兩列互換后，則行列式改變符號(hào)。

類似證明：

性質(zhì)5如果矩陣中有兩列相等，則行列式為0。

性質(zhì)6將矩陣中任一列的c倍加到另一列后，行列式不變。

證明任意指定兩個(gè)不同列號(hào)i,j(i

由性質(zhì)4，det(A)=(-1)det(B); 由性質(zhì)2，det(B)=det(C)。再由性質(zhì)4，det(C)=(-1)·det(D)。所以，det(A)=det(D)。即：矩陣中任一列的c倍加到另一列后，行列式不變。

(4,6,5,7)(9)，其輪換的長(zhǎng)度分別為：3，1，4，1。通常，略去單點(diǎn)輪換后簡(jiǎn)單地表示為π=(1,2,8)(4,6,5,7)，未出現(xiàn)的元素表示自己映射到自己。長(zhǎng)度為2的輪換稱為對(duì)換。

可以驗(yàn)證：一個(gè)長(zhǎng)度為k的輪換可以表示成k-1個(gè)對(duì)換的復(fù)合，且為最小對(duì)換個(gè)數(shù)：

(2)

假定將一個(gè)置換π分解為輪換積形式時(shí)，其輪換的長(zhǎng)度序列為l1,l2,…,lm，則一個(gè)置換π可以表示為(l1-1)+(l2-1)+…+(lm-1)個(gè)對(duì)換的復(fù)合。記

τ#(π)=(l1-1)+(l2-1)+…+(lm-1)

它記錄了置換π表示為對(duì)換的個(gè)數(shù)。

可以證明：

1)τ#(π)是置換π表示為對(duì)換時(shí)的最小對(duì)換個(gè)數(shù)。因此，這個(gè)數(shù)是唯一的，以其奇偶性定義置換π的奇偶性。

2)τ#(π)的奇偶性與整數(shù)序列(π(1),π(2),…,π(n))的逆序數(shù)τ(π)的奇偶性一致。

事實(shí)上，單位置換對(duì)應(yīng)的逆序數(shù)為0。相鄰元素作對(duì)換后，逆序數(shù)的改變相差1。一般對(duì)換使用后，逆序數(shù)的改變相差一個(gè)奇數(shù)。

有了上面的討論，由性質(zhì)4，我們有：

性質(zhì)7給定A=(A1,A2,…,An)，以及列號(hào)集{1,2,…,n}上的一個(gè)置換π，有如下關(guān)系：

det(Aπ(1),Aπ(2),…,Aπ(n))=(-1)τ(π)det(A)

(3)

特別，由公理3，對(duì)于單位矩陣U=(U1,U2,…,Un)，我們有

det(Uπ(1),Uπ(2),…,Uπ(n))=(-1)τ(π)

(4)

其中，τ(π)為整數(shù)序列(π(1),π(2),…,π(n))的逆序數(shù)。

矩陣行列式有一條重要性質(zhì)：矩陣A的行列式與其轉(zhuǎn)置AT的行列式相等。我們現(xiàn)在來看一下這條性質(zhì)如何從公理出發(fā)得到。

首先，我們注意到：對(duì)于單位矩陣U,U=UT。從而，det(UT)=det(U)=1。

(5)

將公理1應(yīng)用到式(5)中，從第1個(gè)列向量開始，依次展開：

(6)

由性質(zhì)5，矩陣中有兩列相等時(shí)，行列式為0。因此，式(6)完全展開后得到

其中，求和是取遍{1,2,…,n}上所有置換。

為了看清楚這一點(diǎn)，讀者可以n=3推導(dǎo)式(7)。

由式(4),我們有

(8)

同樣，對(duì)于

A=(A1,A2,…,An)

行列式展開后得到

我們知道：{1,2,…,n}的所有置換π，在復(fù)合運(yùn)算下構(gòu)成對(duì)稱群Sn，對(duì)于任一個(gè)置換π，逆元π-1與π一一對(duì)應(yīng)，并且τ#(π-1)=τ#(π)。從而，(-1)τ(π-1)=(-1)τ#(π-1)=(-1)τ#(π)=(-1)τ(π)。

改寫式(9)中系數(shù)項(xiàng)中的乘積順序。行、列下標(biāo)自然順序調(diào)整有如下關(guān)系：

aπ(1),1aπ(2),2…aπ(n),ndet(Uπ(1),Uπ(2),…,Uπ(n))

=a1,π-1(1)a2,π-1(2)…an,π-1(n)det(Uπ-1(1),Uπ-1(2),…,

Uπ-1(n))

(10)

因?yàn)?/p>

aπ(1),1aπ(2),2…aπ(n),ndet(Uπ(1),Uπ(2),…,Uπ(n))

=(-1)τ(π)a1,π-1(1)a2,π-1(2)…an,π-1(n)det(U1,U2,…,Un)

=(-1)τ(π)+τ(π-1)a1,π-1(1)a2,π-1(2)…an,π-1(n)·

det(Uπ-1(1),Uπ-1(2),…,Uπ-1(n))

=a1,π-1(1)a2,π-1(2)…an,π-1(n)det(Uπ-1(1),Uπ-1(2),…,Uπ-1(n))

因此，我們有

det(Uπ(1),Uπ(2),…,Uπ(n))

det(Uπ-1(1),Uπ-1(2),…,Uπ-1(n))

(11)

當(dāng)求和取遍所有置換時(shí)，有如下關(guān)系：

Uπ-1(2),…,Uπ-1(n))

Uπ(2),…,Uπ(n))

Uπ(2),…,Uπ(n))

=det(AT)

(12)

由此，有如下性質(zhì)：

性質(zhì)8矩陣A的行列式與其轉(zhuǎn)置AT的行列式相等。即，det(AT)=det(A)。

有了性質(zhì)8，行列式公理定義中，由“列向量”改為“行向量”作變量定義行列式函數(shù)同樣有上述平行性質(zhì)。因此，有關(guān)“列”的性質(zhì)，對(duì)“行”同樣成立。

同時(shí)，由式(11)及公理3，我們可以得到通常行列式的定義公式：

det(Uπ(1),Uπ(2),…,Uπ(n))

det(U1,U2,…,Un)

(13)

其中，求和是取遍{1,2,…,n}的所有置換π，一共有n!項(xiàng)。

由于是對(duì)全體置換求和，我們有

(14)

在上述推導(dǎo)過程中，我們得到如下2個(gè)有用公式：

(15)

(16)

注意：式(15)中，等式左端矩陣部分是取A的轉(zhuǎn)置AT，右端求和系數(shù)項(xiàng)中的乘積是以行標(biāo)為自然順序；而在式(16)中，等式左端矩陣部分是取A，右端求和系數(shù)項(xiàng)中的乘積是以列標(biāo)為自然順序。

用公理方法，可以自然地推出著名的Laplace定理：矩陣乘積(矩陣乘)的行列式等于矩陣行列式的乘積(實(shí)數(shù)乘)。即，det(AB)=det(A)det(B)。

設(shè)有兩個(gè)同階方陣A和B，其矩陣乘法AB形式可以表示為

(17)

仿式(15)(16)，我們有

det(Aπ(1),Aπ(2),…,Aπ(n))

det(Aπ-1(1),Aπ-1(2),…,Aπ-1(n))

det(A1,A2,…,An)

=det(A1,A2,…,An)·

=det(A1,A2,…,An)·

=det(A1,A2,…,An)·det(B1,B2,…,Bn)

=det(A)·det(B)

其中

是依據(jù)求和取遍{1,2,…,n}的所有置換。

3 高階行列式降階與行列式的分解

我們知道：高階行列式的計(jì)算主要是通過降階。

由公理得到的一些主要性質(zhì)(如：行(列)互換行列式變號(hào)，一行(列)的c倍加到另行(列)行列式不變)，以及如下的降階原理可以計(jì)算n階行列式。

對(duì)于n階矩陣A=(A1,…,An)，可以表示成如下形式：

(18)

如果A中第1行全為0，則det(A)=0，否則至少有一個(gè)不為0。于是，通過適當(dāng)?shù)牧谢Q，以及第1(行)列的某個(gè)倍數(shù)加到另一(行)列，矩陣可化為如下形式：

(19)

其中，矩陣B=(B2,…,Bn)為n-1階方陣。對(duì)于形如式(19)的矩陣，利用已經(jīng)由公理推導(dǎo)出的行列式

可以推導(dǎo)出：det(A)=adet(B)。

我們現(xiàn)在觀察n階矩陣A的行列式分解為它的n-1階子矩陣的行列式之間的關(guān)系：

對(duì)于式(18)表示的矩陣A，記A[i,j]為在矩陣A中刪去第i行、第j列后得到的n-1階子矩陣。

展開如下公式：

=a1,1det(A[1,1])+

=a1,1det(A[1,1])+…

我們有

(20)

一般地，對(duì)于固定的i，有如下分解計(jì)算公式：

(21)

對(duì)于式(20),對(duì)于2≤i≤n,由兩行相同行列式為0，我們有

(22)

由此導(dǎo)出一條重要性質(zhì)：如果行列式det(A)≠0，對(duì)于第1個(gè)n維單位向量U1，線性方程組

x1A1+…+xnAn=U1

有解，其中解x(1)的分量

(23)

一般地，類似方法可得到：如果行列式det(A)≠0，對(duì)于第k個(gè)n維單位向量Uk，線性方程組

x1A1+…+xnAn=Uk

有解，其中解x(k)的分量計(jì)算公式為

(24)

請(qǐng)注意：在det(A)≠0條件下，文中線性方程組解的存在性完全由公理及行列式性質(zhì)獨(dú)立推出,并非由det(A)≠0條件按如下路徑得到：A1,…,An線性無關(guān)，由A1,…,An構(gòu)成Rn空間的生成系，再由生成系生成B(如文獻(xiàn)[2])。由生成系生成B要用到線性方程組的解，從邏輯上講，這里出現(xiàn)一個(gè)循環(huán)推導(dǎo)問題。

將式(23)(24)組合在一個(gè)公式中，得到

由此，矩陣A的逆矩陣

ci,j=(-1)i+jdet(A[i,j])

(25)

這就是可逆矩陣的逆矩陣計(jì)算公式(方法)——代數(shù)余子式方法。

4 線性方程組解的存在性及解的表示

行列式之所以重要，最主要的原因之一是它提供了求解線性方程組的一般方法：克萊姆(Cramer)法則。

公理化定義的行列式很容易導(dǎo)出克萊姆法則，這體現(xiàn)了公理化方法的優(yōu)點(diǎn)。

給定一組n維列向量A1,…,Am，稱A1,…,Am線性相關(guān)，指：存在其中一個(gè)列向量Ak，Ak可以表示由其余列向量的線性組合表示。如果不是線性相關(guān)，則稱為線性無關(guān)(或線性獨(dú)立)。

線性相關(guān)等價(jià)于：存在一組不全為0的數(shù)α1,…,αm，使得α1A1+…+αmAm=0。線性無關(guān)等價(jià)于：對(duì)任意一組數(shù)α1,…,αm，如果α1A1+…+αmAm=0，則α1=…=αm=0。

在討論線性方程組解的存在性與系數(shù)矩陣的行列式之間的關(guān)系之前，從邏輯上講，我們應(yīng)該先考慮行列式與列向量線性相關(guān)性質(zhì)之間的關(guān)系，而不是從線性方程組解的性質(zhì)討論這一關(guān)系。

關(guān)于向量組的線性相關(guān)性質(zhì)，容易驗(yàn)證：①在一個(gè)向量組中，如果存在部分向量構(gòu)成的子向量組線性相關(guān)，則該向量線性相關(guān)；②對(duì)一組n維列向量A1,…,Am，B1,…,Bm是由A1,…,Am在同一位置(如第1位置)插入0后得到n+1維向量組，則A1,…,Am與B1,…,Bm的線性相關(guān)性一致。

由第3節(jié)得到的性質(zhì)，我們可以得到如下性質(zhì)：

性質(zhì)9對(duì)于矩陣A=(A1,…,An)，如果det(A1,…,An)=0，并且A第1列都不全為0，則通過如下操作：

1)調(diào)整行的順序使(1,1)位置元素不為0；

2)將第1列的某個(gè)倍數(shù)加到第2,…,n列某一列，得到矩陣

其中det(B2,…,Bn)=0。

請(qǐng)注意：第一類操作只可能改變行列式符號(hào)，第二類操作不改變行列式值，因此

性質(zhì)9就是通常將一個(gè)矩陣劃為三角矩陣的方法：行列式值只改變符號(hào)。對(duì)于行列式為0的原始矩陣，每一步操作后所得到矩陣的行列式仍為0，并且A1,…,An與B1,…,Bn的線性相關(guān)性一致。

請(qǐng)注意：對(duì)于一階矩陣，當(dāng)其行列式為0時(shí)，該矩陣只能為零矩陣。

由上述討論，矩陣A=(A1,…,An)中，A1,…,An的線性相關(guān)性可以由A的行列式值是否為0判定。

性質(zhì)10對(duì)于矩陣A=(A1,…,An)，向量組A1,…,An的線性相關(guān)性當(dāng)且僅當(dāng)det(A)=0。

證明如果A1,…,An線性相關(guān)，則存在一個(gè)列向量Ak，Ak=α1A1+…+αk-1Ak-1+αk+1Ak+1…+αnAn,則由性質(zhì)6及性質(zhì)1，det(A1,…,Ak,…,An)=det(A1,…,0,…,An)=0。

反之，如果det(A)=0，如果A中有一列向量全為0，則A1,…,An中含有零向量，從而線性相關(guān)。否則，取定A1，其中至少有一個(gè)非零分量。按性質(zhì)9，從A可以得到一個(gè)形式如下的矩陣：

其中det(B)=det(B2,…,Bn)=0，且A1,…,An與B1,…,Bn的線性相關(guān)性一致。

對(duì)B重復(fù)上述操作(至多n-1次)，可以在某一步上得到C=(Ck,…,Cn)(1

利用性質(zhì)10，我們得到如下行列式與線性方程解之間的關(guān)系。

性質(zhì)11對(duì)于線性方程x1A1+…+xnAn=B，我們有：

1)如果對(duì)于任意一個(gè)n維列向量B，方程有解，則det(A1,…,An)≠0。

2)如果det(A1,…,An)≠0，則對(duì)于任意給定的一個(gè)n維列向量B，方程有解。從而，如果存在某個(gè)n維列向量B，方程無解，則det(A1,…,An)=0。

證明1) 假定對(duì)于任意的n維列向量B，方程有解，分別取B為單位向量U1,…,Un，記方程x1A1+…+xnAn=Ui的解為Bi，則(A1,…,An)(B1,…,Bn)=(U1,…,Un)。因此

det((A1,…,An)(B1,…,Bn))

=det(U1,…,Un)=1

由Laplace定理：

det((A1,…,An)(B1,…,Bn))

=det(A1,…,An)det(B1,…,Bn)

從而，det(A1,…,An)≠0。

而解x*的分量由如下公式計(jì)算：

可以驗(yàn)證

=det(A1,…,Ak-1,B,Ak+1,…,An)

從而

這就是Cramer法則對(duì)解向量的分量的計(jì)算公式。

反過來，如果det(A1,…,An)≠0，且方程x1A1+…+xnAn=B有解(請(qǐng)注意，這里先假設(shè)有解)，則解向量中分量xj可以表示為

這就是著名的Cramer法則。

Cramer法則的推導(dǎo)：由x1A1+…+xnAn=B，以第1列為例：將B替換A1，計(jì)算行列式det(B,A2,…,An)。由性質(zhì)6，從第2列至第n列，分別以第k列的-xk倍加到第1列(行列式不變)，最后使用公理1得

det(B,A2,…,An)

=det(x1A1+x2A2+…+xnAn,A2,…,An)

=det(x1A1+x3A3+…+xnAn,A2,…,An)

…

=det(x1A1,A2,…,An)

=x1det(A1,A2,…,An)

5 結(jié)語

在行列式的公理定義中，三條公理是相互獨(dú)立的。原因是：公理1只考慮矩陣中任一列上的齊次線性性質(zhì)；公理2是考慮矩陣中兩列比較，不涉及運(yùn)算；公理3是界定行列式函數(shù)的“初始”邊界值，以單位矩陣U的行列式為1作為“種子值”。因此，由其中任意兩個(gè)公理不能導(dǎo)出第三條公理。

其次，根據(jù)第1節(jié)中的討論，行列式的公理定義中的公理2有兩種替換方式。

對(duì)于公理3，可修改det(U1,U2,…,Un)的初始值。如：定義一個(gè)函數(shù)F(A)=F(A1,…,An)滿足公理1和公理2，修改公理3為F(U)=F(U1,…,Un)=c≠0,則可以證明：F(A1,…,An)=c·det(A1,…,An)。

更進(jìn)一步，我們可以將F(A)=F(A1,…,An)的取值不是定義在實(shí)數(shù)集上，而是定義于一個(gè)域上(或者一個(gè)抽象空間上)，通過適當(dāng)規(guī)定公理，可以將矩陣的性質(zhì)映射到相應(yīng)空間上。同樣，可以修改公理，研究滿足公理的函數(shù)的相關(guān)性質(zhì)。

有關(guān)行列式的定義方式，還可以用Valiant給出的圖論方法引入，相關(guān)定義請(qǐng)參見文獻(xiàn)[3-4]。對(duì)于高階稀疏矩陣行列式的計(jì)算，有時(shí)以圖的方式引入的定義更有效。對(duì)于MacMahon主定理、Cayley-Hamilton定理、矩陣樹定理、特征多項(xiàng)式定理等的證明，圖論方法引入行列式體現(xiàn)了新的思路。詳細(xì)介紹請(qǐng)參閱文獻(xiàn)[5-10]。

a1,π(1)a2,π(2)…an,π(n)，其求和是對(duì)集合J={1，2，…,n}上所有置換求和，即在對(duì)稱群Sn中取每個(gè)元素作用于列標(biāo)號(hào)集J。我們知道：對(duì)于定義在J上的任意一個(gè)(有限)群G,都對(duì)應(yīng)產(chǎn)生Sn的一個(gè)子群?；蛘哒f，G同構(gòu)于Sn的一個(gè)子群。