矩陣的秩與非零特征值個(gè)數(shù)差的確定

呂洪斌, 楊忠鵬, 馮曉霞, 陳梅香, 梁小春

(1.北華大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院, 吉林 吉林 132013； 2.莆田學(xué)院 數(shù)學(xué)學(xué)院, 福建 莆田 351100；3.閩南師范大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院, 福建 漳州 363000； 4.福州大學(xué) 數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院, 福州 350108)

研究簡(jiǎn)報(bào)

矩陣的秩與非零特征值個(gè)數(shù)差的確定

呂洪斌1, 楊忠鵬2,3, 馮曉霞3, 陳梅香2, 梁小春4

(1.北華大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院, 吉林 吉林 132013； 2.莆田學(xué)院 數(shù)學(xué)學(xué)院, 福建 莆田 351100；
3.閩南師范大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院, 福建 漳州 363000； 4.福州大學(xué) 數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院, 福州 350108)

以矩陣的Jordan標(biāo)準(zhǔn)形為工具, 給出了用矩陣方冪的秩表示的矩陣的秩和非零特征值個(gè)數(shù)差的確定方法, 其結(jié)果不依賴于矩陣的Jordan標(biāo)準(zhǔn)形.

矩陣秩； 矩陣方冪； 矩陣指數(shù)； 冪零矩陣

秩與非零特征值個(gè)數(shù)相等的矩陣在數(shù)理統(tǒng)計(jì)、試驗(yàn)設(shè)計(jì)、多元統(tǒng)計(jì)分析、金融計(jì)量統(tǒng)計(jì)和經(jīng)濟(jì)統(tǒng)計(jì)分析等領(lǐng)域應(yīng)用廣泛.矩陣的Jordan標(biāo)準(zhǔn)形和其初等因子是互為確定的, 因此目前已有關(guān)于矩陣的秩和非零特征值個(gè)數(shù)關(guān)系的研究結(jié)果都依賴于矩陣的Jordan標(biāo)準(zhǔn)形[1-11].本文給出矩陣的秩和非零特征值個(gè)數(shù)差的確定方法, 得到的結(jié)果不依賴于矩陣的Jordan標(biāo)準(zhǔn)形, 而以矩陣方冪的秩為基本工具.文獻(xiàn)[2]從矩陣Jordan標(biāo)準(zhǔn)形出發(fā), 得到:

命題1[2]設(shè)A∈n×n, 如果A的Jordan標(biāo)準(zhǔn)形中共有t個(gè)階數(shù)分別為n1,n2,…,nt的冪零Jordan塊, 則0≤r(A)-u(A)=n1+n2+…+nt-t≤n-1.

1 預(yù)備知識(shí)

參照文獻(xiàn)[2-3], 當(dāng)A∈n×n時(shí), 總設(shè)

其中:P∈n×n可逆;

Jj=λjEnj+Nnj∈nj×nj,λj≠0,j=t+1,t+2,…,s(3)

約定當(dāng)ni=1時(shí)Ji是1階的, 因此1階冪零Jordan塊可記為N1(=0)∈； 當(dāng)t=0時(shí),ni=0,i=1,2,…,t； 顯然n.由文獻(xiàn)[1]或文獻(xiàn)[2]中引理3知k=n1=indA, 因此式(1)中冪零Jordan塊的不同階數(shù)最多有k=n1=indA個(gè).當(dāng)A∈n×n且k=indA時(shí), 總設(shè)

其中di∈(i=1,2,…,k)為冪零Jordan塊不同階數(shù)的個(gè)數(shù), 顯然d1+d2+…+dk=t.

易得下述引理.

引理1設(shè)A∈n×n, 則:

引理2設(shè)A∈n×n, 則:

1) indA=0 ?A可逆, 即r(A0)=r(Am)=n, ?m∈；

2) 1≤indA=k=n1? 存在1≤k≤m∈, 使得r(Ak-1)>r(Ak)=r(Am).

引理2表明indA即為JA中冪零Jordan塊的最大階數(shù), 從而對(duì)冪零矩陣應(yīng)用引理2即可得文獻(xiàn)[4]中命題2.5; 引理2中的2)表明可不依賴Jordan標(biāo)準(zhǔn)形而只使用矩陣方冪的秩也可以確定矩陣的指數(shù).

引理3設(shè)A∈n×n, 則:

1) 如果存在1≤l≤n-1, 使得r(Al)≠r(Al+1), 則l

2) 當(dāng)r(Al)=r(Al+1)時(shí), 則indA≤l；

3) 當(dāng)1≤k=indA時(shí),r(A)>r(A2)>…>r(Ak-1)>r(Ak)=r(Am),k≤m.

