M矩阵判定定理及证明

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

M 矩阵的性质、判定定理及证明

一、M 矩阵背景介绍:

1、M 矩阵是计算数学中应用极其广泛的矩阵类 。M 矩阵是L 矩阵的一种,M 矩阵要求它自身的逆矩阵为一个非负矩阵。

2、首先,L 矩阵的定义为:若A 一个n*n 的方阵,若0>ii a 而

≤ij a (i ≠j),则称A 为L 矩阵。

3、关于M-矩阵的一篇最早的论文发表于1887年,Stieltje 证明了一个具有非正非对角元的,非奇异对称对角占优矩阵的逆是一个非负矩阵。之后,1937年Ostrowski 提出M 矩阵的定义为:具有非正非对角元,且逆是非负矩阵。近年来,国内外的许多数学工作者对M 矩阵判定方法的研究都极为重视,并开展了深入的研究工作,给出了许多判定方法。但就目前的研究成果来看,所提出的M 矩阵的判定方法仅是、且仅能对M 矩阵作整体判定,这对高阶矩阵来说,在计算上较为困难,判定方法难以实现,因而现有M 矩阵的判定方法存在着相当大的局限性。

二、M 矩阵的概念

定义1 设n n ij a A ⨯=)(,且0≤ij a ,j i ≠,01≥-A ,称A 为M 矩阵。 定义2 设n n ij a A ⨯=)(,且0≥ij a ,若1-A 为M 矩阵,则称A 为逆M

矩阵。

引理1 如果n n ij a A ⨯=)(,且0≤ij a ,j i ≠,A 为M 矩阵的充要条件是A 可做三角分解,R L A •=,其中L 为下三角阵,R 为上三角阵,L 和R 的主对角元都是正值。

三、M 矩阵的判定定理与证明

定理1 若n n ij a A ⨯=)(为M 矩阵,则R L A ⨯=,其中下三角阵L 和上三角阵R 的主对角线元素为正,且其余元素为非正值。

证明 若A 为M 阵,则当j i ≠,0≤ij a ;j i =,0>ij a 。由引理1,A 可做三角分解R L A •=。设

⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=nn n n l l l l l l L 21222111000 , ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣

⎡=nn n n r r r r r r R 00022211211 则⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+++++=nn nn n n n n n n n r l r l r l r l r l l r l r l r l r l r l r l r l A 1122

21211112212122221221112111112111111, 故0,,1111211≤n r l r l 。

因011>l ,故0,,112≤n r r ;因,0,0,,111111121>≤r r l r l n 故0,,121≤n r r ;因

022321231≤+r l r l ,故02221≤r l ,从而021≤l ;因023221321≤+r l r l ,故023≤r 。类

似的有02≤i r ,02≤i l (n i ,,5,4 =)。又因有0343324321421≤++r l r l r l 及

0334323421341≤++r l r l r l 故相应有014≤r ,043≤l 。类似的有03≤i r ,0

3≤i l (n i ,,6,5 =)。

假设k n =时有0≤ik l ,0≤ki r ,(n k i ,,1 +=),当1+=k n 时,由于

02,11,12,,12,22,12,11,1≤++++++++++++++k k k k k k k k k k k k r l r l r l r l ,故02,1≤++k k r 。又由于

01,11,21,11,2≤++++++++k k k k k k r l r l ,故01,2≤++k k l ;类似的可得到0,1≤+i k r ,0

1,≤+k i l (n k i ,,2 +=)。证毕。

定理2 设n n ij a A ⨯=)(,ij a 的代数余子式为ij A ,n j i ,,2,1, =,如果

,,0j i a ij ≠≥则1-A 为M 矩阵的充要条件0,0≤>ij ii A A 。

证明 必要性:如果1-A 为

M 矩阵,由于

))(,0(11n n ij A A A d A d

A ⨯**

-=≠==

,故0,0≤>ij ii A A )(j i ≠。 充分性:由于*-=A d

A 1

1,且0

,0≤>ij ii A A ,

),,2,1,(,0,)(n j i a a A ij n n ij =≥=⨯,就由定义1知1-A 为M 矩阵,证毕。

定义3 设有n 阶矩阵n n ij a A ⨯=)(,如果存在正向量X (即它的分量i x 都是正值),使得),,2,1(n i x a x a i

j j ij i ii =>∑≠成立,则称A 为拟对角

占优。

引理2 设n n ij a A ⨯=)(,满足)(0,0j i a a ij ii ≠≤>,并且矩阵T A A B +=为拟对角占优,则A 为M 矩阵。

定理 3 设n n ij a A ⨯=)(,如果k i n k i a a k i kk ii ≠=ΛΛ>,,,2,1,,4

1

则A 为M 矩阵(其中∑≠+=Λ=>≠≤i

j ji ij i ii ij a a j i a j i a ,0,,0)。

证明 若i ii a Λ>2

1

对n i ,,2,1 = 皆成立,则由定义3 知T

A A

B +=为拟对角占优。由引理2知A 为M 矩阵,为此,只需证明对某个i 有

i ii a Λ≤

21的情形。不失一般性,不妨设12

1

Λ≤ii a 。由 k i kk ii a a ΛΛ>41

,可得

,,,3,2,21n k a k kk =Λ>用111/)2

1

(a Λ乘以矩阵B 的第一列,得新矩阵

)()

1()1(ij b B =,则有

1)

1(11Λ=b ,

相关文档
最新文档