肖麗霞,張俊麗

(內(nèi)蒙古民族大學(xué)數(shù)學(xué)學(xué)院,內(nèi)蒙古 通遼 028043)

非奇異H-矩陣新的含參數(shù)細(xì)分迭代判別法

肖麗霞,張俊麗

(內(nèi)蒙古民族大學(xué)數(shù)學(xué)學(xué)院,內(nèi)蒙古 通遼 028043)

結(jié)合矩陣自身的元素,構(gòu)造了含參數(shù)的迭代公式,進(jìn)而細(xì)分了矩陣非對角占優(yōu)行指標(biāo)集.利用廣義嚴(yán)格α-對角占優(yōu)矩陣與非奇異H-矩陣的關(guān)系,給出了非奇異H-矩陣一組新的細(xì)分迭代判定準(zhǔn)則,推廣和改進(jìn)了已有的結(jié)果,通過數(shù)值算例說明了結(jié)果的優(yōu)越性.

非奇異H-矩陣;α-對角占優(yōu)矩陣;不可約;非零元素鏈

1 引言與符號

非奇異H-矩陣是一類應(yīng)用廣泛的特殊矩陣,在計(jì)算數(shù)學(xué)、矩陣?yán)碚?、控制論等領(lǐng)域發(fā)揮著重要作用.對非奇異H-矩陣判定方法的研究,近年來引起許多數(shù)學(xué)工作者的關(guān)注,并取得了一系列的研究成果[19].本文給出了一類新的含參數(shù)細(xì)分迭代判別法,對文獻(xiàn)[1-5]的結(jié)果進(jìn)行了推廣和改進(jìn).

為敘述方便,引進(jìn)下列記號和定義,設(shè)A=(aij)∈Cn×n為n階復(fù)方陣,N?{1,2,···,n}, α∈[0,1].記

顯然有δk+1,i≤rk+1≤rk≤r1<1(?k∈Z+,i∈N2).

定義 1.1[6]設(shè)A=(aij)∈Cn×n,若|aii|≥(>)Ri(A)(?i∈N),則稱A為(嚴(yán)格)對角占優(yōu)矩陣,記為A∈D0(A∈D);若存在正對角矩陣X使得AX∈D,則稱A為廣義嚴(yán)格對角占優(yōu)矩陣,記為A∈D?.

定義 1.2[6]設(shè)A=(aij)∈Cn×n,若存在α∈[0,1],使得

則稱A為α-(嚴(yán)格)對角占優(yōu)矩陣,記為A∈D0(α)(A∈D(α));若存在正對角矩陣X 使得AX∈D(α),則稱A為廣義嚴(yán)格α-對角占優(yōu)矩陣,記為A∈D?(α).

2 主要結(jié)果

引理 2.1 [7]設(shè)A=(aij)∈Cn×n,若A∈D(α),則A∈D?.

引理2.2 [8]設(shè)A=(aij)∈Cn×n,若存在正對角矩陣X,使AX∈D?,則A∈D?.

定理 2.1設(shè)A=(aij)∈Cn×n,且N′1,N′2/=?,若存在k∈Z+(?i∈N′1,?j∈N′2),滿足:

則A為非奇異H-矩陣.

證明令

因?yàn)?

取正對角矩陣X1=diag(d1,d2,···,dn),并記B=AX1,其中

3) ?i∈N2,可得

綜上所述,有|bii|>αRi(B)+(1?α)Qi(B)(?i∈N)成立,則B∈D(α),由引理2.1知B=AX1∈D?,其中X1為正對角矩陣,根據(jù)引理2.2,則A∈D?,因此矩陣A為非奇異H-矩陣.

引理 2.3 [7]設(shè)A=(aij)∈Cn×n,若A∈D0(α),A不可約,且N2/=?,則A∈D?.

定理 2.2設(shè)A=(aij)∈Cn×n,且A不可約,若存在滿足:

且至少有一嚴(yán)格不等式成立,則A為非奇異H-矩陣.

證明如同定理2.1的證明,記Mi,mj,則因?yàn)?

取正對角矩陣X2=diag(d1,d2,···,dn),并記B=AX2,其中

類似于定理2.1的證明過程,可得|bii|≥αRi(B)+(1?α)Qi(B)(?i∈N),且至少有一個嚴(yán)格不等式成立.由A不可約,可得B不可約,則B為不可約α-對角占優(yōu)矩陣,由引理2.3可知B=AX2∈D?,其中X2為正對角矩陣,根據(jù)引理2.1,則A∈D?,因此矩陣A為非奇異H-矩陣.

記

騰訊董事會主席馬化騰憑借328億美元的身家蟬聯(lián)榜單第二名，但他的財(cái)富縮水了62億美元。去年的首富、中國恒大董事局主席許家印的排名跌至第三名。許家印的身家為308億美元，下降28%，折合約117億美元，他是今年財(cái)富值降低最多的富豪。

