利用文[郑巧娟,李耀堂。p-范数双严格对角占优矩阵与新的特征值包含区域。应用数学进展,2017,6(3):367-375。]中所给矩阵的特征值包含区域获得了实对称矩阵正定性的一种判定方法。另外,给出了一个新的正规矩阵Brauer卵形特征值包含区域,使得每个卵形至少包含矩阵的一个特征值。 The eigenvalue inclusion region in [Qiaojuan Zheng, Yaotang Li. p-Norm DSDD Matrices and New Eigenvalue Localization Region, Applied Mathematical Progress, 2017, 6(3): 367-375.] is used to obtain a method for determining the positive definite property of the real symmetric matrix. In addition, a new normal matrix Brauer oval eigenvalue inclusion region is given, so that each oval contains at least one eigenvalue of the matrix.
郑巧娟
云南大学数学与统计学院,云南 昆明
收稿日期:2017年12月18日;录用日期:2018年1月1日;发布日期:2018年1月8日
利用文[郑巧娟,李耀堂。p-范数双严格对角占优矩阵与新的特征值包含区域。应用数学进展,2017,6(3):367-375。]中所给矩阵的特征值包含区域获得了实对称矩阵正定性的一种判定方法。另外,给出了一个新的正规矩阵Brauer卵形特征值包含区域,使得每个卵形至少包含矩阵的一个特征值。
关键词 :p-范数DSDD矩阵,实对称矩阵,正定性,特征值包含区域
Copyright © 2018 by author and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
正定矩阵及其正定矩阵的判定是矩阵理论研究的重要内容之一,很多学者对其做出了研究。对于正定矩阵的研究,起初出现在实二次型与Hermite型的研究中。而这种正定性的研究只局限于实对称矩阵和Hermite矩阵。随着数学以及矩阵应用等的发展,有不少学者开始研究非对称的较为广义的正定矩阵,并取得了丰富的成果,这些成果在自动控制、系统理论等许多领域都有广泛的应用。本文中,我们还将致力于实对称矩阵的正定性研究。利用文献 [
此外,矩阵特征值的包含区域在动力系统的稳定性分析,控制系统的可控制性研究,线性方程组的算法分析等问题的研究中有着重要的作用,是矩阵应用与分析中的一个重要课题。而我们熟悉的特征值包含区域,有著名的Geršgorin定理、Brauer定理等等,但是无论Geršgorin定理还是Brauer定理都不能保证每个Geršgorin圆盘和每个Brauer卵形都包含矩阵的特征值。对于正规矩阵,本文将给出一种新的Brauer卵形,使其每个卵形至少包含矩阵的一个特征值。
本节给出一些基本符号和基本概念,以备后用。
首先给出文中所用的符号和术语。设n为自然数,记 N = { 1 , 2 , ⋯ , n } 。 C n ( R n ) 为 n 维复(实)列向量空间, C n × n ( R n × n ) 为 n × n 维复(实)矩阵组成的集合。设向量 x = ( x 1 , x 2 , ⋯ , x n ) T ∈ C n , A = ( a i j ) ∈ C n × n , q , p ∈ [ 1 , ∞ ] 。记
‖ x ‖ q = ( ∑ i ∈ N | x i | q ) 1 q , r i p ( A ) = ( ∑ j ∈ N \ { i } | a i j | p ) 1 p , i ∈ N 。
特别地,当 p = 1 时, r i ( A ) = ∑ j ∈ N \ { i } | a i j | 。
定义1 [
x T A x > 0 ( x T A x ≥ 0 ) (2.1)
则称 A 是实对称正定矩阵(实对称半正定矩阵)。
文献 [
定义2 [
| a i i | ⋅ | a j j | ≥ r i ( A ) r j ( A ) , i , j ∈ N ,且 i ≠ j (2.2)
则称 A 为双对角占优矩阵,记为DDD矩阵。如果
| a i i | ⋅ | a j j | > r i ( A ) r j ( A ) , i , j ∈ N ,且 i ≠ j
则称 A 为双严格对角占优矩阵,记为DSDD矩阵。
定理1 [
σ ( A ) ⊆ K ( A ) = ∪ i , j ∈ N i ≠ j K i , j ( A ) (2.3)
其中 σ ( A ) 为 A 的谱,
K i , j ( A ) = { z ∈ C : | z − a i i | ⋅ | z − a j j | ≤ r i ( A ) r j ( A ) } ,
称为矩阵 A 的Cassini卵形区域。
在文献 [
定义3 [
使得
x i x j | a i i | | a j j | > r i p ( A ) r j p ( A ) , i , j ∈ N ; i ≠ j (2.4)
成立,则称 A 为p-范数DSDD矩阵,其中 q 为 p 的Hölder补,即 1 p + 1 q = 1 。
定理3 [
定理4 [
∑ i , j ∈ N i ≠ j ( r i p ( A ) r j p ( A ) | a i i | | a j j | ) p p − 1 < 1 。 (2.5)
定理5 [
σ ( A ) ⊆ Φ p , x ( A ) : = ∩ ‖ x ‖ q ≤ 1 ∪ i , j ∈ N Φ i , j p , x ( A ) 。 (2.6)
其中 Φ i , j p , x ( A ) : = { z ∈ C : x i x j | z − a i i | | z − a j j | ≤ r i p ( A ) r j p ( A ) } ; q 为 p 的Hölder补,即 1 p + 1 q = 1 。
定理6:设 A ∈ C n × n 且对任意 i ∈ N ,存在 j ≠ i ,使 a i j ≠ 0 ,且 p ∈ ( 1 , ∞ ] ,则
σ ( A ) ⊆ Φ p ( A )
其中,
Φ p ( A ) : = { z ∈ C : ∑ i , j ∈ N i ≠ j ( r i p ( A ) r j p ( A ) | z − a i i | ⋅ | z − a j j | ) p p − 1 ≥ 1 } (2.7)
证明:设 A ∈ C n × n , p > 1 , λ ∈ σ ( A ) ,假若 λ ∉ Φ p ( A ) ,则
∑ i , j ∈ N i ≠ j ( r i p ( A ) r j p ( A ) | z − a i i | ⋅ | z − a j j | ) p p − 1 < 1 。(2.8)
令 B = λ I − A = ( b i j ) ,由 λ ∈ σ ( A ) 知, 0 ∈ σ ( B ) ,故 B 为奇异矩阵。由 B = λ I − A = ( b i j ) 知
r i p ( B ) = r i p ( A ) , | b i i | = | λ − a i i | , i ∈ N ,
r j p ( B ) = r j p ( A ) , | b j j | = | λ − a j j | , j ∈ N 。
于是由(2.8)式可得
∑ i , j ∈ N i ≠ j ( r i p ( B ) r j p ( B ) | b i i | ⋅ | b j j | ) p p − 1 < 1 。
再由定理4知矩阵 B 为p-范数DSDD矩阵,于是由定理3知矩阵 B 为非奇异的,这与矩阵 B 是奇异的相矛盾。综上所述, λ ∈ Φ p ( A ) 。
下面给出判定实对称矩阵正定的一些必要条件和充分条件。
定理7 [
定理8 [
定理9 [
本节我们给出一些判定实对称矩阵正(半)定性的充分条件。
定理10:设 A ∈ R n × n 为对称矩阵,且 a i i > 0 ( i ∈ N ) 。如果 A 为p-范数DSDD矩阵,则 A 是正定的。
证明:反证法。设 λ ∈ σ ( A ) ,则 λ 为实数。由定理5知,存在满足 ‖ x ‖ q ≤ 1 的正向量 x = ( x 1 , x 2 , ⋯ , x n ) T > 0 和 i 0 , j 0 ∈ N 且 i 0 ≠ j 0 ,使得
λ ∈ Φ i 0 , j 0 p , x ( A )
即
x i 0 x j 0 | λ − a i 0 i 0 | | λ − a j 0 j 0 | ≤ r i 0 p ( A ) r j 0 p ( A ) (3.1)
假若 λ ≤ 0 ,由于 a i i > 0 , i ∈ N 和 A 为p-范数DSDD矩阵得
x i 0 x j 0 | λ − a i 0 i 0 | | λ − a j 0 j 0 | ≥ x i 0 x j 0 a i 0 i 0 a j 0 j 0 > r i 0 p ( A ) r j 0 p ( A )
这与(3.1)式相矛盾,因此 λ > 0 。再由定理8知, A 是正定的。
引理1 [
定理 11:设 A ∈ R n × n 为对称矩阵,且 a i i > 0 ( i ∈ N ) , p ∈ [ 1 , ∞ ] 。如果
∑ i , j ∈ N i ≠ j ( r i p ( A ) r j p ( A ) | a i i | | a j j | ) p p − 1 ≤ 1 (3.2)
则 A 是正半定的。
证明:若 p = 1 ,则
max i , j ∈ N , i ≠ j r i ( A ) r j ( A ) | a i i | | a j j | ≤ 1 ,
由定义2知, A 为双对角占优矩阵。根据引理1可知, A 是正半定的。若 p > 1 ,设 λ 是 A 的特征值。由定理6知 λ ∈ Φ p ( A ) 。这表明对于
∑ i , j ∈ N i ≠ j ( r i p ( A ) r j p ( A ) | λ − a i i | ⋅ | λ − a j j | ) p p − 1 ≥ 1
假若 λ < 0 ,由于 a i i > 0 , i ∈ N ,故有
∑ i , j ∈ N i ≠ j ( r i p ( A ) r j p ( A ) | a i i | ⋅ | a j j | ) p p − 1 > ∑ i , j ∈ N i ≠ j ( r i p ( A ) r j p ( A ) | λ − a i i | ⋅ | λ − a j j | ) p p − 1 ≥ 1
这与(3.2)式相矛盾,因此 λ ≥ 0 。再由定理9可知, A 是正半定的。
本节,我们给出正规矩阵的新型的特征值Brauer卵形包含区域。
定义4 [
A A ∗ = A ∗ A
则称 A 为正规矩阵。
我们知道,即使著名的Geršgorin’s定理也不能保证每个Geršgorin圆盘都包含矩阵的一个特征值。但在文献 [
定理12 [
E i = { z ∈ C : | z − a i i | ≤ r i 2 ( A ) } , (4.1)
必包含矩阵 A 的一个特征值,其中,
r i 2 ( A ) = ( ∑ i , j ∈ N , j ≠ i | a i j | 2 ) 1 2 。
同样,矩阵 A 的某个Brauer Cassini卵形可能也不包含矩阵 A 的特征值。然而对于正规矩阵,应用定理11,我们可得到下面的结论。
定理13:设 A ∈ C n × n 为正规矩阵,则 A 的每个卵形区域
E i , j = { z ∈ C : | z − a i i | | z − a j j | ≤ r i 2 ( A ) r j 2 ( A ) } (4.2)
必包含矩阵 A 的一个特征值。
证明:反证法。假若存在 E i , j 不包含矩阵 A 的任何特征值,则对任意的 λ ∈ σ ( A ) 有
| λ − a i i | | λ − a j j | > r i 2 ( A ) r j 2 ( A )
即
| λ − a i i | | λ − a j j | r i 2 ( A ) r j 2 ( A ) = | λ − a i i | r i 2 ( A ) ⋅ | λ − a j j | r j 2 ( A ) > 1 。
因此, | λ − a i i | r i 2 ( A ) 与 | λ − a j j | r j 2 ( A ) > 1 至少有一个成立,即
| λ − a i i | > r i 2 ( A ) , | λ − a j j | > r j 2 ( A )
至少有一个成立。这表明 A 的第i个和第j个Geršgorin圆盘中至少有一个不包含矩阵 A 的任何特征值。这与定理12的结论矛盾。
例1:设
A = [ 10 3 7 7 3 4 2 2 7 2 6 5 7 2 5 3 ] 。
易知 A 为正规矩阵且 σ ( A ) = { − 1.4474 , 0.7211 , 3.0056 , 20.7208 } 。见图1~图2。
图1. 例1中矩阵A的新的Brauer卵形Ei,j
图2. 例1中矩阵A的Brauer Cassini卵形Ki,j(A)
定理14:设矩阵 A ∈ C n × n 为正规矩阵,则 E i , j ⊆ K i , j 。
证明:由Cauchy-Schwartz不等式知,
r i 2 ( A ) ≤ r i ( A ) ≤ n − 1 r i 2 ( A ) , i ∈ N 。
因此
r i 2 ( A ) r j 2 ( A ) ≤ r i ( A ) r j ( A ) ≤ ( n − 1 ) r i 2 ( A ) r j 2 ( A ) 。
故 r i 2 ( A ) r j 2 ( A ) ≤ r i ( A ) r j ( A ) ,由此即得, E i , j ⊆ K i , j 。
在例1中,对于矩阵
A = [ 10 3 7 7 3 4 2 2 7 2 6 5 7 2 5 3 ] 。
我们从图3中,可以看出 E i , j ( A ) ⊆ K i , j ( A ) 。但是,正规矩阵 A 的某些特征值可能不在 ∪ i , j ∈ N , i ≠ j E i , j 之内。
例2
B = [ 10 7 8 7 7 5 6 5 8 6 10 9 7 5 9 10 ] 。
矩阵 B 为实对称矩阵,则 B 为正规矩阵。用Matlab计算得
σ ( B ) = { 0.0102 , 0.8431 , 3.8581 , 30.288 } ,
图3. 同例1中的公式、图示
图4. 例2中矩阵B的所有新的Brauer卵形的并集
其Brauer型特征值包含区域如图4所示。显然,
30.288 ⊄ ∪ i , j ∈ N , i ≠ j E i , j , N = { 1 , 2 , 3 , 4 } 。
郑巧娟. 正定矩阵的判定方法和新的Brauer卵形 Determination Method of Positive Definite Matrix and New Brauer Oval[J]. 理论数学, 2018, 08(01): 14-21. http://dx.doi.org/10.12677/PM.2018.81003