仲從磊，李靈曉

(河南科技大學 數(shù)學與統(tǒng)計學院，河南 洛陽 471023)

M矩陣Fan乘積行列式下界估計

仲從磊，李靈曉

(河南科技大學 數(shù)學與統(tǒng)計學院，河南 洛陽 471023)

針對M矩陣Fan乘積行列式的計算問題，結合M矩陣Hadamard乘積行列式的估計式，給出了M矩陣Fan乘積行列式下界的幾個估計。對這些結果進行了比較和推廣，并給出了相應的算例。研究結果表明：采用矩陣的順序主子式替代矩陣元素后行列式多一正項，下界變大，結論更精確，有利于M矩陣Fan乘積行列式的計算。

M矩陣;Hadamard乘積;Fan 乘積;行列式;下界

0 引言

矩陣的Fan乘積是一種特殊的矩陣乘積，廣泛應用于特征函數(shù)以及偏微分方程[1-4]的弱極小原理等方面的研究中。M矩陣是計算數(shù)學中重要的特殊矩陣類[5]，有著廣泛的應用。生物、物理和社會科學中的諸多問題都與M矩陣有著密不可分的聯(lián)系。因此，M矩陣Fan乘積行列式的計算和估計[6]成為特殊矩陣論研究的主要內容之一。其中，關于M矩陣Fan乘積行列式下界的估計是近幾年的研究熱點，并得到了一系列結論。但這些結論中的估計式需要借助于矩陣的所有元素進行計算，所需計算量較大且所得下界的誤差較大。

本文結合M矩陣Hadamard乘積行列式的估計給出M矩陣Fan乘積的一些估計式，并對這些估計式做了比較。對于M矩陣Fan乘積的計算和特征值界的估計問題[6-7]，以及相關的偏微分方程數(shù)值解等有一定參考價值。

1 矩陣Fan 乘積定義

設Rm×n為全體m×n實矩陣的集合，矩陣A=(aij)∈Rm×n，矩陣B=(bij)∈Rm×n,其Hadamard乘積定義為一個m×n矩陣A°B，其中(A°B)ij=aijbij。A和B的Fan乘積定義為一個m×n矩陣A*B，其中：

如果aij≥bij對所有的i,j成立，則稱A≥B。矩陣A∈Rn×n，如果A=sI-B且s>0,B≥0，s>ρ(B),其中,s為某一實數(shù)，I為n階單位矩陣，ρ(B)為B的譜半徑；稱A為非奇異M矩陣，記Mn為全體n×n非奇異M矩陣的集合。 矩陣A∈Rn×n的比較矩陣定義為一個n×n矩陣μ(A)=(mij)，其中:

矩陣A∈Rn×n的全體k×k順序主子矩陣記為Ak。

文獻[6-10]給出了一些關于兩個M矩陣Hadamard乘積行列式下界的估計不等式，經典的結果是：

定理1[6]若A=(aij)∈Mn,B=(bij)∈Mn,則

(1)

本文將上述結論推廣到兩個M矩陣Fan乘積上,并利用文獻[5]的相關結論對其改進。

2 幾個引理

引理1[9]如果A∈Mn,那么Ak∈Mk。

引理2[9]如果A=(aij)∈Mn,那么

det (Ak)≥0,k=1,…,n。

(2)

引理3[10]如果A∈Mn，B∈Mn，那么

(3)

證明 根據(jù)引理1,

因此,

diag(a11,a22,…,ak-1,k-1)≥Ak-1；

并且

從而,

沖砂系統(tǒng)的開挖與常規(guī)的地下隧洞開挖有區(qū)別的，坡積體段管棚灌漿處理形成了封閉裂隙，加強基巖的完整性，達到提高巖體強度和剛度的。在開挖過程中明顯出現(xiàn)砂礫石形成整體現(xiàn)象，加上臨時支護及時，未出現(xiàn)大量塌方現(xiàn)象，說明該隧洞的坡積體段施工方法可靠。

