APP下载

基于同伦方法构造绝对值方程解存在的条件

2018-11-28姜兴武姜舶洋杨雪莹王秀玉

吉林大学学报(理学版) 2018年6期
关键词:方程解有界将式

姜兴武, 姜舶洋, 杨雪莹, 王秀玉

(1. 吉林工商学院 基础部, 长春 130507; 2. 北京工业大学 信息学部, 北京 100022;3. 北京邮电大学 国际学院, 北京 100876; 4. 长春工业大学 基础科学学院, 长春 130012)

1 引言与预备知识

考虑如下形式的绝对值方程:

Ax-|x|=b,

(1)

其中:A∈n×n;b∈n; |x|表示向量x的每个分量取绝对值. 目前, 关于绝对值方程的研究已有很多结果: 文献[1]获得了绝对值方程与水平线性互补问题的等价关系; 文献[2-4]描述了求解绝对值方程的有效算法; 文献[5-7]研究了绝对值方程解存在的条件. 组合同伦方法是求解优化及互补问题的一个有效方法[8-9]. 本文结合绝对值方程与互补问题的关系, 运用组合同伦方法给出绝对值方程解存在的一个条件, 该条件不同于目前常用的正则性条件. 实例分析表明, 该条件不满足正则性, 但满足本文条件的绝对值方程有解且解不唯一. 用xT表示向量x的转置, (xT,yT,zT)T记为(x,y,z).分别表示向量的分量非负(正)的向量集合. 若也记为x≥0(x>0),I表示单位矩阵.

定义1[2]如果对任意的矩阵B∈[A-I,A+I]均为非奇异的, 则方阵区间[A-I,A+I]称为正则的; 否则该矩阵区间称为奇异的.

引理1[1]绝对值方程(1)等价于如下广义线性互补问题:

y=(A-I)x-b≥0,z=(A+I)x-b≥0,yTz=0.

由引理1知, 求解绝对值方程问题可转化为求解广义线性互补问题.

假设:

(H1) 矩阵A可逆, 且有

[(A-I)x]∘[(A+I)x]≤0⟹[(A-I)x]∘[(A+I)x]=0,

即矩阵对(A-I,A+I)为列充分矩阵对, 其中:I为单位矩阵; ∘表示对应分量乘积.

(H2) 绝对值方程Ax-|x|=0只有零解x=0, 即广义线性互补问题

y=(A-I)x≥0,z=(A+I)x≥0,yTz=0

只有零解.

下面举例说明条件(H1),(H2)与区间矩阵[A-I,A+I]的正则性不同.

若(A-I)x∘(A+I)x≤0, 则有

若x1+x2>0, 由(x1+x2)(3x1+x2)≤0及(x1+x2)(x1+3x2)≤0, 必有3x1+x2≤0和x1+3x2≤0, 这与x1+x2<0矛盾. 若x1+x2<0, 由(x1+x2)(3x1+x2)≤0且(x1+x2)(x1+3x2)≤0, 则有3x1+x2≥0和x1+3x2≥0, 这与x1+x2<0矛盾. 因此x1+x2=0, 即假设(H1)成立.

对于

2 主要结果

(2)

其中Y,Y(0)分别为以向量y,y(0)的元素为对角元素的对角矩阵.

其有唯一解x=x(0),y=y(0),z=z(0). 若μ=0, 则式(2)为

证明: 由于w(0)也可视为变量, 因此H的Jacobi矩阵为

非奇异, 其中μ∈(0,1]. 从而0是H(w,w(0),μ)的正则值, 其余证明与文献[9]中命题2.1相同, 故略.

(A-I)x(k)-b-y(k)+μk(A+I)x(k)=0;

(3)

由同伦方程(2)的第二个式子得

(1-μk)[(A+I)x(k)-b]-z(k)+μkz(0)=0;

(4)

由同伦方程(2)的第三个式子知, 对任意的i=1,2,…,n, 有

(5)

情形1)μ*≠1. 将式(3),(4)两边同乘1/‖x(k)‖, 得

(6)

(7)

(A-I)x(*)-y(*)+μ*(A+I)x(*)=0,

(8)

(1-μ*)(A+I)x(*)-z(*)=0.

(9)

由式(8),(9)得

(10)

(11)

由式(5)可知, 对任意的i=1,2,…,n, 有

(12)

若μ*=0, 则‖x(*)‖=1, 式(10)~(12)与条件(H2)矛盾, 因此μ*≠0. 由式(12), 对任意的i=1,2,…,n, 有

(13)

结合式(10),(11),(13)及条件(H1), 可得z(*)=0. 但z(*)=0, ‖x(*)‖=1, 式(10),(11),(13)与条件(H2)矛盾.

情形2)μ*=1.

① {(1-μk)x(k)}仍无界, 即‖(1-μk)x(k)‖→∞,k→∞. 将式(3)两边同乘(1-μk), 得

(1-μk)(A-I)x(k)-(1-μk)b-(1-μk)y(k)+μk(1-μk)(A+I)x(k)=0.

(14)

将式(4),(14)两边同乘1/‖(1-μk)x(k)‖, 得

(15)

(16)

(17)

(A-I)x(*)-y(*)+(A+I)x(*)=0.

(18)

由式(17),(18), 得

(19)

由式(5)得

(20)

由式(20)得

(21)

② {(1-μk)x(k)}有界. 注意到{(1-μk)x(k)}有界, 有界序列必有收敛子列, 其收敛子列仍记为{(1-μk)x(k)}. 令

对式(4),(14)两边取极限, 并注意到μ*=1,

(22)

由式(5)得

(23)

由式(23), 得

(24)

z(k)→z(0).

(25)

由式(5),(25)及μ*=1得

y(k)→y(0).

(26)

由式(3),(26)条件(H1)和μ*=1, 得

x(k)→x(0).

(27)

式(25)~(27)与{(x(k),y(k),z(k))}无界矛盾. 因此,Γw(0)为有界光滑曲线.

证明: 由定理1和定理2知,Γw(0)为有界光滑曲线. 由一维流形分类定理知,Γw(0)或者微分同胚单位圆周, 或者微分同胚单位区间(0,1]. 而

非奇异, 因此Γw(0)只能微分同胚单位区间. 令(w(*),μ*)为Γw(0)的极限点, 则可能出现以下4种情形:

猜你喜欢

方程解有界将式
Navier-Stokes-Coriolis方程解的长时间存在性
平均值不等式的引伸
AKNS方程的三线性型及周期孤立波解
对数Bloch空间上某些算子有界的充要条件
指数有界双连续n阶α次积分C群的次生成元及其性质
单自由度系统
浅谈正项有界周期数列的一些性质
几类可积的Riccati方程解的性质
一类Kirchhoff-Poisson方程解的存在性
基于sub-tile的对称有界DNA结构自组装及应用