求用matlab算矩阵的最大特征根和对应的特征向量

矩阵A=[0,1,1,0,1,1;0,0,1,1,1,1;0,0,0,1,1,1;1,0,0,0,1,1;0,0,0,0,0,1;0,0,0,0,0,0]求这个6阶矩阵的最大特征根和对应的特征向量

[V,D]=eig(A);
%V特征值，D特征向量；
tz=max(D);%最大特征值
[max_column, index_row] = max(D);%最大特征值所在位置
a=V(:,index_row(2))对应特征向量

假设特征值为λ，其对应的特征向量是a=【x1;x2;x3;x4】， 则a*a=λ*a (λ*i-a)*a=0（a!=0）（1） 则det（λ*i-a） =0 即【λ，λ-1，λ-1，λ；λ，λ，λ-1，λ-1；λ，λ，λ，λ-1；λ-1，λ，λ，λ】=0 解此方程组可得λ=1.3953 -0.4604 + 1.1393i -0.4604 - 1.1393i -0.4746（兄弟，你那个算的最大特征根λ=2.232好像是错的呀） 下面是将算的的λ=1.3953带入（1）解出a（有无数个解，只需求一个特殊解即可） 比如单位矩阵为a=【0.6256；0.5516；0.3213；0.4484】 用matlab实现 a=[0 1 1 0;0 0 1 1;0 0 0 1;1 0 0 0]; [m,n]=eig(a) 结果： m = 0.6256 0.6247 0.6247 0.1424 0.5516 0.0100 + 0.4243i 0.0100 - 0.4243i -0.6999 0.3213 -0.2976 + 0.2874i -0.2976 - 0.2874i 0.6323 0.4484 -0.1905 - 0.4714i -0.1905 + 0.4714i -0.3001 n = 1.3953 0 0 0 0 -0.4604 + 1.1393i 0 0 0 0 -0.4604 - 1.1393i 0 0 0 0 -0.4746