(4)

(5)

由式(2)、式(4)和式(5),可得：

引理4[10]若A=(aij)∈Mn,B=(bij)∈Mn,則

3 主要結論

本部分將給出兩個M矩陣A和B的det(A*B)界的估計不等式。

定理2 設A=(aij)∈Mn,B=(bij)∈Mn, 那么

A*B=μ(A°B)。

證明 因為A∈Mn、B∈Mn,易知aii>0,bii>0且aij≤0,bij≤0(i≠j)，從而

故

A*B=μ(A°B)。

根據(jù)定理2和式(1)可得：

定理4 設A=(aij)∈Mn、B=(bij)∈Mn,那么

證明 由引理4可得：

根據(jù)引理3和定理2,易得：

det (A*B) =det [μ(A°B)] ≥

以上結論通過算例加以驗證。

算例：

易知A和B都是M矩陣，且

det(A*B)=1 120；

顯然

det (A*B)= det [μ(A°B) ] ≥

上述算例表明定理4的結論比定理3的結論更加精確。

需要注意的是，矩陣的階數(shù)每提高一階，行列式的順序余子式的個數(shù)會增加一個，估計算法計算量也會相應地增加，這個問題可以借助于計算機得到快速解決。對于偏微分方程數(shù)值解等問題中產生的稀疏矩陣，如果滿足M矩陣的條件，尤其是對角占優(yōu)矩陣，本文中的估計不等式也是成立的。文中的推導方法依據(jù)M矩陣Hadamard 乘積行列式的相關結論得出，今后也會有更加精確的估計式產生，同時也可以把相關的結論推廣到兩個H矩陣[11]的Fan乘積問題上。

[1] 王天軍,殷政偉.Legendre-Gauss-Lobatto節(jié)點的一個注記[J].河南科技大學學報(自然科學版),2012,33(1):71-74.

[2] 王天軍,賈麗蕊.非線性熱傳導方程的 Lagrange 插值逼近[J].河南科技大學學報(自然科學版),2011,32(2):67-71.

[3] 王天軍,李清,殷艷紅.非線性熱傳導方程混合問題插值逼近[J].航空計算技術，2012,42(2):1-3，12.

[4] 陶冬亞,焦琳,王天軍.非線性Klein-Gordon方程的廣義Hermite譜方法[J].河南科技大學學報(自然科學版),2015,35(6):87-90.

[5] 黃廷祝，楊傳勝.特殊矩陣分析及應用[M].北京：科學出版社，2009.

[6] 周平,李耀堂.非負矩陣Hadamard積和M矩陣Fan積的特征值界的估計[J].純粹數(shù)學與應用數(shù)學,2012,28(6):826-833.

[7]王峰.非奇異M矩陣的逆矩陣和M矩陣的Hadamard積的最小特征值下界估計[J].應用數(shù)學,2013,26(2):341-345.

[8]周平,高美平,李艷艷.M矩陣Hadamard積的新估計[J]．河北北方學院學報(自然科學版),2016,32(1):1-7.

[9]LIU J Z,ZHU L.Some improvement of Oppenheim’s inequality forM-matrices[J].SIAM journal on matrix analysis and applications,1997,18(2):305-311.

[10] YANG Z P,LIU J Z.Some results on Oppenheim’s inequalities forM-matrices[J].SIAM journal on matrix analysis and applications,2000,21(3):904-912.

[11]LI Y T,ZHONG C L.Some estimations for determinant of the Hadamard product ofH-matrices[J].Journal of computational mathematics,2005,23:401-407.

國家自然科學基金項目(11301152)；河南省教育廳自然科學基金項目(13B110999)

仲從磊(1979-)，男，山東濟寧人，助教，碩士，主要研究方向為數(shù)值代數(shù)與特殊矩陣論.

2017-01-06

1672-6871(2017)06-0082-04

10.15926/j.cnki.issn1672-6871.2017.06.016

O241.6

A