✅ 操作成功!

次方怎么算

发布时间:2023-06-04 作者:admin 来源:文学

次方怎么算

次方怎么算

-公务员退休金

2023年2月15日发(作者:西洋音乐)

1/16

矩阵n次方几种求法

1.利用定义法

,,

ijkj

snnm

AaBb



则,

ij

sm

Cc

其

1122

...

ijijijinnj

cababab

称为A及B乘积,记为,则由定义可以看出矩阵A及B

乘积C第i行第j列元素等于第一个矩阵A第i行及第二个

矩阵B第j列对应元素乘积之和,且由定义知:第一个矩

阵列数及第二个矩阵行数要相1同。

例1:已知矩阵,,求

解:设CAB=

34

ij

c

,其中1,2,3i;1,2,3,4j

由矩阵乘积定义知:

11

1526533032c

12

1122543231c

13

1321553030c

14

102051305c

21

150623101c

22

110224129c

23

130125107c

24

100021102c

31

c

2/16

32

c

33

c

34

001031403c

将这些值代入矩阵C中得:

CAB=

则矩阵An次方也可利用定义方法来求解。

2.利用矩阵分块来求解

这类方法主要是把一个大矩阵看成是由一些小矩阵

组成,就如矩阵由数组成一样在运算中将这些小矩阵当

做数一样来处理,再由矩阵乘法定义来求解这些小矩阵

乘积所构成矩阵。即设,,

ijkj

snnm

AaBb



把A,B分解成

一些小矩阵:

,,其中

ij

A是

ij

sn小矩阵且1,2...it,1,2...jl,且

12

...

t

ssss,

12

...

l

nnnn;

ij

B是

jk

nm小矩阵且

1,2...jl,1,2...kr;且

12

...

l

nnnn,

12

...

r

mmmm;

令CAB=,其中

ij

C是

ij

sm小矩阵且1,2...it,1,2,...,jr,

12

...

t

ssss,

12

...

r

mmmm;其中

3/16

1122

...

ijijijillj

CABABAB。这里我们应注意:矩阵A列分法

必须及矩阵B行分法一1致。

例2:已知矩阵

45

10025

01013

00128

00006

A













,,求AB

解:将

45

45

10025

10025

01013

01013

00128

00128

00006

00006

A





























,其中

,

22

06A,,

由矩阵乘积法则知:

由矩阵加法和乘积法则1知:

则矩阵An次方求解也可利用以上方法来求解。

3.利用数学归纳法求解

这种方法及矩阵定1义和数学归纳3法相结合,从而找

出规律再求解,但是这种方法比较适合低阶且有规律方

4/16

阵n次方运2算。

例3:已知,求nA

解:当2n时

2cossincossincossin

sincossincossincos













22

22

cos2sin2

cossin2cossin

sin2cos2

2cossincossin

























当3n时

32cossincossincossin

sincossincossincos













cos2cossin2sincos2sinsin2cos

cos2sinsin2coscos2cossin2sin















所以假设nA=

当1k时成立,假设当1kn时成立;则当kn时

1cossincossin

sincossincos

n

nA

















cos1sin1

cossin

sin1cos1

sincos

nn

nn

























由矩阵乘法定及三角函数知:nA=则假设成立。

所以nA=

5/16

4.利用分拆法求解

这类方法主要是将一个矩阵分解成一个单位矩阵和

另外一个矩阵之和再求1解,且另外这个矩阵n次方计算

起来比较简2单。

例4:已知,求nA

解:AEB,其中,矩阵E为单位阵且2EE

EBBEB;故

nA=122+CCCn

nn

nnn

EBEBBB

由2

010010001

001001000

000000000

B













2

3

B













则3n时,nB=0。故122n

nn

AECBCB

由矩阵加法运算法则1知:

nA=

5.利用相似矩阵求解(利用对角矩阵来求)

定义:设矩阵A,B为数域P上两个n级矩阵,如果可

6/16

以找到数域P上n级可逆阵X,使得矩阵1BXAX,就说A

及B相1似。如果矩阵A或B有一个可以化成对角矩阵则计

算比较简便。而判断矩阵A可对角化条件1有:

1)矩阵A可对角化必要条件是矩阵A有n个不同特征值

