From 21d03bc7cab88ceec46beee1f4d9fb44dd409129 Mon Sep 17 00:00:00 2001 From: gummif Date: Wed, 26 Mar 2014 23:45:25 -0400 Subject: [PATCH] removed files --- spectra_error_1d_plotload1.m | 188 ------------------- spectra_error_1d_plotload2.m | 255 -------------------------- spectra_error_1d_plotload3.m | 276 ---------------------------- spectra_error_1d_plotload31.m | 276 ---------------------------- spectra_error_1d_plotload4.m | 335 ---------------------------------- spectra_error_1d_plotload5.m | 140 -------------- spectra_error_2d_plotload4.m | 333 --------------------------------- spectra_error_2d_plotload5.m | 149 --------------- spectra_error_2d_plotload6.m | 217 ---------------------- 9 files changed, 2169 deletions(-) delete mode 100755 spectra_error_1d_plotload1.m delete mode 100755 spectra_error_1d_plotload2.m delete mode 100755 spectra_error_1d_plotload3.m delete mode 100755 spectra_error_1d_plotload31.m delete mode 100755 spectra_error_1d_plotload4.m delete mode 100755 spectra_error_1d_plotload5.m delete mode 100755 spectra_error_2d_plotload4.m delete mode 100755 spectra_error_2d_plotload5.m delete mode 100755 spectra_error_2d_plotload6.m diff --git a/spectra_error_1d_plotload1.m b/spectra_error_1d_plotload1.m deleted file mode 100755 index 3a04e39..0000000 --- a/spectra_error_1d_plotload1.m +++ /dev/null @@ -1,188 +0,0 @@ - -%========== error vs M ========= -set(0,'defaulttextinterpreter','latex'); -FIGCOLOR='-rgb'; % RGB,CMYK -FIG_SIZE(1)=560; %x (*1.8)<-(1,2) -FIG_SIZE(2)=420; %y -FIGLOC='./spectra_fig/'; -SAVEDLOC='./spectra_var/'; -FIGQUAL=''; -CSNAME='SpESO'; % Spectrum Estimation by Sparse Optimization -SAVEFIG=1; - - -COMPTYPE=1; - -if COMPTYPE==1 % '64', [32,16,8,4], [1,2,3,4] - SAVEDLOC='./spectra_var/if/'; - SIMN0t='64'; %8,16,64 - Mratio=8; %16,8,4 - RUNTYPEv=[1,2,3,4]; - %RUNTYPEv=[1]; -end - - -for RUNTYPE = RUNTYPEv; - % 1=w353, 2=w533 wavelet - % 3=f353, 4=f533 fourier - - % f_cs_res1{k0,k}=fhap; - % f_os{k0,k}=f; - % E_cs_res1{k0,k}=fhap; - % E_os{k0,k}=fspecf; - % k_os{1}=kvec; - - if Mratio==16 - FIGQ2='_16'; - M1=num2str(2066); - M2=num2str(1033); - elseif Mratio==8 - FIGQ2='_8'; % ''=m8 - M1=num2str(4132); - M2=num2str(2066); - elseif Mratio==4 - FIGQ2='_4'; - M1=num2str(8263); - M2=num2str(4132); - end - if RUNTYPE==1 - RUNSTR='_w353'; - %do1vec=[8]; (-3,-5/3); qmf=MakeONFilter('Coiflet',3); - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_w353.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_w353.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_w353.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_w353.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_w353_d2if.mat']); - elseif RUNTYPE==2 - RUNSTR='_w533'; - %do1vec=[8]; (-5/3,-3); qmf=MakeONFilter('Coiflet',3); - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_w533.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_w533.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_w533.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_w533.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_w533_d2if.mat']); - elseif RUNTYPE==3 - RUNSTR='_f353'; - %do1vec=[8]; (-5/3,-3); qmf=MakeONFilter('Coiflet',3); - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_f353.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_f353.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_f353.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_f353.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_f353_d2if.mat']); - elseif RUNTYPE==4 - RUNSTR='_f533'; - %do1vec=[8]; (-5/3,-3); qmf=MakeONFilter('Coiflet',3); - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_f533.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_f533.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_f533.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_f533.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_f533_d2if.mat']); - end - - - kvec=var1.k_os{1}; - %plot 1 spectrum - ee=var1.E_os{1,1}; - ee1=var1.E_cs_res1{1,1}; - figure,loglog(kvec,ee),hold all,loglog(kvec,ee1) - ee1b=var1b.E_cs_res1{1,1}; loglog(kvec,ee1b) - ee2=var2.E_cs_res1{1,1}; loglog(kvec,ee2) - ee3=var3.E_cs_res1{1,1}; loglog(kvec,ee3) - ee4=var4.E_cs_res1{1,1}; loglog(kvec,ee4) - - - %============================== - - N=length(var1.f_os{1}); - J=log2(N); - kvec=var1.k_os{1}; - jvec=7:J-1; - - [SIMN0,~]=size(var1.f_os); - err1=zeros(SIMN0,J-1); - err1b=err1;err2=err1; - err3=err1;err4=err1; - for j=jvec - wj=find(kvec>2^(j-1) & kvec<=2^j); - for i=1:SIMN0 - err1(i,j) =norm(log(var1.E_os{i,1}(wj))-log(var1.E_cs_res1{i,1}(wj)))/length(wj); - err1b(i,j)=norm(log(var1b.E_os{i,1}(wj))-log(var1b.E_cs_res1{i,1}(wj)))/length(wj); - err2(i,j) =norm(log(var2.E_os{i,1}(wj))-log(var2.E_cs_res1{i,1}(wj)))/length(wj); - err3(i,j) =norm(log(var3.E_os{i,1}(wj))-log(var3.E_cs_res1{i,1}(wj)))/length(wj); - err4(i,j) =norm(log(var4.E_os{i,1}(wj))-log(var4.E_cs_res1{i,1}(wj)))/length(wj); - end - end - err1=err1(:,jvec); - err1b=err1b(:,jvec); - err2=err2(:,jvec); - err3=err3(:,jvec); - err4=err4(:,jvec); - - %============================ - - figure; - set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); - - - % M/2 best - hh1=semilogy(jvec,mean(err1,1),'--','Color',[0.2 0.2 1],'LineWidth',0.8); - hold all - errorbar(jvec,mean(err1,1),std(err1,1,1),std(err1,1,1),'x','Color',[0.2 0.2 1],'LineWidth',0.8) - - % M best - hh1b=semilogy(jvec,mean(err1b,1),'--','Color',[0.2 0.2 1],'LineWidth',1.3); - errorbar(jvec,mean(err1b,1),std(err1b,1,1),std(err1b,1,1),'x','Color',[0.2 0.2 1],'LineWidth',0.8) - - % shannon - hh2=semilogy(jvec,mean(err2,1),':','Color',[1 0.2 0.2],'LineWidth',1.3); - errorbar(jvec,mean(err2,1),std(err2,1,1),std(err2,1,1),'x','Color',[1 0.2 0.2],'LineWidth',0.8) - - % LOMP - lstd=std(err3,1,1); - temp=mean(err3,1); - lstd(lstd>=mean(err3,1))=temp(lstd>=mean(err3,1))-1e-10; - hh3=semilogy(jvec,mean(err3,1),':','Color',[0 0.6 0],'LineWidth',0.8); - errorbar(jvec,mean(err3,1),lstd,std(err3,1,1),'x','Color',[0 0.6 0],'LineWidth',0.8) - - - % QOM - lstd=std(err4,1,1); - temp=mean(err4,1); - lstd(lstd>=mean(err4,1))=temp(lstd>=mean(err4,1))-1e-10; - hh4=semilogy(jvec,mean(err4,1),'-','Color',[0 0 0],'LineWidth',1.3); - errorbar(jvec,mean(err4,1),lstd,std(err4,1,1),'x','Color',[0 0 0],'LineWidth',0.8) - - % lines - semilogy([11.5,11.5],[1e-5,1e3],'-.','Color',[0 0 0],'LineWidth',0.6); - %semilogy([10.5,10.5],[1e-5,1e3],':','Color',[0 0 0],'LineWidth',0.8); - - grid on - axis tight - % \hspace{-15mm} $ \hspace{10mm} - ylabel('\hspace{-2mm} $ \| \log E-\log E^{\star} \|_{w_j} $','FontSize', 11); - xlabel('$j$','FontSize', 11); - set(gca,'Outerposition',[-0.08,0,1.1,1]) - set(gca,'XTick',jvec) - - ylim([9.5e-5 8e0]) - xlim([jvec(1)-0.2 jvec(end)+0.2]) - - - legend([hh1,hh1b,hh2,hh3,hh4],{'$M/2$-best','$M$-best','Shannon','LOMP',CSNAME}... - ,'interpreter', 'latex','Location','southeast'); - - - if SAVEFIG==1 - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,FIGQ2,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,FIGQ2,RUNSTR], '-eps', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'e1',FIGQUAL,FIGQ2,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'e1',FIGQUAL,FIGQ2,RUNSTR], '-eps', '-transparent', FIGCOLOR) - - end - -end - - - -% Copyright (C) 2014 Gudmundur Adalsteinsson -% See file LICENCE for licence and warranty details diff --git a/spectra_error_1d_plotload2.m b/spectra_error_1d_plotload2.m deleted file mode 100755 index ca18355..0000000 --- a/spectra_error_1d_plotload2.m +++ /dev/null @@ -1,255 +0,0 @@ - -%========== spectrum cases ========= -set(0,'defaulttextinterpreter','latex'); -FIGCOLOR='-rgb'; % RGB,CMYK -FIG_SIZE(1)=560; %x (*1.8)<-(1,2) -FIG_SIZE(2)=420; %y -FIGLOC='./spectra_fig/'; -FIGQUAL=''; -CSNAME='SpESO'; % Spectrum Estimation by Sparse Optimization -SAVEFIG=1; - - - - -COMPTYPEv=[1,10,11]; -%COMPTYPEv=[10,11]; - -for COMPTYPE=COMPTYPEv - - EXTRASTR=''; - if COMPTYPE==1 % '64', [32,16,8,4], [1,2,3,4] - SAVEDLOC='./spectra_var/if/'; - SIMN0t='64'; %8,16,64 - Mratio=8; %16,8,4 - RUNTYPEv=[1,2,3,4]; - %RUNTYPEv=[4]; - elseif COMPTYPE==10 % '16', 8, 4 +s2/s4 - SAVEDLOC='./spectra_var/if/'; - SIMN0t='16'; %8,16,64 - Mratio=8; %16,8,4 - RUNTYPEv=4; - EXTRASTR='s2'; - elseif COMPTYPE==11 - SAVEDLOC='./spectra_var/if/'; - SIMN0t='16'; %8,16,64 - Mratio=8; %16,8,4 - RUNTYPEv=4; - EXTRASTR='s4'; %(av2) - end - - - for RUNTYPE = RUNTYPEv - - % 1=w353, 2=w533 wavelet - % 3=f353, 4=f533 fourier - - % f_cs_res1{k0,k}=fhap; - % f_os{k0,k}=f; - % E_cs_res1{k0,k}=fhap; - % E_os{k0,k}=fspecf; - % k_os{1}=kvec; - - if Mratio==16 - FIGQ2='_16'; - M1=num2str(2066); - M2=num2str(1033); - elseif Mratio==8 - FIGQ2='_8'; % ''=m8 - M1=num2str(4132); - M2=num2str(2066); - elseif Mratio==4 - FIGQ2='_4'; - M1=num2str(8263); - M2=num2str(4132); - end - if RUNTYPE==1 - RUNSTR=['_w353',EXTRASTR]; - %var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_8_synth_284_1_31_w353.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - %var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_8_synth_284_1_1_w353.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); - n_good=41; % for QOM - n_bad=57; %24 - - elseif RUNTYPE==2 - RUNSTR=['_w533',EXTRASTR]; - % var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_8_synth_284_1_31_w533.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - % var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_8_synth_284_1_1_w533.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); - n_good=27; % for QOM - n_bad=34; %36 - - elseif RUNTYPE==3 - RUNSTR=['_f353',EXTRASTR]; - %var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_8_synth_284_1_31_f353.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - %var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_8_synth_284_1_1_f353.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); - n_good=8; % for QOM - n_bad=14; %28 - - elseif RUNTYPE==4 - RUNSTR=['_f533',EXTRASTR]; - % var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_8_synth_284_1_31_f533.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - % var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_8_synth_284_1_1_f533.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); - n_good=5; % for QOM - n_bad=50; %81 - if strcmp(RUNSTR,'_f533s2') - n_good=3; % for QOM - n_bad=4; - elseif strcmp(RUNSTR,'_f533s4') - n_good=2; % for QOM - n_bad=4; - end - end - - if RUNTYPE==1 - ylimit=[8e-16 4e-10]; - C1=1;C2=1;p1=-3;p2=-5/3; - elseif RUNTYPE==2 - ylimit=[8e-16 4e-10]*1000; - C1=1;C2=1e8;p1=-5/3;p2=-3; - elseif RUNTYPE==3 - ylimit=[8e-16 4e-10]*100; - C1=60;C2=60;p1=-3;p2=-5/3; - elseif RUNTYPE==4 - ylimit=[8e-16 4e-10]*100000; - C1=60;C2=60e8;p1=-5/3;p2=-3; - if COMPTYPE==10 - C1=60;C2=15*60e8;p1=-5/3;p2=-3; - elseif COMPTYPE==11 - C1=60;C2=6.5*60e8;p1=-5/3;p2=-3; - end - end - - - %============================== - - N=length(var1b.f_os{1}); - J=log2(N); - kvec=var1b.k_os{1}; - [SIMN0,~]=size(var1b.f_os); - - - % find best/worst cases -% E_av2=var1b.E_cs_res1{1,1}*0; -% E_avo=E_av2*0; -% af=@(x)log(x);afi=@(x)exp(x); %log -% %af=@(x)x;afi=af; % linear -% for i=1:SIMN0 % average spectra -% E_av2=E_av2+af(var4.E_cs_res1{i,1}); -% E_avo=E_avo+af(var1b.E_os{i,1}); -% end -% -% E_av2=afi(E_av2/SIMN0); -% E_avo=afi(E_avo/SIMN0); -% -% -% wj=find(kvec>2^9 & kvec<=2^(J-2)); -% for i=1:SIMN0 -% d{i}=log(var4.E_cs_res1{i,1}(wj)./E_av2(wj)); % E_av{4} -% Vd(i)=var(d{i}); % a smooth spectrum (E_av*C) has 0 error -% -% do{i}=log(var4.E_cs_res1{i,1}(wj)./E_avo(wj)); -% Vdo(i)=norm(do{i}); % compare to org norm -% end -% -% [~,Vdi]=sort(Vd); -% [~,Vdoi]=sort(Vdo); -% -% for i=1:5 -% figure,set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); -% loglog(kvec,var1b.E_os{1,1}) -% hold all -% loglog(kvec,var4.E_cs_res1{Vdi(i),1}) -% loglog(kvec,var4.E_cs_res1{Vdoi(i),1}) -% ylim(ylimit) -% xlim([50 N/2]) -% title(['best ',num2str(i)]) -% end -% for i=1:5 -% figure,set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); -% loglog(kvec,var1b.E_os{1,1}) -% hold all -% loglog(kvec,var4.E_cs_res1{Vdi(end+1-i),1}) -% loglog(kvec,var4.E_cs_res1{Vdoi(end+1-i),1}) -% ylim(ylimit) -% xlim([50 N/2]) -% title(['worst ',num2str(i)]) -% end -% Vdoi(1) -% Vdoi(end) -% -% return - - figure; - set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); - - % original - hh0=loglog(kvec,var1b.E_os{1,1},'-','Color',[0.7 0 0.7],'LineWidth',1.4); - hold all - - % M-best - hh1b=loglog(kvec,var1b.E_cs_res1{1,1},'-','Color',[0.2 0.2 1],'LineWidth',0.7); - - % Shannon - hh2=loglog(kvec,var2.E_cs_res1{1,1},':','Color',[1 0.2 0.2],'LineWidth',1.3); - - % QOM bad - hh4b=loglog(kvec,var4.E_cs_res1{n_bad,1},'--','Color',[0.6 0.6 0.6],'LineWidth',0.8); - - % QOM good - hh4=loglog(kvec,var4.E_cs_res1{n_good,1},'-','Color',[0 0 0],'LineWidth',0.8); - - - if COMPTYPE==1 - loglog([100,1400],C1*3e-3*[100,1400].^p1,'k--','LineWidth',0.8) - loglog([700,1e4],C2*3e-7*[700,1e4].^p2,'k--','LineWidth',0.8) - elseif COMPTYPE==10 - loglog([100,1e3*8*1.2],C1*3e-3*[100,1e3*8*1.2].^p1,'k--','LineWidth',0.8) - loglog([1e3*8/1.2,1.5e4],C2*3e-7*[1e3*8/1.2,1.5e4].^p2,'k--','LineWidth',0.8) - elseif COMPTYPE==11 - loglog([100,1.04e3*4*1.2],C1*3e-3*[100,1.043e3*4*1.2].^p1,'k--','LineWidth',0.8) - loglog([1.03e3*4/1.2,1.5e4],C2*3e-7*[1.03e3*4/1.2,1.5e4].^p2,'k--','LineWidth',0.8) - end - - - axis tight - % \hspace{-15mm} $ \hspace{10mm} - ylabel('\hspace{-2mm} $ E(k)$','FontSize', 11); - xlabel('$k$','FontSize', 11); - %set(gca,'Outerposition',[-0.08,0,1.1,1]) %posx,posy,widx,widy - %set(gca,'XTick',jvec) - - ylim(ylimit) - xlim([50 N/2]) - - - legend([hh0,hh1b,hh2,hh4,hh4b],{'original','$M$-best','Shannon',[CSNAME,' best'],[CSNAME,' worst']}... - ,'interpreter', 'latex','Location','southwest'); - - - if SAVEFIG==1 - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,FIGQ2,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,FIGQ2,RUNSTR], '-eps', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'c1',FIGQUAL,FIGQ2,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'c1',FIGQUAL,FIGQ2,RUNSTR], '-eps', '-transparent', FIGCOLOR) - - end - - - end - -end - - -% Copyright (C) 2014 Gudmundur Adalsteinsson -% See file LICENCE for licence and warranty details diff --git a/spectra_error_1d_plotload3.m b/spectra_error_1d_plotload3.m deleted file mode 100755 index 51c9bf2..0000000 --- a/spectra_error_1d_plotload3.m +++ /dev/null @@ -1,276 +0,0 @@ - -%========== slope errors table ========= -set(0,'defaulttextinterpreter','latex'); -FIGCOLOR='-rgb'; % RGB,CMYK -FIG_SIZE(1)=560; %x (*1.8)<-(1,2) -FIG_SIZE(2)=420; %y -FIGLOC='./spectra_fig/'; -FIGQUAL='_load3'; -CSNAME='SpESO'; % Spectrum Estimation by Sparse Optimization -SAVEFIG=1; -FIGQ2='_mv'; % ''=m8 - - - -COMPTYPEv=1; - -for COMPTYPE=COMPTYPEv - -if COMPTYPE==1 % '64', [32,16,8,4], [1,2,3,4] - SAVEDLOC='./spectra_var/if/'; - SIMN0t='64'; %8,16,64 - Mrvec=[16,8,4]; - RUNTYPEv=[1,2,3,4]; - %RUNTYPEv=[1]; -end - - - -for RUNTYPE = RUNTYPEv; %1,2,3,4 - % 1=w353, 2=w533 wavelet - % 3=f353, 4=f533 fourier - - %Mrvec=8; - rr=1; - for Mratio=Mrvec - - if Mratio==16 - M1=num2str(2066); - M2=num2str(1033); - elseif Mratio==8 - M1=num2str(4132); - M2=num2str(2066); - elseif Mratio==4 - M1=num2str(8263); - M2=num2str(4132); - end - if RUNTYPE==1 - RUNSTR='_w353'; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_w353.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_w353.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_w353.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_w353.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_w353_d2if.mat']); - elseif RUNTYPE==2 - RUNSTR='_w533'; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_w533.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_w533.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_w533.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_w533.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_w533_d2if.mat']); - elseif RUNTYPE==3 - RUNSTR='_f353'; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_f353.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_f353.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_f353.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_f353.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_f353_d2if.mat']); - elseif RUNTYPE==4 - RUNSTR='_f533'; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_f533.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_f533.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_f533.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_f533.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_f533_d2if.mat']); - end - - if RUNTYPE==1 - ylimit=[8e-16 4e-10]; - C1=1;C2=1;p1=-3;p2=-5/3; - elseif RUNTYPE==2 - ylimit=[8e-16 4e-10]*1000; - C1=1;C2=1e8;p1=-5/3;p2=-3; - elseif RUNTYPE==3 - ylimit=[8e-16 4e-10]*100; - C1=60;C2=60;p1=-3;p2=-5/3; - elseif RUNTYPE==4 - ylimit=[8e-16 4e-10]*100000; - C1=60;C2=60e8;p1=-5/3;p2=-3; - end - - - N=length(var1.f_os{1}); - J=log2(N); - kvec=var1.k_os{1}; - % jvec=J-5:J-3; - % kmax=N/4; - jvec=J-4:J-2; - kmin=N/32; - inguess=[5/3,3e-6]; - - - [SIMN0,~]=size(var1.f_os); - err=zeros(4,J-1); - orgslope=err; - % err1b=err1;err2=err1; - % err3=err1;err4=err1; - - - E_av(1:6)={var1.E_cs_res1{1,1}*0}; - %E_av{4}=E_av{4}+1e20; % minimum - af=@(x)log(x);afi=@(x)exp(x); %log - %af=@(x)x;afi=af; % linear - % average spectra - for i=1:SIMN0 - E_av{1}=E_av{1}+af(var1.E_cs_res1{i,1}); - E_av{2}=E_av{2}+af(var1b.E_cs_res1{i,1}); - E_av{3}=E_av{3}+af(var3.E_cs_res1{i,1}); - E_av{4}=E_av{4}+af(var4.E_cs_res1{i,1}); - %E_av{4}=min([E_av{4}';var4.E_cs_res1{i,1}'])'; % minimum - E_av{5}=E_av{5}+af(var2.E_cs_res1{i,1}); - E_av{6}=E_av{6}+af(var1.E_os{i,1}); - end - %for i=[1,2,3,5] % minimum - for i=1:6 - E_av{i}=afi(E_av{i}/SIMN0); - end - - %E_avM{rr}=E_av; - - % plot all - figure - loglog(kvec,var4.E_os{1,1},'r') - hold all - for i=1:SIMN0 - loglog(kvec,var4.E_cs_res1{i,1},'Color',[0.7 0.7 0.7]) - end - loglog(kvec,var4.E_os{1,1},'r') - loglog(kvec,E_av{4},'k') - - % for original - wj=find(kvec>kmin & kvec<2^jvec(end)); - ys=[E_av{6}(wj(1)),E_av{6}(wj(end))]; % a line through range - ks=kvec([wj(1),wj(end)]); - inguesso(1)=-log(ys(1)/ys(2))/log(ks(1)/ks(2)); % exponent - inguesso(2)=ys(1)/ks(1)^(-inguesso(1)); % constant - - for j=jvec - %wj=find(kvec>2^j & kvec<=kmax); - wj=find(kvec>kmin & kvec<2^j); - for i=1:5 - - ys=[E_av{i}(wj(1)),E_av{i}(wj(end))]; - ks=kvec([wj(1),wj(end)]); - inguess(1)=-log(ys(1)/ys(2))/log(ks(1)/ks(2)); - inguess(2)=ys(1)/ks(1)^(-inguess(1)); - %loglog([100,1e4],inguess(2)*[100,1e4].^(-inguess(1)),'k--','LineWidth',0.8) - - %fitscale=@(t) sqrt(t); - fitscale=@(t) log(t); - s1=exp_fit(E_av{i}(wj),kvec(wj)',fitscale,inguess,'log'); - s1o=exp_fit(E_av{6}(wj),kvec(wj)',fitscale,inguesso,'log'); - err(i,j)=s1o(1)-s1(1); - orgslope(i,j)=s1o(1); - if i==4 && j==jvec(end) - loglog([100,1e4],s1(2)*[100,1e4].^(-s1(1)),'b','LineWidth',2) - end - - end - end - iorder=[4,5,1,2,3]; - err=err(iorder,jvec); - orgslope=orgslope(iorder,jvec); - - % Mratio <-> rr - errM{rr}=err; - orgM{rr}=orgslope; - rr=rr+1; - end - - % special [16,8,4] at jvec(end) - em16=errM{1}(:,end); - em8=errM{2}(:,end); - em8(iorder==5)=errM{2}(iorder==5,1); %shannon @ kvec=2048 - em4=errM{3}(:,end); - em4(iorder==5)=errM{3}(iorder==5,2); %shannon @ kvec=4096 - err=[em16,em8,em4]; - orgslope=[orgM{1}(:,end),orgM{2}(:,end),orgM{3}(:,end)]; - - %tablerr=flipud(err')./flipud(orgslope')*100; % relative percentage - tablerr=flipud(err'); %difference - - disp(tablerr) - - - %============================ - -% figure; -% set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); -% -% -% % M/2 best -% hh1=plot(Mrvec,err(1,:),'--','Color',[0.2 0.2 1],'LineWidth',0.8); -% hold all -% % errorbar(jvec,mean(err1,1),std(err1,1,1),std(err1,1,1),'x','Color',[0.2 0.2 1],'LineWidth',0.8) -% -% % M best -% hh1b=plot(Mrvec,err(2,:),'--','Color',[0.2 0.2 1],'LineWidth',1.4); -% % errorbar(jvec,mean(err1b,1),std(err1b,1,1),std(err1b,1,1),'x','Color',[0.2 0.2 1],'LineWidth',0.8) -% % -% % % shannon -% % % hh2=semilogy(jvec,mean(err2,1),':','Color',[1 0.2 0.2],'LineWidth',1.3); -% % % errorbar(jvec,mean(err2,1),std(err2,1,1),std(err2,1,1),'x','Color',[1 0.2 0.2],'LineWidth',0.8) -% % -% % % LOMP -% % lstd=std(err3,1,1); -% % temp=mean(err3,1); -% % lstd(lstd>=mean(err3,1))=temp(lstd>=mean(err3,1))-1e-10; -% hh3=plot(Mrvec,err(3,:),':','Color',[0 0.6 0],'LineWidth',0.8); -% % errorbar(jvec,mean(err3,1),lstd,std(err3,1,1),'x','Color',[0 0.6 0],'LineWidth',0.8) -% % -% % -% % % QOM -% % lstd=std(err4,1,1); -% % temp=mean(err4,1); -% % lstd(lstd>=mean(err4,1))=temp(lstd>=mean(err4,1))-1e-10; -% hh4=plot(Mrvec,err(4,:),'-','Color',[0 0 0],'LineWidth',1.3); -% % errorbar(jvec,mean(err4,1),lstd,std(err4,1,1),'x','Color',[0 0 0],'LineWidth',0.8) -% -% % lines -% %semilogy([11.5,11.5],[1e-5,1e3],'-.','Color',[0 0 0],'LineWidth',0.6); -% %semilogy([10.5,10.5],[1e-5,1e3],':','Color',[0 0 0],'LineWidth',0.8); -% -% grid on -% axis tight -% % \hspace{-15mm} $ \hspace{10mm} -% ylabel('\hspace{-2mm} $ s-s^{\star} $','FontSize', 11); -% xlabel('$N/M$','FontSize', 11); -% %set(gca,'Outerposition',[-0.08,0,1.1,1]) -% set(gca,'XTick',fliplr(Mrvec)) -% -% ylim([-4 4]) -% xlim([Mrvec(end)-0.2 Mrvec(1)+0.2]) -% -% legend([hh1,hh1b,hh3,hh4],{'$M/2$-best','$M$-best','LOMP',CSNAME}... -% ,'interpreter', 'latex','Location','best'); - - - if SAVEFIG==1 - % not used - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,FIGQ2,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,FIGQ2,RUNSTR], '-eps', '-transparent', FIGCOLOR) - % table - latextable(tablerr,... - 'format','%5.2f',... % %4.0f - 'Horiz',{CSNAME,'Shannon','$M/2$-best','$M$-best','LOMP'},... - 'Vert',{' ','4';RUNSTR,'8';' ','16'},... - 'name',[FIGLOC,'table1d',FIGQUAL,FIGQ2,RUNSTR,'.tex']) - - end - - -end - -end - -% ================================================================ - - - - - - - - -% Copyright (C) 2014 Gudmundur Adalsteinsson -% See file LICENCE for licence and warranty details diff --git a/spectra_error_1d_plotload31.m b/spectra_error_1d_plotload31.m deleted file mode 100755 index ca41516..0000000 --- a/spectra_error_1d_plotload31.m +++ /dev/null @@ -1,276 +0,0 @@ - -%========== slope errors table (slope range 1 (128-1000)) ========= -set(0,'defaulttextinterpreter','latex'); -FIGCOLOR='-cmyk'; % RGB,CMYK -FIG_SIZE(1)=560; %x (*1.8)<-(1,2) -FIG_SIZE(2)=420; %y -FIGLOC='./spectra_fig/'; -FIGQUAL='_load3'; -CSNAME='SpESO'; % Spectrum Estimation by Sparse Optimization -SAVEFIG=1; -FIGQ2='_s1'; % ''=m8 - - - -COMPTYPEv=1; - -for COMPTYPE=COMPTYPEv - -if COMPTYPE==1 % '64', [32,16,8,4], [1,2,3,4] - SAVEDLOC='./spectra_var/if/'; - SIMN0t='64'; %8,16,64 - Mrvec=[16,8,4]; - RUNTYPEv=[1,2,3,4]; - %RUNTYPEv=[1]; -end - - - -for RUNTYPE = RUNTYPEv; %1,2,3,4 - % 1=w353, 2=w533 wavelet - % 3=f353, 4=f533 fourier - - %Mrvec=8; - rr=1; - for Mratio=Mrvec - - if Mratio==16 - M1=num2str(2066); - M2=num2str(1033); - elseif Mratio==8 - M1=num2str(4132); - M2=num2str(2066); - elseif Mratio==4 - M1=num2str(8263); - M2=num2str(4132); - end - if RUNTYPE==1 - RUNSTR='_w353'; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_w353.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_w353.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_w353.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_w353.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_w353_d2if.mat']); - elseif RUNTYPE==2 - RUNSTR='_w533'; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_w533.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_w533.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_w533.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_w533.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_w533_d2if.mat']); - elseif RUNTYPE==3 - RUNSTR='_f353'; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_f353.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_f353.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_f353.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_f353.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_f353_d2if.mat']); - elseif RUNTYPE==4 - RUNSTR='_f533'; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31_f533.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312_f533.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32_f533.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1_f533.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82_f533_d2if.mat']); - end - - if RUNTYPE==1 - ylimit=[8e-16 4e-10]; - C1=1;C2=1;p1=-3;p2=-5/3; - elseif RUNTYPE==2 - ylimit=[8e-16 4e-10]*1000; - C1=1;C2=1e8;p1=-5/3;p2=-3; - elseif RUNTYPE==3 - ylimit=[8e-16 4e-10]*100; - C1=60;C2=60;p1=-3;p2=-5/3; - elseif RUNTYPE==4 - ylimit=[8e-16 4e-10]*100000; - C1=60;C2=60e8;p1=-5/3;p2=-3; - end - - - N=length(var1.f_os{1}); - J=log2(N); - kvec=var1.k_os{1}; - % jvec=J-5:J-3; - % kmax=N/4; - jvec=J-5; - kmin=2^7; - %inguess=[5/3,3e-6]; - - - [SIMN0,~]=size(var1.f_os); - err=zeros(4,J-1); - orgslope=err; - % err1b=err1;err2=err1; - % err3=err1;err4=err1; - - - E_av(1:6)={var1.E_cs_res1{1,1}*0}; - %E_av{4}=E_av{4}+1e20; % minimum - af=@(x)log(x);afi=@(x)exp(x); %log - %af=@(x)x;afi=af; % linear - % average spectra - for i=1:SIMN0 - E_av{1}=E_av{1}+af(var1.E_cs_res1{i,1}); - E_av{2}=E_av{2}+af(var1b.E_cs_res1{i,1}); - E_av{3}=E_av{3}+af(var3.E_cs_res1{i,1}); - E_av{4}=E_av{4}+af(var4.E_cs_res1{i,1}); - %E_av{4}=min([E_av{4}';var4.E_cs_res1{i,1}'])'; % minimum - E_av{5}=E_av{5}+af(var2.E_cs_res1{i,1}); - E_av{6}=E_av{6}+af(var1.E_os{i,1}); - end - %for i=[1,2,3,5] % minimum - for i=1:6 - E_av{i}=afi(E_av{i}/SIMN0); - end - - %E_avM{rr}=E_av; - - % plot all - figure - loglog(kvec,var4.E_os{1,1},'r') - hold all - for i=1:SIMN0 - loglog(kvec,var4.E_cs_res1{i,1},'Color',[0.7 0.7 0.7]) - end - loglog(kvec,var4.E_os{1,1},'r') - loglog(kvec,E_av{4},'k') - - % for original - wj=find(kvec>kmin & kvec<2^jvec(end)); - ys=[E_av{6}(wj(1)),E_av{6}(wj(end))]; % a line through range - ks=kvec([wj(1),wj(end)]); - inguesso(1)=-log(ys(1)/ys(2))/log(ks(1)/ks(2)); % exponent - inguesso(2)=ys(1)/ks(1)^(-inguesso(1)); % constant - - for j=jvec - %wj=find(kvec>2^j & kvec<=kmax); - wj=find(kvec>kmin & kvec<2^j); - for i=1:5 - - ys=[E_av{i}(wj(1)),E_av{i}(wj(end))]; - ks=kvec([wj(1),wj(end)]); - inguess(1)=-log(ys(1)/ys(2))/log(ks(1)/ks(2)); - inguess(2)=ys(1)/ks(1)^(-inguess(1)); - %loglog([100,1e4],inguess(2)*[100,1e4].^(-inguess(1)),'k--','LineWidth',0.8) - - %fitscale=@(t) sqrt(t); - fitscale=@(t) log(t); - s1=exp_fit(E_av{i}(wj),kvec(wj)',fitscale,inguess,'log'); - s1o=exp_fit(E_av{6}(wj),kvec(wj)',fitscale,inguesso,'log'); - err(i,j)=s1o(1)-s1(1); - orgslope(i,j)=s1o(1); - if i==4 && j==jvec(end) - loglog([128,1e3],s1(2)*[128,1e3].^(-s1(1)),'b','LineWidth',2) - end - - end - end - iorder=[4,5,1,2,3]; - err=err(iorder,jvec); - orgslope=orgslope(iorder,jvec); - - % Mratio <-> rr - errM{rr}=err; - orgM{rr}=orgslope; - rr=rr+1; - end - - % special [16,8,4] at jvec(end) - em16=errM{1}(:,end); - em8=errM{2}(:,end); - %em8(iorder==5)=errM{2}(iorder==5,1); %shannon @ kvec=2048 - em4=errM{3}(:,end); - %em4(iorder==5)=errM{3}(iorder==5,2); %shannon @ kvec=4096 - err=[em16,em8,em4]; - orgslope=[orgM{1}(:,end),orgM{2}(:,end),orgM{3}(:,end)]; - - %tablerr=flipud(err')./flipud(orgslope')*100; % relative percentage - tablerr=flipud(err'); %difference - - disp(tablerr) - - - %============================ - -% figure; -% set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); -% -% -% % M/2 best -% hh1=plot(Mrvec,err(1,:),'--','Color',[0.2 0.2 1],'LineWidth',0.8); -% hold all -% % errorbar(jvec,mean(err1,1),std(err1,1,1),std(err1,1,1),'x','Color',[0.2 0.2 1],'LineWidth',0.8) -% -% % M best -% hh1b=plot(Mrvec,err(2,:),'--','Color',[0.2 0.2 1],'LineWidth',1.4); -% % errorbar(jvec,mean(err1b,1),std(err1b,1,1),std(err1b,1,1),'x','Color',[0.2 0.2 1],'LineWidth',0.8) -% % -% % % shannon -% % % hh2=semilogy(jvec,mean(err2,1),':','Color',[1 0.2 0.2],'LineWidth',1.3); -% % % errorbar(jvec,mean(err2,1),std(err2,1,1),std(err2,1,1),'x','Color',[1 0.2 0.2],'LineWidth',0.8) -% % -% % % LOMP -% % lstd=std(err3,1,1); -% % temp=mean(err3,1); -% % lstd(lstd>=mean(err3,1))=temp(lstd>=mean(err3,1))-1e-10; -% hh3=plot(Mrvec,err(3,:),':','Color',[0 0.6 0],'LineWidth',0.8); -% % errorbar(jvec,mean(err3,1),lstd,std(err3,1,1),'x','Color',[0 0.6 0],'LineWidth',0.8) -% % -% % -% % % QOM -% % lstd=std(err4,1,1); -% % temp=mean(err4,1); -% % lstd(lstd>=mean(err4,1))=temp(lstd>=mean(err4,1))-1e-10; -% hh4=plot(Mrvec,err(4,:),'-','Color',[0 0 0],'LineWidth',1.3); -% % errorbar(jvec,mean(err4,1),lstd,std(err4,1,1),'x','Color',[0 0 0],'LineWidth',0.8) -% -% % lines -% %semilogy([11.5,11.5],[1e-5,1e3],'-.','Color',[0 0 0],'LineWidth',0.6); -% %semilogy([10.5,10.5],[1e-5,1e3],':','Color',[0 0 0],'LineWidth',0.8); -% -% grid on -% axis tight -% % \hspace{-15mm} $ \hspace{10mm} -% ylabel('\hspace{-2mm} $ s-s^{\star} $','FontSize', 11); -% xlabel('$N/M$','FontSize', 11); -% %set(gca,'Outerposition',[-0.08,0,1.1,1]) -% set(gca,'XTick',fliplr(Mrvec)) -% -% ylim([-4 4]) -% xlim([Mrvec(end)-0.2 Mrvec(1)+0.2]) -% -% legend([hh1,hh1b,hh3,hh4],{'$M/2$-best','$M$-best','LOMP',CSNAME}... -% ,'interpreter', 'latex','Location','best'); - - - if SAVEFIG==1 - % not used - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,FIGQ2,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,FIGQ2,RUNSTR], '-eps', '-transparent', FIGCOLOR) - % table - latextable(tablerr,... - 'format','%5.2f',... % %4.0f - 'Horiz',{CSNAME,'Shannon','$M/2$-best','$M$-best','LOMP'},... - 'Vert',{' ','4';RUNSTR,'8';' ','16'},... - 'name',[FIGLOC,'table1d',FIGQUAL,FIGQ2,RUNSTR,'.tex']) - - end - - -end - -end - -% ================================================================ - - - - - - - - -% Copyright (C) 2014 Gudmundur Adalsteinsson -% See file LICENCE for licence and warranty details diff --git a/spectra_error_1d_plotload4.m b/spectra_error_1d_plotload4.m deleted file mode 100755 index 1c8b2a0..0000000 --- a/spectra_error_1d_plotload4.m +++ /dev/null @@ -1,335 +0,0 @@ - -%========== averaged spectra ========= -set(0,'defaulttextinterpreter','latex'); -FIGCOLOR='-rgb'; % RGB,CMYK -FIG_SIZE(1)=560; %x (*1.8)<-(1,2) -FIG_SIZE(2)=420; %y -FIGLOC='./spectra_fig/'; -FIGQUAL=''; -CSNAME='SpESO'; % Spectrum Estimation by Sparse Optimization -SAVEFIG=1; - -COMPTYPEv=[1,2,3,4,10,11]; -%COMPTYPEv=[4]; - -for COMPTYPE=COMPTYPEv - - EXTRASTR=''; - if COMPTYPE==1 % '64', [32,16,8,4], [1,2,3,4] - SAVEDLOC='./spectra_var/if/'; - SIMN0t='64'; %8,16,64 - Mratio=4; %16,8,4 - RUNTYPEv=[1,2,3,4,10]; - elseif COMPTYPE==2 - SAVEDLOC='./spectra_var/if/'; - SIMN0t='64'; %8,16,64 - Mratio=8; %16,8,4 - RUNTYPEv=[1,2,3,4,10]; - elseif COMPTYPE==3 - SAVEDLOC='./spectra_var/if/'; - SIMN0t='64'; %8,16,64 - Mratio=16; %16,8,4 - RUNTYPEv=[1,2,3,4,10]; - elseif COMPTYPE==4 - SAVEDLOC='./spectra_var/if/'; - SIMN0t='64'; %8,16,64 - Mratio=32; %16,8,4 - RUNTYPEv=[1,2,3,4]; %not with 10 - elseif COMPTYPE==10 % '16', 8, 4 +s2/s4 - SAVEDLOC='./spectra_var/if/'; - SIMN0t='16'; %8,16,64 - Mratio=8; %16,8,4 - RUNTYPEv=4; - EXTRASTR='s2'; - elseif COMPTYPE==11 - SAVEDLOC='./spectra_var/if/'; - SIMN0t='16'; %8,16,64 - Mratio=8; %16,8,4 - RUNTYPEv=4; - EXTRASTR='s4'; %(av2) - end - - %SAVEDLOC='./spectra_var/'; - % SAVEDLOC='./spectra_var/if/'; - % SIMN0t='16'; %8,16,64 - % Mratio=8; %16,8,4 - - - % 1=w353, 2=w533 wavelet - % 3=f353, 4=f533 fourier - % 10=rhh - %RUNTYPEv=10; %1,2,3,4,5 - for RUNTYPE = RUNTYPEv - - if Mratio==32 - AVM='_32'; - M1=num2str(1033); - M2=num2str(517); - elseif Mratio==16 - AVM='_16'; - M1=num2str(2066); - M2=num2str(1033); - elseif Mratio==8 - AVM='_8'; - M1=num2str(4132); - M2=num2str(2066); - elseif Mratio==4 - AVM='_4'; - M1=num2str(8263); - M2=num2str(4132); - end - if RUNTYPE==1 - RUNSTR=['_w353',EXTRASTR]; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); - elseif RUNTYPE==2 - RUNSTR=['_w533',EXTRASTR]; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); - elseif RUNTYPE==3 - RUNSTR=['_f353',EXTRASTR]; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); %d2av2d10 - elseif RUNTYPE==4 - RUNSTR=['_f533',EXTRASTR]; %s4 (av2),s2 - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); %d2av2d10 - elseif RUNTYPE==10 - RUNSTR=['_',EXTRASTR]; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_rhh_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_rhh_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_rhh_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_rhh_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_rhh_284_1_82',RUNSTR,'d2if.mat']); - RUNSTR=['_rhh',EXTRASTR]; - end - - if RUNTYPE==1 - ylimit=[8e-16 4e-10]; - C1=1;C2=1;p1=-3;p2=-5/3; - elseif RUNTYPE==2 - ylimit=[8e-16 4e-10]*1000; - C1=1;C2=1e8;p1=-5/3;p2=-3; - elseif RUNTYPE==3 - ylimit=[8e-16 4e-10]*100; - C1=60;C2=60;p1=-3;p2=-5/3; - elseif RUNTYPE==4 - ylimit=[8e-16 4e-10]*100000; - C1=60;C2=60e8;p1=-5/3;p2=-3; - elseif RUNTYPE==10 - ylimit=[8e-12 4e-6]*100000; - C1=60;C2=60e8;p1=-5/3;p2=-3; - end - - - N=length(var1.f_os{1}); - J=log2(N); - kvec=var1.k_os{1}; - [SIMN0,~]=size(var1.f_os); - - E_av(1:6)={var1.E_cs_res1{1,1}*0}; - E_min=E_av{4}+1e20; % minimum - af=@(x)log(x);afi=@(x)exp(x); %log - %af=@(x)x;afi=af; % linear - % average spectra - for i=1:SIMN0 - E_av{1}=E_av{1}+af(var1.E_cs_res1{i,1}); - E_av{2}=E_av{2}+af(var1b.E_cs_res1{i,1}); - E_av{3}=E_av{3}+af(var3.E_cs_res1{i,1}); - E_av{4}=E_av{4}+af(var4.E_cs_res1{i,1}); - E_min=min([E_min';var4.E_cs_res1{i,1}'])'; % minimum - - E_av{5}=E_av{5}+af(var2.E_cs_res1{i,1}); - E_av{6}=E_av{6}+af(var1.E_os{i,1}); - end - - for i=1:6 - E_av{i}=afi(E_av{i}/SIMN0); - end - - - % logscale difference from average - E_av2=E_av{4}; - nrem_max=round(0.3*SIMN0); - nrem_max=0; % org - nremk=44; % k-means - d={};Vd=[]; - for nrem=1:nrem_max % 24 - wj=find(kvec>2^9 & kvec<=2^(J-2)); - for i=1:SIMN0 - d{i}=log(var4.E_cs_res1{i,1}(wj)./E_av2(wj)); % E_av{4} - Vd(i)=var(d{i}); % a smooth spectrum (E_av*C) has 0 error - %Vd(i)=mean(d{i}.^2); - %d{i}=log(var4.E_cs_res1{i,1}(wj)./E_min(wj)); % min - %Vd(i)=var(d{i}); - %Vd(i)=mean(d{i}.^2); - %d{i}=log(var4.E_cs_res1{i,1}(wj)./E_av{6}(wj)); - %Vd(i)=norm(d{i}); % compare to org norm - end - %figure,plot(ones(SIMN0,1),Vd,'*') - - [~,Vind0]=sort(Vd); - - %====== least Var - Vind=Vind0(1:end-nrem); % remove nrem worst - %====== kmeans Var - % data=Vd(:); %column - % Classes=dcKMeans(data,2,[min(data);median(data)]); - % xvec=1:length(Vd); - % figure - % hold all, - % plot(xvec(Classes==1),data(Classes==1),'b+') - % plot(xvec(Classes==2),data(Classes==2),'r+') - % sum(Classes==1) - % Vind=xvec(Classes==1); - %====== k-means k wise - % addpath('C:\Users\Lenovo\Documents\MATLAB\MAC\clusteringtoolbox','-end'); - % [score,Ek] = kmeans_spectra(var4.E_cs_res1,wj,@log); - % [~,Vind]=sort(score); - % Vind=Vind(1+nremk:end); - % for c=1:SIMN0 - % Ek{c}=exp(Ek{c}); - % end - %%%var4.E_cs_res1=Ek; % bad freq. changed - %====== - - E_av2=E_av{4}*0; - af=@(x)log(x);afi=@(x)exp(x); %log - %af=@(x)x;afi=af; % linear - for i=Vind - E_av2=E_av2+af(var4.E_cs_res1{i,1}); - %E_av2=E_av2+af(Ek{i,1}); - end - E_av2=afi(E_av2/length(Vind)); - - end - % sorted Var - %figure, plot(Vd(Vind0),'*'),hold all,plot(Vd(Vind),'*') - - - - % plot all cases -% figure -% for i=setdiff(1:SIMN0,Vind) -% loglog(kvec,var4.E_cs_res1{i,1},'+') -% hold all -% end -% for i=Vind -% loglog(kvec,var4.E_cs_res1{i,1},'+k') -% end -% loglog(kvec,var4.E_os{1,1},'+r') -% loglog(kvec,E_av2,'g') -% -% -% % different averages types -% figure; -% set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); -% -% % original -% hh0=loglog(kvec,E_av{6},'-','Color',[0.7 0 0.7],'LineWidth',1.4); -% hold all -% % QOM -% hh4=loglog(kvec,E_min,'--','Color',[0.6 0.6 0.6],'LineWidth',0.8); -% hh4=loglog(kvec,E_av{4},'-','Color',[0 0 0],'LineWidth',0.8); -% -% hh4=loglog(kvec,E_av2,'-','Color',[1 0 0.2],'LineWidth',0.8); - - - % hh4=loglog(kvec,var4.E_cs_res1{15,1},'-','Color',[1 0 0],'LineWidth',0.8); - % hh4=loglog(kvec,var4.E_cs_res1{14,1},'-','Color',[0 1 0],'LineWidth',0.8); - - %============================================================ - - - figure; - if COMPTYPE==1 || COMPTYPE==2 || COMPTYPE==3 %three figures. - if COMPTYPE==2 || COMPTYPE==3 - x_size=0.58*FIG_SIZE(1)-16; - else - x_size=0.58*FIG_SIZE(1); - end - if RUNTYPE==1||RUNTYPE==2||RUNTYPE==3 - y_size=0.58*FIG_SIZE(2)-16; - else - y_size=0.58*FIG_SIZE(2); - end - set(gcf, 'Position', [700 300 x_size y_size]); - else - set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); - end - - - % original - hh0=loglog(kvec,E_av{6},'-','Color',[0.7 0 0.7],'LineWidth',1.4); - hold all - - % M-best - hh1=loglog(kvec,E_av{2},'--','Color',[0.2 0.2 1],'LineWidth',1.3); - - % M/2 - hh1b=loglog(kvec,E_av{1},'--','Color',[0.2 0.2 1],'LineWidth',0.8); - - % shannon - hh2=semilogy(kvec,E_av{5},':','Color',[1 0.2 0.2],'LineWidth',1.3); - - % LOMP - hh3=loglog(kvec,E_av{3},':','Color',[0 0.6 0],'LineWidth',0.8); - - % QOM - %hh4=loglog(kvec,E_av{4},'-','Color',[0 0 0],'LineWidth',0.8); - hh4=loglog(kvec,E_av2,'-','Color',[0 0 0],'LineWidth',0.8); - - - % loglog([100,1400],C1*3e-3*[100,1400].^p1,'k--','LineWidth',0.8) - % loglog([700,1e4],C2*3e-7*[700,1e4].^p2,'k--','LineWidth',0.8) - - - axis tight - % \hspace{-15mm} $ \hspace{10mm} - if ~(COMPTYPE==2 || COMPTYPE==3) - ylabel('\hspace{-2mm} $ E(k)$','FontSize', 11); - end - if ~((COMPTYPE==1||COMPTYPE==2||COMPTYPE==3) && (RUNTYPE==1||RUNTYPE==2||RUNTYPE==3)) - xlabel('$k$','FontSize', 11); - set(gca,'XTick',[10^2,10^3,10^4]) - end - %set(gca,'Outerposition',[-0.08,0,1.1,1]) %posx,posy,widx,widy - %set(gca,'XTick',jvec) - - ylim(ylimit) - xlim([50 N/2]) - - - legend([hh0,hh1,hh1b,hh2,hh3,hh4],{'original','$M$-best','$M/2$-best','Shannon','LOMP',CSNAME}... - ,'interpreter', 'latex','Location','southwest'); - - - if SAVEFIG==1 - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,'av',AVM,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - %export_fig([FIGLOC,'spectra_error_1d',FIGQUAL,'av',AVM,RUNSTR], '-eps', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'a1',FIGQUAL,AVM,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'a1',FIGQUAL,AVM,RUNSTR], '-eps', '-transparent', FIGCOLOR) - disp([FIGLOC,'a1',FIGQUAL,AVM,RUNSTR]) - end - - - end - -end - - - -% Copyright (C) 2014 Gudmundur Adalsteinsson -% See file LICENCE for licence and warranty details diff --git a/spectra_error_1d_plotload5.m b/spectra_error_1d_plotload5.m deleted file mode 100755 index d42759a..0000000 --- a/spectra_error_1d_plotload5.m +++ /dev/null @@ -1,140 +0,0 @@ - -%========== signals ========= -set(0,'defaulttextinterpreter','latex'); -FIGCOLOR='-rgb'; % RGB,CMYK -FIG_SIZE(1)=560; %x (*1.8)<-(1,2) -FIG_SIZE(2)=420; %y -FIGLOC='./spectra_fig/'; -FIGQUAL=''; -CSNAME='SpESO'; % Spectrum Estimation by Sparse Optimization -SAVEFIG=1; - -COMPTYPEv=[1]; - -for COMPTYPE=COMPTYPEv - - EXTRASTR=''; - if COMPTYPE==1 % '64', [32,16,8,4], [1,2,3,4] - SAVEDLOC='./spectra_var/if/'; - SIMN0t='64'; %8,16,64 - Mratio=4; %16,8,4 - RUNTYPEv=[1,2,3,4,10]; - end - - - % 1=w353, 2=w533 wavelet - % 3=f353, 4=f533 fourier - % 10=rhh - %RUNTYPEv=10; %1,2,3,4,5 - for RUNTYPE = RUNTYPEv - - if Mratio==4 - AVM='_4'; - M1=num2str(8263); - M2=num2str(4132); - end - if RUNTYPE==1 - RUNSTR=['_w353',EXTRASTR]; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); - elseif RUNTYPE==2 - RUNSTR=['_w533',EXTRASTR]; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); - elseif RUNTYPE==3 - RUNSTR=['_f353',EXTRASTR]; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); %d2av2d10 - elseif RUNTYPE==4 - RUNSTR=['_f533',EXTRASTR]; %s4 (av2),s2 - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_synth_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_synth_284_1_82',RUNSTR,'_d2if.mat']); %d2av2d10 - elseif RUNTYPE==10 - RUNSTR=['_',EXTRASTR]; - var1=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_rhh_284_1_31',RUNSTR,'.mat']); - var1b=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_rhh_284_1_312',RUNSTR,'.mat']); - var2=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_rhh_284_1_32',RUNSTR,'.mat']); - var3=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M1,'_',SIMN0t,'_rhh_284_1_1',RUNSTR,'.mat']); - var4=load([SAVEDLOC,'spectra_error_1d_saved_32768_',M2,'_',SIMN0t,'_rhh_284_1_82',RUNSTR,'d2if.mat']); - RUNSTR=['_rhh',EXTRASTR]; - end - - if RUNTYPE==1 - ylimit=[8e-16 4e-10]; - C1=1;C2=1;p1=-3;p2=-5/3; - elseif RUNTYPE==2 - ylimit=[8e-16 4e-10]*1000; - C1=1;C2=1e8;p1=-5/3;p2=-3; - elseif RUNTYPE==3 - ylimit=[8e-16 4e-10]*100; - C1=60;C2=60;p1=-3;p2=-5/3; - elseif RUNTYPE==4 - ylimit=[8e-16 4e-10]*100000; - C1=60;C2=60e8;p1=-5/3;p2=-3; - elseif RUNTYPE==10 - ylimit=[8e-12 4e-6]*100000; - C1=60;C2=60e8;p1=-5/3;p2=-3; - end - - N=length(var1.f_os{1}); - J=log2(N); - kvec=var1.k_os{1}; - [SIMN0,~]=size(var1.f_os); - - f=var1.f_os{1}; - - ind=N-2200+1:N-100; - - figure; - set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); - - hh=plot(0:length(ind)-1,f(ind),'-','Color',[0 0 0],'LineWidth',0.8); - set(gca,'YTick',[]) - - axis tight - % \hspace{-15mm} $ \hspace{10mm} -% if ~(COMPTYPE==2 || COMPTYPE==3) -% ylabel('\hspace{-2mm} $ E(k)$','FontSize', 11); -% end -% if ~((COMPTYPE==1||COMPTYPE==2||COMPTYPE==3) && (RUNTYPE==1||RUNTYPE==2||RUNTYPE==3)) -% xlabel('$k$','FontSize', 11); -% set(gca,'XTick',[10^2,10^3,10^4]) -% end - %set(gca,'Outerposition',[-0.08,0,1.1,1]) %posx,posy,widx,widy - %set(gca,'XTick',jvec) - -% ylim(ylimit) - % xlim([50 N/2]) - - - %legend([hh0,hh1,hh1b,hh2,hh3,hh4],{'original','$M$-best','$M/2$-best','Shannon','LOMP',CSNAME}... - % ,'interpreter', 'latex','Location','southwest'); - - - if SAVEFIG==1 - export_fig([FIGLOC,'u1',FIGQUAL,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'u1',FIGQUAL,RUNSTR], '-eps', '-transparent', FIGCOLOR) - disp([FIGLOC,'u1',FIGQUAL,RUNSTR]) - end - - - end - -end - - - -% Copyright (C) 2014 Gudmundur Adalsteinsson -% See file LICENCE for licence and warranty details diff --git a/spectra_error_2d_plotload4.m b/spectra_error_2d_plotload4.m deleted file mode 100755 index 5dffdc5..0000000 --- a/spectra_error_2d_plotload4.m +++ /dev/null @@ -1,333 +0,0 @@ - -%========== averaged spectra 2d ========= -set(0,'defaulttextinterpreter','latex'); -FIGCOLOR='-rgb'; % RGB,CMYK -FIG_SIZE(1)=560; %x (*1.8)<-(1,2) -FIG_SIZE(2)=420; %y -FIGLOC='./spectra_fig/'; -FIGQUAL=''; -CSNAME='SpESO'; % Spectrum Estimation by Sparse Optimization -SAVEFIG=1; - -COMPTYPEv=[1,2,10]; -%COMPTYPEv=1; - -for COMPTYPE=COMPTYPEv - -if COMPTYPE==1 % N/M=16 % '8'*2, [4,8], [1,5] (if) - SAVEDLOC='./spectra_var/if/'; - SIMN0t='8'; - SEEDR=1:2; - Mratio=4; - RUNTYPEv=[1,2,3,5,6,7,15,16]; - %RUNTYPEv=[15]; -elseif COMPTYPE==2 % N/M=64 - SAVEDLOC='./spectra_var/if/'; - SIMN0t='8'; - SEEDR=1:2; - Mratio=8; - RUNTYPEv=[1,5]; -elseif COMPTYPE==10 % N/M=16 % '1'*1, 4, 10 - SAVEDLOC='./spectra_var/'; - SIMN0t='1'; - SEEDR=1; - Mratio=4; - RUNTYPEv=10; -elseif COMPTYPE==12 % N/M=16 *testing* less coef - SAVEDLOC='./spectra_var/if/'; - SIMN0t='8'; %2 actually - SEEDR=1; - Mratio=4; - RUNTYPEv=15; -end - - -% 1=w1, 2=w2, 3=w3 wavelet -% 5=f1, 6=f1, 7=f3 fourier -% 10=r1 radial -% 15=jh 16=jhvor JH -%RUNTYPEv = 7 -for RUNTYPE = RUNTYPEv - EE=cell(6,1); - if Mratio==4 - AVM='_4'; - M1=num2str(65536); - M2=num2str(32761); - elseif Mratio==8 - AVM='_8'; - M1=num2str(16384); - M2=num2str(8281); - end - mstr='synth'; - if RUNTYPE==1 - RUNSTR='_w1'; - ylimit=[2e-5 0.2]; - C1=1e40;C2=4e7;p1=-3;p2=-5/3; - elseif RUNTYPE==2 - RUNSTR='_w2'; - ylimit=[6e-8 0.01]; - C1=3*2e0;C2=5*1e7;p1=-5/3;p2=-3; - elseif RUNTYPE==3 - RUNSTR='_w3'; - ylimit=[2e-5 0.2]; - C1=3e6;C2=4e7;p1=-3;p2=-5/3; - elseif RUNTYPE==5 - RUNSTR='_f1'; - ylimit=[8e-6 8e-2]; - C1=1e40;C2=1e7;p1=-3;p2=-5/3; - elseif RUNTYPE==6 - RUNSTR='_f2'; - ylimit=[2e-8 3e-3]; - C1=2e0;C2=1e7;p1=-5/3;p2=-3; - elseif RUNTYPE==7 - RUNSTR='_f3'; - ylimit=[2e-6 6e-2]; - C1=1e6;C2=1e7;p1=-3;p2=-5/3; - elseif RUNTYPE==10 - RUNSTR='_r1'; %r1, r1_s - ylimit=[0.3 2e3]; - C1=1e40;C2=8e11;p1=-3;p2=-2; - elseif RUNTYPE==15 - mstr='jh'; - RUNSTR='_'; %jh Le1,Le2 - ylimit=[2e-2 1e4]; - C1=1;C2=2e12;p1=-3;p2=-5/3; - elseif RUNTYPE==16 %fix - mstr='jh'; - RUNSTR='_vor'; %jh - ylimit=[4e-1 1e2]; - C1=1;C2=10e7;p1=-3;p2=-5/3+2; - end - RUNSTR4=RUNSTR; - SIMN0t4=SIMN0t; - if RUNTYPE==15 - RUNSTR4='_Le3'; %empty or Le1,2,3 - RUNSTR='_'; - %SIMN0t4='2'; - end - - ics=[1,2,4,5,6]; - for seedn=SEEDR - seedstr=['_',num2str(seedn)]; - var1t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M1,'_',SIMN0t,'_',mstr,'_1_2_31',seedstr,RUNSTR,'.mat']); - var2t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M1,'_',SIMN0t,'_',mstr,'_1_2_312',seedstr,RUNSTR,'.mat']); - var5t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M1,'_',SIMN0t,'_',mstr,'_1_2_32',seedstr,RUNSTR,'.mat']); - var4t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M2,'_',SIMN0t4,'_',mstr,'_1_2_82',seedstr,RUNSTR4,'.mat']); - EEt{1}=var1t.E_cs_res1; - EEt{2}=var2t.E_cs_res1; - EEt{4}=var4t.E_cs_res1; - EEt{5}=var5t.E_cs_res1; - EEt{6}=var1t.E_os; - for i=ics - EE{i}=[EE{i};EEt{i}]; - end - - end - if RUNTYPE==15 - RUNSTR='_jh'; %jh - end - - % =========================== - - Ns=size(var1t.f_os{1},1); - J=log2(Ns); - kvec=var1t.k_os{1}; - [SIMN0,~]=size(EE{4}); - - E_av=cell(6,1); - E_av(1:6)={EE{1}{1,1}*0}; - E_min=E_av{4}+1e20; % minimum - af=@(x)log(x);afi=@(x)exp(x); %log - %af=@(x)x;afi=af; % linear - % average spectra - for i=1:SIMN0 - for j=ics - E_av{j}=E_av{j}+af(EE{j}{i,1}); - end - E_min=min([E_min';EE{4}{i,1}'])'; % minimum - end - - for i=1:6 - E_av{i}=afi(E_av{i}/SIMN0); - end - - - % logscale difference from average - E_av2=E_av{4}; - nrem_max=round(0.3*SIMN0); - nrem_max=0; % org - nremk=44; % k-means - Vind=1:SIMN0; - d={};Vd=[]; - for nrem=1:nrem_max % 24 - wj=find(kvec>2^5 & kvec<=2^(J-2)); - for i=1:SIMN0 - d{i}=log(EE{4}{i,1}(wj)./E_av2(wj)); % E_av{4} - Vd(i)=var(d{i}); % a smooth spectrum (E_av*C) has 0 error - %Vd(i)=mean(d{i}.^2); - %d{i}=log(EE{4}{i,1}(wj)./E_min(wj)); % min - %Vd(i)=var(d{i}); - %Vd(i)=mean(d{i}.^2); - %d{i}=log(EE{4}{i,1}(wj)./E_av{6}(wj)); % E_av{4} - %Vd(i)=norm(d{i}); % compare to org norm - end - %figure,plot(ones(SIMN0,1),Vd,'*') - - [~,Vind0]=sort(Vd); - - %====== least Var - Vind=Vind0(1:end-nrem); % remove nrem worst - %====== kmeans Var - % data=Vd(:); %column - % Classes=dcKMeans(data,2,[min(data);median(data)]); - % xvec=1:length(Vd); - % figure - % hold all, - % plot(xvec(Classes==1),data(Classes==1),'b+') - % plot(xvec(Classes==2),data(Classes==2),'r+') - % sum(Classes==1) - % Vind=xvec(Classes==1); - %====== k-means k wise - % addpath('C:\Users\Lenovo\Documents\MATLAB\MAC\clusteringtoolbox','-end'); - % [score,Ek] = kmeans_spectra(var4.E_cs_res1,wj,@log); - % [~,Vind]=sort(score); - % Vind=Vind(1+nremk:end); - % for c=1:SIMN0 - % Ek{c}=exp(Ek{c}); - % end - %%%var4.E_cs_res1=Ek; % bad freq. changed - %====== - - E_av2=E_av{4}*0; - af=@(x)log(x);afi=@(x)exp(x); %log - %af=@(x)x;afi=af; % linear - for i=Vind - E_av2=E_av2+af(EE{4}{i,1}); - %E_av2=E_av2+af(Ek{i,1}); - end - E_av2=afi(E_av2/length(Vind)); - - end - % sorted Var - %figure, plot(Vd(Vind0),'*'),hold all,plot(Vd(Vind),'*') - - - - % plot all cases - figure - for i=setdiff(1:SIMN0,Vind) - loglog(kvec,EE{4}{i,1},'+') - hold all - end - for i=Vind - loglog(kvec,EE{4}{i,1},'+k') - hold all - end - loglog(kvec,EE{6}{1,1},'+r') - loglog(kvec,E_av2,'g') - -% -% % different averages types -% figure; -% set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); -% -% % original -% hh0=loglog(kvec,E_av{6},'-','Color',[0.7 0 0.7],'LineWidth',1.4); -% hold all -% % QOM -% hh4=loglog(kvec,E_min,'--','Color',[0.6 0.6 0.6],'LineWidth',0.8); -% hh4=loglog(kvec,E_av{4},'-','Color',[0 0 0],'LineWidth',0.8); -% -% hh4=loglog(kvec,E_av2,'-','Color',[1 0 0.2],'LineWidth',0.8); - - - % hh4=loglog(kvec,var4.E_cs_res1{15,1},'-','Color',[1 0 0],'LineWidth',0.8); - % hh4=loglog(kvec,var4.E_cs_res1{14,1},'-','Color',[0 1 0],'LineWidth',0.8); - - %============================================================ - - - figure; - set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); - - % original - hh0=loglog(kvec,moving_average(E_av{6},1,5),'-','Color',[0.7 0 0.7],'LineWidth',1.4); - hold all - - % M/2 - hh1b=loglog(kvec,moving_average(E_av{1},1,5),'--','Color',[0.2 0.2 1],'LineWidth',0.8); - - % M-best - hh1=loglog(kvec,moving_average(E_av{2},1,5),'--','Color',[0.2 0.2 1],'LineWidth',1.3); - - % shannon - hh2=semilogy(kvec,moving_average(E_av{5},1,5),':','Color',[1 0.2 0.2],'LineWidth',1.3); - - % LOMP - % hh3=loglog(kvec,E_av{3},':','Color',[0 0.6 0],'LineWidth',0.8); - - % QOM - %hh4=loglog(kvec,E_av{4},'-','Color',[0 0 0],'LineWidth',0.8); - %hh4=loglog(kvec,E_av2,'-','Color',[0 0 0],'LineWidth',0.8); - hh4=loglog(kvec,moving_average(E_av2,1,5),'-','Color',[0 0 0],'LineWidth',0.8); - - loglog([118,Ns/2-20],C1*3e-3*[118,Ns/2-20].^p1,'k--','LineWidth',0.8) - if RUNTYPE==1 || RUNTYPE==5 - loglog([20,Ns/2-20],C2*3e-7*[20,Ns/2-20].^p2,'k--','LineWidth',0.8) - elseif RUNTYPE==16 - loglog([2,40],C2*3e-7*[2,40].^p2,'k--','LineWidth',0.8) - elseif RUNTYPE==10 - loglog([20,Ns/4],C2*3e-7*[20,Ns/4].^p2,'k--','LineWidth',0.8) - else - loglog([10,138],C2*3e-7*[10,138].^p2,'k--','LineWidth',0.8) - end - - axis tight - % \hspace{-15mm} $ \hspace{10mm} - ylabel('\hspace{-2mm} $ E(k)$','FontSize', 11); - xlabel('$k$','FontSize', 11); - %set(gca,'Outerposition',[-0.08,0,1.1,1]) %posx,posy,widx,widy - %set(gca,'XTick',jvec) - - ylim(ylimit) - if RUNTYPE==16 - xlim([1 Ns/2-6]) - else - xlim([7 Ns/2-6]) - end - - legend([hh0,hh1,hh1b,hh2,hh4],{'original','$M$-best','$M/2$-best','Shannon',CSNAME}... - ,'interpreter', 'latex','Location','southwest'); - - - if SAVEFIG==1 - %export_fig([FIGLOC,'spectra_error_2d',FIGQUAL,'av',AVM,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - %export_fig([FIGLOC,'spectra_error_2d',FIGQUAL,'av',AVM,RUNSTR], '-eps', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'a2',FIGQUAL,AVM,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'a2',FIGQUAL,AVM,RUNSTR], '-eps', '-transparent', FIGCOLOR) - end - - if RUNTYPE==10 % radial image - figure; - set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); - imagesc(1:Ns,1:Ns,var1t.f_os{1}(1:2:end,1:2:end)) - axis square; - colormap(flipud(gray)) - set(gca,'XTickLabel',{' '}) - set(gca,'YTickLabel',{' '}) - if SAVEFIG==1 - %export_fig([FIGLOC,'spectra_error_2d',FIGQUAL,'im',AVM,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - %export_fig([FIGLOC,'spectra_error_2d',FIGQUAL,'im',AVM,RUNSTR], '-eps', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'a2',FIGQUAL,'i',AVM,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'a2',FIGQUAL,'i',AVM,RUNSTR], '-eps', '-transparent', FIGCOLOR) - end - end - -end - -end - - - -% Copyright (C) 2014 Gudmundur Adalsteinsson -% See file LICENCE for licence and warranty details diff --git a/spectra_error_2d_plotload5.m b/spectra_error_2d_plotload5.m deleted file mode 100755 index addaf4d..0000000 --- a/spectra_error_2d_plotload5.m +++ /dev/null @@ -1,149 +0,0 @@ - -%========== signals ========= -set(0,'defaulttextinterpreter','latex'); -FIGCOLOR='-rgb'; % RGB,CMYK -FIG_SIZE(1)=560; %x (*1.8)<-(1,2) -FIG_SIZE(2)=420; %y -FIGLOC='./spectra_fig/'; -FIGQUAL=''; -CSNAME='SpESO'; % Spectrum Estimation by Sparse Optimization -SAVEFIG=0; - -COMPTYPEv=1; - -for COMPTYPE=COMPTYPEv - - if COMPTYPE==1 % N/M=16 % '8'*2, [4,8], [1,5] (if) - SAVEDLOC='./spectra_var/if/'; - SIMN0t='8'; - SEEDR=1; - Mratio=4; - RUNTYPEv=[1,2,3,5,6,7,15,16]; - %RUNTYPEv=15; - end - - - % 1=w1, 2=w2, 3=w3 wavelet - % 5=f1, 6=f1, 7=f3 fourier - % 10=r1 radial - % 15=jh 16=jhvor JH - %RUNTYPEv = 7 - for RUNTYPE = RUNTYPEv - EE=cell(6,1); - if Mratio==4 - AVM='_4'; - M1=num2str(65536); - M2=num2str(32761); - elseif Mratio==8 - AVM='_8'; - M1=num2str(16384); - M2=num2str(8281); - end - mstr='synth'; - if RUNTYPE==1 - RUNSTR='_w1'; - ylimit=[2e-5 0.2]; - C1=1e40;C2=4e7;p1=-3;p2=-5/3; - elseif RUNTYPE==2 - RUNSTR='_w2'; - ylimit=[6e-8 0.01]; - C1=3*2e0;C2=5*1e7;p1=-5/3;p2=-3; - elseif RUNTYPE==3 - RUNSTR='_w3'; - ylimit=[2e-5 0.2]; - C1=3e6;C2=4e7;p1=-3;p2=-5/3; - elseif RUNTYPE==5 - RUNSTR='_f1'; - ylimit=[8e-6 8e-2]; - C1=1e40;C2=1e7;p1=-3;p2=-5/3; - elseif RUNTYPE==6 - RUNSTR='_f2'; - ylimit=[2e-8 3e-3]; - C1=2e0;C2=1e7;p1=-5/3;p2=-3; - elseif RUNTYPE==7 - RUNSTR='_f3'; - ylimit=[2e-6 6e-2]; - C1=1e6;C2=1e7;p1=-3;p2=-5/3; - elseif RUNTYPE==10 - RUNSTR='_r1'; %r1, r1_s - ylimit=[0.3 2e3]; - C1=1e40;C2=8e11;p1=-3;p2=-2; - elseif RUNTYPE==15 - mstr='jh'; - RUNSTR='_'; %jh Le1,Le2 - ylimit=[2e-2 1e4]; - C1=1;C2=2e12;p1=-3;p2=-5/3; - elseif RUNTYPE==16 %fix - mstr='jh'; - RUNSTR='_vor'; %jh - ylimit=[4e-1 1e2]; - C1=1;C2=10e7;p1=-3;p2=-5/3+2; - end - RUNSTR4=RUNSTR; - SIMN0t4=SIMN0t; - if RUNTYPE==15 - RUNSTR4='_Le3'; %empty or Le1,2,3 - RUNSTR='_'; - %SIMN0t4='2'; - end - - ics=[1,2,4,5,6]; - for seedn=SEEDR - seedstr=['_',num2str(seedn)]; - var1t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M1,'_',SIMN0t,'_',mstr,'_1_2_31',seedstr,RUNSTR,'.mat']); -% var2t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M1,'_',SIMN0t,'_',mstr,'_1_2_312',seedstr,RUNSTR,'.mat']); -% var5t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M1,'_',SIMN0t,'_',mstr,'_1_2_32',seedstr,RUNSTR,'.mat']); -% var4t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M2,'_',SIMN0t4,'_',mstr,'_1_2_82',seedstr,RUNSTR4,'.mat']); -% EEt{1}=var1t.E_cs_res1; -% EEt{2}=var2t.E_cs_res1; -% EEt{4}=var4t.E_cs_res1; -% EEt{5}=var5t.E_cs_res1; -% EEt{6}=var1t.E_os; -% for i=ics -% EE{i}=[EE{i};EEt{i}]; -% end - - end - if RUNTYPE==15 - RUNSTR='_jh'; %jh - end - - % =========================== - - Ns=size(var1t.f_os{1},1); - J=log2(Ns); - %kvec=var1t.k_os{1}; - %[SIMN0,~]=size(EE{4}); - - f=var1t.f_os{1}; - %ind=1:2:Ns; - - figure - plot(f(:,400)) - - figure; - set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); - imagesc(1:Ns,1:Ns,f(1:2:end,1:2:end)) - axis square; - %colormap(flipud(gray)) - %set(gca,'XTickLabel',{' '}) - %set(gca,'YTickLabel',{' '}) - set(gca,'YTick',[]) - set(gca,'XTick',[]) - if SAVEFIG==1 - %export_fig([FIGLOC,'spectra_error_2d',FIGQUAL,'im',AVM,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - %export_fig([FIGLOC,'spectra_error_2d',FIGQUAL,'im',AVM,RUNSTR], '-eps', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'u2',FIGQUAL,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'u2',FIGQUAL,RUNSTR], '-eps', '-transparent', FIGCOLOR) - end - - - - end - -end - - - -% Copyright (C) 2014 Gudmundur Adalsteinsson -% See file LICENCE for licence and warranty details diff --git a/spectra_error_2d_plotload6.m b/spectra_error_2d_plotload6.m deleted file mode 100755 index 3d52546..0000000 --- a/spectra_error_2d_plotload6.m +++ /dev/null @@ -1,217 +0,0 @@ - -%========== structure function ========= -set(0,'defaulttextinterpreter','latex'); -FIGCOLOR='-rgb'; % RGB,CMYK -FIG_SIZE(1)=560; %x (*1.8)<-(1,2) -FIG_SIZE(2)=420; %y -FIGLOC='./spectra_fig/'; -FIGQUAL=''; -CSNAME='SpESO'; % Spectrum Estimation by Sparse Optimization -SAVEFIG=1; -tic; - -COMPTYPEv=1; - -for COMPTYPE=COMPTYPEv - - if COMPTYPE==1 % N/M=16 % '8'*2, [4,8], [1,5] (if) - SAVEDLOC='./spectra_var/if/'; - SIMN0t='8'; - SEEDR=1:2; - Mratio=4; - RUNTYPEv=[1,5]; %1 - end - - - % 1=w1, 2=w2, 3=w3 wavelet - % 5=f1, 6=f1, 7=f3 fourier - % 10=r1 radial - % 15=jh 16=jhvor JH - %RUNTYPEv = 7 - for RUNTYPE = RUNTYPEv - EE=cell(6,1); - if Mratio==4 - AVM='_4'; - M1=num2str(65536); - M2=num2str(32761); - elseif Mratio==8 - AVM='_8'; - M1=num2str(16384); - M2=num2str(8281); - end - mstr='synth'; - if RUNTYPE==1 - RUNSTR='_w1'; - ylimit=[2e-5 0.2]; - C1=1e40;C2=4e7;p1=-3;p2=-5/3; - elseif RUNTYPE==2 - RUNSTR='_w2'; - ylimit=[6e-8 0.01]; - C1=3*2e0;C2=5*1e7;p1=-5/3;p2=-3; - elseif RUNTYPE==3 - RUNSTR='_w3'; - ylimit=[2e-5 0.2]; - C1=3e6;C2=4e7;p1=-3;p2=-5/3; - elseif RUNTYPE==5 - RUNSTR='_f1'; - ylimit=[8e-6 8e-2]; - C1=1e40;C2=1e7;p1=-3;p2=-5/3; - elseif RUNTYPE==6 - RUNSTR='_f2'; - ylimit=[2e-8 3e-3]; - C1=2e0;C2=1e7;p1=-5/3;p2=-3; - elseif RUNTYPE==7 - RUNSTR='_f3'; - ylimit=[2e-6 6e-2]; - C1=1e6;C2=1e7;p1=-3;p2=-5/3; - elseif RUNTYPE==10 - RUNSTR='_r1'; %r1, r1_s - ylimit=[0.3 2e3]; - C1=1e40;C2=8e11;p1=-3;p2=-2; - elseif RUNTYPE==15 - mstr='jh'; - RUNSTR='_'; %jh Le1,Le2 - ylimit=[2e-2 1e4]; - C1=1;C2=2e12;p1=-3;p2=-5/3; - elseif RUNTYPE==16 %fix - mstr='jh'; - RUNSTR='_vor'; %jh - ylimit=[4e-1 1e2]; - C1=1;C2=10e7;p1=-3;p2=-5/3+2; - end - RUNSTR4=RUNSTR; - SIMN0t4=SIMN0t; - if RUNTYPE==15 - RUNSTR4='_Le3'; %empty or Le1,2,3 - RUNSTR='_'; - %SIMN0t4='2'; - end - - fv=cell(1,1); - ics=[1,2,4,5,6]; - for seedn=SEEDR - seedstr=['_',num2str(seedn)]; - var1t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M1,'_',SIMN0t,'_',mstr,'_1_2_31',seedstr,RUNSTR,'.mat']); - % var2t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M1,'_',SIMN0t,'_',mstr,'_1_2_312',seedstr,RUNSTR,'.mat']); - % var5t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M1,'_',SIMN0t,'_',mstr,'_1_2_32',seedstr,RUNSTR,'.mat']); - % var4t=load([SAVEDLOC,'spectra_error_2d_saved_1048576_',M2,'_',SIMN0t4,'_',mstr,'_1_2_82',seedstr,RUNSTR4,'.mat']); - % EEt{1}=var1t.E_cs_res1; - % EEt{2}=var2t.E_cs_res1; - % EEt{4}=var4t.E_cs_res1; - % EEt{5}=var5t.E_cs_res1; - % EEt{6}=var1t.E_os; - % for i=ics - % EE{i}=[EE{i};EEt{i}]; - % end - if seedn==1 - fv=var1t.f_os; - else - fv=[fv;var1t.f_os]; - end - end - if RUNTYPE==15 - RUNSTR='_jh'; %jh - end - - % =========================== - - Ns=size(var1t.f_os{1},1); - J=log2(Ns); - %kvec=var1t.k_os{1}; - [SIMN0,~]=size(fv); - - for ii=1:SIMN0 - f=fv{ii}; - %ind=1:2:Ns; - - - % structure function - rmax=81; % p/3 law valid for r in inertial range - - % save('./spectra_var/structure2d', 'S'); - for p=1:12 - S{p} = structure_2d_per(f,p,rmax); %direction (1,0) - S{p} = S{p}/2 + structure_2d_per(f',p,rmax)/2; %direction (0,1) - end - - prange=[1,2,4:12]; - for p=prange - - Ssc=abs(S{p})./abs(S{3}); - - C0=0.6*(1/8).^(p/3-1)*Ssc(8); - e0=(p/3-1); - -% figure -% loglog(Ssc) -% hold all -% loglog(8:rmax, C0*(8:rmax).^e0,'k--') -% xlim([1 rmax]) - - inguess=[e0,C0]; - ine=8:80; - if RUNTYPE==5 - ine=5:50; %F - end - rvec=(1:length(Ssc))'; - fitscale=@(t) log(t); - sco=exp_fit(Ssc(ine),rvec(ine),fitscale,inguess,'log'); - expon(p)=sco(1); -% loglog(ine,sco(2)*ine.^(-sco(1)),'r-') - end - expon=-expon; - expon(3)=0; - expon=expon+1; % /S{3} - prange=1:12; - - exc{RUNTYPE,ii}=expon; - - end - - end - - excm{1}=exc{1,ii}*0; - excm{5}=excm{1}; - %excv{1}=excm{1}; - for ii=1:SIMN0 - excm{1}=excm{1}+exc{1,ii}; - % excv{1}=[excv{1};exc{1,ii}]; - excm{5}=excm{5}+exc{5,ii}; - end - excm{1}=excm{1}./SIMN0; - excm{5}=excm{5}./SIMN0; - %excv{1}=var(excv{1},0,1); - - figure; - set(gcf, 'Position', [700 300 0.66*FIG_SIZE(1) 0.66*FIG_SIZE(2)]); - - hhp=plot(prange,prange/3,':','Color',[0.4 0.4 0.4],'LineWidth',0.8); - hold all - hh=plot(prange,excm{1}(prange),'+-','Color',[0 0 0],'LineWidth',0.8); - hhf=plot(prange,excm{5}(prange),'o-','Color',[0 0 0],'LineWidth',0.8); - - axis tight - xlabel('$p$','FontSize', 11); - ylabel('$\zeta_p$','FontSize', 11); - - %set(gca,'Outerposition',[-0.08,0,1.1,1]) %posx,posy,widx,widy - %set(gca,'XTick',jvec) - %set(gca,'YTick',[]) - ylim([1/3-0.1 4.3]) - xlim([0.85 12.15]) - - legend([hh,hhf,hhp],{'Wavelet $(5/3)$','Fourier $(5/3)$','$p/3$'}... - ,'interpreter', 'latex','Location','southeast'); - - if SAVEFIG==1 - export_fig([FIGLOC,'s2',FIGQUAL,RUNSTR], '-pdf', '-transparent', FIGCOLOR) - export_fig([FIGLOC,'s2',FIGQUAL,RUNSTR], '-eps', '-transparent', FIGCOLOR) - disp([FIGLOC,'s2',FIGQUAL,RUNSTR]) - end - -end -toc - - -% Copyright (C) 2014 Gudmundur Adalsteinsson -% See file LICENCE for licence and warranty details