function?gonngcheng2

s=wavread（‘tingli.wav‘）;?

[sfsbits]?=?wavread（‘tingli.wav‘）;

s（39923:65536）=0;
%??sound（sfs）;
?


?
?S=fft（s）;
?N=length（s）;
?n=0:N-1;
?f=n‘*fs/N;
?s1=s（1:4:end）;?S1=fft（s139922）;?
?s2=decimate（s4）;S2=fft（s239922）;?
?
?
?figure（1）;
?subplot（311）;plot（s）;title（‘原始信號波形‘）;
?
?subplot（313）;plot（fabs（S））;title（‘信號頻譜幅值‘）;
?
%?init=2055615866;
%?randn（‘seed‘init）;
%?k1=s+20*randn（size（s））;
%?figure（2）
%?subplot（311）
%?plot（k1）
%?K1=fft（k1）;
%?subplot（312）
%?plot（fK1）


%?sqrt_snr=3;
%?init=126424;
%?[xxn]=wnoise（216sqrt_snrinit）;
%?figure
%?plot（x）
%?title（‘原始信號‘）
%?axis?square
%?figureplot（xn）
%?title（‘含噪信號‘）
%?axis?square
%?noise11=xn-x;
%?figure
%?plot（noise11）
%?title（‘噪聲信號‘）
%?axis?square
%?noise11=noise11.‘;
%?sss=s+noise11;
%?
%?sound（sssfs）

%?[thrsorhkeepapp]=ddencmp（‘den‘‘wv‘sss）;
%?sss1=wdencmp（‘gbl‘sss‘db4‘thrsorhkeepapp）;

count=length（s）;
noiserr=0.05*randn（1count）;
sss=s.‘+noiserr;



figureplot（noiserr）
title（‘噪聲信號‘）
axis?square
figure
plot（sss）
title（‘信號和噪信‘）
axis?square


%?sound（sssfs）






[cl]=wavedec（sss5‘db6‘）;

sigma=wnoisest（cl1）;

alpha=2;

thr=wbmpen（clsigmaalpha）;

keepapp=1;

sss1=wdencmp（‘gbl‘cl‘db6‘5thr‘s‘keepapp）;
sss1=sss1.‘;

%?sound（sss1fs）

%消噪后的情況


xhpp=fft（sss1）;
?N=length（sss1）;
?n=0:N-1;
?f=n‘*fs/N;

?
?
?figure;
?subplot（311）;plot（sss1）;title（‘消除噪聲后信號波形‘）;
?
?subplot（313）;plot（fabs（xhpp））;title（‘消除噪聲后信號頻譜幅值‘）;



















?
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-04-11?21:32??對聲音信號加入隨機噪聲，去除噪聲，觀察去噪前后頻譜圖\
?????文件????????1663??2011-12-26?23:59??對聲音信號加入隨機噪聲，去除噪聲，觀察去噪前后頻譜圖\gonngcheng2.m
?????文件???????79888??2005-12-22?13:00??對聲音信號加入隨機噪聲，去除噪聲，觀察去噪前后頻譜圖\tingli.wav