块H-矩阵新的等价表征及谱分析
New Equivalent Characterizations and Spectral Analysis of Block H-Matrices
DOI: 10.12677/aam.2025.149407, PDF, HTML, XML,    科研立项经费支持
作者: 朱开心, 谢智慧*, 黄 琦:湖南科技学院理学院,湖南 永州
关键词: 块H-矩阵块严格对角占优矩阵谱分析Block H-Matrix Block Strictly Diagonal Dominance Matrix Spectral Analysis
摘要: 本文首先给出了块严格 Sα 对角占优矩阵的一类等价条件,从而得到非奇异块H-矩阵新的判定条件。同时,给出了一类非奇异H-矩阵的特征值范围。最后通过例子说明了判定条件的有效性以及对近期结果的改进。
Abstract: This paper first presents a class of equivalent conditions for block strictly Sα diagonally dominant matrices, thereby obtaining new criteria for nonsingular block H-matrices. Meanwhile, it provides the eigenvalue range of a class of nonsingular H-matrices. Finally, examples are given to illustrate the effectiveness of the criteria and the improvement over recent results.
文章引用:朱开心, 谢智慧, 黄琦. 块H-矩阵新的等价表征及谱分析[J]. 应用数学进展, 2025, 14(9): 130-140. https://doi.org/10.12677/aam.2025.149407

1. 引言

在矩阵理论这一充满深度与广度的研究领域中,矩阵分块技术凭借其独特的优势占据着举足轻重的地位,有着极为广泛的应用。它能够将复杂的高维矩阵问题巧妙分解为若干个低维子矩阵问题,极大地简化了矩阵的运算与分析过程,为众多矩阵理论难题的解决提供了关键思路。1962年,Feingold和Varga首次提出了分块矩阵的对角占优性,吸引众多学者对此进行了颇有价值的推广与改进。如李敏、孙玉祥、

孙吉荣、高会双、Ren Y等[1]-[12]诸多学者,分别利用了 A ii 1 1 > ( R i ) α ( C i ) ( 1α ) A ii 1 1 >α R i +( 1α ) C i

两种等价形式,以及相关的诸多块H-矩阵类,从不同角度、不同层面深入探究,研究了一系列的块H-矩阵等价形式的判定条件。这些条件为块H-矩阵的识别与分析提供了有力的依据。不仅如此,部分学者还在此基础上,在等价形式的框架下进一步开展了矩阵特征值包含域的谱分析。

本文在诸多学者研究的基础上,主要结合文献[3]和文献[9]讨论了一类块H-矩阵等价表征的新的判定条件。同时,针对该类块H-矩阵,对其特征值包含域进行了细致且深入的分析,旨在进一步丰富块H-矩阵的理论体系。最后,通过数值算例直观且有效地验证了新判定条件的可行性与有效性,充分展现了本研究的理论价值与实际意义。

2. 记号与相关定理定义

A=( a ij ) M n ( C ) ,且分块如下

A=( A 11 A 12 A 1m A 21 A 22 A 2m A m1 A m2 A mm ) .(1)

其中 A ij M n i ( C ) 且非奇异, i=1 m n i =n iM={ 1,2,,m } 。对 i,jM 给出定义 R i ( A )= ji A ij = R i C i ( A )= ji A ij = C i ,记 σ( A ) 为矩阵 A 的谱, · 为任意的矩阵诱导范数。定义分块矩阵 A 的块比较矩阵 T( A )= ( t ij ) m×m R m×m ,其中