文獻(xiàn)[9]證明了當(dāng)A∈n×n為n-冪零矩陣時(shí),r(Al)=n-l(l=1,2,…,n).實(shí)際上, 由An=0且An-1≠0即A∈Wn和引理2中2)知,n-1r(A2)>…>r(An-1)>r(An)=0, 即得文獻(xiàn)[9]的結(jié)果.

引理4設(shè)A∈n×n,mA(x)=xkh(x),h(x)∈[x]且h(0)≠0, 則indA=k.

引理5設(shè)A∈n×n, 則A∈Wn?u(A)=0.

引理6設(shè)A∈Wn, 則A為k-冪零的 ? indA=k.

2 主要結(jié)果

定理1設(shè)A∈n×n, 如果indA=k, 則:

證明： 當(dāng)indA=k=0時(shí), 由引理2中1)知式(7)成立, 因此只證明r(A)≤n-1的情形.

由引理2及r(A)≤n-1知, 1≤t且r(Am)=r(Ak)

又由式(1),(4)知d1k+d2(k-1)+…+2dk-1+dk=n1+n2+…+nt, 再由式(5),(6),(9)知

從而由命題3可知式(7)成立.

當(dāng)1≤l≤k=indA時(shí), 由式(1)～(6),(9)有

即

再由式(4),(10)有

即

又由r(A0)=r(E)=n=kd1+(k-1)d2+…+2dk-1+dk+r(Ak)和式(7)有

r(A0)-r(A)=n-r(A)=d1+d2+…+dk-1+dk=t,

可見(jiàn)式(11)對(duì)l=1也成立, 即

令k-l+1=f, 則l=k-f+1, 進(jìn)而可得式(12)的等價(jià)形式

又由式(13)有dl=(r(Ak-l)-r(Ak-l+1))-(r(Ak-l+1)-r(Ak-l+2)), 即式(8)成立.證畢.

式(8)表明只用矩陣方冪的秩就可得到由式(7)確定的矩陣的秩和非零特征值個(gè)數(shù)的差.

如果indA=k,A∈n×n, 則有唯一的X滿足XAX=X,AX=XA,Ak+1X=Ak,X=AD稱為A的Drazin逆[5].由A的Drazin逆的唯一性知下述結(jié)論成立：

推論1設(shè)A∈n×n, 當(dāng)存在k(≥1)∈, 使得r(Ak-1)>r(Ak)=r(Am),k≤m時(shí), 則在式(1)～(3)下, 有

推論1在應(yīng)用上比文獻(xiàn)[1]中定理2.5及文獻(xiàn)[6]中定理2.1.2更方便些.

定理2設(shè)A∈n×n, 則下述條件等價(jià)：

1)r(A)=u(A)；

2)A可逆或?qū)儆诹闾卣髦档腏ordan塊的階數(shù)都等于1；

3) indA=0或1且d1=n-r(A)；

4)mA(x)=h(x)或mA(x)=xh(x)且h(0)≠0；

5)r(A)=r(A2)；

6) 對(duì)任意的m(≥2)∈, 都有r(A)=r(Am).

證明: 1)?2).由式(7),(9)知indA=k, 即r(A)=r(Ak).由引理3中2)知indA≤1； 當(dāng)indA=0時(shí), 由引理2中1)知A可逆； 當(dāng)indA=1時(shí), 由引理2中2)得n1=indA=1, 即A屬于零特征值的Jordan塊的階數(shù)都等于1.

2)?3).當(dāng)A可逆時(shí), 由引理2中1)知indA=0且d1=0=n-r(A)； 當(dāng)A屬于零特征值的Jordan塊的階數(shù)都等于1時(shí), 由引理2中2)得indA=n1=1, 再由式(8)得,d1=r(A0)-2r(A)+r(A2)=n-r(A).

3)?4).由引理2和引理4可得結(jié)論.

4)?5).由引理2和引理4知indA=n1≤1, 再由引理2得r(A)=r(A2).

5)?6).由r(A)=r(A2)和引理3中2)得indA≤1, 于是由引理2知, 對(duì)任意的m(≥2)∈, 都有r(A)=r(Am).

6)?1).由m≥2任意性知r(A)=r(A2), 從而由引理3中2)得indA=k=0或1, 進(jìn)而由式(7),(9)得r(A)-u(A)=r(A)-r(Ak)=0.證畢.

定理2簡(jiǎn)化了文獻(xiàn)[2]中定理3的證明.

