当前位置:文档之家› 严格对角占优矩阵与SOR迭代法的收敛性定理

严格对角占优矩阵与SOR迭代法的收敛性定理

第34卷第1期 2011年3月 长春理工大学学报(自然科学版) JournalofChangchunUniversity ofScience andTechnology(Natural ScienceEdition) VO}.34 No.1 Mar.201l 

严格对角占优矩阵与SOR迭代法的收敛性定理 

宋岱才,敬长红,陈德艳 

(辽宁石油化工大学理学院,抚顺113001) 

摘要:针对线性方程组的系数矩阵为 链严格对角占优矩阵和双严格对角占优矩阵的情况,讨论了线性方程组求 

解时常用的SOR迭代方法的收敛性,给出了迭代法收敛性定理,解决了以往估计迭代矩阵谱半径的问题。结果不 

仅适用于这两类矩阵.还适用于广义严格对角占优矩阵类,最后举例说明了所给结果的优越性。 关键词:a.链严格对角占优矩阵;双严格对角占优矩阵;迭代法;收敛性 

中图分类号:O241.6;O151.2 文献标识码:A 文章编号:1672—9870(2011)01-0170-03 

Diagonal Strictly Dominance Matrix and Convergence 

Theorem 0f SoR Iteration Method 

SONG Daicai,JING Changhong,CHEN Deyan 

(School ofSciences,LiaoningUniversityofPetroleum&ChemicalTechnology,Fushun 113001) 

Abstract:InthispaperConvergencetheoremofSORiterationme.odOrsolvingl earsystemis studied.whencoefficientma- 

trix is 6c—chain diagonal strictly dominance or doubly diagonal strictly dominance,and some convergence theorems are given. 

which solves the problem ofspectral radius ofiterative matrices.Results obtained are applicable for。[一chain diagonal strictly domi— 

nancematrixordoublydiagonal strictlydominancematrix,andimprovetheknown resultsandareapplicableforgeneralizeddi— 

agonal strictly dominance matrices.Finally,a numerical example is given for illustrating advantage ofthe results in this paper. Key words:口一chain diagonal strictly dominance matrix;doubly diagonal strictly dominance matrix;iteration method; 

convergence theorem 

1基本概念及引理 

给定线性方程组Ax=b,其中A∈C 为非奇异 

矩阵,b为,2维列向量。文献[1.7]对于迭代矩阵 

为严格对角占优矩阵、 .链严格对角占优矩阵和双 

“.严格对角占优矩阵等情形分别讨论了常用的几种 

迭代法的谱半径的上界估计问题。文献[8-9]针对 

系数矩阵为以上对角占优的情况讨论了Jacobi迭代 

法、Gauss.Seidel迭代法以及JOR迭代法的收敛性 

问题。本文的主要工作是:针对方程组的系数矩阵 

4为 一链严格对角占优矩阵以及双严格对角占优矩 

阵,讨论了常用到的SOR迭代法的收敛性,得到 

几个重要的结论,解决了估计迭代矩阵谱半径的 界限问题。最后举例说明这一结果的适用性。 

设方程组的系数矩阵A分解为A=D—L—U,其 

中D=diag(alj'a2:,…, ),一 是矩阵A的严格 

下三角矩阵,一魄矩 的严格上三角矩阵。SOR 

迭代法的公式为: 

“’= +厂,k=0,1,2… (1) 

其中L。=(D-coL) [(1一∞)D+∞ 称为SOR迭代法的 

迭代矩阵,厂=(D-coL)~b, 称为松弛因子。SOR 

迭代法收敛的必要条件为O<w<2。我们在文献 

[8-9]基础上讨论SOR迭代法的收敛性。 

设 =( )∈C ,记R ):∑ I. (,4)=∑{ l, l≠|l≠) i∈Ⅳ {1,2,…, )。 

定义i¨叫 设 =( 1∈c” ,若对任意的i∈ 

收稿日期:2010—09—26 基金项目:国家自然科学基金项目(20273028);辽宁省教育厅高校科研项目(2004F100) 作者简介:宋岱才(1954一),男,教授,主要从事数值代数的研究,E-mail:sdc1@163.corn。

 第1期 宋岱才,等:严格对角占优矩阵与SOR迭代法的收敛性定理 l7l 

Ⅳ,存在瑾∈[0,1],皆有la l> ) ),则称 为 

6c.链严格对角占优矩阵,记为A∈D 。若存在正对 

角矩阵d=diag , ,…, )使得AdE D ,则称 

为广义口一链严格对角占优矩阵,记为A∈GD 。 

定义2¨ 设 ∈C ,若 ∽,V 

