function simsensors cd Z:/Education/cvs/uio/postdoc/luminy/matlab set(gca,'FontSize',24); c =0.5; K=8; Nvalues = [64 512]; for N=Nvalues sprintf('N=%i\n',N) R = zeros(N,N); for (k=1:K) R(k,k) = 1; end momsR = (K/N) .* ones(1,8); L = N/c; X = randn(N,L); eigenvals = eig((1/L) .* X * X' * R); hist(eigenvals,50); title(''); xlabel(''); ylabel(''); ax = axis; axis([0 2 0 10]); switch N case 64 print -deps ../images/simsensors1.eps; case 512 print -deps ../images/simsensors2.eps; end end