黎 穩(wěn), 陳艷美, 莫榮華

(1. 華南師范大學數學科學學院, 廣州 510631; 2. 廣東技術師范學院計算機科學學院, 廣州 510665;3. 廣東工業(yè)大學應用數學學院, 廣州 510006)

?

可對角化矩陣的特征值與特征空間的擾動

黎 穩(wěn)1*, 陳艷美1,2*, 莫榮華3

(1. 華南師范大學數學科學學院, 廣州 510631; 2. 廣東技術師范學院計算機科學學院, 廣州 510665;3. 廣東工業(yè)大學應用數學學院, 廣州 510006)

矩陣特征值和特征空間的計算是數值代數的重要課題之一,在科學工程計算等領域有重要的作用. 而特征值與特征空間的擾動分析是有關特征值數值分析的一個重要研究方向,它的經典結果分別是特征值擾動的Hoffman-Wielandt定理和特征空間的sinθ定理. 文中所考慮的是可對角化矩陣的乘法與加法擾動下的特征值與特征空間的組合擾動分析，給出了組合擾動界,所得到的結果推廣了Hermite矩陣的組合擾動的相關結果. 另一方面，從新得到的結果可以分別導出有關特征值和特征空間的擾動界.

可對角化矩陣； 加法擾動界； 乘法擾動界

許多科學與工程計算都涉及矩陣特征值計算問題，而特征值的敏感性分析是特征值的數值分析中的重要方向. 不少學者研究了特征系統的擾動界[1-8],提出了矩陣分解的一些組合擾動界[6-7，9-10].特別地,給出了Hermitian矩陣的特征系統的一些組合擾動界[6-7]. 為了對矩陣的特征系統給出組合擾動界,首先介紹一些符號.

AX=X，

(1)

(2)

(3)

(4)

眾所周知,矩陣的擾動有2種不同的方式：

定義

(5)

(6)

(7)

界(7)的1個漸近形式是

(8)

本文首先給出了加法擾動的組合擾動界,改進了界(7)的結果;然后得到了乘法擾動的組合擾動界,推廣了式(8)的結果.

1 加法擾動界

本節(jié)主要得到特征值與特征空間在加法擾動下的組合擾動界. 首先給出一些引理.

引理1[12]52假定X=(X1,X2)n×n是一個非奇異矩陣,其中X1n×m,那么對于2-或者F-范數‖·‖和任意的列滿秩矩陣n×m,有 ‖‖≤‖‖2‖‖2‖‖,

(9)

其中用M?表示矩陣M的Moore-Penrose逆.

