HM30 and 31

HM30 and 31 - on plot(m,b'x 30.7 n=100 xv= rand(n-1,1...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
31.4 %%%%%%%%%%%%%%%%%%%%%%%%%%%%% m=5; A=triu(ones(m)); for i=1:1:m A(i,i)=0.1; end %%%%%%%%%%%%%%%%%%%%%%%%% min(svd(A)); %%%%%%%%%%%%%%%%%%%%%% (min(eig(A*A)))^(1/2); %%%%%%%%%%%%%%%%%%% n=30; a=zeros(n,1); b=zeros(n,1); for m=1:1:30 A=triu(ones(m)); for i=1:1:m A(i,i)=0.1; end a(m)=min(svd(A)); b(m)=(min(eig(A*A)))^(1/2); end m=1:1:30; plot(m,a, 's' ); hold
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 2
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: on plot(m,b, 'x' ); 30.7 n=100; xv= rand(n-1,1); A=diag(rand(n,1))+diag(xv,1)+diag(xv,-1) [V,D] = eig(A); m=1:1:30 plot(m,abs(V(1:100,m))) % T is the discrete Laplacian with entries 1,-2,1 n=100; xv= ones(n-1,1); A=diag(ones(n,1))+diag(xv,1)+diag(xv,-1); [V,D] = eig(A); m=1:1:30 plot(m,abs(V(1:100,m)))...
View Full Document

Page1 / 2

HM30 and 31 - on plot(m,b'x 30.7 n=100 xv= rand(n-1,1...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online