-
线性互补问题在众多领域具有广泛的应用,如数学规划、最优停止、期权定价、市场均衡、自由边界和弹性接触等[1-3].线性互补问题的数学模型为:求z∈Rn,满足:
记作LCP(A,q),其中A=(aij)∈Rn×n和q∈Rn为给定的矩阵和向量.
线性互补问题不论是解的存在性、唯一性,还是算法的收敛性,都与矩阵A的结构有着密切的关系.如当矩阵A是P-矩阵(即它的所有主子式都是正的[4])时,线性互补问题有唯一解[2].文献[5]给出了如下结果:设A是P-矩阵,则
其中x*是LCP(A,q)的解,r(x)=min{x,Ax+q},D=diag(d1,…,dn)(0≤di≤1).近年来,众多学者对特殊结构矩阵线性互补问题解的误差界、扰动界进行了分析,给出了一系列估计式[6-13].本文在文献[14]的基础上,继续讨论弱链对角占优B-矩阵线性互补问题解的误差界,获得了新估计式.同时,证明了新估计式在一定条件下优于文献[12-13]的结果.最后,用数值算例表明了所得结果的有效性.
全文HTML
-
用
$\mathbb{N}$ +表示全体正整数的集合.记A=(aij)∈Rn×n,若则称A为严格对角占优的[15];若aij < 0(i≠j),则称A为Z-矩阵[4];若A为Z-矩阵且A-1≥0,则称A为M-矩阵[15];若|aii|≥ri(A)(∀i∈
$\mathbb{N}$ +),且对每个i∉J(A)={i∈$\mathbb{N}$ +:|aii|>ri(A)}≠Ø,存在非零元素链aii1,ai1i2,…,airj,满足j∈J(A),则称A为弱链对角占优矩阵[15].定义1[4] 设A=(aij)∈Rn×n,若对任意的i,j∈
$\mathbb{N}$ +,有:则称A为B-矩阵.
定义2[13] 设M=(mij)∈Rn×n,若M可表示为M=B++C,其中:
ri+=max{0,mij:j≠i},B+是弱链对角占优矩阵且其主对角元素为正,则称M为弱链对角占优B-矩阵.
引理1[14] 设A=(aij)∈Rn×n是弱链对角占优的M-矩阵,则
其中
$\prod\limits_{j = 1}^{i - 1} {\left({1 + \frac{{{u_\mathit{i}}\left(\mathit{\boldsymbol{A}} \right)}}{{1 - {u_\mathit{j}}\left(\mathit{\boldsymbol{A}} \right){h_\mathit{j}}\left(\mathit{\boldsymbol{A}} \right)}}} \right)} = 1$ ,uj(A),hj(A)见后面的记号说明.引理2[12] 设γ>0且η≥0,则对任意的x∈[0, 1],有:
为叙述方便,给出如下记号:设A=(aij)∈Rn×n,i,j,k∈
$\mathbb{N}$ +,i≠j,i+1≤k≤n,记:
-
文献[9]给出结果:设M=(mij)∈Rn×n是B-矩阵,可表示为M=B++C,其中B+=(bij)形如(1)式,则
其中
$\mathit{\beta = }\mathop {{\rm{min}}}\limits_{i \in {{\rm{\mathbb{N}}}_ + }} \left\{ {{\mathit{\beta }_\mathit{i}}} \right\}$ 和${\mathit{\beta }_\mathit{i}} = {b_{\mathit{ii}}} - \sum\limits_{j = 1, j \ne i}^n {\left| {{b_{ij}}} \right|} $ .文献[12]给出了优于(2)式的估计式:设M=(mij)∈Rn×n是B-矩阵,可表示为M=B++C,其中B+=(bij)形如(1)式,则
其中
${\mathit{\bar \beta }_\mathit{i}} = {b_{ii}} - \sum\limits_{j = i + 1}^n {\left| {{b_{ij}}} \right|} {l_i}\left({{\mathit{\boldsymbol{B}}^\mathit{\boldsymbol{ + }}}} \right) = \mathop {\max }\limits_{k \le i \le n} \left\{ {\frac{1}{{{b_{ii}}}}\sum\limits_{j = k, j \ne i}^n {\left| {{b_{ij}}} \right|} } \right\}$ 且$\prod\limits_{j = 1}^{i - 1} {\left({1 + \frac{1}{{{{\mathit{\bar \beta }}_\mathit{j}}}}\sum\limits_{k = j + 1}^n {\left| {{b_{jk}}} \right|} } \right)} = 1$ .文献[13]给出了新估计式:设M=(mij)∈Rn×n是弱链对角占优B-矩阵,可表示为M=B++C,其中B+=(bij)形如(1)式,则
其中
${\mathit{\tilde \beta }_\mathit{i}} = {b_{ii}} - \sum\limits_{j = i + 1}^n {\left| {{b_{\mathit{ij}}}} \right|} > 0$ 且$\prod\limits_{j = 1}^{i - 1} {\frac{{{b_{jj}}}}{{{{\mathit{\tilde \beta }}_\mathit{j}}}}} = 1$ .下面,给出弱链对角占优B-矩阵线性互补问题解的新误差界.
定理1 设M=(mij)∈Rn×n是弱链对角占优B-矩阵,可表示为M=B++C,其中B+=(bij)形如(1)式,则
其中
${\mathit{\hat \beta }_\mathit{i}} = {b_{ii}} - \sum\limits_{j = i + 1}^n {\left| {{b_{ij}}} \right|} {h_i}\left({{\mathit{\boldsymbol{B}}^ + }} \right)$ ,且$\prod\limits_{j = 1}^{i - 1} {\left({1 + \frac{1}{{{{\mathit{\hat \beta }}_\mathit{j}}}}\sum\limits_{k = j + 1}^n {\left| {{b_{jk}}} \right|} } \right) = 1} $ .证 令MD=I-D+DM,则
其中BD+=I-D+DB+,CD=DC.类似于文献[13]中定理2的证明,可得BD+是弱链对角占优M-矩阵,主对角元为正,且
由引理1知
由引理2知,对任意的i,j,k∈
$\mathbb{N}$ +,有:且
进而,由引理2知
由(6)式和(7)式知
因此,由(6)式及(9)式知(5)式成立.
推论1 设M=(mij)∈Rn×n是B-矩阵,可表示为M=B++C,其中B+=(bij)形如(1)式,则
其中
${\mathit{\hat \beta }_i}$ 如定理1所定义.接下来,对(4)式及(5)式进行比较.
定理2 设M=(mij)∈Rn×n是弱链对角占优B-矩阵,可表示为M=B++C,其中B+=(bij)形如(1)式,且hi(B+)≤li(B+)(∀i∈
$\mathbb{N}$ +),则其中βi,
${\mathit{\tilde \beta }_\mathit{i}}$ 及${\mathit{\hat \beta }_i}$ 分别如(3)式、(4)式及(5)式所定义.证 因B+是弱链对角占优的,且主对角元为正,故对任意的i∈
$\mathbb{N}$ +,有由(11)式可知,对任意的i∈
$\mathbb{N}$ +,有由(12)式,可得对任意的i∈
$\mathbb{N}$ +,有且对任意的1≤j≤n-1,有
由(13)式及(14)式,可得
另外,注意到:
且lki(B+D)≤lk(B+) < 1.因此,对任意的i∈
$\mathbb{N}$ +,有${\mathit{\tilde \beta }_\mathit{i}} \le {\mathit{\bar \beta }_\mathit{i}}$ ,且同时,对任意的1≤j≤n-1,有
由(16)式及(17)式,可得
综上所述,由(15)式及(18)式知结论成立.