2)矩阵A可对角化充要条件是矩阵A有个n线性无关特

征向量

3)在复数域上矩阵A没有重根

而求矩阵A特征值和特征向量方法1有:

1)求矩阵A特征多项式EA在数域P中全部根,这

些根是矩阵A全部特征值。把这些所求特征值逐个代入方

程组0EAX中,对于每一个特征值,解方程组

0EAX,求出一组基础解系,那么这个基础解系就

是属于这个特征值特征向量。

再利用判别法判断矩阵A是否可对角化。

例5:已知矩阵,求nA

解:易知矩阵A特征多项式EA=

由行列式计算方法知:

EA=213113

7/16

所以矩阵A特征值为1,1,3。

当特征值为1时,解方程0EAX,由齐次线性方

程组计算方法知:0EAX基础解系为

1

a=111

;

所以矩阵A属于特征值1全部特征向量为

1

111k

,其

1

k0。

当特征值为1时,解方程0EAX,由齐次线性

方程组计算方法知:0EAX基础解系为

2

a=110

;所以矩阵A属于特征值1全部特征向量为



2

110k

,其中

2

k0。

当特征值为3时,解方程30EAX,由齐次线性方

程组计算方法知:30EAX基础解系为

3

a=011

,

所以矩阵A属于特征值3全部特征向量为

3

011k

,其

3

k0。

则由矩阵A可对角化条件知:矩阵A可对角化且对

角阵为

B

令=,由求逆矩阵方法知:

因为线性变换在不同基下所对应矩阵是相似知:

1CACB

8/16

所以11

n

nnCACCACB,则



33

33

100

100

010010

003003

n

n

n

n

B





















由1nnACBC,由矩阵乘法运算法则知:





33

11111

3113111

13131

nn

nn

nnn

nn

A



















2)对方阵A,设

1

FEA

,对1n

FE做初

等变换,化成DP其中D为上三角阵,则矩

阵D主对角线上元素乘积多项式根即为A特征根

i

。对矩阵A任一特征根

i

,代入DP中,若

i

D

中非零向量构成一满秩矩阵,则

i

D行向量所对应



i

P中行向量

i

即为

i

特征向量;否则,继续施行初等

行变换,使得

i

D中非零向量构成一满秩矩阵,则

i

D

中零向量所对应

i

P中行向量

i

即为

i

特征向8量。

这类问题所涉及定理是:对任意方阵A特征矩阵F经

过行变换,可化为上三角矩阵G,且G主对角线上

元素乘积多项式根即为矩阵A特征值。

9/16

例6:已知矩阵A,求nA

解:

3

211100

121010

112001

FE



















作初等行变换

121010

211100

112001

















2

121010

01020

43

011011



























33

121010

011011

00113

41

DP





















由上述定理知:矩阵A特征值为1(二重),4。

当1时,

111010

11000011

000112

DP















,由2)

中判别法知:矩阵特征向量为:

1

011

,



2

112

。

当4时,

121010

44033011

000111

DP















,由2)

中判别法知:矩阵A特征向量为:

3

111

。

10/16

则由相似矩阵条件知:矩阵及对角矩阵相似且对角矩阵

则存在可逆阵使得

由求可逆阵方法知:

由1

100

010

004

n

nATT





















知:

nA=







111

424141

333

111

414241

333

111

434142

333

nnn

nnn

nnn

























6.利用若当形矩阵求解

这类方法主要是运用任何一个n级复矩阵都相似一

个若当形矩阵和利用相似矩阵相关定理及化若当形矩阵

方1法。

例7:已知矩阵A,求nA

11/16

解:

126

13

114

EA

















,由求初等因子方法知:

EA初等因子为1,21;所以矩阵A若当标

准形为:

则存在可逆阵P,使得1PAPJ,则APPJ。

设,其中

1

a

,

2

a

,

3

a

为列向量

将矩阵P代入APPJ得

11

Aaa



,

223

Aaaa



,

33

Aaa



由齐次线性方程组:0AEX,即

1

2

3

2260

1130

1130

x

x

x

















,则

1

301a



,

3

211a



是

齐次线性方程组解且

1

a

,

3

a

是线性无关,则

1

a

,

3

a

是由齐次线性方程组:0AEX基础解系。

由:

3

AEXa

有解

2

100a



且

1

a

,

2

a

,

3

a

线

性无关。由数学归纳法知:

由求可逆阵方法知:

由1PAPJ知:1APJP

则1nnAPJP=

33

1226

13

31

nnn

nnn

nnn

















12/16

7.利用多项式求解

主要运用带余除法即:对于数域Px中任意两个多项式

fx和gx,其中gx0,一定有Px中多项式qx,rx

存在使得fxqxgxrx成立,其中rxgx或

rx=0,并且这样qx和rx是唯一1的。

7.1特征多项式无重根

例8:已知矩阵A,求nA

解:设f为矩阵A特征多项式,则fEA

由计算行列式方法知:

211f

由带余除法及辗转相除法则1知:设

nfqr,其中rxfx;由

3fx,所以设2rabc。将特征多项式

0f根代入nfqr中得:

解得,0b,;

所以2

11

2142

33

nnnfq

由哈密顿—凯莱定1理:A是数域P上一个级矩阵,

13/16

fEA是矩阵A特征多项式则0fA。

将A代入2

11

2142

33

nnnfq中得:

2

11

2142

33

nnnAAE

由矩阵乘法定义知:,

所以由矩阵加法运算法则知:



1

33

1220

020

5

0211

3

n

nn

n

A

















7.2特征多项式有重根

例9:已知矩阵A,求nA

解:设f为矩阵A特征多项式,则fEA

由行列式计算方法知:221f

由带余除法及辗转相除法知:nfqr,

其中rxfx;由3fx,所以设

2rabc。将特征多项式0f根代入

nfqr中得:

14/16

因为1是0f2重根。

由定理:如果不可约多项式Px是fxk重因式

(1k),则它微商fx

是1重因式.则1是f

3根。

则由导数定义及性质:对nfqr等号

两边同时求导得:1nnfqfqr



则将1代入1nnfqfqr

中得:

2abn;则由

解得:21nan,1322nbn,22ncn。

由哈密顿—凯莱定理知:0fA。

则将矩阵A代入nfqr中得:

212132222nnnnAnAnAnE

由矩阵乘法运算法则知:

由矩阵加法运算法则知:

33

120

4210

212212

n

nnn

nn

Ann

nn















8.总结

上述七种方法求解矩阵n次方乘积适用于求低阶矩

15/16

阵n次方乘积适用于求低阶矩阵n次方计算,而对于高阶

矩阵求解则比较困难。利用方块、拆项、数学归纳法和

相似矩阵方法求解适用于比较特殊一些矩阵求解;利用

定义、若尔当形矩阵和多项式方法对于普通矩阵都适用,

但利用定义方法对于求矩阵n次方计算比较复杂;而利用

多项式和若尔当形矩阵方法有利于对所学知识及时巩

固、能加深对所知识理解,而这两种方法提供了解这类

问题行之有效方法且容易掌握。

参考文献

[1]同济大学应用数学系,高等代学,高等教育出

版社,2008.

[2]钱吉林.高等代数解题精粹.北京:中央民族大学

出版社,2002.

[3]华东师范大学数学系.数学分析(第二版).高

等教育出版社.

[4]刘嘉.矩阵相似及其应用.中国西部科

16/16

技,2010,(26)

[5]袁进.特征值及特征向量.高等数学研

究,2004,(02)

[6]张斌斌.矩阵特征值及特征向量研究.才

智,2010,(08)

[7]施劲松,刘剑平.矩阵特征值、特征向量确定.

大学数学,2003,(06).第19卷第6期

[8]汪庆丽.用矩阵初等变换求矩阵特征值及特征

向量.岳阳师范学院学报(自然科学

版),2001,(03)

[9]刘学鹏,王文省.关于实对称矩阵对角化问题[J]

聊城师院学报(自然科学版),2003,(03)

[10]李佩贞.矩阵对角化及相似变换矩阵.中山大

学学报论丛,2000,(04)

[11]朱靖红,朱永生.矩阵对角化相关问题[J]辽宁

师范大学学报(自然科学版),2005,(03).

[12]熊凯俊,李丽萍.矩阵多项式特征值、特征向

量简单求法.科协论坛(下半月),2008,(02)1007-3973

👁️ 阅读量:0