t ij ={ A ii 1 1 ,i=j A ij ,ij i,jM .

故可定义其范数矩阵为

A =( A 11 1 1 A 12 A 1m A 21 A 22 1 1 A 2m A m1 A m2 A mm 1 1 ) .

S M 的任意非空子集, S ¯ M S 的补集,即 S, S ¯ M S S ¯ = S S ¯ =M

iS ,设 R i S ( A )= tS,ti A it = R i S R i S ¯ ( A )= t S ¯ A it = R i S ¯ C i S ( A )= tS,ti A ti = C i S C i S ¯ ( A )= t S ¯ A ti = C i S ¯ ,故有 R i = R i S + R i S ¯ C i = C i S + C i S ¯ 。另设 r i S ( A )= tS,ti | a it | r i S ¯ ( A )= t S ¯ | a it | c i S ( A )= tS,ti | a ti | c i S ¯ ( A )= t S ¯ | a ti | ,同样有 r i ( A )= r i S ( A )+ r i S ¯ ( A ) c i ( A )= c i S ( A )+ c i S ¯ ( A )

对矩阵指标集进行划分,具体如下:

M 1 ={ i| R i S < A ii 1 1 R i S ¯ < C i S } , M 2 ={ i| C i S < A ii 1 1 R i S ¯ < R i S } ,

M 3 ={ i| A ii 1 1 R i S ¯ C i S > R i S } , M 4 ={ i| A ii 1 1 R i S ¯ R i S > C i S } ,

M 5 ={ i| A ii 1 1 R i S ¯ > C i S = R i S } , M 6 ={ i| A ii 1 1 R i S ¯ R i S , A ii 1 1 R i S ¯ C i S } ,

显然, M= M 1 M 2 M 3 M 4 M 5 M 6

定义1 [2] 存在 α[ 0,1 ] kN ,使得

| a ii |> ( r i S ) α ( c i S ) 1α + r i S ¯ ,SN,| S |=k ,(2)

则称 A 为严格 Sα 对角占优矩阵。

定义2 [5] A=( a ij ) M n ( C ) 且分块如式(1),若存在 α[ 0,1 ] ,使得

A ii 1 1 > ( R i ) α ( C i ) ( 1α ) iM ,

则称 A 为严格 α 2 型块对角占优矩阵,记为 A α 2 -BSD

定义3 [6] A=( a ij ) M n ( C ) 且分块如式(1),若存在 α[ 0,1 ] ,使得

A ii 1 1 >α R i +( 1α ) C i iM ,

则称 A 为严格 α 1 型块对角占优矩阵,记为 A α 1 -BSD

引理1 [6] A=( a ij ) M n ( C ) 且分块如式(1),若矩阵 A 满足 A α 1 -BSD A α 2 -BSD ,则 A 是一个块H-矩阵。

引理2 [7] τ,σ 为两个非负实数,则对 α[ 0,1 ]

ατ+( 1α )σ τ α σ ( 1α ) ,

当且仅当 α=0 α=1 时,有 τ=σ

引理3 [4] 假设 f( t )=at+b( 1t ) ,对 t( 0,1 )

a>b>0 ,则 f( t ) 是一个单调递增的函数。

b>a>0 ,则 f( t ) 是一个单调递减的函数。

引理4 [4] ε>0( ε=o( x ),x0 ) ,定义

E={ zC| | za |b+ε } , F={ zC| | za |b } ,

E=F

引理5 [8] 设矩阵 B=( b ij ) C m×m ,其中

b ij ={ A ii 1 1 ,i=j A ij ,ij ,i,jM ,

若矩阵 B 满足(2)式为严格 Sα 对角占优矩阵,则有

[ A ii 1 1 R i S ¯ ]> ( R i S ) α ( C i S ) 1α , SN,| S |=k ,

成立,并称矩阵 A 为块严格 Sα 对角占优矩阵。

证明 由定义1可知,若矩阵 B 满足(2)式为严格 Sα 对角占优矩阵,则存在正数 x 1 , x 2 ,, x m >0 ,使得

x i | b ii |> x j ji | b ij | ,i,jM .

作正对角矩阵 X=diag{ x 1 , x 2 ,, x m } ,即可得 BX 为严格对角占优矩阵。于是作 X 1 =diag{ x 1 I n 1 , x 2 I n 2 ,, x m I n m } ,则 A X 1 为块严格对角占优矩阵,故有引理5结论成立。

引理6 [3] A=( a ij ) M n ( C ) 且分块如式(1),若 M 6 = ,对 i M 1 j M 2

max j M 2 R j A jj 1 1 R j C j < min i M 1 A ii 1 1 R i C i R i ,

A 为非奇异块H-矩阵。

引理7 [3] A=( a ij ) M n ( C ) 且分块如式(1),若矩阵 A α 1 -BSD ,那么

σ( A )G( A )= i=1 m σ( A ii ) i M 1 M 3 G i j M 2 M 4 M 5 G j ,

其中

G M = iM { zC| ( A ii z I n i ) 1 1 R i S ¯ ( A ) } ,

G i ={ zC| ( A ii z I n i ) 1 1 ( min t M 2 A tt 1 1 R t C t R t ) R i +( max t M 2 C t A tt 1 1 C t R t ) C i } , i M 1 M 3 ;

G j ={ zC| ( A jj z I n j ) 1 1 ( max v M 1 R v A vv 1 1 R v C v ) R j +( min v M 1 A vv 1 1 C v R v C v ) C j } , j M 2 M 4 M 5 .

引理6和引理7为文献[3]中的主要结果,本文在此基础上利用不等式的放缩技巧,扩大的块H-矩阵的判定范围,给出一类新的等价表征,并对其进行谱分析。从而对文献中的结果进行改进。

3. 主要结果

3.1. 定理1

A=( a ij ) M n ( C ) 且分块如式(1),若存在 S, S ¯ M S S ¯ = S S ¯ =M 。对 iS ,有 A ii 1 1 > R i S ,对 i M 1 j M 2

max j M 2 R j S [ A jj 1 1 R j S ¯ ] R j S C j S < min i M 1 [ A ii 1 1 R i S ¯ ] R i S C i S R i S , (3)

当且仅当 M 6 = ,则 A 为非奇异块H-矩阵。

证明 充分性,对 i M 1 j M 2

R i S < A ii 1 1 R i S ¯ < C i S C j S < A jj 1 1 R j S ¯ < R j S ,

则有 0< [ A ii 1 1 R i S ¯ ] R i S C i S R i S <1 0< R j S [ A jj 1 1 R j S ¯ ] R j S C j S <1

由(3)式和引理2可知,存在 α( 0,1 ) ,使得

max j M 2 R j S [ A jj 1 1 R j S ¯ ] R j S C j S <1α< min i M 1 [ A ii 1 1 R i S ¯ ] R i S C i S R i S , (4)

故而可得

[ A jj 1 1 R j S ¯ ]>α R j S +( 1α ) C j S ( R j S ) α ( C j S ) 1α ,

[ A ii 1 1 R i S ¯ ]>α R i S +( 1α ) C i S ( R i S ) α ( C i S ) 1α .

i M 3 M 4 M 5 以及 α( 0,1 ) ,显然有

[ A ii 1 1 R i S ¯ ]>α R i S +( 1α ) C i S ( R i S ) α ( C i S ) 1α ,

由条件 M 6 = ,对 i M 1 M 2 M 3 M 4 M 5 M 6 α( 0,1 ) ,使得

[ A ii 1 1 R i S ¯ ]>α R i S +( 1α ) C i S ( R i S ) α ( C i S ) 1α ,

A 为非奇异块H-矩阵。

必要性,因为 A 为非奇异块H-矩阵,显然有 M 6 = ,且对 i M 1

[ A ii 1 1 R i S ¯ ]>α R i S +( 1α ) C i S ( R i S ) α ( C i S ) 1α ,

1α< [ A ii 1 1 R i S ¯ ] R i S C i S R i S . (5)

而对 j M 2 ,有 [ A jj 1 1 R j S ¯ ]>α R j S +( 1α ) C j S ( R j S ) α ( C j S ) 1α ,故有

R j S [ A jj 1 1 R j S ¯ ] R j S C j S <1α . (6)

由(5)式和(6)式可知对 i M 1 j M 2

max j M 2 R j S [ A jj 1 1 R j S ¯ ] R j S C j S < min i M 1 [ A ii 1 1 R i S ¯ ] R i S C i S R i S .

综上,矩阵 A 为非奇异块H-矩阵。该定理对文献中的引理6进行改进,相较于文献中的结果,此定理的判定范围更广。最后由数值算例2进行相关证明。

3.2. 定理2

A=( a ij ) M n ( C ) 且分块如式(1),若矩阵 A 为非奇异块H-矩阵,那么

σ( A )G( A )= i=1 m σ( A ii ) G M i M 1 M 3 G i j M 2 M 4 M 5 G j ,

其中

G M = iM { zC| ( A ii z I n i ) 1 1 R i S ¯ ( A ) } ,

G i ={ zC| [ ( A ii z I n i ) 1 1 R i S ¯ ( A ) ]( min t M 2 [ A tt 1 1 R t S ¯ ] R t S C t S R t S ) R i S +( max t M 2 C t S [ A tt 1 1 R t S ¯ ] C t S R t S ) C i S }

i M 1 M 3 ;

G j ={ zC| [ ( A jj z I n j ) 1 1 R j S ¯ ( A ) ]( max v M 1 R v S [ A vv 1 1 R v S ¯ ] R v S C v S ) R j S +( min v M 1 [ A vv 1 1 R v S ¯ ] C v S R v S C v S ) C j S }

j M 2 M 4 M 5 .

证明 因 [ A ii 1 1 R i S ¯ ]>α R i S +( 1α ) C i S αK ,其中

K=( max j M 2 R j S [ A jj 1 1 R j S ¯ ] R j S C j S , min i M 1 [ A ii 1 1 R i S ¯ ] R i S C i S R i S )( 0,1 ) .

λσ( A ii ) ,则 iM ,有 [ ( A ii z I n i ) 1 1 R i S ¯ ( A ) ]α R i S +( 1α ) C i S αK

接着根据引理3考虑以下五种情况:

i M 1 ,则 0< R i S < C i S ,因为 f( t )= R i S t+ C i S ( 1t ) 是一个单调递减的函数,故对 ε>0 ,令 ε 1 = ε C i S R i S ,可以得到 ε 1 >0

[ ( A ii z I n i ) 1 1 R i S ¯ ( A ) ]( min t M 2 [ A tt 1 1 R t S ¯ ] R t S C t S R t S ε 1 ) R i S +( 1( min t M 2 [ A tt 1 1 R t S ¯ ] R t S C t S R t S ε 1 ) ) C i S =( min t M 2 [ A tt 1 1 R t S ¯ ] R t S C t S R t S ε 1 ) R i S +( 1( 1 max t M 2 C t S [ A tt 1 1 R t S ¯ ] C t S R t S )+ ε 1 ) C i S =( min t M 2 [ A tt 1 1 R t S ¯ ] R t S C t S R t S ) R i S +( max t M 2 C t S [ A tt 1 1 R t S ¯ ] C t S R t S ) C i S + ε 1 ( C i S R i S ) =( min t M 2 [ A tt 1 1 R t S ¯ ] R t S C t S R t S ) R i S +( max t M 2 C t S [ A tt 1 1 R t S ¯ ] C t S R t S ) C i S +ε.

因此,我们可以得到

λ{ zC| [ ( A ii z I n i ) 1 1 R i S ¯ ( A ) ]( min t M 2 [ A tt 1 1 R t S ¯ ] R t S C t S R t S ) R i S +( max t M 2 C t S [ A tt 1 1 R t S ¯ ] C t S R t S ) C i S +ε }

λ G i , i M 1 .

2) 若 j M 2 ,则 0< C j S < R j S ,因为 f( t )= R j S t+ C j S ( 1t ) 是一个单调递增的函数,故对 ε>0 ,令 ε 2 = ε R j S C j S ,可以得到 ε 2 >0

[ ( A jj z I n j ) 1 1 R j S ¯ ( A ) ]( max v M 1 R v S [ A vv 1 1 R v S ¯ ] R v S C v S + ε 2 ) R j S +( 1( max v M 1 R v S [ A vv 1 1 R v S ¯ ] R v S C v S + ε 2 ) ) C j S =( max v M 1 R v S [ A vv 1 1 R v S ¯ ] R v S C v S + ε 2 ) R j S +( 1( 1 min v M 1 [ A vv 1 1 R v S ¯ ] C v S R v S C v S ) ε 2 ) C j S =( max v M 1 R v S [ A vv 1 1 R v S ¯ ] R v S C v S ) R i S +( min v M 1 [ A vv 1 1 R v S ¯ ] C v S R v S C v S ) C j S + ε 2 ( R j S C j S ) =( max v M 1 R v S [ A vv 1 1 R v S ¯ ] R v S C v S ) R i S +( min v M 1 [ A vv 1 1 R v S ¯ ] C v S R v S C v S ) C j S +ε.

因此,我们可以得到

λ{ zC| [ ( A jj z I n j ) 1 1 R j S ¯ ( A ) ]( max v M 1 R v S [ A vv 1 1 R v S ¯ ] R v S C v S ) R j S +( min v M 1 [ A vv 1 1 R v S ¯ ] C v S R v S C v S ) C j S +ε }

λ G j , j M 2 .

3) 若 i M 3 ,则 0< R i S < C i S ,因为 f( t )= R i S t+ C i S ( 1t ) 是一个单调递减的函数,这与1)类似,故可得 λ G i i M 3

