数字信号处理的课程设计——IIR数字滤波器冲激响应不变法设计.docVIP

课 程 设 计 任 务 书 学 院 班 级 姓 名 设 计 起 止 日 期 设计题目：语音信号的数字滤波—IIR数字滤波器冲激响应不变法设计 设计任务（主要技术参数）： 语音信号的采样录音（、N）； 语音信号的频谱分析（） 语音信号的加噪（）； 数字滤波器的设计（，滤波器的阶数）； 滤波结果输出及效果分析（） 指导教师评语： 成绩： 签字： 年 月 日 语音信号的数字滤波 ——IIR数字滤波器的冲激响应不变换法设计 设计题目：语音信号的数字滤波——IIR数字滤波器的冲激相应不变法设计 一.课程设计的目的: 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；掌握利用数字滤波器对信号进行滤波的方法。并能够对设计结果加以分析。 二.设计步骤: 2.1．语音信号的采样录音: 利用windows对语音信号进行采集，按“开始－＞程序－＞附件－＞娱乐－＞录音机”的顺序操作打开Windows系统中的录音机软件,开始录音。 图1 利用Window打开录音机 （1）用麦克风录入自己的声音信号并保存成qiufang.wav文件； （2）记录语音信号文件保存的文件名qiufang.wav, 格式为PCM格式，属性为语音信号的采样频率fs等于8000Hz，8位，单声道，7KB/秒的语音信号，语音信号的时间长度为4.66s。 2语音信号的频谱分析 （1）将上一步骤中保存下来的语音信号文件“qiufang.wav”复制到计算机装有Matlab软件的磁盘F中相应Matlab目录中的“work”文件夹中； （2）双击桌面上Matlab软件的快捷图标，打开Matlab软件； （3）在菜单栏中选择“File－＞new－＞M-File”或是点击快捷按钮，打开m文件编辑器； （4）在m文件编辑器中输入相应的指令将自己的语音信号导入Matlab工作台。 图3 原信号频谱图 代码如下： close all clear all [x,fs]=wavread(qiufang.wav); %将语音文件“qiufang.wav”中的波形信号值赋值给数组x %将采样速率赋值给变量fs subplot(211) %把一个图形窗口中分2行一列，在第一个分区中画图 plot(x) %绘制原始语音信号的时间域波形 x1=x(4000:4000+fs-1); %取出被处理信号x1 wavwrite(x1,fs,x1.wav) %将被处理信号x1输出为语音文件“x1.wav” subplot(212) %把一个图形窗口中分2行一列，在第二个分区中画图 plot(x1) %绘制被处理的语音信号的时间域波形 运行得到如图4所示2个信号图形，第一个为原始语音信号“qiufang.wav”的时间域波形，第二个为被处理的语音信号“x1.wav”的时间域波形。舍弃语音信号“zhanghan.wav”起首空白部分，取被处理信号的长度为L=8000，将被处理信号x1 输出为语音文件“x1.wav”。 图4 语音信号的时间域波形图 X1=fft(x1); %对被处理语音信号x1求取FFT，其频谱赋值给数组X1 figure plot(abs(X1)) %绘制被处理语音信号的FFT幅频特性 运行得如图4所示波形，该波形为被处理的语音信号FFT幅频特性，从图中可以看出3bd截止频率fs为8000Hz，通带截止频率fh为1700Hz，阻带截止频率fc为3500Hz 根据公式wp=2πfh/fs ,wst=2πfc/fs,wc=1/2(wp+wst)分别计算出wp=0.425π，wst=0.875π,wc=0.65。 图5 被处理信号的幅频特性 2.3 迭加噪声 x2=awgn(x1,20,0); %在被处理信号x1上迭加20dB的噪声，并赋值给数组x2 wavwrite(x2,fs,x2.wav) %将含噪声的信号x2输出为语音文件“x2.wav” figure subplot(211) plot(x2) %绘制叠加噪声语音信号的时间域波形 X2=fft(x2); %对叠加噪