陳恒新
(華僑大學(xué)數(shù)學(xué)科學(xué)學(xué)院,福建泉州 362021)
可正定化矩陣的判別定理
陳恒新
(華僑大學(xué)數(shù)學(xué)科學(xué)學(xué)院,福建泉州 362021)
對(duì)有關(guān)可正定化矩陣的理論做進(jìn)一步的研究,給出有關(guān)可正定化矩陣的充分必要性定理.有關(guān)可正定化矩陣的主要判別定理是構(gòu)造性的,即相關(guān)的對(duì)角陣D0,D*是可由矩陣A的元素確定構(gòu)造的.數(shù)值例子表明,定理具有較好的實(shí)用性.
可正定化矩陣;判別定理;充分必要性;構(gòu)造性
對(duì)于解線性方程組Ax=f的許多迭代法,當(dāng)系數(shù)矩陣A正定時(shí)的收斂性定理可直接推廣到A為可正定化矩陣[1-4].關(guān)于可正定化矩陣,文獻(xiàn)[1-2]從理論上進(jìn)行研究,給出了一些相關(guān)定理及判定方法.然而,目前關(guān)于這一類問(wèn)題的研究尚不夠深入.為此本文對(duì)有關(guān)可正定化矩陣的理論做進(jìn)一步的研究,給出了一些可正定化矩陣的充分必要性定理.
為敘述簡(jiǎn)便,先引入如下記號(hào).n階實(shí)矩陣A=[ai,j];對(duì)角陣 D=diag(d1,d2,…,dn),若di>0,i= 1,2,…,n,則稱D為正對(duì)角陣;集合N={1,2,…,n}.
定義1 若存在對(duì)角陣 P和Q使PAQ為正定矩陣,則稱A為可正定化矩陣.
由文獻(xiàn)[1]的定理1,2可知有如下引理.
引理1n階矩陣A是可正定化矩陣的充分必要條件是存在對(duì)角陣D,使DA為正定矩陣.
引理2n階矩陣A是可正定化矩陣的充分必要條件是存在正對(duì)角陣D,使DA+為正定矩陣.
顯然,由引理1,2可知有如下引理.
引理3n階矩陣A是可正定化矩陣的充分必要條件是A+為可正定化矩陣.
由文獻(xiàn)[3]的定理5,6可知有如下引理.
引理4 實(shí)對(duì)稱矩陣A正定的充分必要條件是,它的所有順序主子式det Ak>0,k=1,2,…,n.
對(duì)任正對(duì)角陣D,DA+非對(duì)稱;若不然,DA+為對(duì)稱矩陣.因?yàn)?DA+=[di~ai,j],所以有=
顯然,若的某一對(duì)對(duì)稱元素不同時(shí)為零(即一個(gè)為零,另一個(gè)非零),則對(duì)任正對(duì)角陣D,DA+非對(duì)稱.
若矩陣A有零對(duì)角元,則DA亦有零對(duì)角元.綜上所述,由引理5及引理1,2可得如下定理.
定理1 若n階矩陣A有零對(duì)角元或其A+的某一對(duì)對(duì)稱元素異號(hào)或不同時(shí)為零,則A不是可正定化矩陣.
由引理3及定理1可知,判別n階矩陣A是否為可正定化矩陣,只需考察相應(yīng)的主對(duì)角元全為正的,且有如下所定義的矩陣.
定義2 若n階矩陣A=[ai,j]的主對(duì)角元全為正,即ai,i>1,i=1,2,…,n,且A的對(duì)稱元素ai,j與aj,i(i≠j),i,j=1,2,…,n同號(hào)或同時(shí)為零.即ai,j aj,i≥0且ai,j aj,i=0,當(dāng)且僅當(dāng)ai,j=aj,i=0,則記為A∈S+.
引理6 設(shè)n階矩陣A=[ai,j]中存在某一個(gè)l∈N,使得al,j=aj,l≠0,j=1,2,…,l-1,l+1,…,n.若A不是對(duì)稱矩陣,則對(duì)任非奇異對(duì)角陣D,DA不是對(duì)稱矩陣.
證明 因?yàn)锳是非對(duì)稱矩陣,則必存在i0≠j0,使≠
于是有 D=diag(d,d,…,d),d≠0.因此,DA=[dai,j].由于 DA為對(duì)稱陣,則有=故可得=與≠矛盾.
所以可知,不存在非奇異對(duì)角陣D,使DA為對(duì)稱矩陣.證畢.又由已知有d1d2…dk>0,可知引理8成立.證畢.
定理2 設(shè)n階矩陣A∈S+,且A=[ai,j]中存在al,j≠0,aj,l≠0,j=1,2,…,n(l為 N中的某一個(gè)數(shù)),則A為可正定化矩陣的充分必要條件是D0A為正定矩陣,而D0與引理7中 D0=diag(d1,d2,…, dn)所確定的對(duì)角陣相同.
證明 因A∈S+,可知D0為正對(duì)角陣,且矩陣A相應(yīng)的A+=A.
充分性.若D0A為正定矩陣,由引理2可知A是可正定化矩陣.
必要性.若A是可正定化矩陣.
反證 假設(shè)D0A不是正定矩陣,則有下列2種情況.
(2)D0A為對(duì)稱矩陣,但D0A不為正定矩陣.
由引理4可知,必存在某一個(gè)k,使det(D0A)k≤0,1,≤k≤n.因 D0為正對(duì)角陣,由引理8可知有det Ak≤0,1≤k≤n.但因A是可正定化矩陣,由引理2可知,必存在正對(duì)角陣D使DA+=DA為正定矩陣.于是,由引理4可知有det(DA)k>0,k=1,2,…,n.又由引理8可知有det Ak>0,k=1,2,…,n.這與det Ak≤0,1≤k≤n矛盾.即情況(2)亦不可能出現(xiàn).因此,可知D0A為正定矩陣.證畢.
引理9 假設(shè)n階矩陣A的三對(duì)角矩陣Atrd為兩條斜對(duì)角線元素皆不為零的對(duì)稱陣.即ai,i+1=ai+1,i≠0,i=1,2,…,n-1.若A不是對(duì)稱矩陣,則對(duì)任非奇異對(duì)角陣D,DA不是對(duì)稱矩陣.
證明 因A為非對(duì)稱矩陣,則必存在i0≠j0,使ai0j0≠aj0i0.
于是有 D=diag(d,d,…,d),d≠0,故 DA=[dai,j].因?yàn)?DA為對(duì)稱矩陣,則有dai0,j0=daj0,i0,由此可得ai0,j0=aj0,i0與ai0,j0≠aj0,i0矛盾.所以可知,不存在非奇異對(duì)角陣 D使DA為對(duì)稱矩陣.證畢.
其次,若ai,i+1≠0,ai+1,i≠0,i=1,2,…,n-1中ai,i+1·ai+1,i>0,則 D*為正對(duì)角陣.
據(jù)引理9,10及相關(guān)引理,完全類似定理2的證明,便有下列定理.
定理3 設(shè)n階矩陣A∈S+,且A的三對(duì)角矩陣Atrd中兩條斜對(duì)角線元素皆不為零,即ai,i+1≠0, ai+1,i≠0,i=1,2,…,n-1.則A是可正定化矩陣的充分必要條件,D*A為正定矩陣.其中:D*同引理10中 D*=diag(d1,d2,…,dn)所確定的對(duì)角陣.
[1]王偉賢,王志偉.關(guān)于可正定化矩陣的判定[J].數(shù)值計(jì)算與計(jì)算機(jī)應(yīng)用,1999,20(3):215-222.
[2]胡家贛,劉興平.EPEk方法和可正定化矩陣[J].數(shù)值計(jì)算與計(jì)算機(jī)應(yīng)用,1997,18(1):30-39.
[3]蔣爾雄,高坤敏,吳景琨.線性代數(shù)[M].北京:人民教育出版社,1978.
[4]陳恒新.關(guān)于非負(fù)矩陣Perron特征值的上、下界[J].應(yīng)用數(shù)學(xué)與計(jì)算數(shù)學(xué)學(xué)報(bào),2007,21(1):1-8
(責(zé)任編輯:陳志賢英文審校:張金順,黃心中)
Criteria theorem of Positive-Definable Matrix
CHEN Heng-xin
(School of Mathematical Sciences,Huaqiao University,Quanzhou 362021,China)
Theory of the positive-definable matrix be further studied.The necessary and sufficient theo rem s for the positive-definable matrix are given in the paper.These theorem s have better uses than the present critical theorems,that is, the co rrelated diagonal matrix D0,D*in our theorem s can be structured determinately by the elements of the matrix A. For this reason,it has good practical value.Four numericial examples are given here,that show s these theorems had better practical uses.
positive-definable matrix;critical theorem;necessity and sufficiency;structure
O 241.6
A
1000-5013(2011)03-0356-05
2009-12-16
陳恒新(1956-),男,副教授,主要從事計(jì)算數(shù)學(xué)和數(shù)值代數(shù)的研究.E-mail:chenhx@hqu.edu.cn.
福建省自然科學(xué)基金計(jì)劃資助項(xiàng)目(S0650018)