杨泰山,姜兴武,王秀玉
(1.吉林大学 数学学院,长春130012;2.吉林工商学院 基础部,长春130062;3.长春工业大学 基础科学学院,长春130012)
线性互补问题是:求x≥0,使得y=Mx+q≥0且有xTy=0,其中M为一个n阶方阵,该互补问题记为LCP(M,q),称为非齐次互补问题.当q=0时,记为LCP(M,0),称为齐次互补问题.线性互补问题是互补问题的一个重要组成部分,并且二次规划的K-K-T方程也为线性互补问题.
同伦方法由于具有大范围收敛性,目前已成为求解数学问题的一个重要工具.文献[1]构造了一类同伦方程求解互补问题;文献[2-6]对文献[1]的同伦方程给出了不同条件下同伦路径的存在性;文献[7]建立了与文献[1]完全不同的同伦方程,获得了互补问题的可解性,但文献[7]的同伦方程当参数为零时不能回到原互补问题,也未给出具体条件;文献[8]改进了文献[7]的结果,利用同伦方法对互补问题进行求解.本文运用文献[1]的同伦方程给出半单调线性互补问题同伦路径的存在性、有界性及收敛性,并给出LCP(M,q)有解与LCP(M,0)只有零解的关系.
本文用x≥0或x∈ℝn+(x>0或x∈ℝn++)表示向量x的每个分量为非负(正数),用w=(x,y)表示向量w=(xT,yT)T.
定义1[9]如果对任意的x≥0,x≠0,存在一个分量xk>0,使得(Mx)k≥0,则n阶方阵M∈ℝn×n称为半单调矩阵.
当M为半单调矩阵时,互补问题LCP(M,q)称为半单调线性互补问题.
假设条件:
(H1)M∈ℝn×n是半单调矩阵;
(H2)y=Mx,x≥0,y≥0,xTy=0只有零解;
(H3)存在u∈ℝn++,使得v=Mu+q>0;
(H4)存在常数τ≥0,α≥0,1<β<2,使得对任意的x∈ℝn,y∈ℝn,下式成立:
例1
则M显然为半单调矩阵.令
则M不是P0矩阵.
只有零解.M 满足假设条件(H1),(H2).
记w=(x,y),任取x(0)>0,y(0)>0及w=(x(0),y(0)),构造如下同伦方程:
其中X=diag(x).同伦方程(1)也记为Hw(0)(w,μ),并记
证明:Γw(0)的存在性由定理1可得.Γw(0)⊆ℝn+×ℝn+×(0,1]显然成立.若Γw(0)无界,则必存在子列(x(k),y(k),μk)∈Γw(0),使得当k→∞时,有‖(x(k),y(k),μk)‖→∞,由同伦方程(1)的第二个等式得
由式(2)可知
由同伦方程(1)的第一个等式得
情形1)μ*∈[0,1).
由式(2)及x(*),y(*)的定义,对任意的i=1,2,…,n,得
又由式(6)知(x(*),y(*))为齐次互补问题LCP(M,0)的非零解,与假设(H2)矛盾.
情形2)μ*=1.分两种情形论证.
令
而由式(2)可知,对任意的i=1,2,…,n,有
因而由式(2),(9)得
式(10)与{x(k)}无界性矛盾.
式(4)两边取极限得
再由式(2)得
由式(4)-式(13),得
将式(2)代入式(14)的分量形式,对所有的i=1,2,…,n,得
由式(15)可知
由式(16)可知,当x(k)i→∞时,有
且对所有的i=1,2,…,n,有
由于{x(k)}为无穷序列,因而存在整数s,p,使得
显然有
由式(18),(19)可知,对充分大的k,有
由条件(H3)知下式成立:
整理式(21)得
证明:由定理1~定理3易知Γw(0)为有界曲线.由一维流形分类定理知,Γw(0)微分同胚于单位圆周或单位区间(0,1](证明与文献[7]的定理2.1类似).注意到
是非奇异的,得Γw(0)不能微分同胚于单位圆周,而只能微分同胚于单位区间.记(w(*),μ*)为Γw(0)的极限点,则只有下列4种情形可能发生:
[1]Kojima M,Megiddo N,Mizuno M.A General Framework of Continuation Methods for Complementarity Problems[J].Math Oper Res,1993,18(4):945-963.
[2]XU Qing,DANG Chang-yin.A New Homotopy Method for Solving Non-linear Complementarity Problems[J].Optimization,2008,57:681-689.
[3]YU Qian,HUANG Chong-chao,WANG Xian-jia.A Combined Homotopy Interior Point Method for the Linear Complementarity Problem [J].Applied Mathematics and Computation,2006,179(2):696-701.
[4]ZHAO Yun-bin,LI Gong-nong.Properties of a Homotopy Solution Path for Complementarity Problems with Quasi-monotone Mappings[J].Applied Mathematics and Computation,2004,148:93-104.
[5]LI Gong-nong.Analysis for a Homotopy Path of Complementarity Problems Based onμ-Exceptional Family[J].Applied Mathematics and Computation,2005,169(1):657-670.
[6]WANG Xiu-yu,JIANG Xing-wu, LIU Qing-huai.The Combined Homotopy Method for Nonlinear Complementarity Problems[J].Acta Mathematicae Applicatae Sinica,2012,35(3):430-440.(王秀玉,姜兴武,刘庆怀.非线性互补问题的组合同伦算法 [J].应用数学学报,2012,35(3):430-440.)
[7]DING Jun-di,YIN Hong-you.A New Homotopy Method for Nonlinear Complementarity Problems [J].Numericla Mathematics,A Journal of Chinese Universities:English Series,2007,16(2):155-163.
[8]WANG Xiu-yu, JIANG Xing-wu, LIU Qing-huai.New Homotopy Method for Solving Nonlinear Complementarity Problems[J].Journal of Jilin University:Science Edition,2012,50(3):494-498.(王秀玉,姜兴武,刘庆怀.求解互补问题的新同伦算法 [J].吉林大学学报:理学版,2012,50(3):494-498.)
[9]韩继业,修乃华,戚厚铎.非线性互补理论与算法 [M].上海:上海科学技术出版社,2006.