肖麗霞,張俊麗
(內(nèi)蒙古民族大學數(shù)學學院,內(nèi)蒙古 通遼 028043)
非奇異H-矩陣新的含參數(shù)細分迭代判別法
肖麗霞,張俊麗
(內(nèi)蒙古民族大學數(shù)學學院,內(nèi)蒙古 通遼 028043)
結(jié)合矩陣自身的元素,構(gòu)造了含參數(shù)的迭代公式,進而細分了矩陣非對角占優(yōu)行指標集.利用廣義嚴格α-對角占優(yōu)矩陣與非奇異H-矩陣的關(guān)系,給出了非奇異H-矩陣一組新的細分迭代判定準則,推廣和改進了已有的結(jié)果,通過數(shù)值算例說明了結(jié)果的優(yōu)越性.
非奇異H-矩陣;α-對角占優(yōu)矩陣;不可約;非零元素鏈
非奇異H-矩陣是一類應用廣泛的特殊矩陣,在計算數(shù)學、矩陣理論、控制論等領(lǐng)域發(fā)揮著重要作用.對非奇異H-矩陣判定方法的研究,近年來引起許多數(shù)學工作者的關(guān)注,并取得了一系列的研究成果[19].本文給出了一類新的含參數(shù)細分迭代判別法,對文獻[1-5]的結(jié)果進行了推廣和改進.
為敘述方便,引進下列記號和定義,設(shè)A=(aij)∈Cn×n為n階復方陣,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ōu)矩陣,記為A∈D0(A∈D);若存在正對角矩陣X使得AX∈D,則稱A為廣義嚴格對角占優(yōu)矩陣,記為A∈D?.
定義 1.2[6]設(shè)A=(aij)∈Cn×n,若存在α∈[0,1],使得
則稱A為α-(嚴格)對角占優(yōu)矩陣,記為A∈D0(α)(A∈D(α));若存在正對角矩陣X 使得AX∈D(α),則稱A為廣義嚴格α-對角占優(yōu)矩陣,記為A∈D?(α).
引理 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-矩陣.
證明令
因為0 取正對角矩陣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不可約,若存在滿足: 且至少有一嚴格不等式成立,則A為非奇異H-矩陣. 證明如同定理2.1的證明,記Mi,mj,則因為0 取正對角矩陣X2=diag(d1,d2,···,dn),并記B=AX2,其中 類似于定理2.1的證明過程,可得|bii|≥αRi(B)+(1?α)Qi(B)(?i∈N),且至少有一個嚴格不等式成立.由A不可約,可得B不可約,則B為不可約α-對角占優(yōu)矩陣,由引理2.3可知B=AX2∈D?,其中X2為正對角矩陣,根據(jù)引理2.1,則A∈D?,因此矩陣A為非奇異H-矩陣. 記 騰訊董事會主席馬化騰憑借328億美元的身家蟬聯(lián)榜單第二名,但他的財富縮水了62億美元。去年的首富、中國恒大董事局主席許家印的排名跌至第三名。許家印的身家為308億美元,下降28%,折合約117億美元,他是今年財富值降低最多的富豪。 引理 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-矩陣. 輸入:已知矩陣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.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不能由文獻[1]中定理2判定.又因為 所以矩陣A不能由文獻[2]中定理1判定.又因為 所以矩陣A不能由文獻[3]中定理1判定.又因為 所以矩陣A不能由文獻[4]中定理1判定.經(jīng)驗證,對α=0.1,0.2,···,1.0,矩陣A不能由文獻[5]中定理1判定. [1]黃廷祝.非奇H矩陣的簡捷判據(jù)[J].計算數(shù)學,1993,15(3):318-328. [2]高中喜,黃廷祝,王廣彬.非奇H-矩陣的充分條件[J].數(shù)學物理學報:A輯,2005,25(3):409-413. [3]黃澤軍,劉建州.非奇異H矩陣的一類新迭代判別法[J].工程數(shù)學學報,2008,25(5):939-942. [4]孫德淑.非奇異H-矩陣的判定準則[J].溫州大學學報,2009,30(3):18-21. [5]尹如軍,徐仲,陸全.非奇H-矩陣的細分迭代判別準則[J].工程數(shù)學學報,2013,30(3):433-441. [6]孫玉祥.廣義對角占優(yōu)矩陣的充分條件[J].高等學校計算數(shù)學學報,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-矩陣的實用判定 [J].純粹數(shù)學與應用數(shù)學, 2013,29(6):601-608. New subdivided and iterative criteria with parameter for nonsingular H-matrices Xiao Lixia,Zhang Junli 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ū)高等學??茖W技術(shù)研究項目(NJZY13159). 肖麗霞(1980-),碩士,講師,研究方向:數(shù)值代數(shù). 2010 MSC:15A573 數(shù)值算法
4 數(shù)值算例
(School of Mathematics,Inner Mongolia University for the Nationalities,Tongliao 028043,China)