maxeigenvalue.m

用牛顿法解方程组

function Ray=maxeigenvalue(A,X0,tol,nmax)%请输入合适的初始值，本题使用相对误差。
X0=X0';
X1=X0+[1 10 100]';
niter=0;
while abs((R(A,X0)-R(A,X1))/R(A,X0)) >= tol && niter < nmax%R为内嵌函数。
    niter=niter+1;
    X=A*X0/max(X0);
    X1=X0;
    X0=X;
end
Ray=R(A,X0);%采用瑞利商来计算最大特征值。
niter
X0
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%下为计算结果
%maxeigenvalue([4 2 2;2 6 1;2 1 6],[1 1 1],1.e-6,1000)
%niter =5
%X0 =

%    7.4059
%    8.7030
%    8.7030
%ans =8.7016