dsp课程设计讲述正文.doc

语音信号滤波去噪——使用脉冲响应不变法设计的巴特沃斯滤波器 摘 要 本课程设计主要运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用一个满足指标的，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。MATLAB 1 引 言 本课程设计主要利用麦克风采集一段语音信号，1.2 课程设计的要求 （1）滤波器指标必须符合工程设计。 （2）设计完后应检查其频率响应曲线是否满足指标。 （3）处理结果和分析结论应该一致，而且应符合理论。 （4）独立完成课程设计并按要求编写课程设计报告。 2 设计原理 用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用一个满足指标的，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析 3.2语言信号的采集 利用PC 机上的录音软件进行采集。将话筒输入计算机的语音输入插口上,启动录音机。按下录音按钮,接着对话筒说话“大萝卜”,说完后停止录音,屏幕左侧将显示所录声音的长度。点击放音按钮,可以实现所录音的重现。以文件名“speech”保存入g :\ MATLAB \ work 中。可以看到,文件存储器的后缀默认为. wav ,这是WINDOWS 操作系统规定的声音文件存的标准。然后，通过文件中的属性，将此录音转换成8000Hz，8k的单声道语音信号。 图3.2 语音信号设置 3.3语音信号的频谱分析 利用MATLAB中的wavread”命令来读入（采集）语音信号，将它赋值给某一向量。 [y,fs,nbits]=wavered (‘speech.wav’); 然后，画出语音信号的时域波形，再对语音信号进行频谱分析。MATLAB提供了快速傅里叶变换算法FFT计算DFT的函数fft,其调用格式如下： Xk=fft(xn,N) 域序列向量，N是DFT变换区间长度，当N大于xn的长度时，fft函数自动在xn后面补零。，当N小于xn的长度时，fft函数计算xn的前N个元素，忽略其后面的元素。 在本次设计中，我们利用fft对语音信号进行快速傅里叶变换，就可以得到信号的频谱特性。其程序如下： [y,fs,nbits]=wavread(44.wav); n=length(y); Y=fft(y,n); subplot(2,1,1);plot(y);title(原始信号波形); subplot(2,1,2);plot(abs(Y));title(原始信号频谱); 图3.3 时域波形 图3.4 频域幅度谱 图3.5 频域幅度谱和相位谱 3.4滤波器设计 将数字滤波器的设计指标设为通带截止频率fb=1100HZ,阻带频率fc=1200HZ,通带波纹Ap=1dB，阻带波纹As=20dB,要求确定H(z)。设计步骤如下： （1）确定所需类型数字滤波器的技术指标。 （2）将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为Ω=2/T tan(0.5ω) (3)将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标。 （4）设计模拟低通滤波器。 （5）通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。 （6）采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。 程序： Fp=1200;%阻带截止频率 Fs=1100;%通带截止频率 Ft=8000;%采集频率 As=20;%通带波纹Ap=1 Ap=1;%阻带波纹As=20 wp=2*pi*Fp/Ft; ws=2*pi*Fs/Ft; fp=2*Ft*tan(wp/2); fs=2*Ft*tan(ws/2); [n,wn]=buttord(wp,ws,Ap,As,s); %求低通滤波器的阶数和截止频率 [b,a]=butter(n,wn,s); %求S域的频率响应的参数 [num,den]=bilinear(b,a,1); %双线性变换实现S域到Z域的变换 [h,w]=freqz(num,den); %根据参数求出频率响应 plot(w*fs/(2*pi),abs(h)); grid; xlabel(频率／Hz); ylabel(频率响应幅度); title(IIR低通滤波器); legend(用butter设计); 图3.6 利用脉冲响应不变法设计的数字巴特沃斯滤波器（w(单位：)） 3.5 信号滤波信号 图3.7 滤波前后的时域对比波形 滤波前后频域频谱图3.8所示 图3.8