数字信号处理课程设计报告
noise=0.01*randn(n,2); %随机函数产生噪声 s=y+noise; %语音信号加入噪声 S=fft(s); %傅里叶变换 z12=filter(num12,den12,s); sound(z12);
m12=fft(z12); %求滤波后的信号 subplot(2,2,1); plot(abs(S),'g');
title('滤波前信号的频谱'); subplot(2,2,2); plot(abs(m12),'r'); title('滤波后信号的频谱'); subplot(2,2,3); plot(s);
title('滤波前信号的波形'); subplot(2,2,4); plot(z12);
title('滤波后的信号波形');
图10 双线性法带通滤波
%*****************************高通滤波器*****************************
[y,fs,nbits]=wavread ('speech');
14
数字信号处理课程设计报告
n = length (y) ; %求出语音信号的长度 noise=0.01*randn(n,2); %随机函数产生噪声 s=y+noise; %语音信号加入噪声 S=fft(s); %傅里叶变换 z13=filter(num13,den13,s); sound(z13);
m13=fft(z13); %求滤波后的信号 subplot(2,2,1); plot(abs(S),'g');
title('滤波前信号的频谱'); subplot(2,2,2); plot(abs(m13),'r'); title('滤波后信号的频谱'); subplot(2,2,3); plot(s); title('滤波前信号的波形'); subplot(2,2,4); plot(z13);
title('滤波后的信号波形');
图11 双线性法高通滤波
%========================窗函数法==========================
15
数字信号处理课程设计报告
%**************************低通滤波器****************************
[y,fs,nbits]=wavread ('speech');
n = length (y) ; %求出语音信号的长度 noise=0.01*randn(n,2); %随机函数产生噪声 s=y+noise; %语音信号加入噪声 S=fft(s); %傅里叶变换 z21=fftfilt(b21,s); sound(z21);
m21=fft(z21); %求滤波后的信号 subplot(2,2,1); plot(abs(S),'g');
title('滤波前信号的频谱'); subplot(2,2,2);
plot(abs(m21),'r'); title('滤波后信号的频谱'); subplot(2,2,3);
plot(s); title('滤波前信号的波形'); subplot(2,2,4);
plot(z21); title('滤波后的信号波形');
图12 窗函数法低通滤波
16
数字信号处理课程设计报告
%****************************带通滤波器***************************
[y,fs,nbits]=wavread ('speech');
n = length (y) ; %求出语音信号的长度 noise=0.01*randn(n,2); %随机函数产生噪声 s=y+noise; %语音信号加入噪声 S=fft(s); %傅里叶变换 z22=fftfilt(b22,s); sound(z22);
m22=fft(z22); %求滤波后的信号 subplot(2,2,1); plot(abs(S),'g');
title('滤波前信号的频谱'); subplot(2,2,2); plot(abs(m22),'r'); title('滤波后信号的频谱'); subplot(2,2,3); plot(s);
title('滤波前信号的波形'); subplot(2,2,4); plot(z22);
title('滤波后的信号波形');
图13 窗函数法带通滤波
17
数字信号处理课程设计报告
%*************************高通滤波器*****************************
[y,fs,nbits]=wavread ('speech');
n = length (y) ; %求出语音信号的长度 noise=0.01*randn(n,2); %随机函数产生噪声 s=y+noise; %语音信号加入噪声 S=fft(s); %傅里叶变换 z23=fftfilt(b23,s); sound(z23);
m23=fft(z23); %求滤波后的信号 subplot(2,2,1); plot(abs(S),'g'); subplot(2,2,2); plot(abs(m23),'r'); title('滤波后信号的频谱'); subplot(2,2,3); plot(s);
title('滤波前信号的波形'); subplot(2,2,4); plot(z23);
title('滤波后的信号波形');
图14 窗函数法高通滤波
18