• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      求解對稱不定線性系統(tǒng)的吉爾-默里強(qiáng)迫正定方法

      2014-10-09 03:39:02
      關(guān)鍵詞:徐成默里吉爾

      程 軍

      (曲靖師范學(xué)院教師教育學(xué)院,云南曲靖655011)

      考慮如下2×2塊狀線性系統(tǒng)

      其中,A∈Rn×n是對稱不定矩陣,B∈Rm×n(m≤n)滿秩,即秩(B)=m,令BT表示B的轉(zhuǎn)置.向量x,f∈Rn,y,g∈Rm.在此假設(shè)條件下易知線性方程組(1)的解是存在且唯一的,并且方程組的系數(shù)矩陣是非奇異的.具有形如方程組(1)的線性系統(tǒng)有許多實(shí)際應(yīng)用背景,如計算流體力學(xué)[1-2]、電磁計算[3]、Stokes方程和二階橢圓形的混合有限元方法,以及帶約束的優(yōu)化問題等[4-13].線性系統(tǒng)(1)中的A矩陣為對稱正定或?qū)ΨQ半正定的,有許多不同的迭代方法來求解這類問題[8-9],但是當(dāng)(1,1)塊矩陣A是不定矩陣的研究工作相對來說則少很多.本文針對系數(shù)矩陣(1,1)塊矩陣A是不定矩陣,運(yùn)用吉爾-默里強(qiáng)迫正定分裂方法[14]使分解成一個對稱正定矩陣和一個對角矩陣,構(gòu)造一個新的迭代方法,并給出該算法的收斂條件.

      1 吉爾 -默里強(qiáng)迫正定分解算法

      2 吉爾-默里強(qiáng)迫正定迭代方法的收斂性分析

      3 數(shù)值算例

      表1 吉爾-默里強(qiáng)迫正定迭代方法的迭代數(shù)及運(yùn)行時間Table 1 Number of iterations and running time of Gill-Murry forced positive definite splitting methods

      表1列出了迭代矩陣G的譜半徑的值以及迭代格式(5)收斂所需要的時間.由結(jié)果可知迭代格式(5)收斂,故此算法是有效的.

      [1]Cliffe K A,Garratt T J,Spence A.Eigenvalues of block matrices arising from problems in fluid mechanics[J].SIAM J Matrix Analy Appl,1994,15:1310-1318.

      [2]Glowinski R.Finite element methods for incompressible viscous flow[C]//Handbook of Num Anal.Amsterdam:North-Holland,2003.

      [3]Arbenz P,Geus R.Multilevel preconditioned iterative eigensolvers for Maxwell eigenvalue problems[J].Appl Num Math,2005,54:107-121.

      [4]Zhou Y Y,Zhang G F.A generalization of parameterized inexact Uzawa methods for generalized saddle point problems[J].Appl Math Comput,2009,215:599-607.

      [5]Ling X F,Hu X.On the iterative algorithm for large sparse saddle point problems[J].Appl Math Comput,2006,178:372-379.

      [6]Jiang M Q,Cao Y.On local Hermitian and skew-Hermitian splitting iteration methods for generatized saddle point problems[J].J Comput Appl Math,2009,231:973-982.

      [7]Bai Z Z,Wang Z Q.Restrictive preconditioners for conjugate gradient methods for symmetric positive definite linear systems[J].J Comput Appl Math,2006,187:202-226.

      [8]Cao Z H.Fast Uzawa algorithm for generalized saddle point problems[J].Appl Num Math,2003,46:157-171.

      [9]Chen F,Jiang Y L.A generalization of the inexact parameterized Uzawa methods for saddle point problems[J].Appl Math Comput,2008,206:765-771.

      [10]Cao Z H.Constraint Schur complement preconditioners for nonsymmetric saddle point problems[J].Appl Num Math,2009,59:151-169.

      [11]Bai Z Z.Structured preconditioners for nonsingular matrices of block two-by-two structures[J].Math Comput,2006,75:791-815.

      [12]Bai Z Z,Parlett B N,Wang Z Q.On generalized successive overrelaxation methods for augmented linear systems[J].Num Math,2005,102:1-38.

      [13]Cui M R.Analysis of iterative algorithms of Uzawa type for saddle point problems[J].Appl Num Math,2004,56:133-146.

      [14]徐成賢,陳志平,李乃成.近代優(yōu)化方法[M].北京:科學(xué)出版社,2002:62-67.

      [15]Cao Z H.Positive stable block triangular precondetioners for symmetric saddle point problems[J].Appl Num Math,2007,57:899-910.

      [16]Bai Z Z,Pan J Y,Golub G H.Preconditioned Hermitian and skew-Hermitian splitting methods for non-Hermitian positive semidefinite linear systems[J].Num Math,2004,98:1-32.

      [17]Bai Z Z,Pan J Y,Ng M K.New preconditioners for saddle point problems[J].Appl Math Comput,2006,172:762-771.

      [18]Bai Z Z,Golub G H,Pan J Y.Preconditioned Hermitian and skew-Hermitian splitting methods for non-Hermitian positive semidefinite linear systems[J].Num Math,2004,98:1-32.

      [19]程云鵬.矩陣?yán)碚揫M].西安:西北工業(yè)大學(xué)出版社,2005:266-271.

      猜你喜歡
      徐成默里吉爾
      一點(diǎn)小忙
      尋找寶藏
      歡迎來到恐龍園
      笑對人生
      城市住宅小區(qū)園林規(guī)劃設(shè)計探討
      拿開以后
      讀者(2018年21期)2018-10-17 03:31:08
      教授猛如虎
      桃之夭夭B(2018年6期)2018-09-14 10:55:38
      吉爾伽美什,尋找永生的故事
      兒童時代(2017年10期)2017-06-21 10:00:12
      人類進(jìn)化得真快
      西湖(2017年4期)2017-04-07 07:46:56
      照相日
      万全县| 松溪县| 香河县| 建平县| 彰化市| 大理市| 安顺市| 卓资县| 合水县| 临沭县| 乃东县| 临漳县| 漯河市| 永康市| 东兰县| 江北区| 陆良县| 临沧市| 安宁市| 纳雍县| 平塘县| 海晏县| 吐鲁番市| 太原市| 荔浦县| 土默特右旗| 常熟市| 阿尔山市| 平山县| 西城区| 天水市| 荣昌县| 东乡族自治县| 辉南县| 财经| 双桥区| 乐都县| 渭源县| 伊宁市| 百色市| 武威市|