广义鞍点问题的GSOR方法及其收敛性-应用数学专业论文.docxVIP

  • 5
  • 0
  • 约1.82万字
  • 约 24页
  • 2019-04-05 发布于上海
  • 举报

广义鞍点问题的GSOR方法及其收敛性-应用数学专业论文.docx

太原理T大学硕士研究生学位论文 如果令 -nD =I I.L = I _, I.U =I I -n O Q) αQF OβQ) 川 )=0(;)+叫巾ν[((Iυ1 一川 则含参数的修正 SOR-like 13)去; ,简记为 MPSOR-Jike,它的法代形式 定义如下, 这理脱矩阳际=竹 r J 3 拟高斯一赛德尔迭代方法 12) 出约束优化 [; S:](:J = ( ) n其小 A E R时为非对称矩 阵,B E Rrn.n(m n) 且 rank(B) = m 记 A = r s:ì为 n t O} 非奇异 阵,当 A x1.fjj\ 正 定或 A 的对称部分 正 定时,对于 I I(x ì= (fì I B 0 I\Y) 飞 g ) 用拟高斯-赛德尔法代方法求解时,其对 应的炖阵分裂为 : 叫:;]-[:11 其小取 Q 为对称jE定矩阵,对应的法代计算 (13) 为, ( 礼Xn.卢1 卢==--A-1B yη叶 = Yn -Q-1Bxl ( A 13 ,(x (b4 ( A 13 ,(x (b 由于问题|二 : 11 1 = 1 1 的系数矩阵是大型稀疏矩阵,很多学者将:ì!51 \.li U)飞y) 飞q) 方法应用到求解鞍点问题上得到了许多有效地方法如 SOR-Jike 方法以及这些 方法的推广等,这些方法都要涉及到求逆矩阵 A-1 这样使得方法的有效性在很 大程度上依赖于矩阵A 的性质和结构,本文构造了 一种新的算法,该算法将通常 4 、BEEEt 、 / B) BT B) BTB ) r_αB 1B) r _αB 1B) l(B o )l(y);)千lB -αB1B X)lq; 记 B,:= BTA 记 B,:= BTA -αBTBBY,Bz:= BTB ,纠=:b ,ql := BTb 一αBTBq , 可以得到如下线性系统 、 ‘、 ‘、 B 、1A:.:= A +αB 、 1 A:.:= A +αBBT -B( 81B)-l BTA , B) :=BT A- αB T ?B T , 乌句 - B(BTBfqp 这早 α 这早 α 是一个非零的实参数 I 足 m 阶单位矩阵注意到 _αBγB Jl BT 0 J l BTA -αBTBB7 B ,何 , B t(lm O)( t (lm O)(A B)=(A  Yn I Y ,,‘ ,, PU . BBE,Et ‘、? , 、 / 统 、m????,, 、 m 一、nuJLU b q 一、 nu J Z 2 AU Qa 系 首先将上述鞍点问题转化为下面的问题性V fIll--、 / / t 首先将上述鞍点问题转化为下面的问题 性 tIll-k q2 = q2 = ql 、SEEE,, \ fe--l = 、线-- 、 线 -- / r \--BE r \ 1 11ll Et l 、j 、 = x Y 、 / 、 ‘、x y x ‘、 、11 F J / / \ 的 III-E J FIll- OE-- FBBEEt 1JSll- 1 J 61llt 价0ll 价 吁 \ J 飞 / 、1 / f \ lIll- B 鸟 B n吨 O 乓 、ill-t 、 / 1 1 7 Y O B A 鸟 A D叫 4 Ad 等 ι、Ill-l ι、 J u oa α tlt、 / 、 / / \ FEltt‘1 、1SEE-, 、1 lt1·、 FE--EEE = ‘、, 、飞 ,/ ,7 DHM , I B T hEEEE ,x y , /plEE飞\ / \ B I 、 , Ill- J BT.Ill BT. ‘,飞 B(BTo B ( BT B O Ar ttt‘ \太原理 正大学硕士研究生学位论文算法中计算 大型矩阵 \ 太原理 正大学硕士研究生学位论文 算法中计算 大型矩阵 A 的逆转化为求一个较小维数的矩阵 BTB 的逆,在 mn 时这种算法是有放的 5 IF ,,1 I / 气 与 It- - 由卜式得 由卜式得: 从而有这里不z gt 此 因 从而有 这里 太原理 r火学硕士研究生学位论文 5 广义 SAOR 方法 117,18,19凡21) 极点问题的基本形式为 (A3X=b2 B3x+B2Y=Q2 (:r )(:)=(:) 其中知阵 A ε R mX fn 是正定 对称矩阵 ,矩阵B E Rmxm 是列满秩的且 m 三 n ,即 rank( 的= n ,向量b l.? Rm ,q 仨 R 是给定的向量. 首先将(; )(:)=( J改写为 (_;r )(:)=(

文档评论(0)

1亿VIP精品文档

相关文档