矩阵行列式和代数多项式根的计算问题
电子科技大学
硕士学位论文
矩阵行列式和代数多项式根的计算问题
姓名:刘伟
申请学位级别:硕士
专业:计算数学
指导教师:蒋泽云
20080501
摘要
摘要
矩阵行列式和代数多项式根的计算问题,实际上是复杂而又很经典的数学问题之一,很早人们就对其进行了研究。因此,对其进行研究具有很高的理论和应用价值。本文结合最近出现的一些相关结论,对此类问题,在进行深入研究的基础之上,给出了几种简单估计。所得结果推广或改进了一些经典结果。
本文的研究主要分为两大部分:
1.主对角占优矩阵类的行列式估计问题:主要证明了当满足一定条件的时候,通过研究逆元素的估计,获得了一些有趣的主对角占优矩阵类行列式的上下界的几个简单估计,如:
(1)对一般对角占优矩阵,我们获得了如下几个结果:
定理2.1如果A=‰】eD,那么对任意指标k∈N,有
]detAI习akk兀(]ajji-c,Tj)兀(1%一哆‘),
其中哆=m叫ax{v,。},哆=龄{%)和%=鱼篙≠,J∈Ⅳ.
rj=o(.jf=l,…,k—1)和0=o(jf=七+1,…,嚣)。
IdetAI-qakk
..,并且等号成立当且仅当定理2.2如果A=【%】∈D,那么对任意指标keN,有ILI(1a∥l+g芬)兀(Ia//I+哆‘),歹2I歹=七+1
其中D,,访,同定理2.1(a)。并且等号成立当且仅当rj=0(歹=l,...,七一1)和‘=0(歹=k+l,..卅)。
定理2.3如果A=‰】∈D,那么A~=嗡】,并且对任意后∈M有
I口舫I兀(I%I-∑I%f%)n(I口ffl_∑I乃1w,j)---IdetA]<
la胜l兀(J口lfI+∑Ia,j1%)兀(1%l+∑Ia,j
其中%和嘞(i,jef%),N)同定理2.1(a).并且上式等号成立当且仅当i<k,all=o(/>D和,l≥f>|i},口ff=oU<O.
摘要
(2)对双对角占优矩阵,获得了如F结论:
定理2.4如果A=[av]eDD,且Ⅳo=UIla口I<Rj(A),jeN}≠o,那么
(I%I-YJa材I乃)…如[-<ldctAI-<(1a欣I+∑IalFIpj)IdetA灶I,keNo.
2.关于经典代数多项式的根模估计问题,获得了一系列有趣的结果,推广或改进了一些经典的相关结论。如:
定理3.1对多项式(3.2),其任意实根兄将满足:…鲕ax{・,F啊,F硒’,
lA其中口=min{o,一ajI2≤/≤刀},b=min{0,q2≤jf≤刀}。定理3.3对多项式(3.2),以及任意1≤S≤n,o->0,其任意根旯将满足:(1)I<max{Pl,P2),
其中届=黔…小…竹,pz<--…m—ax.n掣,掣}.
(2)
I兄阵max{・,三[IqI+√iI面],去[Iq¨l+√:二面]}.关键词:行列式,主对角占优矩阵类,代数多项式,根Ⅱ
ABSTRACT
Theproblems・----——thedetermentof
oneamarxandthepolynomialinoneindeterminateevaluation・——・——areinfactofcomplexandclassical
onMathematicalaproblems.Manyfamousmathematiciansdolotsofinvestigationthemsincelong
longago.Therefore,itisveryimportantandusefultoestimate
onthem.Inthispaper,wegivelotsofevaluationsforthese
resultsgeneralizeproblems,basingones.somenewinvestigations.Ourandmodifysomeclassical
Thispapermainlyincludestwoparts:
1.Boundsfordeterminantswitll
newupperdominantprincipaldiagonal:Inthispart,someandlowerbounds
arefordeterminantswiⅡldominantprincipaldiagonalarepresented.Theseboundssome
firstlyimprovementsofsomegivethat
allsimilarbounds.Forexample:(1)InChaptertwo,weTheorem2.1.IfA=略]eD,then,for
k-Iarbitraryindexk∈N,Wehave
-DTJ)兀(1a口一哆‘),detAl4a船兀(I%I
j=lj=k+1
where111e唧1ali够D,:x{vA,仍f:燃{wⅣ)册JwⅣ:冬生牟,_,∈Ⅳ.11le唧lali够holdsg2。,哆2龄{~)册J坳22盖f’,∈Mif髓donlyifrj=o(/=l,…,J|}一1)and‘=o(_,=七+l,…,以).
Theorem2.2.IfA=【aF】∈D,then,for
k-Iallarbitraryindexk∈Ⅳ,^
detAl4aaIH(1a∥I+t5fj)兀(1%I+谚/A,
Theorem2.1.Theequalitiesholdifandonlyifwhere巧and髟(,∈N)arethesame船thosein
吩=o(,=l,…,J|}一1)and‘=o(/=七+1,…,咒)・
Theorem2.3.If
k-IA一[吩]eD,thenA一=[bo.]exists,and,foranarbitraryindexk∈Ⅳ,“
口肚皿(1a。I-Y,lauI%)I-I(1auI-∑}嘞I嘞)sldctAI<_
l=lj>ii=k+lksj<i
k-1月
口船I兀(1口fjI+∑l%I%)兀(I魄l+∑laF
i=lI嘞),j>ii=k+lIs,《fIII
where%and~(i,ja忉asinTheorem2.1.Theequalitiesholdifandonlywheni<k,a/j=0(pf)andwhenn≥i>k,嘞20(歹<f).
alsoshowthat(2)Fordoubly(row)diagonallydominant,we
Theorem2.4.If
(IaaA=【嘞】∈DD,andNo={川a//峰Rj(A),.,∈册≠0,thenl-∑1%Ipj)tdetA艟I<_1detAl<(Ia娃I+∑j%I&)ldetA特j,七∈Ⅳ0.
,靠户I
inone2.Evaluationforthepolynomialindeterminate:Wegivesomesimple
bounds,whichgeneralizeandmodifysomeclassicalones.Forexample:
Theorem3.1.Forthepolynomial(3-2),itsanyrealeigenvalue允satisfiesthat:
那一{-,
where口=min0,一aj2≤/≤刀),b=min0,哆12_<j-<玎).
of1≤s≤n,盯>0,itsTheorem3.3.Forthepolynomial(3—2)andanyany
eigenvalue五satisfiesthat:
(1)
五|<maX{岛,岛),
砒ere亿m华a州x{1+I吼…+仃)∥州m华ax川{-+掣,料
(动
兄降m觚L1—2q+
一¨一
Keywords:determent,dominantprincipal一H小
IVdiagonalmatrix,polynomial,value
独创性声明
本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。据我所知,除了文中特别加以标注和致谢的地方外,论文中不包含其他人已经发表或撰写过的研究成果,也不包含为获得电子科技大学或其它教育机构的学位或证书而使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。
签名:.一生堑笙——日期:≯。。6年厂月p日
关于论文使用授权的说明
本学位论文作者完全了解电子科技大学有关保留、使用学位论文的规定,有权保留并向国家有关部门或机构送交论文的复印件和磁盘,允许论文被查阅和借阅。本人授权电子科技大学可以将学位论文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存、汇编学位论文。
(保密的学位论文在解密后应遵守此规定)
签名:导师签名:蝎缉么导师签名:筠主聋亟
日期:≯矿略年厂月f,)日
第一章引言
第一章引言
1.1选题背景
随着科学技术的迅猛发展,特别是计算机的广泛应用,矩阵理论越来越受到数学工作者,科技和工程人员的重视。它不仅是一门重要的数学分支,而且在数值分析,最优化方法,微分方程,控制理论,数学模型等分支及各种工程科学领域有者极其重要的应用。
矩阵行列式和代数多项式根的计算问题,实际上是复杂的而又很经典的数学问题之一,很早人们就对其进行了研究。因此,对其进行研究具有很高的理论和应用价值。如在线性代数一些问题研究中,如线性方程组、矩阵秩和特征值等问题,常要用到行列式作为工具。
关于矩阵行列式的计算问题,自从矩阵出现以来就是人们研究的热点之一。至今在几乎所有的矩阵或线性代数教材中对其皆有非常详细的介绍【1‘51。有关它的相关理论也几乎家喻户晓,如二、三阶行列式的计算,拉普拉斯展开定理以及求解线性方程组的克莱默法则等等。然而,对高阶行列式计算来说,虽然目前随着高性能计算机的发展有所改善,但对于某些较病态的行列式来说,采用目前的计算机算法来计算,数值并不稳定。因此,目前对矩阵行列式的估计问题仍然是国内外研究的热点之一。
代数方程是最古老的方程,有着重要的应用。代数多项式根的计算问题,已经有着更加悠久的历史。比方说对一元二次多项式根的计算问题,很久以前就有了根的判别式和求根公式(韦达定理)。到了16世纪,意大利一个靠自学成才的数学家塔尔塔利亚发现了一般一元三次方程的求根公式‘2,4】:
1.给定方程z3—1=0,则其三个根为:
而:l,恐:国:—-1丁+以i,x3=c02-—-1丁-’vf3i.
且玉+恐+毛=1+∞+∞2=O,jcl・恐・x3=1・国・∞2=0)3=1.
2.给定方程≯+∥+bx+c=o.令x=y一:a,代入得Y3+PY+q=0.设其根为Yl,Y2,Y3,则有:
(1)根:炉丽+丽。
J,22电子科技大学硕士学位论文
式中国:_-1+x/3i.一:+一
Yl+此+乃:O,1+1+1:一旦,乃耽咒----q.
f>0,有一个实根与一对共轭复根.
A{=0,有三个实根,其中有两个相等.
I<0,有三个不相等的实根.
Y3一秒2+(耐一4力y—b2P+4钟一d2=0(3,笋n别式:△=(兰]2+(詈)3.路[2’41:给定方程p+W+“:+dx+c=0.先求出方程
的任一实根Y。,当by.一2c>0时,再解下列两个方程:
,+圭(6±乒i石"芝1帆±-4c)=o.
当砒一2c<0时,再解下列两个方程:
x2+三(6±4—b2-4c—+4yo)z+(Yo干√;F二石)=o.
随后在他们工作的鼓舞下,人们开始继续寻求更高次数的多项式的求根公式。然而,三百多年过去了,仍没有太大的进展。1824年,天才的挪威青年数学家阿
定理1.1(阿贝尔定理)四次以上的文字系数的代数方程,没有一般的由方程2
第一章引言
的系数经有限次四则运算和开方运算求根的方法1。
随后,法国数学家伽罗瓦(EvaristeGalois,1811.1832)在其20岁时,彻底解决了这一问题一给出了一般n次方程能用根式解出的充分和必要条件。其研究论文“关于用根式解方程的可解性条件",于1846年发表在“数学杂志"上。其发表后,并没被人们所理解,直到1870年法国数学家若尔当详细解释了这篇论文后,才被人们所广泛接受。其解决这一问题所用的工具“伽罗瓦群”,现己被公认为19世纪数学的最杰出的成就之一。
然而,由代数基本定理知道,每个次数不小于1的复系数多项式,在复数域上至少有一个根。然而要把这些根全部求解出来,一般来说还是很难的。目前,一般都采用代数方程的牛顿算法【4】:
给定多项式
f(x)=aox”+axx4_1+…+an_iX+a。.(1-1)
式中的q(f=0,1,..卅)均为实数。
利用迭代公式:
‰2而一器舷-l,2,…)・
计算(1.1)的根。
进行求解。(1-2)显然,这需计算厂(讫)及厂仇)。对此一般地采用多项式求值的秦九韶算法,来
然而,正如上面所看到的,牛顿算法对初值区间的选择要求是很高的,在某些情况下,对初值也十分敏感。因此,若能给出初值区间的一个较好的选择,必将对牛顿算法十分有益。
1.2对角占优矩阵与友矩阵简介
由前所述,对矩阵行列式和代数多项式根的估计问题,有着重要而广泛的应用,对其进行研究具有很高的理论和应用价值。在本论文中,我们将针对某些有着广泛应用背景的矩阵类一一般对角优势矩阵类和与代数多项式求根密切相关的友矩阵进行展开研究。在此我们先对本文所涉及的一些矩阵类和基本概念,作以简单介绍。
1注:一般的一元n次方程的根,可以用福克斯函数表示。五次方程的根可用椭圆函数表示。3
电子科技大学硕士学位论文
众所周知,许多实际问题最后常常归结为一个或一些大型系数矩阵为特殊矩阵的线性方程组的求解问题【3】,如在控制论及神经网络大系统的稳定性、线性时滞系统的稳定性研究中常常需要H.矩阵和正稳定矩阵;在生物学、物理学、数学和社会科学、偏微分方程中的有限元方法、经济学中的投入一产出分析和增长模型、运筹学中的线性互补问题、概率统计中的马尔科夫过程中需要的M.矩阵等。这些特殊矩阵在矩阵分析和科学计算中占有十分重要的地位,它们在计算数学、经济学、物理学、生物学、应用数学等方面都有着广泛的应用。本文所涉及的一般对角优势矩阵类就是与这些特殊矩阵密切相关的一类矩阵。
为了叙述上的方便,让4是一个nXn阶复矩阵,并记作A=【aO】∈r”。定义N:{1,2,..棚),并且在此文中我们定义如下两个常用符号:
墨(彳)=∑Iaol,c;(4)=∑f%I,f∈Ⅳ.
j≠ij{i
定义1.1设彳=【口{『】∈C…,若IauP-置(么),i∈<n>,则称么为(行)对角占优矩阵,记为A∈Do。若上式中均为严格不等号,则称彳为(行)严格对角占优矩阵,记为A∈D。若存在正对角矩阵x,使B=似为严格对角占优矩阵,则称彳为广义严格对角占优矩阵【3,6】(也就是H.矩阵),记为A∈西。
定义1.2【7】设彳=‰】∈C~,若l嘞怕f序R,(A)Rj似),vi≠je<刀>,则称彳为(行)双对角占优矩阵,记为A∈DoDo。若上式中均为严格不等号,则称彳为(行)双严格对角占优矩阵,记为A∈DD。
显然,(行)对角占优矩阵一定是(行)双对角占优矩阵;(行)严格对角占优矩阵也一定是(行)双严格对角占优矩阵。但反之,不一定成立。然而,他们都是H.矩阵。在本文中,我们通称他们为主对角占优矩阵类。
另外,这些对角占优矩阵类也与下半强对角占优矩阵,半强对角占优矩阵,非零元素链对角占优矩阵和口一对角占优矩阵等有着天然和密切的联系,有兴趣的读者可参见文献【s】,在此不再详述。
行列式是矩阵理论中,最早出现的概念之一,一般定义为:
定义1.3Eun阶行列式4
第一章引言
口lla21
q2
a22
●
…口ln…口2^
●
●
%%~
.
=
●
‘‘’
ailafl
●
ain
∑(一1)7‘^,止’…・^’口l^口2矗…口峨.(1-3)
(^,J2,…,^)
%~
●
anlan2
~,
…%。
即取自不同行、不同列的n个元素乘积的代数和。
从行列式的定义可以看到,若直接采用(1.3)式来计算行列式的值计算工作量是
很大的。
下面我们介绍,一元n次多项式的根与矩阵特征值之间的关系。定义1.4【61对于多项式以(z)=r+alx剃+...+%一。x+an,称11阶矩阵
-aI
-a2
10l
…一口^一l…O…
0
-a^
O0
l0
…
’.
OO1O
0O01
-aI
.
00
q=
O
或
0
0
…1
(1-4)
OO
…
1O
一口n—口^一1。。。-a2
为多项式Pn(x)的友矩阵,利用数学归纳法可以证明:对于,z≥2,恒有
(一1)”pax)=det(Cv—x,)或Pn(X)=det(x/一c:).
(1—5)
其中J表示n阶单位矩阵。
从上面(1.5)式可知,多项式见@)=∥+alx”1+...+口州x+a。的根实际上就是其所对应的友矩阵的特征值。因此,目前关于多项式根的估计问题都是从其所对应的友矩阵的角度来研究的。另外,关于友矩阵的研究常常用到如下几个性质
性质1.1【6】每个首一多项式既是它的友矩阵的极小多项式又是其特征多项式。性质1.2【6】矩阵彳相似于特征多项式的友矩阵,当且仅当么的极小多项式与特征多项式恒等。
性质1.3【6】矩阵彳相似于特征多项式的友矩阵,当且仅当彳是非减次矩阵。
1.3本文主要工作
已有的文献已经表明,关于主对角占优矩阵行列式的估计,以及一元n次多项式根的估计问题,仍然是一个十分活跃的研究领域。目前已取得的很多优美的结果。在本文中,我们将在这些研究的基础上,进一步开拓思路,得出一些更好
5
电子科技大学硕士学位论文
的结果。我们的安排如下:
第二章主要介绍主对角占优矩阵行列式的估计问题,通过研究逆元素的估计,获得了一些有趣的主对角占优矩阵类行列式的上下界的几个简单估计。
第三章关于经典代数多项式的根模估计问题,提出一系列有趣的结果,推广或改进了一些经典的相关结论。
第四章展望未来,总结所得结果。
6
第二章主对角占优矩阵类的行列式的界
第二章
主对角占优矩阵类的行列式的界
2.1引言与记号
让彳是一个nXn阶复矩阵,记作彳=k】∈C““。定义Ⅳ={l,2,...,刀),并且在此章中我们定义如下几个常用符号:
州)=驯,cl(铘=渺№=等,c:f=筲'f∈Ⅳ.
(2.1)
众所周知,对矩阵彳来说,如果对任意i∈N,皆有届≤1,那么矩阵彳被说成是行对角占优矩阵,记作A∈或;若对任意i∈N,皆有岛<1,则被称为行严格对角占优矩阵,记作A∈D。同样地,对参数CI(f∈N)类似地可以定义列对角占优矩阵和列严格对角占优矩阵。显然,行对角占优矩阵的任意阶主子矩阵亦是行对角占
优的。
近几十年来,对主对角占优矩阵行列式的研究一直是矩阵代数与计算数学领域里的热门课题之一,无论是主对角占优的等价表征,判定条件,还是更深入的性质研究,都出现了一系列既有重要的理论价值,又在实际计算中有广泛用途的
重要成果。如:
对严格对角占优矩阵类一些优美的行列式不等式已被获得:首先,国际著名数值代数专家A.M.Ostrowski教授在文献‘9】中,首次给出了如下著名的不等式:
[dctAI≥兀[1a豇I一足(铆
(2-2)
这个不等式提供了一个简单的下界。随后,G.B.Price[101和Ocder_【111分别采用如下不
同的表达式
H
i-1
e
,;=∑hI,‘=∑l嘞I,ij≠t+l
1=1
N.
G.B.Price在假设彳∈D情况下,表吲1川):
兀[I
i=l
ait
h】习dotAl>n【|aiiI-r,].
i=1
(2—3)
Ocdcr在假设么∈域情况下,获得了【ll】:
7
电子科技大学硕士学位论文
IdetA巨]aa。I丌(hl_ff).(2-4)
然而,在文献【12】中,所有的上述结果又可以以这样的方式做进一步的改进,即对任意指标k∈N,有:
I%I讨(%I+pr;)卉(I口ff+P1,)->1detA辟1口肚I兀k-I(1%I—p,;)n(I
其中P=m衅{岛j.
ati
I一矾).(2-5)
受这些工作的鼓舞,下面我们将对此类问题作迸一步的研究。改进上述结果,并获得一些新的行列式上下界估计。
2.2一般对角占优矩阵的情形
首先,给出如F几个定义和引理:
定义2.1【13】矩阵B=【龟,】∈C…是非负矩阵(记为B≥0),若对任意口≥pp),其中p(B)表示矩阵B的谱半径,矩阵彳=k】∈C“”满足
A=al—B,口>0,
则称矩阵彳为M.矩阵,这里J表示单位矩阵。
引理2.1t71设矩阵么=‰】∈∥“,记4∈c‘¨脚哪是位置为(ij)处的余子式,
则对任意f∈N,有
de%.detA-Y.(一1)‘H力嘞de%,或detA-∑(一1)‘H力aij
产l
J21
(2-6)
引理2.2171若矩阵A=【%】∈C”,相对于M.矩阵B=【%】来说,若满足玩<1%I,
且对任意f≠,,有l口f,I-<1%I,那么
ldetA除detB≥0.(2-7)
引理2.3如果么=‰】∈D,那么A一1_【%】存在,并且(a)(见文献【14】中的引理2.2)
1%匿乃l玩I<1%l,对任意f≠歹・1%I<-p;1%H%I,对任意f≠.,.
(2.8)(2・9)
∽(见文献【15】中的引理4)
8
第二章主对角占优矩阵类的行列式的界
耶丽面1砺,对某~/“
下面为了叙述上的方便,现在对任意f∈Ⅳ,我们引入下列记号:
(2-10)
吩2融1咖j窆=i+lh旧2等//l妒掣t4
产
I“
lⅣl,
其中4.,(f<_,)表示矩阵彳的由集合{i,i+l….,j)所决定的主子矩阵。
定理2.1(a)如果彳=[aU]ED,那么对任意指标keN,有
IdetAI到a肚兀(I%I-g乃)兀(I%一嘭‘),
其中哆2
(2.11)
z学{%),哆=m…ax,{w.}和wj,=鱼篙半,_,∈Ⅳ.
并且等号成立当且仅当
乃=o(/=1,…,Jj}一1)和‘=o(,=Ji}+1,…,n)。
(b)若A7∈D,那么对任意k∈N,有
IdetAI>_la船II-[(I%I-乃乃)兀(I%I-r。/:u^
(2-12)
其中膨2蛩{心),彩=畏焉{旦高半),,∈Ⅳ-并且等号
嘭=o(,=1,...,七一1)和“,=o
U=k+l,...,,1)。
成立当且仅当
证明:我们的方法是简单的,注意到
因此,从引理2.30),我们获得对任意ieN有
%=警'f,_,∈M…=黜≤丽1,对某些删纠
I拙Al>I晋(1a“I-/≈(A)p:)ldet4,I,ieN.
(2・13)
即
IdetAI_>(J%I嵋(A)pj)Idet4-,-J,对某些,eN,,≠f.
既然对任意ieN,和某些歹eN,,≠f,上面的不等式总是成立,因此
(2—14)
(2—15)
现在我们从f-1开始,对Idet
Anl在使用卢2,如此继续下去,到徘l为止。然
后,从i=n开始,类似上述过程到到f=甜1为止。这样就不难获得不等式(2.11)。
9
电子科技大学硕士学位论文
(b)既然A7∈D,那么对∥重复以上的证明,那么不等式(2.12)也不难得到。
证毕。
注2.1.-有时候为了获得更好的结果,(2.11)和(2.12)可以混合使用,见后面的数值试验部分。
另外,利用下面的引理2.4,我们也可获得一个相似的行列式上界。引理2.4如果4=【嘞]∈D,那么A叫=【6:f『】存在,并且
I%巨西南,对某些歹≠£
j=l
(2。6)
证明:由引理2.3(a),我们知道I%巨啦蟹{‰|}.因此,我们可以假设
1%剐%圈%l,对任意f≠-,≠七.
由彳爿’1:,,可知
∑%%=l,对任意ieN,∑吩%=0,对任意j≠k.
这样,对任意i≠,≠k,有
『I%lI6flI+写(么)lbfi险1,
【1%II%|_弓(4)I%|<0.
解上面的不等式,我们获得
悱鬲‰2雨面1砺,对某些川・
(2-17)
这样证明被完成。
现在,由上面的引理2.4,类似于上面的定理2.1的证明,不难得到:定理2.2(a)如果A=[%】∈D,那么对任意指标七∈N,有
Idet彳匡I%IⅡ(I%I+哆弓)n(1ajjI+哆‘),
其中哆,髟同定理2.1(a)。并且等号成立当且仅当乃=0(.,=l,...,k-1)和‘=0(,=七+1,..卅)。
(b)若Ar∈D,那么对任意k∈N,有
10
第二章主对角占优矩阵类的行列式的界
[detA[4a肚I兀(1a∥I+/西jaj)H(I%I+帝juj),
j=l
』=t+l
(2—18)
其中廖,彩同定理2.1(b)。并且等号成立当且仅当嘭=o(_『=1,...,七一1)和
uj=0(,=七+l,…,玎)。
定理2.3如果A=[口{,】∈D,那么A一=【6:f『】,并且对任意后∈Ⅳ,有
la肚I兀(I%I-∑1%I%)兀(I%I_∑l嘞1w。j)-<ldetAI<
i=1
j>i
i=k+l
k<j<i
k一1
H
I+∑I%1w,j),I口齄In(I%I+∑I%f%)兀(I口if
j>i
i=l
i=k+l
k£j<i
其中%和%(f,J∈忉同定理2.1(a).并且上式等号成立当且仅当i<七,口{『=o
(_,>i)和以≥f>k,de=o(-,<O.
证明:首先,我们考虑上述不等式的左手边,由引理2.2,我们可以假设矩阵彳是个M.矩阵。下面由引理2.1,可得
detA---Ia.I
detA.一∑(一1)‘“叫吻Idet4i.
j-'{i
(2・19)
注意到A是一个M.矩阵,因此对任意f,jeN,b0≥o和det4i>oB31.这样由
(2.13),可知
(一1)‘H7’det4≥0.
现在,由引理2.3(a),可得
Ideta口I-<pjdet4,对任意i≠/.(2-20)
因此,由(2.19)推出
IdetAl2(Iau—ZEaF
』耐
IP,)ldet%jI,对任意f∈Ⅳ.
(2-20
以下的证明相似于定理2.1,在此不再详述。现在我们考虑上述不等式的右手边,由引理2.1,有
IdetAI_<)-?1au
j=l
IIdet
I,对任意i∈Ⅳ.
注意至lJ(2.20),那么有
IdetA[_<(Ia.I+∑l勺tP,)ldet4iI,对任意f∈Ⅳ.(2-22)
j=1.≠f
重复前面定理2.1(a)的相应的证明过程,易知结论成立.
注2.2:特征值是关于元素的连续函数,因此行列式也是关于元素的连续函数。
电子科技大学硕士学位论文
这样对上述定理2.1—2.3来说,在更弱的条件~一般主对角占优下也是成立的。
2.3双对角占优矩阵的情形
现在我们将上述结果推广到双对角占优矩阵的情形。为此我们先定义双对角
占优矩阵如下:
定义2.2【71设矩阵A=[av]EC一,那么矩阵A被称作双对角占优的(记为AE坟坟),如果对任意is,,有
I%11%巨g(a)Rj(铘.
如果上述不等式是严格的,则称之为严格双对角占优矩阵,记为彳∈DD。
引理2.5若么:f鬈】∈肋,并且鹄=0"11%Rj(A),J∈奶≠0,那么A~=晦]满
足
(a)
%巨型等掣I吮匿以Mf∈Ⅳ0,七∈Ⅳ\鼢
%笛肛IkI,i∈Ⅳo,k∈Ⅳ\{i);
(2—23)(2-24)
‰巨型铲㈧,f∈w雕砌.
4胜I
(2-25)
其中d是一个满足
蹴
证明:(a)既然No≠0,因此彳仨D。不失一般性,我们可以假设第i行非严格
对角占优(即,f∈Ⅳo),记
而:承l砺乳阵志,‘涎Ⅳo小Ⅳ\(。,止踊揣小揣.
L<d<,≤1.
12
嵩刚
习
%
正数
(2.26)
既然么=魄】∈DD,这里存在~个正数d使得
现在,构造一个正对角矩阵A=diag(人,,人:,...,A。),其中
第二章主对角占优矩阵类的行列式的界
中{L嚣
让c=【cF】_么人,那么我们有c{,=~人J,i,jeN。容易证明CeD。下面对于矩阵
c,由引理2.3(a),我if].-7pA获得
(1)Ik
(2)I%I<型型竺酱塑掣lb.|<以I%l,f∈Ⅳ0,后∈Ⅳ\{班
I<竽.掣:辟Iba胁Ⅳ0,七∈Ⅳ\{f};
aii
l
口
(3)I%B型尘掣16删悱Ⅳo,,,m∈Ⅳ\”
a
ao
因此,结论(a)成立。
(b)相似地,我们可以获得c=b】-从∈D.由引理2.3(b)和引理2.4,对于矩阵
c和i∈Ⅳo,/∈Ⅳ\{f>,我们有
la。I+皿(4)形
其中
她|<丽丽1
ajj
,
(2—27)
d?:垫丝±!!二型丝
。d
I
对任意i∈Ⅳ0,,eN\{玑我们又有
蟛型等型≤等嘶
’
l口口l
I口Fl
叫
(2-28)
、
。
这样由(2.27)和(2.28),可知结论成立。…
定理2.4如果彳=[口{,]eDD,且%={Jlla//|<吩(栅,jeN}≠o,那么
(|口肚I一∑la茸I乃)Idet如I_<IdetAI<(Ia艟I+Xl
a眵Ipj)IdetA,,l,k
e
N。.
(2-29)
证明:由引理2.5(a),相似于定理2.3的证明,结论不难获得。证明完毕。注2.3:注意到在定理2.4中如eD,因此我们对如可以连续执行定理2.1—2.3的界,获得矩阵A的行列式的界,见下一节中的例4.2.
2.4数值例子
例2.4.1令
13
电子科技大学硕士学位论文
73
A=:
—22
O-3
70
48
O6
—10
1—2
应用上面的结果到它,其中当k=l我们可以获得如下结果:
表2—1在不同定理下行列式的界
相关结果
(2—2)
行列式的界IdetA巨15
6336>ldetA|>480ldetAI>315
6423.7>IdetA巨471.184128>l
detA
(2—3)(2_4)(2-5)
(2.11)和(2.17)定理2.3
I'z-936
3387.9>ldetA陲1389.9
另外,在同样的假设k=-I下,若我们首先应用(2.11)然后应用(2.12),那么我们获
得IdetA险980.57,这皆好于直接使用(2.11)或(2—12)的结果。
例2.4.2.考虑双对角占优矩阵
彳=㈠].
现在,应用定理2.4,我们获得
[detAN(4+警+扣et如鲫6,
这接近于真实值detA=98。因此,定理2.4是文献。乒12】的推广和改进。
14
第三章代数多项式的根的几个简单估计
第三章代数多项式的根的几个简单估计
3.1引言与记号
代数方程是最古老的方程,有着重要的应用。正如前面第一章所述,古代就已知一次、二次代数方程的解法。然而,直到16世纪意大利一个靠自学成才的数学家塔尔塔利亚发现了一般一元三次方程的求根公式,在他的启发下,
意大利数
学家费拉里也很快给出了一元四次方程的求根公式。随后在他们工作的鼓舞下,人们开始继续寻求更高次数的多项式的求根公式。然而,三百多年过去了,仍没有太大的进展。直到1870年法国数学家若尔当介绍了法国数学家伽罗瓦的论文“关于用根式解方程的可解性条件”后这一问题才得以解决。以此文为基础的研究已成为近世代数所研究的最重要的课题之一。
由代数基本定理知道,每个次数不小于1的复系数多项式,在复数域上至少有一个根。然而要把这些根全部求解出来,一般来说还是很难的。目前,一般都采用代数方程的牛顿法算法,然而,正如第一章所述,牛顿算法对初值区间的选择要求是很高的,在某些情况下,对初值也十分敏感。因此,若能给出初值区间的一个较好的选择,必将对牛顿算法十分有益。
一般说来,实数域上的代数方程可以设为
只(z)=aoX”+口l,一1+…+an_iX+an,
(3・1)
式中的at(f=0,1,..棚)均为实数,且a。≠o。在本节中为了研究的方便,且不失一般
性,我们可以假设:
只(功=X“+aix”-1+…+an—lx+巳,
且a。≠0。
(3—2)
在很多数学理论与实际应用中需要对(3—2)的根的模进行估计。如齐次常系数线性差分方程
厂(DUk=Ut+H+aluI+月一l+…+%一luk+l+anUk
20
其所对应的特征方程就是上述的n次方程(3.2)。
关于上述n次方程(3.2)的根模的估计,最早由Cauchy给出【16】:
15
一一一_一———————————————————————————————————————————————————————————————一
(3・3)1名|<max{la.1,1+IaII,…,1+J口。一。0.
其中1名I表示上述多项式(3.2)的任意根的模。近来的主要结果有Carmichael-Mason的估计【7】:
电子科技大学硕士学位论文
IA阵、fl+∑I口i12,(3-4)
’
i=1
厂—i—一
Farmer.Loizou的估计【17】:
[3.1<2m。;。a;x。[at
Il
7‘,
(3-5)
还有Montel的上界估计:
r
。
1
lA
max{l,∑la,l}
L
i=1
J
(3-6)
等等。这些界形式简单,却有着很广泛的应用。
在本节中,我们将利用有关实矩阵特征值估计的几个最新结果,结合已有的研究技巧,改进上述的几个结论,这些结果对某些矩阵来说,有着更好的上界。
3.2主要结果
为了书写方便,我们将采取下边的记号:对任意A=[ag]ER脒“和f∈N,分别
定义矿=max{0,a0.1j『≠f),r-=vain{0,a/jIJ*f)和哆=aI口,一口川,歹∈N(a州2o)。
首先我们介绍几个有趣的引理如下:
引理3.1t17】设矩阵么=[嘞]∈R棚,名为矩阵A的任意实特征值,那么
旯∈G垒U[%一矿-Y.1r,+一%l,att--rl一寸∑I,;一一%I】.(3-7)
引理3.2【18】让1≤s≤n,仃>0,设
qlcr=r—lalIxPl一…一Ias—llx一(1asl+仃)
(3_8)(3-9)
%彳’一掣厂1…一掣
分别取其唯一正根岛。,,02矿。则多项式(3—2)的根满足:
I兄阵nlax以仃,岛口).
引理3.3【18】对多项式(3.2),其任意根允将满足:
(3-10)
16
第三章代数多项式的根的几个简单估计
……㈦hI+盯面]).
首先,我们先给出多项式(3.2)的任意实根的几个估计。定理3.1对多项式(3.2),其任意实根力将满足:
m阵叫LF可,F硒,,B,・,
-aI
-a2
其中口=min{0,一aI2≤,≤,l},b=min{O,勺l2≤/≤,z}。
证明由第一章可知,多项式f3.21的根即为其友矩阵
…一口n一1…O…
0
一an
1O1
OO
q=
0
(3—12)
0O
…
1O
所对应的特征值。那么允2将是矩阵q‘191:
6l
-aI
62
-a2
一
2
吃一。
0
吃
0
屯
玎、J
-a^一1—口n
q=
1O
(3・13)
OO
址%o;●
O0
对上述矩阵,显然有
min{O,哆I2≤歹≤珠f=1;min{O,-ayI2≤歹≤以},f=2;
0,
i≠1,2.
因此,由引理3.1不难得到
咫懋∞m驴刮,-a:-a+Z脚I…。I),
从而,易知结果成立。证毕。
注3.1:由定理3.1的证明可知,当多项式(3.2)的系数相差不大时。在此条件下的定理3.1的估计结果是很好的,见下面数值试验部分。
17
电子科技大学硕士学位论文
下面,我们用“广义友阵”的定义,进一步优化上述结果。
定义3.1C201令只(x)是一型如(3.2)的多项式,那么这里存在复数c1,c2,o..gCn和不
为零的复数dl,畋,...,dn,且满足
-a1.=CI,
_%≥c2西’
:
(3.14)
一an暑qd。一l…以西.
使得矩阵彳:
0以一,
:
・.
‘
O
●●
●
…
A=
:
:
●●
●
●●
●
(3-15)
0
・・・
…
0
qq一1…c2
的特征多项式就是多项式只(功。在本文中,我们称矩阵4为多项式£(功的广义友
阵。
利用广义友阵和引理3.1,我们可以类似地得到如下定理3.2,其在一定程度上改进Y(3.5)式。
…一p乒可,Fi礴丽,,p旧
定理3.2对多项式(3.2),其任意实根五将满足:
l
其中扛m:鲻ax。l%li,q=劳,H,2,…咒,e=min
证明类似于定理3.1的证明,可知
oo
42=
:
●
0,dcj
l
l<_j<n)和
f=min{O,ClCj+dcj+。I1≤歹≤咒}(艺:。=o)。
00
:
●
以一,吒一2
O
。.
●
…
O…
’.
●
●
0o
:
以一:吒一3
‘.
吐乞
clc;
匾巳一。qcol+以一lCn
……
盔c3qc3+以cj
磊c2
qc2+d2c3
西clqq+吐c2
令
第三章代数多项式的根的几个简单估计
dn_|-,o.-西乖戮㈣÷),
q=寿,
推出结论成立。证毕。
下面,对所有根的模进行估计。
定理3.3对多项式(3.2),以及任意1≤s≤n,盯>0,其任意根兄将满足:
(1)
拈l,2,…见
由文献【20l可知,上述选择满足定义3.1的相关条件。因此,由引理3.1,不难
I旯I_<max{,ol,岛},
(3—18)
其中n:罂m儿…+仃}∥,器。ft+掣,掣}o
(2)
l见I≤max{・,三[I
01
00
q
l+√l:『面],三≥[I%¨l+√l≤面]>.
(3—19)
q
证明根据前面第一章中,关于多项式友阵的定义,可知引理3.2中多项式g。,和q:。所对应的友阵,可分别写成如下形式‘181:
0……I
0……
l+盯
;
O1
O……∑型
盯
Iq一。l
,E=
0……
1:j:
;i
(3-20)
1’.
互=
‘.
’.
:
‘.0
1
Ia2l
al
I
1也J
仃
对上述矩阵,分别采用行和范数(即州l。),易得
胪。m蚓ax川吼…竹"02"-<s班m“ax・+掣,掣>
因此,由引理3.2易知,结论(1)成立。
类似地,对(3.20)分别采用引理3.3的结论,易知,结论(2)成立。证毕。推论3.4对多项式(3.2),其任意根名将满足:
19
电子科技大学硕士学位论文
I五I<max{[口。l,1+IqI,…,1+I口。一。I).
证明只要令定理3.3(1)中,J=n,仃=1,即可。证毕。
注3.2:在定理3.3中,对J和仃取不同的值就会得出许多不同的根的估计,如:推论3.4。因此,定理3.3可以说是对一些多项式根模估计式(如(3.3)等)的进一步推广或改进。
3.3数值例子
例3.3.1考虑一个简单多项式只(x)=工3+x2—6x一6,估计此多项式根模的最大值。
解:易知其友阵为
q=㈧.
利用Matlab7.1直接计算,它的3个根分别为:
它们皆为实根。显然这些根的模的最大值为:2.4495。
应用上面的结果到它,我们可以获得如下结果:
表3-1在不同定理下多项式根模的界
相关结果
Cauchy(3-3)
多项式根模的界
。
7.00008.60234.899013.00004.3589
Carmichael-Mason(3-4)
Farmer-Loizou(3—5)
Montel(3-6)
定理3.1(3.11)
显然,我们的结果皆好于这些经典结果。
例3.3.2.考虑圪(x)=∥一,一2x4—3x3—4x2—5x一6,估计此多项式根模的最大值。
解:易知其友阵为
20
第三章代数多项式的根的几个简单估计
ll0
c=
00O
000
10O
0l0
O0l
00O
201
30O
40O
50O
6OO
利用Matlab7.1直接计算,它的6个根分别为:
毛=0.8317+1.1719i,
x2=0.8317-1.1719i.x3=-0.2117+1.3341i,x4=一0.2117-1.3341i.x5=一1.1200+0.5813i.
黾=一1.1200-0.5813i.
这些根的模的最大值为:1.4370。
既然,它没有实根,下面我们利用定理3.3对根模的最大值估计如下:
表3-2在不同s,仃取值下定理33(1)的结果
SS=lS=2s=3J=4s=5s=6
仃=0.87.57.57.57.57.57.57.5
仃=1.06.O6.O6.O6.06.06.06.0
盯=2.03.54.O5.06.07.O8.0
仃=3.04.O5.O6.O7.O8.09.04.O
最小值瓠S
21
电子科技大学硕士学位论文
表3-3在不同S,盯取值下定理3.3(2)的结果
Js=lS=2S=3s=4s=5s=6
仃=0.86.15546.593611.524719.524729.524741.52476.1554
盯=1.05.35895.653311.618019.618029.618041.61805.3589
盯=2.03.54146.012.020.030。O42.03.5414
仃=3.02.80546.302812.302820.302830。302842.30282.8054
最小值
在这些估计中,最小值是2.8054,其已非常接近最大模根1.4370。这表明在
某些情况下,定理3.3的估计也是很精确的。
第四章结论和展望
第四章结论和展望
4.1本论文研究总结
本文的研究主要分为两大部分:
1.在第二章中,主要证明了当满足一定条件的时候,主对角占优矩阵行列式界的问题。通过研究逆元素的估计,获得了一些有趣的行列式估计式。如:
(1)对一般对角占优矩阵,我们获得了如下几个结果:定理2.1如果A=【口f『】∈D,那么对任意指标k∈N,有
IdetAI到a肚H(1a∥I-嘭乃)H(Ia//一哆‘),
其中g=t学{%},哆=嬲{%)和%=鱼篙≠,_,∈Ⅳ-
乃=0(_,=1,...,k-1)和,,=0(,=k+1,...,以)。
k-I
H
并且等号成立当且仅当
定理2.2如果A=k】∈D,那么对任意指标k∈N,有
detAIqaa
ILI(I%I+OjS)兀(Ia//I+哆‘),j=lj=k+l
其中岛,髟同定理2.1(a)。并且等号成立当且仅当,:『=o(.,=1,...,k-1)和‘=0(J=k+1,..・,n)。
-
定理2.3如果彳=lay]∈D,那么A~=【%】,并且对任意尼∈Ⅳ'有
I口肚f兀(1口ffI-∑IavI%)n(I%J_∑l吩1wo)-<ldetAI<
i=l
j>l
i=k+l
Is,可
k-1
n
I%I兀(1q;I+∑l嘞I%)兀(I%I+∑I口{『I~),
1=1
j>i
i=k+l
k≤j<i
其中%和嘞(f,-,∈加同定理2.1(a).并且上式等号成立当且仅当i<后,%=0
(歹>力和以≥i>k,口J,=O(j<f).
(2)对双对角占优矩阵,获得了如下结论:
定理2.4如果彳=心】∈肋,且Ⅳo=Ulla口|<吩(彳),,∈奶≠o,那么
电子科技大学硕士学位论文
(Iart
I一∑1%Ipj)IdetAn
J}k
I-<1detAI<(Iaa
J+∑1%Ipj)ldeta瓣I,keNo.
i{k
2.在第三章中,我们研究了经典代数多项式的根模估计问题,获得了一系列有趣的结果,推广或改进了一些经典的相关结论。如:
定理3.1对多项式(3—2),其任意实根五将满足:
…一tL
(1)
F研,F砸j,
I五喀max{Pl,P2},
其中a=min{0,-ay2≤jf≤以),b=rain{0,岛l2s/≤,z}。
定理3.3对多项式(3-2),以及任意1≤J≤n,仃>0,其任意根兄将满足:
其中Em华a¨xm^hm舻川m轳ax州{・+掣,掣)o
(2)
4.2前景展望
…≤一㈦旧l+一Ⅲ‰l+阿]).
矩阵的行列式和代数多项式的根模估计问题,在科学研究中已有着悠久的历史和广泛的应用背景。国内外很多著名学者都曾对其进行过深入而细致的研究,如著名法国数学家伽罗瓦等。对其研究的成果,也在一定程度上影响了数学的发展,如群论的出现。
如今,矩阵的行列式和代数多项式的根模估计问题也仍备受国内外学者的青睐。总的来说,随着研究的越来越深入,对这些量的刻画也越来越细微,所用知识也越来越复杂。本文所论结果,只能说是“管中窥豹”,还有待继续对其深入研究。望同行多加斧正。
致谢
致谢
在短暂的几年多的学习和生活中,我得到了来自导师,亲人,朋友及同事对我多方面的关怀及帮助,使我得以顺利地满怀自信地投入到工作和学习中,且不断进步和提高,在此我将深深地感谢他们,并以此不断激励自己不断进取!
首先,感谢我的导师蒋泽云副教授,在我生活上,学习上给予我细心的关怀和极大的精神鼓励。每当我科研陷入困境时,是他给予我最诚挚的帮助和热情洋溢的鼓励。可以说我每一点小小的成绩都凝聚了他无数的心血!且他对工作和科研的一丝不苟,追求完美的精神也使我终身受用!他严谨的科学态度,渊博的专业知识,开阔而活跃的思维方法,高度的敬业精神以及对科学问题敏锐的洞察力都深深的影响了我。在此,谨向他表示衷心的感谢和敬意!
最后,我要深深感谢我的父母和家人,是他们的关爱和理解始终激励着我刻苦学习,不断进取!衷心感谢所有帮助或支持过我的人们!
电子科技大学硕士学位论文
参考文献
1J黄廷祝,成孝予.线性代数与空间解析几何(第二版)一匕京:高等教育出版社,2003张顺燕.数学的源与流.北京:高等教育出版社,2003
n口p,
1J陈景良,陈向晖.特殊矩阵_匕京:清华大学出版社,2001沈永欢,梁在中,许覆瑚,蔡蓓蓓。实用数学手册.北京:科学出版社,1992
D.E.Crabtrcc.ApplicationofM-matricesto1J
1Jnonnegmivematrices.DukeMath-J.,1966,197-208R.HornandC.R.Johnson.MatrixAnalysis.CambridgeUniversityPress,NewYork,1985p陋№口1JR.B.Bapat.T.E.S.Raghavan.NonnegafiveMatricesand
Press,1997Applications.CambridgeUniversity
【8】
【9】黄廷祝,杨传胜.特殊矩阵分析及应用一匕京:科学出版社,2007M.I.Gil.StabilityofFiniteandInf'miteDimensionalSystems.KluwerAcademic
Publishers,Dordrecht,1998
[10]J.LiuandYHunag.SomepropertiesOnSchurcomplementsofH—matricesanddiagonally
dominantmatrices.LinearAlgebraAppl.,2004,389:365-380
GB.Price.Boundsfordeterminants、)l,ithdominantprincipaldiagonal.ProceedingsoftheAmer.Math.Soc.,1951,2:497-502
[12】
【13】R.Oeder.ProblemE.949.Alller.Math.Monthly,1951,58:31A.M.Ostrowski.NoteonboundsfordeterminantswiⅡldominantprincipaldiagonal.
ProceedingsoftheAlner.Math.Soc.,1952,3:26—30
【14】B.LiandM.J.Tsatsomeros.Doublydiagonallydominantmatrices.LinearAlgebra
Appl.,1997,261:221-235
[15】KyFan.InequalitiesforsumoftheM-matrices,InequalitiesI.AcademicPress,NewYorkand
London,1967,105-117
[16】X.R.Yong.Proofof
2000,320:167—171aconjectureofHeAlerandMarkham.LinearAlgebraAppl.,
【17】Ch.L.WangandGJ.Zhang.Somesimpleestimatesforthesingularvalues
Math.Appl.Sinica,2002,18:117—122ofmatrices.Acta
【18】M.FujiiandEkubo.Operatornoil31s丛boundsfor
Acad.,1973,49:805—808rootsofalgebraicequations.Proc.Japan.
26
参考文献
[19】J.M.Pena.AclassofP啦atriccSwithapplicationstothelocalizationoftheeigenvaluesofa
realmamx.SIAMJ.MatrixAnal.Appl.,2001,22:1027—1037【20】赵维加.多项式零点的界的一种估计.青岛大学学报.2000,13:14.19
电子科技大学硕士学位论文
攻硕期间取得的研究成果
[1】刘伟,主对角占优矩阵行列式的估计与应用,已投《工程数学学报》[2]刘伟,有关代数多项式根的几个简单估计问题研究,已投《大学数学》
矩阵行列式和代数多项式根的计算问题作者:
学位授予单位:刘伟电子科技大学
相似文献(0条)
本文链接:http://d.g.wanfangdata.com.cn/Thesis_Y1314414.aspx
下载时间:2010年4月2日