用MATLAB求矩阵特征值
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
用matlab求矩阵的特征值和特征向量
我要计算的矩阵：
1 3 5
1/3 1 3
1/5 1/3 1
[v,d]=eig(A);
A为你的矩阵，V为特征向量矩阵，D为特征值矩阵，然后对D求最大值即可得最大特征根！
[V,D] = EIG(X) produces a diagonal matrix D of eigenvalues and a full matrix V whose columns are the corresponding eigenvectors so that X*V = V*D.
V是特征向量，D是特征值
实例：
矩阵：
1 2/3 7/3 7/3
3/2 1 3/2 3/2
3/7 2/3 1 3/2
3/7 2/3 2/3 1
>> format rat
>> A=[1 2/3 7/3 7/3
3/2 1 3/2 3/2
3/7 2/3 1 3/2
3/7 2/3 2/3 1]
A =
1 2/3 7/3 7/3
3/2 1 3/2 3/2
3/7 2/3 1 3/2
3/7 2/3 2/3 1
>> [V,D]=eig(A)
V =
1793/2855 504/3235 - 146/235i 504/3235 + 146/235i 1990/4773
670/1079 -3527/5220 -3527/5220 -509/959
4350/11989 1160/4499 + 287/3868i 1160/4499 - 287/3868i -350/647
838/2819 181/3874 + 1179/4852i 181/3874 - 1179/4852i 1238/2467
D =
810/197 0 0 0
0 -93/4229 + 455/674i 0 0
0 0 -93/4229 - 455/674i 0
0 0 0 -149/2201
***************************************************************************************** 如何归一化求权重呢？
>> a=[1 3 5;1/3 1 3; 1/5 1/3 1]
a =
1.0000 3.0000 5.0000
0.3333 1.0000 3.0000
0.2000 0.3333 1.0000
>> [V,D]=eig(a)
V =
0.9161 0.9161 0.9161
0.3715 -0.1857 + 0.3217i -0.1857 - 0.3217i
0.1506 -0.0753 - 0.1304i -0.0753 + 0.1304i
D =
3.0385 0 0
0 -0.0193 + 0.3415i 0
0 0 -0.0193 - 0.3415i
**************************************************************************
>> a=[1 2 4 8 6 6 8;1/2 1 2 6 4 4 8;1/4 1/2 1 4 2 4 6;1/8 1/6 1/4 1 2 2 4;1/6 1/4 1/2 1/2 1 1 4;1/6 1/4 1/4 1/2 1 1 2;1/8 1/8 1/6 1/4 1/4 1/2 1]
a =
1.0000
2.0000 4.0000 8.0000 6.0000 6.0000 8.0000
0.5000 1.0000 2.0000 6.0000 4.0000 4.0000 8.0000
0.2500 0.5000 1.0000 4.0000 2.0000 4.0000 6.0000
0.1250 0.1667 0.2500 1.0000 2.0000 2.0000 4.0000
0.1667 0.2500 0.5000 0.5000 1.0000 1.0000 4.0000
0.1667 0.2500 0.2500 0.5000 1.0000 1.0000 2.0000
0.1250 0.1250 0.1667 0.2500 0.2500 0.5000 1.0000
>> rats(a)
ans =
1 2 4 8 6 6 8
1/2 1 2 6 4 4 8
1/4 1/2 1 4 2 4 6
1/8 1/6 1/4 1 2 2 4
1/6 1/4 1/2 1/2 1 1 4
1/6 1/4 1/4 1/2 1 1 2
1/8 1/8 1/6 1/4 1/4 1/2 1
>> [V,D]=eig(a)
V =
0.7884 0.8327 0.8327 0.8083 0.8083 -0.5119 + 0.3865i -0.5119 - 0.3865i
0.4894 0.3216 + 0.2636i 0.3216 - 0.2636i 0.1760 + 0.0792i 0.1760 - 0.0792i 0.6783 0.6783
0.3038 0.0883 + 0.2728i 0.0883 - 0.2728i -0.4630 + 0.1038i -0.4630 - 0.1038i -0.2011 - 0.2400i -0.2011 + 0.2400i