%exm0531_3.m

%指令eig與eigs的比較。
rand('state',1),A=rand(100,100)-0.5;
t0=clock;[V,D]=eig(A);T_full=etime(clock,t0)	%指令eig的運(yùn)作時間。
options.tol=1e-8;									%為eigs設(shè)定計算精度。
options.disp=0;									%使中間迭代結(jié)果不顯示。
t0=clock;[v,d]=eigs(A,1,'lr',options);		%計算最大實(shí)部特征值和特征向量。
T_part=etime(clock,t0)							%指令eigs的運(yùn)作時間。
[Dmr,k]=max(real(diag(D)));						%在eig求得的全部特征值中找最大實(shí)部的那個。
d,D(1,1)  

vk1=V(:,k+1);										%與d相同的特征向量應(yīng)是V的第k+1列。
vk1=vk1/norm(vk1);v=v/norm(v);					%向量長度歸一。
V_err=acos(norm(vk1'*v))*180/pi				%求復(fù)數(shù)向量之間的夾角（度）。
D_err=abs(D(k+1,k+1)-d)/abs(d)					%求兩個特征值間的相對誤差。