23.12.2012 Views

MATLAB Function Reference Volume 1: A - E - Bad Request

MATLAB Function Reference Volume 1: A - E - Bad Request

MATLAB Function Reference Volume 1: A - E - Bad Request

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

eigs<br />

2-478<br />

subplot(2,1,1)<br />

plot(dlm,'k+')<br />

hold on<br />

plot(d(ind(end:-1:end-5)),'ks')<br />

hold off<br />

legend('eigs(A)','eig(full(A))',3)<br />

set(gca,'XLim',[0.5 6.5])<br />

subplot(2,1,2)<br />

plot(dsm,'k+')<br />

hold on<br />

plot(d(ind(1:6)),'ks')<br />

hold off<br />

legend('eigs(A,6,''sm'')','eig(full(A))',2)<br />

set(gca,'XLim',[0.5 6.5])<br />

8<br />

7.95<br />

7.9<br />

7.85<br />

7.8<br />

0.2<br />

0.15<br />

0.1<br />

0.05<br />

0<br />

eigs(A)<br />

eig(full(A))<br />

1 2 3 4 5 6<br />

eigs(A,6,’SM’)<br />

eig(full(A))<br />

1 2 3 4 5 6<br />

However, the repeated eigenvalue at 4 must be handled more carefully. The<br />

call eigs(A,18,4.0) to compute 18 eigenvalues near 4.0 tries to find<br />

eigenvalues of A - 4.0*I. This involves divisions of the form 1/

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!