第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—