4) 若 j M 4 ,则 0< C j S < R j S ,因为 f( t )= R j S t+ C j S ( 1t ) 是一个单调递增的函数,这与2)类似,故可得 λ G j j M 4

5) 若 j M 5 ,则 0< C j S = R j S <[ ( A jj z I n j ) 1 1 R j S ¯ ( A ) ] ,故有

[ A jj 1 1 R j S ¯ ]>α R j S +( 1α ) C j S , α( 0,1 ) ,

若令这 α= max v M 1 R v S [ A vv 1 1 R v S ¯ ] R v S C v S ,则有

[ ( A jj z I n j ) 1 1 R j S ¯ ( A ) ]( max v M 1 R v S [ A vv 1 1 R v S ¯ ] R v S C v S ) R j S +( min v M 1 [ A vv 1 1 R v S ¯ ] C v S R v S C v S ) C j S ,

即可得出 λ G j j M 5

综上所述,矩阵 A 为非奇异块H-矩阵且有 M= M 1 M 2 M 3 M 4 M 5 时,有

σ( A )G( A )= i=1 m σ( A ii ) iM G M i M 1 M 3 G i j M 2 M 4 M 5 G j .

4. 论数值算例

4.1. 算例1

考虑矩阵 A

A=( 8 1 1 3 1 1 1 0 2 13 2 5 1 2 0 1 3 1 13 8 1 1 1 3 0 1 7 31 1 1 2 1 1 4 1 1 14 2 1 2 1 1 1 1 1 35 1 1 1 0 1.5 0.5 0 0.5 7 2 0 1 0 1 0.5 0 3 11 ) .

