csdms-contrib/slepian_alpha

Switch branches/tags
Nothing to show
Fetching contributors…
Cannot retrieve contributors at this time
88 lines (71 sloc) 1.88 KB
 function sdwvals % SDWVALS % % Eigenvalue structure for various concentration regions. % Mixed-order ranking. % % Last modified by fjsimons-at-alum.mit.edu, 04/24/2009 TH=[10 20 30 40]; L=18; clf [ah,ha,H]=krijetem(subnum(2,2)); yls=[-0.1 1.1]; symbs={'o','x','s','+','v','*','^','d','<','>','p','h',... 'o','x','s','+','v','*','^','d'}; xmax=60; more off for index=1:length(TH) legsi{index}=sprintf('%s = %i%s ','\Theta',TH(index),str2mat(176)); SN(index)=floor(L/180*TH(index)); Nall(index)=(L+1)^2*(1-cos(TH(index)*pi/180))/2; [lrnk,mrnk,lval,VV,Vsum]=sdwelm(TH(index),L); ldubs=gamini(lval,~~mrnk+1); mrnk=gamini(mrnk,~~mrnk+1); axes(ah(index)) for ondi=1:length(ldubs) p(ondi,index)=plot(ondi,ldubs(ondi),symbs{mrnk(ondi)+1}); hold on end hold on plot(round([Nall(index) Nall(index)]),yls,'k:') plot([0 xmax],[0.5 0.5],'k:') plot([0 xmax],[0 0],'k:') plot([0 xmax],[1 1],'k:') set(ah(index),'xlim',[0 xmax],'ylim',yls,'xgrid','off','ygrid','off',... 'xtickl',[1 10:10:xmax],'xtick',[1 10:10:xmax],... 'ytick',[0:0.25:1]) [bh(index),th(index)]=boxtex('ur',ah(index),legsi{index},12); drawnow end axes(ah(4)) fb=fillbox([2 18 0.88 -0.05],'w'); for ondi=1:12 ypo=0+0.075*(ondi-1); pl(ondi,1)=plot(4,ypo,symbs{ondi}); hold on tl(ondi,1)=text(7,ypo,sprintf('m = %s %i','\pm',ondi-1),'FontS',8); end % Now make the plot beautiful movev(th,-.00) longticks(ah) set([p(~~p(:)) ; pl(~~pl(:))],'MarkerS',4,'MarkerF',grey,'MarkerE','k') axes(ha(1)) al(1)=ylabel('eigenvalue \lambda'); axes(ha(2)) al(2)=ylabel('eigenvalue \lambda'); axes(ha(2)) xl(1)=xlabel('rank'); axes(ha(4)) xl(2)=xlabel('rank'); nolabels(ha(3:4),2) nolabels(ah(1:2),1) serre(H',1/2,'down') serre(H,1/2,'across') for ind=1:4 xx(ind)=xtraxis(ah(ind),round(Nall(ind)),... {sprintf('N = %i',round(Nall(ind)))}); end longticks(xx) set([xl al],'FontS',13) set([ ah],'FontS',12) figdisp