[P,D]=eig(A)
——计算出A的全部特征值和对应的特征向量.
其中,
D是对角矩阵,保存矩阵A的全部特征值;
P是满阵,
P的列向量构成对应于D的特征向量组。
用eig就可以了.特征多项式是poly.转为一般形式是poly2sym.
a=magic(5),%产生一个魔方阵;
[v,d]=eig(a),计算特征向量与特征值;
ap=poly(a),计算特征多项式;
aps=poly2sym(ap),计算一般形式;
结果:
a
=
17
24
1
8
15
23
5
7
14
16
4
6
13
20
22
10
12
19
21
3
11
18
25
2
9
v
=
-0.4472
0.0976
-0.6330
0.6780
-0.2619
-0.4472
0.3525
0.5895
0.3223
-0.1732
-0.4472
0.5501
-0.3915
-0.5501
0.3915
-0.4472
-0.3223
0.1732
-0.3525
-0.5895
-0.4472
-0.6780
0.2619
-0.0976
0.6330
d
=
65.0000
0
0
0
0
0
-21.2768
0
0
0
0
0
-13.1263
0
0
0
0
0
21.2768
0
0
0
0
0
13.1263
ap
=
1.0e+006
*
0.0000
-0.0001
-0.0006
0.0406
0.0780
-5.0700
aps
=
x^5
-
65*x^4
-
625*x^3
+
40625*x^2
+
78000*x
-
2721935523840001/536870912
用
[d,v]
=
eig(A)
就可以了
如:
>>
A=[1,2;3,4]
A
=
1
2
3
4
>>
[d,v]=eig(A)
d
=
-4216/5113
-250/601
671/1186
-1736/1909
v
=
-736/1977
0
0
1977/368
v
中是特征值,
d中是对应的特征向量
满意请采纳^_^