留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

弱链对角占优B矩阵线性互补问题的误差界新估计

上一篇

下一篇

徐玉梅, 王峰. 弱链对角占优B矩阵线性互补问题的误差界新估计[J]. 西南师范大学学报(自然科学版), 2019, 44(2): 18-24. doi: 10.13718/j.cnki.xsxb.2019.02.005
引用本文: 徐玉梅, 王峰. 弱链对角占优B矩阵线性互补问题的误差界新估计[J]. 西南师范大学学报(自然科学版), 2019, 44(2): 18-24. doi: 10.13718/j.cnki.xsxb.2019.02.005
Yu-mei XU, Feng WANG. New Error Bounds for Linear Complementarity Problems of Weakly Chained Diagonally Dominant Matrices[J]. Journal of Southwest China Normal University(Natural Science Edition), 2019, 44(2): 18-24. doi: 10.13718/j.cnki.xsxb.2019.02.005
Citation: Yu-mei XU, Feng WANG. New Error Bounds for Linear Complementarity Problems of Weakly Chained Diagonally Dominant Matrices[J]. Journal of Southwest China Normal University(Natural Science Edition), 2019, 44(2): 18-24. doi: 10.13718/j.cnki.xsxb.2019.02.005

弱链对角占优B矩阵线性互补问题的误差界新估计

  • 基金项目: 国家自然科学基金项目(11601473);贵州省科学技术基金项目(20181079);贵州省教育厅自然科学基金项目(2015420)
详细信息
    作者简介:

    徐玉梅(1992-), 女, 硕士研究生, 主要从事数值计算和数学模型的研究 .

    通讯作者: 王峰, 博士, 教授
  • 中图分类号: O151.21

New Error Bounds for Linear Complementarity Problems of Weakly Chained Diagonally Dominant Matrices

计量
  • 文章访问数:  827
  • HTML全文浏览数:  548
  • PDF下载数:  55
  • 施引文献:  0
出版历程
  • 收稿日期:  2018-03-05
  • 刊出日期:  2019-02-20

弱链对角占优B矩阵线性互补问题的误差界新估计

    通讯作者: 王峰, 博士, 教授
    作者简介: 徐玉梅(1992-), 女, 硕士研究生, 主要从事数值计算和数学模型的研究
  • 贵州民族大学 数据科学与信息工程学院, 贵阳 550025
基金项目:  国家自然科学基金项目(11601473);贵州省科学技术基金项目(20181079);贵州省教育厅自然科学基金项目(2015420)

摘要: 利用弱链对角占优M-矩阵的逆矩阵的无穷大范数的范围,结合不等式放缩技术,给出了弱链对角占优B-矩阵线性互补问题误差界新的估计式,进而给出B-矩阵线性互补问题误差界的估计式.新估计式改进了已有文献的结果.

English Abstract

  • 线性互补问题在众多领域具有广泛的应用,如数学规划、最优停止、期权定价、市场均衡、自由边界和弹性接触等[1-3].线性互补问题的数学模型为:求zRn,满足:

    记作LCP(Aq),其中A=(aij)∈Rn×nqRn为给定的矩阵和向量.

    线性互补问题不论是解的存在性、唯一性,还是算法的收敛性,都与矩阵A的结构有着密切的关系.如当矩阵AP-矩阵(即它的所有主子式都是正的[4])时,线性互补问题有唯一解[2].文献[5]给出了如下结果:设AP-矩阵,则

    其中x*LCP(Aq)的解,r(x)=min{xAx+q},D=diag(d1,…,dn)(0≤di≤1).近年来,众多学者对特殊结构矩阵线性互补问题解的误差界、扰动界进行了分析,给出了一系列估计式[6-13].本文在文献[14]的基础上,继续讨论弱链对角占优B-矩阵线性互补问题解的误差界,获得了新估计式.同时,证明了新估计式在一定条件下优于文献[12-13]的结果.最后,用数值算例表明了所得结果的有效性.

  • $\mathbb{N}$+表示全体正整数的集合.记A=(aij)∈Rn×n,若

    则称A为严格对角占优的[15];若aij < 0(ij),则称AZ-矩阵[4];若AZ-矩阵且A-1≥0,则称AM-矩阵[15];若|aii|≥ri(A)(∀i$\mathbb{N}$+),且对每个iJ(A)={i$\mathbb{N}$+:|aii|>ri(A)}≠Ø,存在非零元素链aii1ai1i2,…,airj,满足jJ(A),则称A为弱链对角占优矩阵[15].

    定义1[4]  设A=(aij)∈Rn×n,若对任意的ij$\mathbb{N}$+,有:

    则称AB-矩阵.

    定义2[13]  设M=(mij)∈Rn×n,若M可表示为M=B++C,其中:

    ri+=max{0,mijji},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×nijk$\mathbb{N}$+iji+1≤kn,记:

  • 文献[9]给出结果:设M=(mij)∈Rn×nB-矩阵,可表示为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×nB-矩阵,可表示为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知,对任意的ijk$\mathbb{N}$+,有:

    进而,由引理2知

    由(6)式和(7)式知

    因此,由(6)式及(9)式知(5)式成立.

    推论1  设M=(mij)∈Rn×nB-矩阵,可表示为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≤jn-1,有

    由(13)式及(14)式,可得

    另外,注意到:

    lki(B+D)≤lk(B+) < 1.因此,对任意的i$\mathbb{N}$+,有${\mathit{\tilde \beta }_\mathit{i}} \le {\mathit{\bar \beta }_\mathit{i}}$,且

    同时,对任意的1≤jn-1,有

    由(16)式及(17)式,可得

    综上所述,由(15)式及(18)式知结论成立.

  • 例1  考虑B-矩阵[12]

    Mk=Bk++Ck,其中

    由(2)式得

    易见30(k+1)→+∞(k→+∞).

    由(3)式得

    由(4)式得

    由(5)式得

    可见(5)式优于(2)-(4)式.

    例2  考虑弱链对角占优B-矩阵[13]

    M=B++C,其中

    由(4)式得

    由(5)式得

    可见(5)式优于(4)式.

  • 文中讨论了弱链对角占优B-矩阵线性互补问题解的误差界估计问题,得到了几个新的上界估计式,理论证明及数值算例均表明了新估计式在一定条件下优于文献[9, 12-13]的结果.

参考文献 (15)

目录

/

返回文章
返回