将矩阵  A 作如下分块,并得其范数矩阵

A=( A 11 A 12 A 13 A 14 A 21 A 22 A 23 A 24 A 31 A 32 A 33 A 34 A 41 A 42 A 43 A 44 ) , A =( A 11 1 1 A 12 A 13 A 14 A 21 A 22 1 1 A 23 A 24 A 31 A 32 A 33 1 1 A 34 A 41 A 42 A 43 A 44 1 1 ) .

即将矩阵 A 划分成每个子块均为 2×2 阶的分块矩阵。令 S={ 1,2 } ,则 S ¯ ={ 3,4 } 取范数为矩阵的2-范数。

通过计算可以得

A 11 1 1 =7.5964 , A 12 =6.2429 , A 13 =2.3028 , A 14 =1 ;

A 21 =3.1796 , A 22 1 1 =10.2903 , A 23 =1.4142 , A 24 =3.1926 ;

A 31 =4.1926 , A 32 =1.4142 , A 33 1 1 =13.8985 , A 34 =2.3028 ;

A 41 =1 , A 42 =1.6283 , A 43 =0.5 , A 44 1 1 =5.8123 .

对于定理1有以下结果:

第一行: R 1 S =6.2429>[ A 11 1 1 R 1 S ¯ ]=7.79642.30281=4.2936> C 1 S =3.1796

