-
引理:r[ACOB]=r[AOOB]的充要條件爲存在矩陣X,Y使得C=XA+BY
證:
充分性:r[ACOB]=r[AXA+BYOB]=行變換r[ABYOB]=列變換r[AOOB]
必要性:由等價標準形定理,存在可逆矩陣P1,Q1,P2,Q2使得P1AQ1=[ItOOO]且P2BQ2=[IsOOO],其中t=r(A),s=r(B)。計算可得[P1OOP2][AOOB][Q1OOQ2]=⎣⎢⎢⎡ItOIsO⎦⎥⎥⎤以及[P1OOP2][ACOB][Q1OOQ2]=⎣⎢⎢⎡ItD1D3OD2D4IsO⎦⎥⎥⎤其中[D1D3D2D4]=P2CQ1。由於左乘/右乘可逆矩陣不改變矩陣的秩,上面兩個式子等號右端的分塊矩陣有相同的秩。第二式等號右端的矩陣可通過分塊行初等變換消去子塊D1,D3,再通過分塊列初等變換消去子塊D2,於是r⎣⎢⎢⎡ItOIsO⎦⎥⎥⎤=r⎣⎢⎢⎡ItD1D3OD2D4IsO⎦⎥⎥⎤=r⎣⎢⎢⎡ItOD4IsO⎦⎥⎥⎤故t+s=r⎣⎢⎢⎡ItOIsO⎦⎥⎥⎤=r⎣⎢⎢⎡ItOD4IsO⎦⎥⎥⎤=t+s+r(D4)這就得到r(D4)=0,D4=O,P2CQ1=[D1D3D2O]。令S=[D1D3OO],T=[OOD2O],則有SP1AQ1+P2BQ2T=S+T=P2CQ1。兩端左乘P2−1,右乘Q1−1得P2−1SP1A+BQ2TQ1−1=C,令X=P2−1SP1,Y=Q2TQ1−1就有XA+BY=C。
【注】這個引理的證明可以說是充分利用了分塊矩陣的技巧。結論的意義不僅在於給出Frobenius不等式取等號的充要條件(如下),更重要的是這個引理給出了一類矩陣方程有解的充要條件,這類矩陣方程具有形式XA+BY=C,稱爲Lyapunov方程,是一種雙邊多項式矩陣方程(bilateral polynominal matrix equation),其中X,Y是未知量。這種方程經常在自動控制理論、線性系統理論中遇到。
-
公式4(Frobenius不等式):設A∈Fm×n,B∈Fn×s,C∈Fs×t,則r(ABC)⩾r(AB)+r(BC)−r(B)。取等號的充要條件爲存在矩陣X,Y使得B=XAB+BCY。
證:先證明不等式,再考慮取等號的條件。
法1:分塊矩陣初等變換
r(AB)+r(BC)=r[ABOOBC]⩽r[ABBOBC]=行變換r[OB−ABCBC]=列變換r[OB−ABCO]=r(B)+r(ABC)或將上述初等變換過程用分塊初等矩陣寫出,即[ImO−AIn][ABBOBC][IsO−CIt]=[OB−ABCO]於是r(ABC)+r(B)=r[OB−ABCO]=r[ABBOBC]⩾r(AB)+r(BC)。
【注】令B=I,可見Sylvester不等式只是Frobenius不等式的一個特例,因此還可從推廣Sylvester不等式的角度出發證明Frobenius不等式:
法2:利用Sylvester不等式
若r(B)=0,顯然結論成立;若r(B)=0,考慮B的一個滿秩分解B=KL,其中K,L分別爲列滿秩矩陣(有r(B)列)和行滿秩矩陣(有r(B)行),則r(ABC)=r(AKLC)⩾r(AK)+r(LC)−r(B)=r(AKL)+r(KLC)−r(B)=r(AB)+r(BC)−r(B)。
.
【取等號的充要條件】
從法1的證明過程容易看出,取等號的充要條件爲r[ABOOBC]=r[ABBOBC]於是由引理立即可得,取等號的充要條件爲存在X,Y使得B=XAB+BCY。
-
公式5:r(I−AB)⩽r(I−A)+r(I−B)
證:利用公式1和2
r(I−AB)=r(I−A+A−AB)⩽r(I−A)+r(A(I−B))⩽r(I−A)+r(I−B)
【注】這個感覺更像是一個課後習題,不過這個公式從形式上來看挺有趣的~
-
公式6:設A和B爲同階方陣,滿足AB=BA,則r(A+B)⩽r(A)+r(B)−r(AB)
證:分塊矩陣的初等變換+分塊矩陣恆等式的構造
注意到下述分塊矩陣恆等式[−(A+B)OAI][AAAA+B]=[−BAAOA+B]於是r(A)+r(B)=r[AOOB]=行變換r[AAOB]=列變換r[AAAA+B]⩾r[−BAAOA+B](公式2)⩾r[−BAOOA+B]=r(BA)+r(A+B)=r(AB)+r(A+B)【注】這個公式可以看做是在A與B可交換(也就是AB=BA)的條件下,公式1的一個加強版。與公式1相比,這個公式顯然更加緊湊。
-
公式7:r(AHA)=r(AAH)=r(A)
證:只需證r(AHA)=r(A),因爲由r(AHA)=r(A)可推出r(AAH)=r(AH)=r(A)。
法1:線性方程組
考慮兩個齊次線性方程組AHAx=0和Ax=0,顯然Ax=0⇒AHAx=0,現證明AHAx=0⇒Ax=0:用xH左乘AHAx=0的兩端得到xHAHAx=(Ax)H(Ax)=∑i∣(Ax)i∣2=0,其中(Ax)i是Ax的第i個分量,故Ax=0。這說明AHAx=0和Ax=0是同解方程組,故它們的基礎解系所含解向量個數相等,即n−r(AHA)=n−r(A),故r(AHA)=r(A)。
【題外話】通過齊次線性方程組的基礎解系分析矩陣的秩在公式7的證明上不得不說是十分簡潔、巧妙的,網上大多都是這種解法。然而,這並不意味着公式7沒有其他證明思路,下面舉兩例。
法2:向量組
取A列向量組的一個極大無關組a1,a2,...,ar,其中r=r(A),並令K=[a1a2⋯ar],考察AHa1,⋯,AHar是否是AHA列向量組的一個極大無關組即可。考慮齊次線性方程組KHKx=0的解,兩端左乘xH得到xHKHKx=∣∣Kx∣∣22=0,從而Kx=0。根據K的列向量組線性無關,得知x=0,即只有零解。根據線性無關的定義知KHK的列向量組KHa1,KHa2,⋯,KHar是線性無關的,從而延伸組AHa1,AHa2,⋯,AHar也線性無關。假設AHA有多於r個線性無關的列向量,可設爲AHai1,⋯,AHais,s>r,並設ki1ai1+⋯+kisais=0,兩端左乘AH得ki1AHai1+⋯+kisAHais=0,從而ki1=⋯=kis=0,這說明ai1,⋯,ais是線性無關的,而這與r(A)=r矛盾,因此AHA至多有r個線性無關的列向量。這就證明了AHa1,⋯,AHar是AHA列向量組的一個極大無關組,從而r(AHA)=r=r(A)。
法3:滿秩分解+特徵值
先證明Am×n列滿秩的情形。由A列滿秩,A有一個n階非零子式,其對應的子矩陣是可逆的,不失一般性,設A=[BC],其中B是可逆子矩陣。則AHA=BHB+CHC=BH(I+(CB−1)HCB−1)B,注意到(CB−1)HCB−1的特徵值是非負的,從而I+(CB−1)HCB−1的特徵值都不小於1,I+(CB−1)HCB−1可逆,進而BH(I+(CB−1)HCB−1)B即AHA也是可逆的,於是r(AHA)=r(A)。
接下來證明一般情形。若A=O則結論顯然成立,否則可設A=KL是A的一個滿秩分解,則AHA=LH(KHK)L,因爲L行滿秩,LH列滿秩,故r(AHA)=r((KHK)L)=r(KHK)。由已證情形的結論知KHK是可逆的r(A)階方陣,因此r(AHA)=r(A)。
【注】證明中用到的特徵值相關結論見矩陣論(零):線性代數基礎知識整理(5)——特徵值與相似。
-
公式8:設A∈Fm×n,X∈Fn×k,Y∈Fm×k。若W=YTAX非奇異,則r(A−AXW−1YTA)=r(A)−r(AXW−1YTA)
證:
法1:注意到這個公式不過是公式1取等號情形的一個特殊情況,因此可以利用公式1取等號的充要條件。
記B=AXW−1YTA,原公式即r(A−B)+r(B)=r(A)。只要證明存在矩陣C,D使得C(A−B)=A−B,CB=O,(A−B)D=O,BD=B,結論便成立。取C=I−AXW−1YT,D=XW−1YTA,證畢。
法2:使用分塊矩陣初等變換技術直接證明
記B=AXW−1YTA。實際上是要證r(A−B)+r(B)=r(A),故考慮下述分塊初等變換過程:r(A−B)+r(B)=r[A−BOOB]=列變換r[A−BBOB]=行變換r[ABBB]=列變換r[ABOB(I−XW−1YTA)]=行變換r[AOOB(I−XW−1YTA)]注意到右下角子塊B(I−XW−1YTA)=B−AXW−1(YTAX)W−1YTA=B−B=O,故r(A−B)+r(B)=r(A)。證畢。
-
推論(Wedderburn秩1化簡公式,Wedderburn rank-one reduction formula):設A∈Fm×n,x∈Fn,y∈Fm。若w=yTAx=0,則r(A−w−1AxyTA)=r(A)−1
證:
由公式8可得r(A−w−1AxyTA)=r(A)−r(w−1AxyTA),所以這裏的關鍵是證明r(w−1AxyTA)=r(AxyTA)=1。由公式2可得r(AxyTA)⩽r(x)=1(注意r(x)=0,因爲如果這樣的話就有x=0,則AxyTA=O與r(AxyTA)=1矛盾),故只需證明r(AxyTA)⩾1。實際上,假設r(AxyTA)=0,即AxyTA=O,則w=yT(w−1AxyTA)x=0,與已知矛盾,故假設不成立。得證。
【注】這個公式是說,給矩陣A一個特定形式的秩1擾動,則A的秩的改變量恰好也是1。而這對於一般的秩1擾動並不成立。
-
公式9:設A是n階方陣,則r(An)=r(An+1)=r(An+2)=...
【不算是題外話】根據公式2,我們知道任意k=1,2,...都有r(Ak+1)⩽r(Ak),也就是說,A的冪的秩隨着冪指數增加是單調減的。如果我們把r(Ak+1)<r(Ak)(也就是當冪指數增加1,A的冪的秩減少了)稱作是一次“衰減”,那麼這種“衰減”最多可以發生n次(n=A的階數),否則A的冪的秩就會變成負的。那麼問題來了,既然衰減是有限次的,什麼時候衰減會結束?也就是說,當冪指數增加到多大以後,就再也不會發生衰減?公式9將告訴我們,衰減全部在指數達到n之前發生,指數達到n以後不可能再發生衰減。
證:
法1:線性方程組
考慮齊次線性方程組Anx=0和An+1x=0,顯然前者的解都是後者的解,現在證明後者的解也都是前者的解。假設An+1x=0有一解x0不是Anx=0的解,即Anx0=0,那麼向量組An−1x0,An−2x0,...,Ax0,x0中每個向量一定都非零(否則就可推出Anx0=0,矛盾)。設k0x0+k1Ax0+...+knAnx0=0,其中k0,k1,...,kn爲常數。該式兩端左乘An,可得k0Anx0=0,根據Anx0=0知k0=0,等式變爲k1Ax0+k2A2x0+...+knAnx0=0。同理,該式兩端左乘An−1,可推出k1=0。採用類似的做法,最終可得到k0=k1=...=kn=0,這就證明了Anx0,An−2x0,...,Ax0,x0是線性無關的。一個包含n+1個向量的n維向量組是線性無關的,這是矛盾的,因此假設不成立,An+1x=0的解一定都是Anx=0的解。兩方程組同解,則它們的基礎解系所含解向量個數相等,即n−r(An+1)=n−r(An),故r(An)=r(An+1)。同理可證任意m>n有r(Am)=r(Am+1),因此結論成立。證畢。
法2:分塊矩陣初等變換
如果r(An)=0或者r(A)=n,則結論顯然成立。現考慮r(An)>0且r(A)<n的情形。考慮不等式鏈r(An)⩽r(An−1)⩽...⩽r(A)<n,式中的⩽不可能全部取<,否則將有r(An)⩽0與r(An)>0矛盾。因此,存在正整數1⩽p⩽n−1使得r(Ap)=r(Ap+1)。r(Ap+1)=r(Ap)=r[OAp]=列變換r[Ap+1Ap]根據定理21(見後文),存在矩陣Y使得Ap+1Y=Ap,兩端左乘An−p得到An+1Y=An。r(An)=r[OAn]=列變換r[An+1An]=r[An+1An+1Y]=列變換r[An+1O]=r(An+1)如果式Ap+1Y=Ap兩端左乘Am−p並運用同樣的初等變換,就得證任意m⩾n有r(Am)=r(Am+1),因此結論成立。證畢。