?? exm04031_3.m
字號:
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) %求兩個特征值間的相對誤差。
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -