
次方怎么算
-公务员退休金
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
k0。
当特征值为1时,解方程0EAX,由齐次线性
方程组计算方法知:0EAX基础解系为
2
a=110
;所以矩阵A属于特征值1全部特征向量为
2
110k
,其中
2
k0。
当特征值为3时,解方程30EAX,由齐次线性方
程组计算方法知:30EAX基础解系为
3
a=011
,
所以矩阵A属于特征值3全部特征向量为
3
011k
,其
中
3
k0。
则由矩阵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,其中gx0,一定有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是0f2重根。
由定理:如果不可约多项式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