第5卷第l期
2002年6月
沙洲职业工学院学报
Journal of Shazhou Polytechnical Institute of Technology
V01.5.No.1
Jun.,2002
关于单调变分不等式的一种邻近
增广拉格朗El方法
王治华
(沙洲职业工学院,张家港215600)
摘要:对单调变分不等式的一种新的拉格朗日方法(AL)进行讨论。这种方法只需要解一系列强
单调变分不等式子问题。允许参数从一个叠代点到另一个叠代点变化。在比较弱的条件下算法的
收敛性得到了证明。
关键词:单调变分不等式;邻近点方法;增广拉格朗日方法
中图分类号:0224
文献标识码:A
文章编号:1009—8429(2002)01—0001—06
1 引言
设Q是Rn的一个非空闭的凸集,F是一个从R“到自身的一个单调连续映射,变分不等
式问题是寻找一个向量u+E Q,使得
(U一“+)rF(“。)≥0
V “E n
(1.1)
用VI(Q,F)来表示。
在算子研究,经济平衡,非线性分析中,变分不等式方法是一种经常使用的方法。对于变
分不等式,现在有许多叠代方法n¨2】[3】㈨1 3。。
在这篇文章里,我们关心的是具有下面这种特殊结构的变分不等式:求z。E n
(z—z+)Tf(z。)≥0
n={z I触=b,z E X}
V
M∈0
(1.2)
(1.3)
其中X是Rn的非空闭凸子集,A∈R“”,b∈R”,f:X—X是一个给定的单调算子,
问题(1.2)一(1.3)等价于下面一个问题:
求z’∈X,Y’∈R”
(z—z。)丁(f(z+)一A0’)≥0
Ax+=b
(1.4)
(1.5)
这个问题的解集用Q。表示。对于问题(1.4)一(1.5),可以用增广拉格朗日方法来
解吲㈨Ⅲ。增广拉格朗日方法的思路如下,从给定的一对x,Y)开始,它的校正(文,歹)是通过下
面的计算得到
(z’一i)r{,(童)一A丁[Y—fl(AX—b)]}≥0
(1.6)
收稿日期:2002—04—10
作者简介:王治华(1969一),男,沙洲职业工学院基础科学系讲师,南京大学硕士生。
万方数据
·
1
·
歹=Y—p(Aj—b)
(1.7)
由于求解的次数明显地依赖于惩罚因子0,Wang和Liao在文献[9]提出了一个修整参数
p的办法,在他们的算法中允许参数|3随叠代点而变化,在文献[11]里,Kontogiorgis和Meyer
给出了另一个算法,允许惩罚因子是一个正定矩阵。B.S.H发展了这方法。川。最近B.S.He
对单调变分不等式提出了一种新的不精确交替方向法,只要不精确的解子问题就可以了馏】。
受上述这些算法的启示,我们讨论一种临近的增广拉格朗日方法:
给定点(x,y),我们通过下面的方法得到新的点(曼,歹):
X’一i)丁{厂(王)一A丁[y—fl(Ai一6)]+),(童一z)}≥0
b:Y一卢(A奎一6)
其中口,y是给定的正常数。,
(1.8)
如果.厂是单调的,那么(1.8)变成了强单调。在下面的算法中,允许卢,y被~系列正定
矩阵{H。},{R。}代替。
2预备知识
用II x|l=~/zTz来表示欧几里德范数,设n是R”的一个非空闭凸子集。欧几里德范
数下的投影用Pn(.)表示,即:
Pn(副)=argmin{{f可一“fl f M∈o i}
设F是一个Q上的一个映射,如果
(乱一口)T(F(“)一F(勘))≥0
V“,“∈n
(2.1)
那么称F在Q上单调。
由于我们用一个正数口乘变分不等式VI(Q,F)是不变的,有下面众所周知的引理:
弓l理1: 设Q是R“的一个非空闭凸子集,那么U。是VI(Q,F)的一个解的充要条件
U。=P“U’一肛(U。)]
因此,解一个变分不等式等价于找残量函数
e(“,卢)=甜一Pn(“~犀F(U))
的一个零点。一个基本的投影性质是
(u—Pn(口))r(Pn(u)一甜)≥0
V可∈R”
V“∈n
可以推出
l|Pn(口)一“lI 2≤l|口一“lI 2~l|勘一Pn(铆)0 2
II Pn(甜)一Pn(口)I|2≤≤Il“一口l|2
V“,可∈R“
V M∈n V u∈R”
设u=(X,y)T∈W=X×R“,问题(1.4)一(1.5)等价于求x’∈X,y。∈R“,满足
{(∥一z’)T(f(z’)一ATy’)≥0
((y—y’)T(Ax’一6)≥0
所以解问题(1.4)一(1.5)相当于找出
V z7∈X
V y∈R”
g(“)=-Px
t
x-If(z’-Ary]}Ax
b
的零点。在下面的讨论中作一个标准的假设:
/、
.
万方数据
(2.2)
(2,3)
(2.4)
(2.5)
(2.6)
(2.7)
(2.8)
假设A:x是非空闭凸集,Q+是非空。
3临近点拉格朗日方法
为了保证算法的收敛性,给出下列条件,
条件c:设{讯}孑是一个非负序列,且∑孑Tlko,帮^≤地+·≤(1+r/k)Hk
R>0,R≤R川≤(1+仇)R^
1:
2:
由数学分析可知,∑仉<。。等价于Ⅱ(1+啦)(≥)H是指H,一H是正定(半正定)。
注:如果R6=0,H。=口,那么这个算法就是(1.6)一(1.7)。由R。>0(指R是正定
的),只要算子f(z)单调的,就保证了算子^的强单调,这样就保证了子问题(3.2)的可解
性,并且解是唯一的。
为了方便起见,记
“=G)柚=f‘可)
引理2: 在假设A下,算法1产生的序列满足
(M‘一“+)_D^(“6一M6十1)≥||“‘一M6+1|J D^
证明:根据(3.2)
(z+一z‘+1)一厂(z‘+1)≥0
在由(1.4)
(z’一.717‘+1)了’(一f(z+)+A7≥’)≥o
(3.6)+(3.7),再利用f的单调性和(3.3)
万方数据
(3.5)
(3.6)
(3.7)
3 ·
(夕一Y。)1了i1(矿一yk“)+ (z6一z+)TRt(z‘ 一z‘+1)≥ ||Ax‘”一b||备。 +
^
z—z
^+l
2
lI
|l R5
引理得证。
定理1: 在假设A下,如果条件C满足,算法1产生的序列具有:
lI“H1一“。II乞川≤(1+仇)II越6一“+l|毛。一II“6一“¨1 l|乞。
(3.8)
证明,由条件C可知
II“¨1一乱’l|乞川≤(1+仉)lI“¨1一M’||乞。
≤(1+仉)Il(U‘一“+)+(“‘+1一“‘)ll乞;
≤(1+玑)|J“6一“’IJ乞。+Jl甜‘+1一“。JI乞。
一2(“6一“’)rD^(U‘一M6+1)
≤(1+仇I|f“。一“’ff乞。一if“。一2‘¨1 f|乞。
根据引理2,得到
4收敛性分析
引理3: 在假设A下,如果条件C满足,算法1产生的序列具有:
lira
Z—+∞
II M‘~M针1 Il 2D.=0
证明:由定理1可得
l|“‘+1一M。l|乞。+。≤(1+叩。)I|“‘一“。II乞。≤…
(4.1)
k
≤Ⅱ(1+17;)lI让。一“。Il乞。≤C,lI“。一M+Il乞。
(4.2)
所以…“‘一“。II乞。}有界,存在常数L>O。使得I|乱‘一“+Il乞。≤L
根据定理1(3.8),两边相加
“
M一
2
D
‘
“
一
比
/一
“
一
甜
2
D
2
D
。量
此因
lira
Il“‘一“¨1 I|2DI=0
+
∑仉II“‘一“+lI乞。
^=O
+
L∑仇=||“。一M’I‰+LC,
(4.3)
定理2: 在假设A下,如果条件C满足,算法1产生的序列具有:
lira e(Ҭ1)=0
其中e(u)是(2.8)所定义
证明:
忆∽+1)俨=扩1。刚奢∥+1卜的¨1’112
(4.4)
根据(3.2),得到
万方数据
ll
∥z¨1。搿≥吼¨L盯bH。“y+1¨I
触针1—6
I
利用(2.6)
≤I
l
一
砭圻
¨z∥
一
缸“y
,0
=||“‘一M¨1 Il 2D^2
(4.5)
由范数的等价性可知,存在常数p>0,使得
|l U4一U¨1 JI乞。:≤矿f|U‘一“¨1 l|乞。
所以
P(甜针1)≤卢2 0“‘一U¨1 I|乞。
再根据引理3,有lira e(uk+1)=0
^—+∞
定理3: 在假设A下,如果条件C满足,{u。}是算法1产生的序列,则{u。}收敛于问题
(1.4)一(1.5)的一个解。
证明:根据(4.2),{uk}有界,所以{uk}至少有一聚点,设U‘是其一聚点,且一子列{uk,l收
敛于u’,所以
V愚>k。
V e>0,]ko>0,使得V k。≥k。,(“~∈{扎‘-}),均有
||“‘t—U。||<£/cp
|I M‘一“+II≤Ⅱ(1+琅)I|“‘。一“’||≤Cp l|M%一“’||<£
i 2
ko
反以lim“‘=“‘,又因为e(u)是连续的,故可得
^—-∞
e(“。)=lira e(“‘)=0
k一∞
根据引理1,u+是问题(1.4)一(1.5)的解。
在实际中,{H。}和{R。}如何选取值得进一步研究。
参考文献:
[1]Noor,m.A.Some recent advances in variational inequalities,part I:Basic concepts[J].New Zealand Journal of
Mathematics,1997,26:53—80.
[2]Dafermos,S.An iterative scheme for variational inqualities[J].Mathematical Programming,1983,26:40—47.
[3]Harker,P.T,Pang,J.S.Finite—dimensional variational inequality and nonoelinear complementarity problems,a
survey of theory,algorithms and applications[J].Mathematical Programming,1990,25:247—262.
[4]B.S.He.A projection and contraction method for a class of linear complementarity problems and its application in
convex quadratic programmingJ[J].AppHed Mathematics and Optimization,1992,25:247—262.
[5]Gabay,D.“Application of the method of multipliers to variational inequalities,”in Augmented Lagrange
Methods:Applieatons to the Solution of boundary—valued Problems.M.Fortin and R.Glowinsky,eds.,North
Holland,Amsterdam,The Netherlands,1983,299—331.
[6]D.Gabay and B.Mercier,“A dual algorithm for the solution of nonlinear variational problems via finite element
approximations.”Computer and Mathematics with Applications 1976,2:17—40.
[7]B.S.He,L.z.Liao and H.Yang,“A decomposition method for a class of monotone variational inequality
problems.”JOTA,1999,103(3):603—622.
[8]B.S.He,L.Z.Liao,D.R.Han and H.yang,A new inexact alternating directions method for monotone
variational inequahties,Mathematical Programming,2002,6:1.
万方数据
· 5 ·
[9Is.L.Wang and L.一Z.Liao,“Decomposition method with a variable parameter for a class of monotone
variational inequality problems.”JOTA 2000,106:349—368.
[10]D.P.Bermekas and J.N.Tsitsiklis,ParaUed and Distributed Computation:Numerical methods,Prentice—Hall
Englewood Cliffs,NJ,1989.
[11]Kontogiorgis,S.and Meyer,R.R,A variable—penaky alternating directions method for convex optimization,
MathmaticaI Programming 1998,83:29—53.
[12]Tseng,P.,Alternating projection—proximal methods for convex programming and variational inequilties,
SIAM J.Optimization 7,951—965.
[13]B.S.He and Yang Hai,Some convergence properties of a method of multipliers for linearly constrained
monotone variational inequalities operations Research Letters,1998,23:151—161.
A Proximal Augmented Lagrange Method for Monotone
Varational Inequailities
WANG Zhi—hua
(Shazhou Polytechnical Institute of Technology)
Abstract:This paper addresses a new augmented Lagrange method(AL)for mono—
tone variational inequality(VI),which needs only to solve one stongly monotone sub
—VI problem.The paramaters are allowed to vary from iteration to iteration.Its con.
vergence of the method is proved under mild assumptions.
Key words:Monotone variational inequalities;Proximal point method;Augmented
Lagrange method
·6 ·
万方数据
关于单调变分不等式的一种邻近增广拉格朗日方法
作者:
王治华, WANG Zhi-hua
作者单位:
刊名:
沙洲职业工学院,张家港,215600
沙洲职业工学院学报
英文刊名:
JOURNAL OF SHAZHOU POLYTECHNICAL INSTITUTE OF TECHNOLOGY
2002,5(1)
年,卷(期):
参考文献(13条)
1.Noor, m. A Some recent advances in variational inequalities, part Ⅰ: Basic concepts 1997
2.Dafermos, S An iterative scheme for variational inqualities 1983
3.Harker, P. T;Pang, J. S Finite - dimensional variational inequality and non0elinear
complementarity problems, a survey of theory , algorithms and applications 1990
4.B. S. He A projection and contraction method for a class of linear complementarity problems and
its application in convex quadratic programmingJ 1992
5.Gabay, D Application of the method of multipliers to variational inequalities 1983
6.D. Gabay;B. Mercier A dual algorithm for the solution of nonlinear variational problems via finite
element approximations 1976
7.B. S. He;L.z. Liao;H. Yang A decomposition method for a class of monotone variational inequality
problems[外文期刊] 1999(03)
8.B. S. He;L. Z. Liao;D. R. Han;H. Yang A new inexact alternating directions method for monotone
variational inequalities[外文期刊] 2002
9.S. L. Wang;L. - Z. Liao Decomposition method with a variable parameter for a class of monotone
variational inequality problems[外文期刊] 2000
10.D. P. Bertsekas;J. N. Tsitsiklis Paralled and Distributed Computation: Numerical methods 1989
11.Kontogiorgis, S;Meyer, R. R A variable - penalty alternating directions method for convex
optimization[外文期刊] 1998(1)
12.Tseng, P Alternating projection- proximal methods for convex programming and variational
inequilties[外文期刊]
13.B S He;Yang Hai Some convergence properties of a method of multipliers for linearly constrained
monotone variational inequalities operations[外文期刊] 1998(3/5)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_szzygxyxb200201001.aspx