%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% %%% showing Fig. 10 & 15 %%% clear; load('eccsc.mat') load('eccsp.mat') load('eccnsp.mat') load('eccnsc.mat') load('pri_noise.mat') %% fs=4096; rs=size(eccsc); rsn=size(eccnsc); T=0:rs(2)-1; t=T/fs; J_eccsc=zeros(rs); J_eccsp=zeros(rs); Tn=0:rsn(2)-1; tn=Tn/fs; J_eccnsc=zeros(rsn); J_eccnsp=zeros(rsn); %% for k=1:rs(1) J_eccsc(k,:)=10*log10(eccsc(k,:).^2)-20*log10(2*10^(-5)); J_eccsp(k,:)=10*log10(eccsp(k,:).^2)-20*log10(2*10^(-5)); end q = 0.995; Jeccsc=filter((1-q),[1 -q],(sum(J_eccsc)/rs(1))); Jeccsp=filter((1-q),[1 -q],(sum(J_eccsp)/rs(1))); for k=1:rsn(1) J_eccnsc(k,:)=10*log10(eccnsc(k,:).^2)-20*log10(2*10^(-5)); J_eccnsp(k,:)=10*log10(eccnsp(k,:).^2)-20*log10(2*10^(-5)); end J_pri=10*log10(pri_noise.^2)-20*log10(2*10^(-5)); q = 0.995; Jeccnsc=filter((1-q),[1 -q],(sum(J_eccnsc)/rsn(1))); Jeccnsp=filter((1-q),[1 -q],(sum(J_eccnsp)/rsn(1))); Jpri=filter((1-q),[1 -q],(J_pri)); %% figure %.......figure 10 hold on plot(t,Jeccsc,'color',[.5 .5 .5],'LineWidth',2) plot(t,Jeccsp,':','color',[0 0 0],'LineWidth',2) legend('Conventional Algorithm','Proposed Algorithm') xlim([0.25,20]); ylim([78,86]); grid on box on set(gca,'fontsize',20,'fontname','Times'); xlabel('Time/s','Fontname', 'Times New Roman','FontSize',24); ylabel('Sound Pressure Level/dB','Fontname',... 'Times New Roman','FontSize',24); figure %......figure 15 hold on plot(tn,Jpri,'color',[.8 .8 .8],'LineWidth',2)% plot(tn,Jeccnsc,'color',[.5 .5 .5],'LineWidth',2)% plot(tn,Jeccnsp,'color',[0 0 0],'LineWidth',2)% legend('Primary Noise','Conventional Algorithm',... 'Proposed Algorithm','Location','best') xlim([0.25,30]); ylim([74,90]); grid on box on set(gca,'fontsize',20,'fontname','Times'); xlabel('Time/s','Fontname', 'Times New Roman','FontSize',24); ylabel('Sound Pressure Level/dB','Fontname',... 'Times New Roman','FontSize',24);