小波去除音频信号中的噪声，观察频谱

对声音信号加入随机噪声，用小波去除噪声，观察去噪前后的时域和频谱图，可听前后声音的变化。

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