(E)用双线性变换法设计的低通滤波器及对加入噪声的信号的滤波处理。 figure(5);
%用双线性变换法设计的低通滤波器 fp=1000;fc=1200;As=100;Ap=1;fs1=22050; wc = 2*fc/fs1;wp = 2*fp/fs1; [n,wn]=ellipord(wp,wc,Ap,As); [b,a]=ellip(n,Ap,As,wn); freqz(b,a,512,fs1);
x = filter(b,a,y2); %对加噪后的语音信号进行低通滤波 X = fft(x,40000); %对滤波后的信号进行FFT figure(6);
subplot(211);plot(x);title('低通滤波后信号波形'); subplot(212);plot(abs(X));title('低通滤波后信号频谱'); sound(x,40000); %播放滤波后语音
图5 滤波器的频率响应
图6 滤波后的信号和频谱
五、实验总结
通过本次实验,进一步加深DFT算法原理和基本性质的理解,熟悉了FFT算法原理,理解和掌握了音频信号各参数的意义,通过对音频信号加入噪声而去掉噪声的过程学会了设计低通滤波器的方法,巩固了课程理论知识。