第二行: R 2 S =3.1796<[ A 22 1 1 R 2 S ¯ ]=10.29031.41423.1926=5.6835< C 2 S =6.2429

第三行: C 3 S ¯ =0.5< R 3 S ¯ =2.3028<[ A 33 1 1 R 3 S ]=13.89854.19261.4142=8.2917

第四行: R 4 S ¯ =0.5< C 4 S ¯ =2.3028<[ A 44 1 1 R 4 S ]=5.812311.6283=3.184

综上, M 1 ={ 2 }, M 2 ={ 1 }, M 3 ={ 4 }, M 4 ={ 3 }, M 5 = M 6 = 。将以上数据代入(4)式计算可得

0.6363 R 1 S [ A 11 1 1 R 1 S ¯ ] R 1 S C 1 S < [ A 22 1 1 R 2 S ¯ ] R 2 S C 2 S R 2 S 0.8174 .

根据定理1可知矩阵 A 为非奇异块H-矩阵。

4.2. 算例2

考虑矩阵 A

A=( 9 1 1 3 1 1 1 0 2 13 2 5 1 2 0 1 3 1 13 8 1 1 1 3 0 1 7 31 1 1 2 1 1 4 1 1 14 2 1 2 1 1 1 1 1 35 1 1 1 0 1.5 0.5 0 0.5 7 2 0 1 0 1 0.5 0 3 11 ) .

