满意答案
ukobdy
2017.05.12
采纳率:54% 等级:10
已帮助:822人
A=[1 2 3 4 5 6 ; 1/2 1 2 3 4 5 ; 1/3 1/2 1 2 3 4 ; 1/4 1/3 1/2 1 2 3 ; 1/5 1/4 1/3 1/2 1 2 ; 1/6 1/5 1/4 1/3 1/2 1] ;
%A为矩阵;ep为精度要求;N为最大迭代次数;m为绝对值最大的特征值;u为对应最大特征值的特征向量.
N=0;
ep=1e-8;
n=length(A);
u=ones(n,1);
index=0;
k=0;m1=0;
while k<=N
v=A*u;
m=max(abs(v));
u=v/m
if abs(m-m1)
index=1;
break;
end
m1=m;
k=k+1;
end
m %特征值
u/norm(u) %特征向量
[vv,ll]=eig(A); %matlab求解的特征值和特征向量
[mm,ii]=max(abs(diag(ll)));
m_matlab=mm
v_matlab=vv(:,ii)
00分享举报