如果AA#A=A,A#AA#=A#,AA#=A#A, 則稱A∈n×n是群可逆的,A#為A的群逆, 此時(shí)也稱A為GP矩陣.由文獻(xiàn)[6-8]知, 不是每個(gè)矩陣都為群可逆的, 當(dāng)然GP矩陣的群逆是唯一的.由文獻(xiàn)[7]中推論5.5.9知,A為GP矩陣 ? indA≤1, 從而定理2給出了GP矩陣的等價(jià)刻畫(huà).

定理3設(shè)A∈n×n,n≥2, 則下述結(jié)論等價(jià)：

1)r(A)-u(A)=n-1；

2)r(A)=n-1,A∈Wn；

3) indA=n,JA=Nn；

4)mA(x)=xn；

5)An-1≠0,An=0；

6)r(Al)=n-l, 1≤l≤n-1且r(Al)=0,n≤l；

7) 對(duì)任意的1≤p,l≤n(p≠l), 總有r(Ap)≠r(Al).

證明: 1)?2)顯然.

2)?3).由A∈Wn和引理5知u(A)=0且n1+n2+…+nt=n, 再由式(7)得

r(A)-u(A)=r(A)-r(Ak)=r(A)=n-t=n-1,

知t=1, 進(jìn)而由引理2知indA=k=n1=n, 即JA=Nn.

3)?4).由引理2和引理4知此時(shí)mA(x)=xn.

4)?5).由mA(x)=xn和引理4知indA=n1=n, 進(jìn)而JA=Nn, 再由式(1),(5)得An-1≠0,An=0.

5)?6).由An-1≠0,An=0, 引理2和引理5知indA=n,u(A)=0； 又由引理3中3)知n-1≥r(A)>r(A2)>…>r(An-1)≥1且r(An)=0, 表明r(Al)=n-l, 1≤l≤n-1且r(Al)=0,n≤l.

6)?7)顯然.

7)?1).由對(duì)1≤p,l≤n(p≠l)總有r(Ap)≠r(Al)知,n-1≥r(A)>r(A2)>…>r(An-1)>r(An)≥0, 表明r(Al)=n-l, 1≤l≤n-1且r(An)=0； 又由引理5知u(A)=0, 進(jìn)而由式(7)知n-1=r(A)-u(A)=r(A). 證畢.

定理3推廣并簡(jiǎn)化了文獻(xiàn)[2]中定理5的相關(guān)結(jié)論.由定理1～定理3可得:

定理4設(shè)A∈n×n, 則:

1)r(A)=u(A) ?r(A)=r(A2)；

2)r(A)-u(A)=n-1 ?r(A)=n-1,r(An)=0 ? indA=n,JA=Nn?A為n-冪零的.

3 冪零矩陣Jordan標(biāo)準(zhǔn)形的確定

定理5設(shè)A∈Wn, 則:

1) 如果存在k-1∈, 使得r(Ak-1)=1, 則indA=k, 即A為k-冪零的；

2)dl=r(Ak-l)-2r(Ak-l+1)+r(Ak-l+2),d1=r(Ak-1), 1≤l≤k=indA.

證明: 當(dāng)r(Ak-1)=1時(shí), 如果0≠r(Ak)≤r(Ak-1)=1, 則r(Ak)=r(Ak-1)=1.由引理3中2)知indA≤k-1, 再由引理6得Ak-1=0, 矛盾； 表明r(Ak)=0, 因此由引理6知indA=k, 即A為k-冪零的. 由式(8)可得2).證畢.

定理6設(shè)A∈Wn, 若存在1≤k∈, 使得r(Ak+1)=r(Ak)

由矩陣相似的充要條件為它們有相同的Jordan標(biāo)準(zhǔn)形及定理5和定理6得:

定理7設(shè)A,B∈Wn, 則:

1)A,B相似 ?r(Al)=r(Bl), 1≤l≤n-1；

2) 如果1≤k=indA, 則A,B相似 ?r(Al)=r(Bl), 1≤l≤k.

文獻(xiàn)[9-11]討論了3-冪零矩陣, 于是由引理6、定理5～定理7, 有:

推論2設(shè)A,B∈Wn都是3-冪零的, 則:

1)A,B相似 ?r(A)=r(B),r(A2)=r(B2)；

2) 當(dāng)A,B相似時(shí),d1=r(A2),d2=r(A)-2r(A2),d3=n-2r(A)+r(A2), 且

[1]Nikuie M, Mirnia M K, Mahmoudi Y.Some Results about the Index of Matrix and Drazin Inverse [J].Mathematical Sciences, 2010, 4(3): 283-294.

