wavelet.rar_声音去噪_声音小波去噪_小波去噪_小波去噪程序_小波声音

clear clc sound=wavread('a.wav'); count1=length(sound); noise=0.05*randn(1,count1); for i=1:count1 signal(i)=sound(i); end for i=1:count1 y(i)=signal(i)+noise(i); end [coefs1,coefs2]=dwt(y,'db3'); %[low frequency, high frequency] count2=length(coefs1); count3=length(coefs2); energy1=sum((abs(coefs1)).^2); energy2=sum((abs(coefs2)).^2); energy3=energy1+energy2; for i=1:count2 recoefs1(i)=coefs1(i)/energy3; end for i=1:count3 recoefs2(i)=coefs2(i)/energy3; end %low frequency parameter estimate the sound signal is unvoice or voice zhen=160; count4=fix(count2/zhen); for i=1:count4 n=160*(i-1)+1:160+160*(i-1); s=sound(n); w=hamming(160); sw=s.*w; a=aryule(sw,10); sw=filter(a,1,sw); sw=sw/sum(sw); %normalization r=xcorr(sw,'biased'); corr=max(r); %if unvoiced, output 1, if voice, output 0; if corr>=0.8 output1(i)=0; elseif corr<=0.1 output1(i)=1; end end for i=1:count4 n=160*(i-1)+1:160+160*(i-1); if output1(i)==1 switch abs(recoefs1(i)) case abs(recoefs1(i))<=0.002 recoefs1(i)=0; case abs(recoefs1(i))>0.002 & abs(recoefs1(i))<=0.003 recoefs1(i)=sgn(recoefs1(i))*(0.003*abs(recoefs1(i))-0.000003)/0.002; otherwise recoefs1(i)=recoefs1(i); end elseif output1(i)==0 recoefs1(i)=recoefs1(i); end end count5=fix(count3/zhen); for i=1:count5 n=160*(i-1)+1:160+160*(i-1); s=sound(n); w=hamming(160); sw=s.*w; a=aryule(sw,10); sw=filter(a,1,sw); sw=sw/sum(sw); r=xcorr(sw,'biased'); corr=max(r); %if unvoiced, output 1, if voice, output 0; if corr>=0.8 output2(i)=0; elseif corr<=0.1 output2(i)=1; end end for i=1:count5 n=160*(i-1)+1:160+160*(i-1); if output2(i)==1 switch abs(recoefs2(i)) case abs(recoefs2(i))<=0.002 recoefs2(i)=0; case abs(recoefs2(i))>0.002 & abs(recoefs2(i))<=0.003 recoefs2(i)=sgn(recoefs2(i))*(0.003*abs(recoefs2(i))-0.000003)/0.002; otherwise recoefs2(i)=recoefs2(i); end elseif output2(i)==0 recoefs2(i)=recoefs2(i); end end %wavelet transfer output3=idwt(recoefs1, recoefs2,'db3'); %normalization maxdata=max(output3); output4=output3/maxdata; wavwrite(y,22050,32,'aa'); wavwrite(output4,22050,32,'ab'); display(' done============');