蔣建新,李艷艷
(文山學(xué)院 數(shù)學(xué)學(xué)院,云南 文山663000)
嚴格對角占優(yōu)矩陣是一類在數(shù)值代數(shù)、數(shù)學(xué)物理和控制論等領(lǐng)域有著廣泛應(yīng)用的特殊矩陣,例如:線性方程組Ax=b,當系數(shù)矩陣A為嚴格對角占優(yōu)矩陣時,許多經(jīng)典的迭代算法均是收斂的,同時對目前提出的一些修正算法也是收斂的,所以在理論探討和實際工作中常要估計矩陣逆的無窮范數(shù),尤其是對大型矩陣的判別,還存在許多困難.經(jīng)過國內(nèi)外許多學(xué)者不懈努力,已獲得一些重要結(jié)果[1]. 本文繼續(xù)研究嚴格對角占優(yōu)M-矩陣A的‖A-1‖∞的上界估計問題,給出其新的收斂的上界序列.
定義1[2]記,稱Zn×n中的矩陣A為Z-矩陣(簡記為;設(shè),如果A可表示為A=αI-P,其中P≥0(即,?i,j∈N),α≥ρ(P),則稱A為M-矩陣(ρ(P)是非負矩陣P的譜半徑).特別,當α >ρ(P)時,稱A為非奇異M-矩陣;當α=ρ(P)時,稱A為奇異M-矩陣.用Mn表示非奇異M-矩陣的集合,q(A)表示非奇異M-矩陣A的最小特征值.
定義2[3]設(shè)且滿足下列條件:
(i,存在非零元素序列aii1ai1i2…airk,其中i≠i1,i1≠i2,…,ir≠k,k∈J(A),則稱A為弱鏈對角占優(yōu)矩陣.
定義3[3]設(shè)A=(aij)∈Rn×n,若J(A)=N,則稱A為嚴格對角占優(yōu)矩陣.
注1:由定義2 和定義3 知,若A為嚴格對角占優(yōu)矩陣,則A為弱鏈對角占優(yōu)矩陣.
引理1[4]設(shè)A=(aij)∈Rn×n是弱鏈對角占優(yōu)M-矩陣,則B=A(2,n)∈R(n-1)×(n-1)也是弱鏈對角占優(yōu)M-矩陣,且B-1=(βij)存在,βij≥0(i,j=2,3,…,n).
引理2[5]設(shè)A=(aij)∈Rn×n是弱鏈對角占優(yōu)M-矩陣,B=A(2,n)∈R(n-1)×(n-1),A-1=(αij),B-1=(βij),則:
其中
若J(A)=N,則Δ≥a11(1-d1l1)≥a11(1-d1).
引理3[4]設(shè)A=(aij)是嚴格對角占優(yōu)的M-矩陣,則Δ≥a11(1-d1l1)>a11(1-d1)>0.
引理4[6]設(shè)A=(aij)∈Rn×n是行嚴格對角占優(yōu)的M-矩陣,則A-1=(αij)滿足:
引理5[6]設(shè)A=(aij)∈Rn×n是行嚴格對角占優(yōu)的M-矩陣,則A-1=(αij)滿足:
引理6[7]設(shè)A=(aij)∈Rn×n是弱鏈對角占優(yōu)M-矩陣,A-1=(αij),令q=q(A),
定理1 設(shè)A=(aij)∈Rn×n是嚴格對角占優(yōu)M-矩陣,則?t=0,1,…,有
由式(5)知:當2≤j≤n時,故當2≤i≤n時
則定理得證.
定理2 設(shè)A=(aij)∈Rn×n是嚴格對角占優(yōu)M-矩陣,則:‖A-1‖∞<Ωt,
定理3 由定理2 得到的‖A-1‖∞的上界序列是單調(diào)遞減的且以‖A-1‖∞為下界,所以該序列是收斂的.
定理5 由定理4 得到的q(A)的下界序列是單調(diào)遞增的且以q(A)為上界,所以該序列收斂.
應(yīng)用文獻[6]中定理3.7,當取迭代總數(shù)T=10 時得‖A-1‖∞≤0.537,應(yīng)用文章[7]中定理3.2 的估計式得‖A-1‖∞≤0.699;應(yīng)用本文定理2,當取迭代總數(shù)T=10 時得‖A-1‖∞≤0.3952;事實上,應(yīng)用Matlab 7.1 計算得
上例表明,在某些條件下本文所得的‖A-1‖∞的上界序列優(yōu)于現(xiàn)有的一些結(jié)果.
[1] Johnson C R.A Hadamard product involvingM-matrix[J].Linear and Multilinear Algebra,1997,4:261-264.
[2] 陳景良,陳向暉.特殊矩陣[M].北京:清華大學(xué)出版社,2000.
[3] 黃廷祝,楊傳勝.特殊矩陣分析及應(yīng)用[M].北京:科學(xué)出版社,2007.
[4] Shivakumar P N,Chew K H.A sufficient condition for nonvanishing of determinants[J].Proc Amer Math Soc,1974,43:63-66.
[5] Cheng G H,Huang T Z.An upper bound for ‖A-1‖∞of strictly diagonally dominantM-matrices[J]. Linear Algebra Appl,2009,43:511-517.
[6] 趙建興.M-矩陣最小特征值估計及其相關(guān)問題研究[D].昆明:云南大學(xué),2014.
[7] 王亞強,李耀堂,孫小軍,等.A new upper bound for ||A-1||∞of strictly diagonally dominantM-matrix[J].山東大學(xué)學(xué)報:理學(xué)版,2010,45(4):43-48.