引理2[13]167設TCn×n和i=diag(n×n(i=1,2,3,4),則一定存在集合〈n〉的一個置換,使得‖1T2-3T4‖,

(10)

其中σmin(T)是矩陣T的最小奇異值.

(11)

證明 由式(1)可得

再由式(4)可得

(12)

顯然有

(13)

綜合式(4)和式(5)可得

(14)

再由引理1 和引理2,存在集合〈r〉的轉置,使得

從而證明了式(11).

(15)

其中κ(X)=‖X‖2‖X-1‖2是矩陣X的譜條件數. 不等式(15)可見文獻[14]335和文獻[15]107.

由式(11)可得sinθ型定理[16]. 事實上，有

(16)

2 乘法擾動的情形

引理3[5]475-476已知Ωs×s和Γt×t是2個Hermite矩陣,而且E,Fs×t. 如果(Ω)(Γ)=?,那么ΩX-XΓ=ΩE+FΓ有唯一解Xs×t,而且

引理4[6]245已知=diag(1,…,r)和(1,…,r)滿足，則一定存在集合〈r〉的一個置換,使得

下面給出可對角化矩陣的特征值和特征空間的組合擾動界.

(17)

(18)

從而得到

(19)

和

(20)

把引理3用到式(19)可得

(21)

又由引理1可得

(22)

由式(21)和式(22)可得

(23)

由引理4 和式(20)可得存在〈r〉上的置換,使得

再結合式(23)可以得到

由此證明了定理2.

[1] CHEN X X,LI W. A note on the perturbation bounds of eigenspaces for Hermitian matrices[J]. Journal of Computational and Applied Mathematics,2006,196(1):338-346.

[2] IPSEN I C F. Relative perturbation results for the matrix eigenvalues and singular values[J]. Acta Numerica,1998,7(1):151-201.

[3] IPSEN I C F. A note on unifying absolute and relative perturbation bounds[J]. Linear Algebra and Its Applications,2003,358(2):239-253.

[4] LI R C. Relative perturbation theory:(I) Eigenvalue and singular value variations[J]. SIAM Journal on Matrix Analysis and Applications,1998,19(4):956-982.

[5] LI R C. Relative perturbation theory:(II) Eigenspace and singular subspace variations[J]. SIAM Journal on Matrix Analysis and Applications,1998,20(2):471-492.

[6] LI W. Multiplicative perturbation bounds for spectral and singular value decompositions[J]. Journal of Computational and Applied Mathematics,2008,217(1):243-251.

[7] LI W,SUN W W. Combined perturbation bounds I:Eigensystems and singular value decomposition[J]. SIAM Journal on Matrix Analysis and Applications,2007,29(2):643-655.

[8] DOPICO F M,MORO J,MOLERA J M. Weyl-type relative perturbation bounds for eigensystems of Hermitian matrices[J]. Linear Algebra and Its Applications,2000,309(1):3-18.

[9] CHEN Y M,CHEN X S,LI W. On perturbation bounds for orthogonal projections[J]. Numerical Algorithms,2016,doi:10. 1007/s11075-016-0102-2.

[10]LI W,SUN W W. Combined perturbation bounds:II. Polar decompositions[J]. Science in China:Series A,2007,50(9):1339-1346.

[11]STEWART G,SUN J G. Matrix perturbation theory[M]. Boston:Academic Press,1990.

[12]LI W,CHEN X X. Some residual bounds for approximate eigenvalues and approximate eigenspaces[J]. Journal of Computational Mathematics,2012,30(1):47-58.

[13]ElSNER L,FRIEDLAND S. Singular values,doubly stochastic matrices and applications[J]. Linear Algebra and Its Applications,1995,220(1):161-169.

[14] 孫繼廣. 關于正規(guī)矩陣特征值的擾動[J]. 計算數學,1984(3):334-336.

SUN J G. On the perturbation of the eigenvalues of a normal matrix[J]. Mathematica Numerica Sinica,1984(3):334-336.

[15] 張振躍. 關于非虧損矩陣特征值的擾動[J]. 計算數學,1986(1):106-108.

ZHANG Z Y. On the perturbation of the eigenvalues of a non-defective matrix[J]. Mathematica Numerica Sinica,1986(1):106-108.

[16] DAVIS C,KAHAN W. The rotation of eigenvectors by a perturbation:III[J]. SIAM Journal on Numerical Analysis,1970,7(1):1-46.

【中文責編：莊曉瓊 英文責編：肖菁】

Perturbation Bounds of Eigenvalues and Eigenspaces for Diagonalizable Matrices

LI Wen1*, CHEN Yanmei1,2*, MO Ronghua3

(1. School of Mathematical Sciences, South China Normal University, Guangzhou 510631, China; 2. School of Computer Science, Guangdong Polytechnic Normal University, Guangzhou 510665, China; 3. School of Applied Mathematics, Guangdong University of Technology, Guangzhou 510006, China)

The computation for eigenvalues and eigenspaces is one of the important research fields in numerical algebra, and it plays an important role in scientific Engineering computations. The perturbation analysis for eigenvalues and eigenspaces is one of the significant topics in the numerical analysis in eigenvalue computing. The classical perturbation results are the Hoffman-Wielandt theorem for eigenvalue perturbation and the sinθfor the eigenspace. In this paper, the combined perturbation analysis for the matrix and it’s perturbed matrix being diagonalizable matrices under the additive and the multiplicative perturbation are considered, and their combined perturbation bounds are given, respectively, which extend the corresponding combined perturbation result for the Hermitian matrix case. On the other hand, the eigenvalue perturbation bound and the eigenspace perturbation bound can be derived from the new bounds, respectively.

diagonalizable matrix; additive perturbation bound; multiplicative perturbation bound

2016-05-12 《華南師范大學學報(自然科學版)》網址：http://journal.scnu.edu.cn/n

國家自然科學基金項目(11271144,11571124,11671158,11601340); 中山大學廣東省計算科學重點實驗室開放基金項目(2016016); 廣東技術師范學院青年科研項目

O241.6

A

1000-5463(2016)05-0082-04

*通訊作者:黎穩(wěn),教授,Email:liwen@scnu.edu.cn;陳艷美,講師,Email:chch1980@163.com.