数字信号滤波器设计 - 图文

山东财经大学学士学位论文

图5-6 理想低通滤波器波形

图5-7 汉明窗滤波器波形

19

山东财经大学学士学位论文

5.4.4 利用布莱克曼窗设计 FIR低通滤波器 利用布莱克曼窗设计的滤波器的主要程序如下： wp=0.2*pi;

ws=0.3*pi; wc=1/2*(wp+ws);

wdelta=ws-wp;%过渡带宽(数字频率) N=ceil(11*pi/wdelta); %取整 a=(N-1)/2 n=0:N-1 c=n+0.0000001

hd=sin(wc*(c-a))./(pi*(c-a)) ;%;理想低通滤波器 x_label=linspace(-pi,pi,N); figure(1) subplot(211) plot(x_label,hd) xlabel('时域') ylabel('幅度')

title('理想低通滤波器时域波形') subplot(212)

plot(x_label,fftshift(abs((fft(hd))))) xlabel('频率(rad/s)') ylabel('频率响应幅度')

title('理想低通滤波器频域波形')

w=0.42-0.5*cos((2*pi*c)./(N-1))+0.08*cos((4*pi*c)./(N-1)); figure(2) subplot(211) plot(x_label,w) xlabel('时域') ylabel('幅度')

title('布莱克曼窗时域波形') subplot(212)

plot(x_label,fftshift(abs((fft(w))))) xlabel('频率(rad/s)') ylabel('频率响应幅度') title('布莱克曼窗频域波形') for i=1:N

%选择布莱克曼窗

山东财经大学学士学位论文

h(i)=hd(i).*w(i) ; %生成窗函数 end figure(3) subplot(211) plot(x_label,h) xlabel('时域') ylabel('幅度')

title('窗函数时域波形') subplot(212)

plot(x_label,fftshift(abs((fft(h))))) xlabel('频率(rad/s)') ylabel('频率响应幅度') title('窗函数频域波形')

波形如下图所示：

图5-8理想低通滤波器波形

21

山东财经大学学士学位论文

图5-8布莱克曼窗滤波器波形

5.4.5 用滤波器对加噪语音信号进行滤波

用设计的汉明窗滤波器和布莱克曼滤波器分别对加噪的语音信号进行滤波，在Matlab中，利用函数conv对加噪后的信号和构造出的滤波器进行卷积，从而实现滤波效果。

比较汉明窗滤波器对原始信号滤波前后的波形，如下图：

图5-9 输入原始信号波形