第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx

上传人:b****6 文档编号:21220679 上传时间:2023-01-28 格式:DOCX 页数:22 大小:25.85KB
下载 相关 举报
第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx_第1页
第1页 / 共22页
第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx_第2页
第2页 / 共22页
第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx_第3页
第3页 / 共22页
第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx_第4页
第4页 / 共22页
第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx_第5页
第5页 / 共22页
点击查看更多>>
下载资源
资源描述

第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx

《第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx》由会员分享,可在线阅读,更多相关《第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx(22页珍藏版)》请在冰豆网上搜索。

第五专题矩阵的数值特征行列式范数条件数迹秩相对特征根讲解Word文档格式.docx

A=0;

8.AB(即AB0)

则tr(A)

tr(B),且等号成

立的充要条件是A=B(A

B

i(A)i(B));

9.对于n阶方阵A,若存在正整数k,使得Ak=0,

则tr(A)=0(从Schur定理或Jordan标准形证明)。

若干基本不等式

对于两个m×

n复矩阵A和B,tr(AHB)是m×

维酉空间上的内积,也就是将它们按列依次排成的两

个mn维列向量的内积,利用Cauchy-schwarz不等式

[x,y]2≤[x,x]﹒[y,y]

定理:

对任意两个m×

n复矩阵A和B

|tr(A

H

2

B)