i。jEN成立,则称 为双严格对角占优矩阵,记为 

A∈D。若存在正对角矩阵d=diag , ,…, ), 

使得 ∈D,则称 为广义双严格对角占优矩阵, 

记为A∈GD。 

引理1[1叫 设 ∈C ,若A∈Do或A∈ 

G ,则A为非奇异矩阵。 

引理2¨。。 设 ==( ∈C ,若A∈D或A∈ 

GD,则A为非奇异矩阵。 

引理3设 是一个常数,0<09<1,则当I ̄1--1 

时,总有12-1+091>__12I.o9>09。 

证明首先证明不等式的左端成立。当4=1 

时,不等式显然成立。 

当2>1时, 一1>0,又因为 1,所以有 

一1 ∞ -1),即得到 一1+09>, ̄zo,即得到 一1+091>-l 

l。 

当 一1时,有1-4>0,同理由0< 1得,1一 

∞(1-4),即1-2-09>_-209,又由于1-2-09>0和 

一209>0,所以I1-2-091>_1-209l,即 一1+∞l ∞l。 

其次,由于I,tl>-1及0<∞ 1,所以有不等式的 

右端自然成立。 

2主要结论 

定理1 EDo鲥∈J[),则当松弛因子09满足 

0<∞ 1时,对任意初始向量解线性方程组Ax=b的 

SOR迭代法都收敛。 

证明首先证明当AEDo时,SOR迭代法收 

敛。 

由条件知,对角矩阵D=-diag(a a::,…, 

an )≠0。设 为SOR迭代法的迭代矩阵 = 一n ) 

[(1—09)D+缈 的任意特征值,则det(2/-L )=det 

{ 一 一co/;)~I(1—09) +_CO )=0,U为单位矩阵), 

即 

det(D一∞ ) {2(D—coL)一[(1—09)D+09 )=0, 

从而得,det(D一∞ )~・det{2(D一砒)一[(1一∞)D+ 

∞ ):0,由于det(D—coL) ≠0,所以 

必有 

det{ 一1+09)D一209L-co =O (2) 

假设 =(D一09L) [(1_∞)D+∞ 至少存在一个 特征值141>-1。由于A∈D ,所以有 

I I> ) +∽ +∽ (3) 

对f∈ 皖立。又由矩阵 的分解以及矩阵 , 

的结构知, 

( )= (上+to= 肛)十R , 

)= +to=& )斗 (∽ (4) 

从而得到,对任意的iEN, ∈[0,1],皆有 

Ia。4>尺 ) )=[尺旺)+ (∽] ・ )+S(to] 

成立。 

由引理3知,此时 一1+091>-12I.09,上式左边 

乘以14-1+o9I,右边乘以 ∞l,得 

一1+09I.1a,,l>lA091 肛)般l(∽] ・ ∞ IS, )+ 

(∽] [ 09L)+R ∞ ]。・ 砒 ∞to] 

由于141>1,由上式得, 

l 一1+o91.Ia.I>[R,(209L)+R (09∽] 

[S,(209L)+S,(2COto] 

 ̄(209L+09to・ (2mL+m∽= 

彤[ 一1+∞)D— ∞L一∞ ・ 

[ 一l+∞ 一 ∞£一∞ 

此说明 一1+09)D— 础一 ∈Da,由引理1得到 

一l+09)D-209L一∞ 乍奇异,与(2)式矛盾。所 

以121>_1不 =(D一095) [(1一CO CO 的特征值。 

从而得知p(三 )<l,所以对任意初始向量解线性方程 

组Ax=b的SOR迭代法都收敛。 

其次,证明当A∈D时,SOR迭代法收敛。 

由于A∈D,则知,对Vf,JEN,I口胁I>R ) 

母(,4)成立,由以上证明得知,对Vf,J∈N, 

胁I>[R肛)+尺 ( ].[ ) ∽] (5) 

同上,假设L = 一coL) [(1一∞ 09 至少存 

在一个特征值141>_1。注意到引理3的结论,(5) 

式左边乘以14-1+091 ,右边乘以12091 得, 

J 一l+0911[a ,a. ̄I>I2091 ER )+ (∽]‘ 

l ∞I )斗R(to-=JR (209L)+R ∞to]. 

[弓 ∞ )十弓 ∞to] 

同样由于141>_1,上式变为 

一1+09] la, l> (209L)+R (09to]. 

[R ̄(209L)+ 『<∞to]= 

R (209L+09to・R ̄(209L+09∽ (6) 

又由于对Vi∈Ⅳ,总有R (209L+09tO=R 一1+09)D一 

209L—CO 成立。所以(6)式说明 一1+09)

D—

相关主题