国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

H—矩陣的預(yù)條件AOR迭代法收斂性

2012-12-31 00:00:00張仕光
考試周刊 2012年67期


  摘 要: 近年來(lái),許多預(yù)條件子被運(yùn)用于線性系統(tǒng).討論了新的多參數(shù)一般下三角預(yù)條件子的AOR迭代法的收斂性.當(dāng)線性系統(tǒng)的系數(shù)矩陣為H-矩陣時(shí),得到了該預(yù)條件子下的AOR迭代法的收斂性定理.
  關(guān)鍵詞: AOR迭代法 預(yù)條件子 H-矩陣 收斂性
  1.引言
  考慮線性系統(tǒng)Ax=b,其中A=(a)∈R,b∈R是已知的,x∈R是未知的.
  不失一般性,令A(yù)=I-L-U,其中I是單位矩陣,-L和-U分別是A的嚴(yán)格下三角部分和嚴(yán)格下三角部分.考慮預(yù)條件線性系統(tǒng)PAx=Pb,其中P是非奇異矩陣.
  文獻(xiàn)[1]考慮具有一般上三角形式的預(yù)條件子,給出當(dāng)線性系統(tǒng)的系數(shù)矩陣為M-矩陣時(shí)預(yù)條件SOR型迭代法與經(jīng)典SOR迭代法的收斂性比較定理.
  考慮一般下三角形式的預(yù)條件子P=I+S,記D(β)=diag(1,β,β,β,…,β),β (i=2,…,n)是非負(fù)實(shí)數(shù),S=D(β)S,m  解決線性系統(tǒng)的經(jīng)典AOR迭代法的迭代矩陣[2]為
  T=(I-rL)((1-ω)I+(ω-r)L+ωU),(1)
  這里ω和r為實(shí)數(shù),且ω≠0.
  定義1.1[3] 如果一個(gè)n×n的矩陣A=(a)滿足i≠j時(shí),a≤0,稱A為Z-矩陣;如果A是Z-矩陣且a>0,稱A為L(zhǎng)-矩陣;如果A是L-矩陣且A≥0,稱A為M-矩陣.
  定義1.2[4] 如果A=(a)是n×n的矩陣,稱〈A〉=()是A的比較矩陣,其中i=j時(shí),=|a|,i≠j時(shí),=-|a|.如果〈A〉是一個(gè)非奇異M-矩陣,則稱為H-矩陣.
  引理1.1[5] 設(shè)是Z-矩陣,則A是M-矩陣當(dāng)且僅當(dāng)存在正向量u=(u,…,u)>0,使得Au>0.
  引理1.2[6] 令A(yù)是H-矩陣,則ρ(T)<1.
  2.主要結(jié)論
  令D,-L,-U分別是A的對(duì)角部分﹑嚴(yán)格下三角部分和嚴(yán)格下三角部分,則
  A=(I+S)A=D-L-U.(2)
  對(duì)應(yīng)的預(yù)條件AOR迭代法的迭代矩陣為
  =(D-γL)[(1-ω)D+(ω-γ)L+ωU].(3)
  定理2.1 令A(yù)是對(duì)角元為1的H-矩陣,a≠0(i=2,…,n),則
  β′=1+>1,i=2,…,n
  證明:因?yàn)锳是H-矩陣,由定義1.2知〈A〉是非奇異M-矩陣,且〈A〉=I-|L|-|U|≤I,得〈A〉≥I≥0,即||〈A〉||≥1,則β′>1,i=2,…,n.
  定理2.2 令A(yù)是對(duì)角元為1的H-矩陣,是(5)中給出的迭代矩陣,a≠0(i=2,…,n),0≤β≤β′,0≤γ≤ω≤1,ω≠0,則A是H-矩陣,且ρ()<1.
  證明:記A=(),當(dāng)i=1時(shí),=a;當(dāng)i=2,…,n時(shí),=a-βaa.令r=〈A〉e,其中e=(1,1,…,1).因A是H-矩陣,由定義1.1和1.2得r≥0,〈A〉r=〈A〉〈A〉e=e,記r=(r,r,…,r),由文獻(xiàn)中的引理3.1得(〈A〉r)=r-|a|r=1>0,則(〈A〉r)=r-|a|r=1>0.令(〈A〉r)(i=2,…,n)是向量〈A〉r的第i個(gè)元素,則
 ?。ā碅〉r)=|1-βaa|r-|a-βaa|r
  ≥r-β|aa|r-|1-β||ar|-|a|r-β|aa|r
  =(〈A〉r)+|a|r-β|a|[-(〈A〉r)m+r]-|1-β||a|r
  =1+β|a|+[(1-β)-|1-β|]|a|r
  當(dāng)0≤β≤1(i=2,…,n)時(shí),有(〈A〉r)≥1+β|a|>0.
  當(dāng)1<β<β′(i=2,…,n)時(shí),有
 ?。ā碅〉r)≥1+β|a|+[(1-β)-|1-β|]|a|r
  =1+β|a|+[(1-β)-|β-1|]|a|r=1-|a|r-β|a|(2r-1)
  >1+2|a|r-(1+)|a|(2r-1)
  ≥(1+|a|)-(2||〈A〉||-1)|a|
  =(1+|a|)-(1+|a|)=0
  因此〈A〉r>0.由引理1.1知〈A〉是非奇異M-矩陣,因此A是H-矩陣,由引理1.2得ρ()<1.
  參考文獻(xiàn):
  [1]蔣小鳳,袁東錦,孫霞,李凱.預(yù)條件SOR迭代法的收斂性[J].江南大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,9(3):339-342.
  [2]A.Hadjimos.Accelerated overrelaxation method[M].Math.Comp,1978,3:149-157.
  [3]D.M.Yong.Iterative solution of large linear systems[M]Academic Press,New York,1971.
  [4]R.S.Varga.Matrix Iterative Analysis[M].Prentice-Hall,Englewood Cliffs,NJ,1981.
  [5]K.Y.Fan.Topological proofs for certain theorems on matrices with non-negative elements[M].Monatsh.Math,1958,62:219-237.
  [6]Yao-Tang Li,Shun-feng Yang.A multi-parameters preconditioned AOR iterative method for linear systems[J].Applied Mathematics and Computation,2008,206:465-473.
  [7]Wu Mei-jun,Wang Li,Song Yong-zhong.Preconditioned AOR iterative method for linear systems[J].Applied Numerical Mathematics,2007,57(5-7):672-685.
  基金項(xiàng)目:2011年衡水學(xué)院科學(xué)研究項(xiàng)目(201102