将矩阵 A 作如下分块,并得其范数矩阵

A=( A 11 A 12 A 13 A 14 A 21 A 22 A 23 A 24 A 31 A 32 A 33 A 34 A 41 A 42 A 43 A 44 ) , A =( A 11 1 1 A 12 A 13 A 14 A 21 A 22 1 1 A 23 A 24 A 31 A 32 A 33 1 1 A 34 A 41 A 42 A 43 A 44 1 1 ) .

即将矩阵 A 划分成每个子块均为 2×2 阶的分块矩阵。令 S={ 1,2 } ,则 S ¯ ={ 3,4 } 取范数为矩阵的2-范数。

通过计算可以得

A 11 1 1 =8.5114 , A 12 =6.2429 , A 13 =2.3028 , A 14 =1 ;

A 21 =3.1796 , A 22 1 1 =10.2903 , A 23 =1.4142 , A 24 =3.1926 ;

A 31 =4.1926 , A 32 =1.4142 , A 33 1 1 =13.8985 , A 34 =2.3028 ;

A 41 =1 , A 42 =1.6283 , A 43 =0.5 , A 44 1 1 =5.8123 .

对于引理2有:

第一行: R 1 =9.5457> A 11 1 1 =8.5114> C 1 =8.3722

第二行: R 2 =7.7864< C 2 =9.2854< A 22 1 1 =10.2903

第三行: C 3 =4.217< R 3 =7.9096< A 33 1 1 =13.8985

第四行: R 4 =3.1283< A 44 1 1 =5.8123< C 4 =6.2206

综上, M 1 ={ 4 }, M 2 ={ 1 }, M 3 ={ 2 }, M 4 ={ 3 }, M 5 = M 6 = 。将以上数据代入(4)式计算可得

0.8814 R 1 A 11 1 1 R 1 C 1 > A 44 1 1 R 4 C 4 R 4 0.8680 .

这不满足引理2的条件。

对于本文定理1有以下结果:

第一行: R 1 S =6.2429>[ A 11 1 1 R 1 S ¯ ]=8.51142.30281=5.2086> C 1 S =3.1796

