a=[1 1/5 1/7 3;5 1 1/7 7;7 7 1 9;1/3 1/7 1/9 1];[v,d]=eig(a);m=d(1,1)%最大特征值n=v(:,1)%最大特征值对应的特征矢量结果：m = 4.4420n = 0.1026 0.2958 0.9483 0.0528
