语音信号的去噪与噪声的定位.docVIP

clc; clear all; nbits=8; %nbit是采样精度 %第一步，读取原始语音信号和噪声 [yuyin,fs,nbits]=wavread(D:\\原始话音.wav); %读取原始语音信号 [yuyin_noise,fs,nbits]=wavread(D:\original_voice_Noise_by_SineWaves.wav); %读取原始语音加噪声信号信号 N=length(yuyin); %第二步，显示原始语音信号和加入噪声的 频谱%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% figure(1) FFT_Signal1=abs(fft(yuyin)); i=1:N/2; x=(i-1)*fs/N; %将横轴改为对应的频率值。 y(i)=FFT_Signal1(i)/(N/2); y(1)=FFT_Signal1(1)/N; subplot(2,1,1); plot(x,y); title(原始语音信号FFT);xlabel(频率(Hz));ylabel(幅值(V)); axis([0 fs/2 0 0.5]); FFT_Signal2=abs(fft(yuyin_noise)); y2(i)=FFT_Signal2(i)/(N/2); y2(1)=FFT_Signal2(1)/N; subplot(2,1,2); plot(x,y2); title(原始带噪声语音信号FFT);xlabel(频率(Hz));ylabel(幅值(V)); axis([0 fs/2 0 0.5]); %第三步，设计带阻滤波器的h(n) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%滤波器%%%%%%%%%%%%%%%%% wc0=2*pi*9950/fs; wc1=2*pi*950/fs; wc2=2*pi*3450/fs; wc3=2*pi*4950/fs; wc4=2*pi*7450/fs; wc5=2*pi*9450/fs; wc6=2*pi*2450/fs; wh1=2*pi*1050/fs; wh2=2*pi*3550/fs; wh3=2*pi*5050/fs; wh0=2*pi*10050/fs; wh4=2*pi*7550/fs; wh5=2*pi*9550/fs; wh6=2*pi*2550/fs; L=3001; window_B=blackman(L); M=L-1; for i=1:L if i==M/2 h1b_p(i)=1-wh1/pi+wc1/pi-wh2/pi+wc2/pi-wh3/pi+wc3/pi; % 当i==M/2时的极限值 else h1b_p(i)=(sin(pi*(i-M))-sin((i-M/2)*wh1)+sin((i-M/2)*wc1)-sin((i-M/2)*wh2)+sin((i-M/2)*wc2)-sin((i-M/2)*wh3)+sin((i-M/2)*wc3)-sin((i-M/2)*wh0)+sin((i-M/2)*wc0)-sin((i-M/2)*wh4)+sin((i-M/2)*wc4)-sin((i-M/2)*wh5)+sin((i-M/2)*wc5)-sin((i-M/2)*wh6)+sin((i-M/2)*wc6)) /((i-M/2)*pi); % 带通滤波器的h(n)表达式 end end h1b_p=h1b_p.*window_B; figure(2) freqz(h1b_p) freqz( h1b_p);title(未加加窗函数后，滤波器的频谱) %第四步，滤波%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% delete_noise=conv( h1b_p,yuyin_noise); %滤出噪声后的的信号 for i=1:L delete_noise(i)=0; end %第五步，显示滤波后的 频谱%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% L_DN=length(delete_noise)