一類反對稱特征值問題的向后誤差

龐澄澄,王衛(wèi)國

(中國海洋大學數學科學學院,山東青島266100)

一類反對稱特征值問題的向后誤差

龐澄澄,王衛(wèi)國

(中國海洋大學數學科學學院,山東青島266100)

本文研究一類反對稱特征值問題在實擾動情形下的向后誤差。該問題的系數矩陣同時具有反對稱性和零塊結構。本文將給出范數型雙結構向后誤差的表達式,并將所得結果與一般的反對稱結構向后誤差進行比較。數值實驗結果表明,雙結構向后誤差和反對稱結構向后誤差的比值有時會非常大。

結構特征值問題;雙結構;向后誤差

0 引言

考慮特征值問題

在研究Stiefel流形上的測地線計算[1]時,導出上述問題。

Higham和Higham[2]指出,向后誤差及條件數在現代數值代數中是2個重要的概念。前者反映數值方法的向后穩(wěn)定性;后者揭示問題的解關于數據擾動的敏感性。在一階近似下,二者結合可給出計算解的誤差估計:

計算解誤差?向后誤差×條件數。

對于問題(1),系數矩陣具有反對稱性及零塊結構。如文獻[2]所指出的,對于結構問題,應該考慮結構向后誤差與結構條件數,一階近似下,計算解的誤差估計為:

計算解誤差?結構向后誤差×結構條件數。

關于結構向后誤差分析,已有很多研究。Sun研究了KKT方程組的結構向后誤差[3]。Tissuer歸納了包括反對稱結構在內的約90種結構特征值問題的向后誤差分析的結果,其中一部分是作為待解決問題提出的[4]。對于系數矩陣為Hermite陣的特征值問題的結構向后誤差,Higham對廣義特征值問題進行了相關研究[2],Tisseur考慮了多項式特征值問題[5], Hachstenbach和Plestenjak研究了多參數特征值問題[6]。劉新國等在近些年研究了多種特征值問題的向后誤差,其中包括Hermite陣,酉陣,一類周期辛矩陣等[7-10]。有關特征值問題(1)的向后誤差還未見到研究。

本文主要研究該特征值問題的結構向后誤差。對于實矩陣,求解相應特征值問題的時候一般盡量規(guī)避復運算(比如雙位移QR方法),因此,本文主要研究實擾動情形。

1 預備

1.1 一般反對稱結構向后誤差

問題(1)可簡寫為

1.2 四個引理

在本文以下推導過程中,將用到以下引理。

此結果給出了矩陣方程XB=C的解X的明顯表達式。

引理1[11]設B∈Rn×p,C∈Rm×p已給定。定義集合

這里Z∈Rm×n,B+表示B的廣義逆。

下述結果給出了方程HB=C的反對稱解H的明顯表達式。

引理2[4]設B,C∈Rn×l已給定。定義集合

即,實部x(1)和虛部x(2)正交,并且等長。

2 雙結構向后誤差

這里

3 數值實驗

所有數值實驗的計算都是使用Matlab7.0運行的,相對機器精度約為2.2202×10-16。將(2)定義的一般反對稱結構向后誤差η()簡記為η,將(4)定義的雙結構向后誤差η)簡記為ηD。

取m=3,n=4,隨機生成H和A,要求HT=-H∈R3×3,A∈R4×3。使用order[V,D]=eig()得到近似特征值和特征向量。選擇第一個特征值λ,以及相應的特征向量z ,來計算ηD/η的值。

將100次實驗得到的結果用圖1表示。可以看出,大部分ηD/η的值比較小,在[1.4,10]之間擾動,只有很少幾個ηD/η的值相差很大。例如,取

可得到ηD/η=239.983 3。

圖1 計算100次的結果Fig.1 The results of computing one hundred times

圖2 計算10 000次的結果Fig.2 The results of computing ten thousand times

將10 000次實驗得到的結果用圖2表示。可以看出,大部分ηD/η的值在[1.4,300]之間,只有幾個ηD/η的值會相差比較大。例如,取

可得到ηD/η=3.325 7×103。但這些少量的ηD/η的值表明,雙結構向后誤差在某些情形比單結構向后誤差可以大很多。

[1] Edelman A,Arias T A,Smith S T.The geometry of algorithms with orthogonally constraints[J].SIAM J Matrix anal Appl, 1999,20(2):303-353.

[2] Higham D J,Higham N J.Structured backward error and condition of generalized eigenvalue problems[J].SIAM J Matrix Anal Appl,1998,20(2):493-512.

[3] Sun J G.Structured backward errors for KKT systems[J].Linear Alg Appl,1999,288:75-88.

[4] Tisseur F.A chart of backward errors for singly and doubly structured eigenvalue problems[J].SIAM J Matrix Anal Appl,2003, 24(3):877-897.

[5] Tisseur F.Backward error and condition of polynomial eigenvalue problems[J].Linear Alg Appl,2000,309(1-3):339-361.

[6] Hochstenback M E,Plestenjak B.Backward error,condition numbers,and pseudospectra for the multiparameter eigenvalue problem [J].Linear Alg Appl,2003,375:63-81.

[7] Xin-Guo Liu,Ze-Xi Wang.A note on the backward errors for Hermite eigenvalue problems[J].Appl Math Comput,2005,165 (2):405-417.

[8] 劉新國,欒世寶.一類周期辛矩陣對特征值問題的向后誤差分析[J].山東大學學報:工學版,2007,37(1):94-98.

[9] 司玉琨.結構酉陣特征值問題向后擾動分析[D].青島:中國海洋大學,2005.

[10] 王澤璽.兩類系統的向后誤差及條件數[D].青島:中國海洋大學,2005.

[11] 孫繼廣.矩陣擾動分析第二版[M].北京:科學出版社,2001.

Abstract: This paper deals with backward errors for a class of skew-symmetric eigenvalue problems.For the doubly structure,both skew-symmetric structure and zeroth structure are considered.We derive explicit formula for the doubly structured normwise backward error,and compare the results with the single structured case.Numerical results demonstrate that the ratio of doubly structured backward error and the single structured one can be large.

Key words: structured eigenvalue problem;doubly structure;backward error

AMS Subject Classification: 15A18

責任編輯 朱寶象

Backward Errors for a Class of Skew-Symmetric Eigenvalue Problems

PANG Cheng-Cheng,WANG Wei-Guo
(School of Mathematics Science,Ocean University of China,Qingdao 266100,China)

O241.1

A

1672-5174(2011)06-130-05

國家自然科學基金資助項目(10971204);山東省自然科學基金資助項目(Y2008A07);中央高校基本科研業(yè)務費專項(201013048)資助

2010-11-20;

2011-03-24

龐澄澄(1985-),女,碩士生。E-mail:pcc192295@163.com