趙 英 霞 王 峰
(貴州民族大學(xué)數(shù)據(jù)科學(xué)與信息工程學(xué)院,貴州 貴陽 550025)
M-矩陣是應(yīng)用廣泛的一類重要矩陣,與許多學(xué)科領(lǐng)域有著密切的關(guān)聯(lián),經(jīng)濟(jì)學(xué)、生物學(xué)及社會(huì)科學(xué)中有大量問題和M-矩陣有著必不可少的聯(lián)系. 且嚴(yán)格對(duì)角占優(yōu)M-矩陣A的逆矩陣的‖A-1‖∞上界估計(jì)在數(shù)值代數(shù)中有著重要的應(yīng)用,在代數(shù)方程組的收斂性條件及條件數(shù)需計(jì)算‖A-1‖∞,可當(dāng)M-矩陣的階較大時(shí),其逆矩陣求解復(fù)雜,因此對(duì)‖A-1‖∞上界進(jìn)行估計(jì)是十分重要的.近年來,許多學(xué)者對(duì)M-矩陣的逆矩陣的無窮大范數(shù)上界進(jìn)行了估計(jì),得到了很多好的結(jié)果[1-12].另外,特殊結(jié)構(gòu)矩陣線性互補(bǔ)問題是一類在工程學(xué)、經(jīng)濟(jì)學(xué)、控制論等領(lǐng)域具有重要應(yīng)用價(jià)值的優(yōu)化問題[13-16].對(duì)于給定的n階矩陣A=(aij)∈Rn×n,n維實(shí)向量q∈Rn,用LCP(A,q)來表示矩陣A的線性互補(bǔ)問題,尋找x∈Rn,使其滿足
為敘述方便,給出一些記號(hào).設(shè)A=(aij)∈Rn×n,aii≠0,m≤i,j,k≤n,ε>0,且
定義1[2]設(shè)A=(aij)∈Rn×n,如果對(duì)任意的j,i∈N,i≠j,都有aij≤0,則稱A為Z-矩陣,記A∈Zn.設(shè)A∈Zn,則A可表示為A=sI-B,其中B≥0.當(dāng)s≥ρ(B)時(shí),稱A為M-矩陣;當(dāng)s>ρ(B)時(shí),稱A為非奇異M-矩陣.
定義2[3]設(shè)A=(aij)∈Rn×n,如果滿足下面條件
(c)對(duì)于任意i∈N,i?J(A),存在i1,i2,…,ik,使得aii1ai1i2,...,aik-1ik≠0,ik∈J(A).則稱A為弱鏈對(duì)角占優(yōu)矩陣.
定義3[3]設(shè)A=(aij)∈Rn×n,若J(A)=N,則稱A為嚴(yán)格對(duì)角占優(yōu)矩陣.
定義4[25]設(shè)A=(aij)∈Rn×n,若
則稱A為B-矩陣.
首先給出一些引理.
故定理4改進(jìn)了文獻(xiàn)[8]中的定理3.4,進(jìn)而優(yōu)于文獻(xiàn)[2]中的定理3.3和文獻(xiàn)[7]的定理3.2.
2009年García-Esnaola等[13]給出如下結(jié)果:設(shè)A=(aij)∈Rn×n為B-矩陣,將A表示為A=B++C的形式,其中
綜合上述知(14)式成立.
下面對(duì)估計(jì)式(11)式與(14)式進(jìn)行比較.
定理7 設(shè)A=(aij)∈Rn×n是B-矩陣,令A(yù)=B++C且B+=(bij)形如式(9),則
綜上可得(14)式優(yōu)于(11)式.
下面用數(shù)值例子說明新估計(jì)式比已有的一些結(jié)果更加精確.
例1 設(shè)
顯然A是嚴(yán)格對(duì)角占優(yōu)的M-矩陣.應(yīng)用文獻(xiàn)[2]中的定理3.3,文獻(xiàn)[7]中的定理3.4及文獻(xiàn)[8]中的定理3.2,分別得
表1 的上界
表1 的上界
n(11)式(12)式(13)式(14)式10 104.9084 90.8111 83.6275 78.5563 20 415.8320 361.6667 332.8539 306.0755 50 2.8859e+3 2.1573e+3 1.9584e+3 1.1876e+3 100 9.5301e+3 7.1042e+3 5.3078e+3 3.2824e+3
本文給出了嚴(yán)格對(duì)角占優(yōu)M-矩陣及逆矩陣之間的元素關(guān)系式,通過迭代法獲得了嚴(yán)格對(duì)角占優(yōu)M-矩陣的逆矩陣的無窮大范數(shù)的新上界.同時(shí),利用新上界與兩個(gè)重要不等式放縮技巧得出B-矩陣線性互補(bǔ)問題誤差界的新估計(jì)式,理論證明及數(shù)值算例表明了新估計(jì)式的有效性.