第二行: R 2 S =3.1796<[ A 22 1 1 R 2 S ¯ ]=10.29031.41423.1926=5.6835< C 2 S =6.2429

第三行: C 3 S ¯ =0.5< R 3 S ¯ =2.3028<[ A 33 1 1 R 3 S ]=13.89854.19261.4142=8.2917

第四行: R 4 S ¯ =0.5< C 4 S ¯ =2.3028<[ A 44 1 1 R 4 S ]=5.812311.6283=3.184

综上, M 1 ={ 2 }, M 2 ={ 1 }, M 3 ={ 4 }, M 4 ={ 3 }, M 5 = M 6 = 。将以上数据代入(4)式计算可得

0.3376 R 1 S [ A 11 1 1 R 1 S ¯ ] R 1 S C 1 S < [ A 22 1 1 R 2 S ¯ ] R 2 S C 2 S R 2 S 0.8174 .

根据定理1可知矩阵 A 为非奇异块H-矩阵。

5. 总结

本文给出了块严格 Sα 对角占优矩阵的一类等价表征并理论证明其为非奇异块H-矩阵,描述了其特征值包含域。数值算例1中可以看出定理1是成立的。数值算例2可以看出其判定范围相比于文献[3]中的结果是更加广泛的。因此,本文结果块H-矩阵中可判定范围更加广泛,特征值包含域更加精确。

致 谢

感谢导师、师姐和同门对本项目的悉心指导和帮助,在此表示由衷的感谢!

基金项目

湖南省科研创新一般项目《广义严格对角占优矩阵的高效迭代算法研究(23C0342)》。

NOTES

*通讯作者。

参考文献

[1] Feingold, D.G. and Varga, R. (1962) Block Diagonally Dominant Matrices and Generalizations of the Gerschgorin Circle Theorem. Pacific Journal of Mathematics, 12, 1241-1250.
https://doi.org/10.2140/pjm.1962.12.1241
[2] Cvetković, L. (2006) H-Matrix Theory vs. Eigenvalue Localization. Numerical Algorithms, 42, 229-245.
https://doi.org/10.1007/s11075-006-9029-3
[3] Gao, H., Han, G., Sun, Y., Sun, F. and Ren, Y. (2020) Block H-Matrices and Spectrum of Block Matrices. Journal of Physics: Conference Series, 1575, Article 012114.
https://doi.org/10.1088/1742-6596/1575/1/012114
[4] 李敏, 孙玉翔. H-矩阵的实用判定及谱分析[J]. 高等学校计算数学学报, 2007(2): 117-125.
[5] 李庆春, 刘磊, 等. 矩阵对角占优性的推广[J]. 吉林师范学院学报, 1996(5): 4-7.
[6] 孙玉祥. 广义对角占优矩阵的充分条件[J]. 高等学校计算数学学报, 1997(3): 216-223.
[7] Xu, C. and Xu, Z. (1991) Matrix Analysis. Northwestern Polytechnical University Press.
[8] 黄廷祝, 蒲和平. (块)H-矩阵与亚正定矩阵[J]. 电子科技大学学报, 1998(2): 106-108.
[9] 李敏, 孙玉翔. 块α-对角占优矩阵的等价表征及应用[J]. 北华大学学报, 2008, 9(6): 485-489.
[10] 孙吉荣, 吴建勇, 王焕东. 块α-对角占优矩阵的等价表征及应用[J]. 武汉理工大学学报, 2009, 31(16): 200-204.
[11] 高会双, 韩贵春. 块α-对角占优矩阵的充要条件[J]. 湖北民族学院学报, 2017, 35(2): 131-133.
[12] Ren, Y. and Gao, H. (2022) Criteria for Block H-Matrices. 2022 8th Annual International Conference on Network and Information Systems for Computers (ICNISC), Hangzhou, 16-19 September 2022, 960-962.
https://doi.org/10.1109/icnisc57059.2022.00193