[2]梁小春, 陳梅香, 楊忠鵬, 等.矩陣的秩和非零特征值個(gè)數(shù)關(guān)系的進(jìn)一步討論 [J].閩南師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 2014(2): 1-6.(LIANG Xiaochun, CHEN Meixiang, YANG Zhongpeng, et al.Some Further Discussions on the Relationship between Ranks and Numbers of Non-zero Eigenvalues of Matrices [J].Journal of Minnan Normal University: Natural Science, 2014(2): 1-6.)

[3]CHEN Meixiang, Lü Hongbin, FENG Xiaoxia, et al.The Essential (m,l)-Idempotent Matrix and Its Minimal Polynomial [J].International Journal of Applied Mathematics and Statistics, 2013, 41(11): 31-41.

[4]WU Yan, Linder D F.On the Eigenstructures of Functionalk-Potent Matrices and Their Integral Forms [J].WSEAS Transactions on Mathematics, 2010, 9(4)： 244-253.

[5]Bernstein D S.Matrix Mathematics Theory, Facts, and Formulas [M].2nd ed.Princeton: Princeton University Press, 2009.

[6]WANG Guorong, WEI Yimin, QIAO Sanzheng.Generalized Inverse: Theory and Computations [M].Beijing: Science Press, 2004.

[7]ZHANG Fuzhen.Matrix Theory: Basic Results and Techniques [M].2nd ed.New York: Springer-Verlag, 1999.

[8]Baksalary O M, Trenkler G.Onk-Potent Matrices [J].Electronic Journal of Linear Algebra, 2013, 26: 446-470.

[9]李殿龍.3-冪零矩陣的Jordan標(biāo)準(zhǔn)型在GLn(F)共軛作用下的迷向子群 [J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí), 2008, 38(8): 143-146.(LI Dianlong.Isotropy Subgroups of Jordan Normal Form of 3-Nilpotent Matrices under the Conjugate Action ofGLn(F) [J].Mathematics in Practice and Theory, 2008, 38(8): 143-146.)

[10]林榮珍, 江飛.3-冪零矩陣的相似等價(jià)類的計(jì)數(shù) [J].數(shù)學(xué)研究, 2006, 39(4): 394-400.(LIN Rongzhen, JIANG Fei.The Enumeration of Equivalence Class of 3-Nilpotent Matrices [J].Journal of Mathematical Study, 2006, 39(4): 394-400.)

[11]曾月迪, 林麗芳.3-冪零矩陣Jordan規(guī)范型的計(jì)數(shù) [J].莆田學(xué)院學(xué)報(bào), 2013, 20(5): 5-8.(ZENG Yuedi, LIN Lifang.The Enumeration of Jordan’s Normal Form of 3-Nilpotent Matrices [J].Journal of Putian University, 2013, 20(5): 5-8.)

DeterminationofDifferencebetweentheRankandtheNumberofNon-zeroEigenvaluesofaMatrix

Lü Hongbin1, YANG Zhongpeng2,3, FENG Xiaoxia3, CHEN Meixiang2, LIANG Xiaochun4
(1.SchoolofMathematicsandStatistics,BeihuaUniversity,Jilin132013,JilinProvince,China；
2.SchoolofMathematics,PutianUniversity,Putian351100,FujianProvince,China;
3.SchoolofMathematicsandStatistics,MinnanNormalUniversity,Zhangzhou363000,FujianProvince,China;
4.CollegeofMathematicsandComputerScience,FuzhouUniversity,Fuzhou350108,China)

Taking the Jordan canonical form as a basic tool, we showed the determination method of the difference between the rank of matrix’s power and the number of non-zero eigenvalues, with its result independent of the Jordan canonical form.

rank of a matrix; matrix power; matrix index; nilpotent matrix

2014-05-21.

呂洪斌(1964—), 男, 漢族, 博士, 教授, 從事數(shù)值代數(shù)和矩陣?yán)碚摷捌鋺?yīng)用的研究, E-mail: lhb1964@126.com.通信作者： 楊忠鵬(1947—), 男, 漢族, 教授, 從事矩陣及其應(yīng)用的研究, E-mail: yangzhongpeng@126.com.

吉林省教育廳科學(xué)技術(shù)研究項(xiàng)目(批準(zhǔn)號(hào)： 2013192)、福建省教育廳科研基金(批準(zhǔn)號(hào)： JA12286; JA08196)、福建省自然科學(xué)基金(批準(zhǔn)號(hào)： 2010J01018)和2008年福建省高校服務(wù)海西建設(shè)重點(diǎn)項(xiàng)目(批準(zhǔn)號(hào)： 2008HX03).

O151.21

A

1671-5489(2014)06-1210-05

10.13413/j.cnki.jdxblxb.2014.06.20

趙立芹)