引理 2.4 [7]設(shè) A=(aij)∈Cn×n,若 A ∈D0(α),并且 ?i∈N3,都有非零元素鏈aik1ak1k2...akpj,使得j∈N2,則A∈D?.

定理 2.3設(shè)A=(aij)∈Cn×n,且N′1,N′2/=?,若存在k∈Z+(?i∈N′1,?j∈N′2),滿足:

證明如同定理2.2的證明,可得|bii|≥αRi(B)+(1?α)Qi(B)(?i∈N).其中

成立;

成立,則B為非零元素鏈對角占優(yōu)矩陣.根據(jù)引理2.4,B=AX2∈D?,其中X2為正對角矩陣.根據(jù)引理2.1,則A∈D?,因此矩陣A為非奇異H-矩陣.

3 數(shù)值算法

輸入:已知矩陣A,參數(shù)α,迭代次數(shù)k.

輸出:正對角矩陣X.

1)若aii=0(?i∈N)或N2=?,輸出:矩陣A不是非奇異H-矩陣”,停止;若N1=?,輸出:矩陣A是非奇異H-矩陣”,停止;否則執(zhí)行2);

2)若N1/=?且N2/=?,

成立,則輸出“矩陣A是非奇異H-矩陣”,停止;否則輸出“不確定矩陣A是否為非奇異H-矩陣”,停止.

4 數(shù)值算例

例4.1設(shè)

取α=0.5,則N1={1,2,3},N2={4,5},令k=1,可得r0=1,r1=0.5,δ2,4=0.45, δ2,5=0.2375,r2=0.45,則取i=3,j=1,則

取i=3,j=2,則

可見矩陣A滿足定理2.1的條件,因此矩陣A為非奇異H-矩陣.但

所以矩陣A不能由文獻(xiàn)[1]中定理2判定.又因?yàn)?/p>

所以矩陣A不能由文獻(xiàn)[2]中定理1判定.又因?yàn)?/p>

所以矩陣A不能由文獻(xiàn)[3]中定理1判定.又因?yàn)?/p>

所以矩陣A不能由文獻(xiàn)[4]中定理1判定.經(jīng)驗(yàn)證,對α=0.1,0.2,···,1.0,矩陣A不能由文獻(xiàn)[5]中定理1判定.

[1]黃廷祝.非奇H矩陣的簡捷判據(jù)[J].計(jì)算數(shù)學(xué),1993,15(3):318-328.

[2]高中喜,黃廷祝,王廣彬.非奇H-矩陣的充分條件[J].數(shù)學(xué)物理學(xué)報(bào):A輯,2005,25(3):409-413.

[3]黃澤軍,劉建州.非奇異H矩陣的一類新迭代判別法[J].工程數(shù)學(xué)學(xué)報(bào),2008,25(5):939-942.

[4]孫德淑.非奇異H-矩陣的判定準(zhǔn)則[J].溫州大學(xué)學(xué)報(bào),2009,30(3):18-21.

[5]尹如軍,徐仲,陸全.非奇H-矩陣的細(xì)分迭代判別準(zhǔn)則[J].工程數(shù)學(xué)學(xué)報(bào),2013,30(3):433-441.

[6]孫玉祥.廣義對角占優(yōu)矩陣的充分條件[J].高等學(xué)校計(jì)算數(shù)學(xué)學(xué)報(bào),1997,19(3):216-223.

[7]Sun Yuxiang.An improvement on a theorem by Ostrowski and its applications[J].Northeastern Math.J., 1991,7(4):497-502.

[8]Berman A,Plemmons R J.Nonnegative Matrices in the Mathematical Sciences[M].Philadephia:SIAM Press,1994.

[9]韓貴春,錢茜,張俊麗.Ostrowski定理的推廣與非奇異 H-矩陣的實(shí)用判定 [J].純粹數(shù)學(xué)與應(yīng)用數(shù)學(xué), 2013,29(6):601-608.

New subdivided and iterative criteria with parameter for nonsingular H-matrices

Xiao Lixia,Zhang Junli
(School of Mathematics,Inner Mongolia University for the Nationalities,Tongliao 028043,China)

Associating the elements of the matrix,the iterative formulas with parameter are constructed,and then the index set of non diagonally dominant rows in a square matrix is subdivided.According to the relations between generalized α-diagonally dominant matrices and nonsingular H-matrices,a set of new subdividing and iterative criteria for nonsingular H-matrices is obtained,which extend and improve some related results.A numerical example is used to show the advantages of the results.

nonsingular H-matrix,α-diagonally dominant matrix,irreducible,non-zero elements chain

O151.21

A

1008-5513(2014)004-0386-07

10.3969/j.issn.1008-5513.2014.04.008

2014-05-25.

內(nèi)蒙古自治區(qū)高等學(xué)?？茖W(xué)技術(shù)研究項(xiàng)目(NJZY13159).

肖麗霞(1980-),碩士,講師,研究方向：數(shù)值代數(shù).

2010 MSC:15A57