exm04031_3.m

MATLAB 6.5 从基础讲解其功能

rand('state',1),A=rand(100,100)-0.5;
t0=clock;[V,D]=eig(A);T_full=etime(clock,t0)%指令eig的运作时间。
options.tol=1e-8;				%为eigs设定计算精度。
options.disp=0;					%使中间迭代结果不显示。
t0=clock;[v,d]=eigs(A,1,'lr',options);%计算最大实部特征值和特征向量。
T_part=etime(clock,t0)			%指令eigs的运作时间。
[Dmr,k]=max(real(diag(D)));	%在eig求得的全部特征值中找最大实部的那个。
d,D(1,1)

vk1=V(:,k);							%与d相同的特征向量应是V的第k列。
vk1=vk1/norm(vk1);v=v/norm(v);	%向量长度归一。
V_err=acos(norm(vk1'*v))*180/pi	%求复数向量之间的夹角（度）。
D_err=abs(D(k,k)-d)/abs(d)	%求两个特征值间的相对误差。