B)|≤tr(AA)﹒tr(B

这里等号成立的充要条件是A=cB,c为一常数。

特别当A和B为实对称阵或Hermit矩阵时

0≤|tr(AB)|≤tr(A2)tr(B2)

设A和B为两个n阶Hermite阵,且A≥0,B≥0,则

0≤tr(AB)≤λ1(B)tr(A)≤tr(A)﹒tr(B)

λ1(B)表示B的最大特征值。

证明:

tr(AB)=tr(A1/2BA1/2)≥0,又因为

A1/2[λ1(B)I-B]A1/2≥0,所以λ1(B)tr(A)≥A1/2BA1/2,得

tr(AB)=tr(A1/2BA1/2)≤tr(λ1(B)A)

=λ1(B)tr(A)≤tr(A)﹒tr(B)

推论:

设A为Hermite矩阵,且A>

0,则

tr(A)tr(A-1)≥n

另外,关于矩阵的迹的不等式还有很多,请参考

《矩阵论中不等式》。

三、矩阵的秩

矩阵的秩的概念是由Sylvester于1861年引进的。

它是矩阵的最重要的数字特征之一。

下面讨论有关矩阵秩的一些性质和不等式。

矩阵A的秩定义为它的行(或列)向量的最大无关组所包含的向量的个数。

记为rank(A)

1.rank(AB)min(rank(A),rank(B));

2.rank(AB)rank(A,B)rank(A)rank(B);

3.rank(AAH)rank(AH)rank(A);

4.rank(A)rank(XA)rank(AY)rank(XAY),

其中X列满秩,Y行满秩(消去法则)。

定理(Sylvester):

设A和B分别为m×

n和n×

l

矩阵,则

rank(A)rank(B)nrank(AB)

min(rank(A),ra

Sylveste定理是关于两个矩阵乘积的秩的不等式。

其等号成立的充要条件请参考王松桂编写的《矩阵论中不等式》,三个矩阵乘积的秩的不等式也一并参考上述文献。

四、相对特征根

设A和B均为P阶实对称阵,B>

0,方程

|A-λB|=0的根称为A相对于B的特征根。

|A-λB|=0等价于|B-1/2AB-1/2-λI|=0

(因为B>

0,所以B1/2>

0)

注:

求A相对于B的特征根问题转化为求B-1/2AB-1/2的特征根问题或AB-1的特征根。

因B-1/2AB-1/2是实对称阵,所以特征根为实数。

使(A-λiB)li=0的非零向量li称为对应于λi

的A相对于B的特征向量。

①设l是相对于λ的AB-1的特征向量,则

-1-1-1

ABl=λl或A(Bl)=λB(Bl)

-1

Bl为对应λ的A相对于B的特征向量

(转化为求AB-1的特征向量问题)。

②设l是相对于λ的B-1/2AB-1/2的特征向量,则

B-1/2AB-1/2l=λl

可得

A(B-1/2l)=λB(B-1/2l)

则B-1/2l为对应λ的A相对于B的特征向量

(转化为求B-1/2AB-1/2对称阵的特征向量问题)。

五、向量范数与矩阵范数

向量与矩阵的范数是描述向量和矩阵“大小”的

一种度量。

先讨论向量范数。

1.向量范数定义:

设V为数域F上的线性空间,若对于V的任一向量x,对应一个实值函数x,并满

足以下三个条件:

(1)非负性x0,等号当且仅当x=0时成立;

()齐次性x

x,

k,xV;

(3)三角不等式xy

x

y,x,yV。

则称x为V中向量x的范数,简称为向量范数。

定义了范数的线性空间定义称为赋范线性空间。

例1.x

Cn,它可表示成x

12

T,i

C,

1

x2

就是一种范数,称为欧氏范数或2-范数。

i

(i)非负性

0,

当且仅当i

0i

1,2,

n

时,即x=0时,x

2=0

(ii)齐次性

(iii)三角不等式

y

T

C

,i

y2

2Reii

2i

2x

根据H?

lder不等式:

p

q

1,1

aibi

aip

biq

,p,q

1,ai,bi0

2.常用的向量范数(设向量为x12

T)

1-范数:

x1

∞-范数:

x

max

i;

1in

P-范数:

2-范数:

i(p>

1,p=1,2,,∞,);

xHx

2;

椭圆范数(2-范数的推广):

xA

xHAx2,A为Hermite正定阵.

加权范数:

xwwii

当AWdiagw1w2wn,wi0

xp显然满足非负性和齐次性

(iii)y

xp

,yp

,xyp

p1

yp

应用H?

lder不等式

p1q

1q

3.向量范数的等价性

定理设、为Cn的两种向量范数,则必定存

在正数m、M,使得mxxMx,(m、M与x

无关),称此为向量范数的等价性。

同时有

M

m

(1)对某一向量X而言,如果它的某一种范数小(或大),那么它的其它范数也小(或大)。

(2)不同的向量范数可能大小不同,但在考虑向量序列的收敛性问题时,却表现出明显的一致性。

4、矩阵范数

向量范数的概念推广到矩阵情况。

因为一个m×

n阶矩阵可以看成一个mn维向量,所以Cmn中任何一种向量范数都可以认为是m×

n阶矩阵的矩阵范数。

1.矩阵范数定义:

设Cmn表示数域C上全体mn

阶矩阵的集合。

若对于Cmn中任一矩阵A,均对应一个实值函数A,并满足以下四个条件:

(1)非负性:

A

0,等号当且仅当A=0时成立;

(2)齐次性:

A,

C;

(3)三角不等式:

AB

AB,A,BCmn,则称

A为广义矩阵范数;

(4)相容性:

ABAB,则称A为矩阵范数。

5.常用的矩阵范数

(1)Frobenius范数(F-范数)

F-范数:

=

A)=

AF

aij

trace(A

i,j

i,j1

矩阵和向量之间常以乘积的形式出现,因而需要

考虑矩阵范数与向量范数的协调性。

如果矩阵范数A和向量范数x满足

AxAx

则称这两种范数是相容的。

给一种向量范数后,我们总可以找到一个矩阵范

数与之相容。

(2)诱导范数

设A∈Cm×

n,x∈Cn,x为x的某种向量范数,

AmaxAx

x1

则A是矩阵A的且与x相容的矩阵范数,也称之为

A的诱导范数或算子范数。

(3)p-范数:

Apmax

Axp,

A

aijmn,x为所有可能的向量,x

2n

T,

xp,Axp

1Ax

maxAx

xp1

A1

1,x1

i1,Ax1

aijj

x11

j1

可以证明下列矩阵范数都是诱导范数:

(1)A

列(和)范数;

1jn

(2)A2

i(AHA)谱范数;

1i

AHA的最大特征值称为AHA的谱半径。

当A是Hermite矩阵时,A

maxi(A)

是A的谱

半径。

谱范数有许多良好的性质,因而经常用到。

AH

A2;

AHA2

A2

(3)

Amax

行(和)范数

1im

i,x

定理矩阵A的任意一种范数A是A的元素的连

续函数;

矩阵A的任意两种范数是等价的。

定理设A∈Cn×

n,x∈Cn,则AF和x2是相容的

Ax2A

由于Ax2A2x

Fx2

2AFx2成立。

n,则AF是酉不变的,即对于任意酉矩阵U,V∈Cn×

n,有

AFUAVF

UAVFtr[(UAV)H(UAV)]

tr[VHAHUHUAV]tr[VHAHAV]

tr[AHAVVH]tr(AHA)

定义设A∈Cn×

n,A的所有不同特征值组成的集合

称为A的谱;

特征值的模的最大值称为

A的谱半径,

记为ρ(A)。

定理ρ(A)不大于A的任何一种诱导范数,即

ρ(A)≤A

设λ是A的任意特征值,x是相应的特征向

量,即

Ax=λx

|λ|·

||x||=||Ax||≤||A||·

||x||,||x|≠0

|λ|≤||A||

试证:

设A是n阶方阵,||A||是诱导范数,当||A||<

1时,I-A可逆,且有

||(I-A)-1||≤-||A||)(1-1

若I-A不可逆,则齐次线性方程组

(I-A)x=0

有非零解x,即x=Ax,因而有

||x||=||Ax||||﹒≤||x||<

||x|||A但这是不可能的,故I-A可逆。

于是(I-A)-1=[(I-A)+A](I-A)-1=I+A(I-A)-1因此||(I-A)-1||≤||I||+||A(I-A)-1||=1+||A(I-A)-1||

≤1+||A||﹒||(I-A)-1||

即证

补充证明||I||=1:

由相容性可知:

||A||﹒||A-1||≥||AA-1||=||I||

xIxIxI1

对于诱导范数(AmaxAx)

x1

ImaxIx1

六、条件数

条件数对研究方程的性态起着重要的作用。

设矩阵A是可逆方阵,称||A||﹒||A-1||为矩

阵A的条件数,记为cond(A),即cond(A)=||A||﹒||A-1||

(1)cond(A)≥1,并且A的条件数与所取的诱导范数的类型有关。

因cond(A)=||A||﹒||A||

≥||AA

||=||I||=1

(2)cond(kA)=cond(A)=cond(A-1),这里k为任

意非零常数。

当选用不用的范数时,就得到不同的条件数,如:

cond1(A)=||A||1﹒||A-1||1

cond∞(A)=||A||∞﹒||A-1||∞

cond2(A)=||A||2﹒||A-1||2=1,其中1,n分别

为AHA的特征值的模的最大值和最小值。

谱条件数特别地,如果A为可逆的Hermite矩阵,则有

cond2(A)=

这里1,n分别为A的特征值的模的最大值和最小值。

如果A为酉阵,则cond2(A)=1

例求矩阵A的条件数cond1(A),cond∞(A)

152

A210

382

解:

||A||1=max{6;

14;

4}=14;

||A||∞=max{8;

3;

13}=14;

6

4

8

13

23

11

||A-1||1=17/4;

||A-1||∞=47/4;

cond1(A)=||A||1﹒||A-1||1=14×

17/4=259/2;

cond∞(A)=||A||∞﹒||A-1||∞=611/4。

例设线性方程组Ax=b的系数矩阵A可逆。

讨论当b有误差δb时,解的相对误差δx的大小。

因矩阵A可逆,所以Ax=b有唯一解x=A-1b,设解的误差为δx,由

A(x+δx)=b+δb

Aδx=δb或δx=A-1δ

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 考试认证 > 从业资格考试

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1