IIR数字滤波的冲击响应不变法设计.docVIP

课程设计说明书 NO. PAGE 1 沈 阳 大 学 语音信号的数字滤波 —IIR数字滤波器的冲激响应不变法设计 设计题目：语音信号的数字滤波—IIR数字滤波器的冲激响应不变法设计 课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；掌握利用数字滤波器对信号进行滤波的方法。并能够对设计结果加以分析。 设计步骤 总体方案 首先需要录制一段音频信号，利用Matlab软件测出原始信号的频率fs、语音信号的采样频率，通过Matlab对原始信号进行频谱分析，记录原始信号的最高频率fH、最低频率fL和频率分辨率Δf，最后对其语音信号加载上噪声生成最终的原始语音信号保存下来。 接下来就是数字滤波器的设计过程，对于本次利用冲激响应不变法设计的IIR数字滤波器，大体可以总结为以下的过程：先利用冲激响应不变法，把模拟低通原型数字化成低通滤波器，然后利用数字频带变换法，将它变换成所需要的低通滤波器，如图1所示。 数字-数字频带变换 数字-数字频带变换 模拟归一化原型Ωc=1 数字低通 数字低通 数字化 图1 设计数字滤波器的过程图 最后根据数字滤波器的性能要求转换为与之相应的作为“样本”的模拟滤波器的性能要求，根据此性能要求设计模拟滤波器，这既可以用查表的方法，也可以用解析的方法。然后通过冲激响应不变法将此“样本”模拟低通滤波器Ha(s)数字化为所需的数字低通滤波器H（z）。 2.2 设计的具体步骤 2.2.1 语音信号的采样录音 （1）按“开始-程序-附件-娱乐-录音机”的顺序操作打开Windows系统中的录音软件； （2）用麦克风录入自己的声音信号并保存成文件； （3）记录以下内容，如下表1所示。 文件名、格式 采样速率 时间长度 luyin.wav 8000Hz 3.16秒 表1 录音语音信号的相关参数 2.2.2 （1）将上一步骤中保存下来的语音信号文件“luyin.wav”复制到计算机装有Matlab软件的磁盘中相应Matlab目录中的“work”文件夹中； （2）双击桌面上Matlab软件的快捷图标，打开Matlab软件； （3）在菜单中选贼“File-new-M-File”，打开m文件编辑器； （4）在m文件编辑器中输入相应的指令将自己的语音信号导入Matlab工作台。 （5）编写程序段，利用Matlab中已有的FFT程序函数和画图程序函数对导入的语音信号进行分析，并根据自己的语音信号频谱中主要能量的分布特征确定自己的语音信号所在的频段； 程序代码如下： [s, fs]=wavread(luyin.wav); %将声音文件“luyin.wav”中的波形信号值赋值给s %将采样速率赋值给变量fs N=25000; %取被处理信号的长度为小于25000 subplot(311),plot(s),title(声音信号的时域波形); %绘制声音信号的时域波形 s1=fft(s); %取被处理信号s的FFT，其频谱赋值给s1 subplot(312),plot(abs(s1)),title(声音信号的频谱); %绘制声音信号的频谱 detf=fs/N; %计算频率分辨率，并赋值给变量detf subplot(313),plot(k(1:N/2)*detf,abs(s1(1:N/2))),title(声音信号的单边谱); %绘制声音信号的单边谱 （6）记录数据，声音信号的波形图，如图2。 图2 频谱分析图 上图是程序执行结果，其中第一个图是声音信号在时域范围内的波形图，从图中可以看出声音的信号主要集中在0-25000的区域内，所以声音的信号长度取25000；第二个图是声音信号的频谱图，其横坐标的具体值是X（N）中的序号N；第三个图是确定滤波频率范围的参考图，其横坐标是通过频率分辨率求得的，频率分辨率，通过第三个图还能看出声音信号的频率主要集中在200Hz-1700Hz之间，选定1100Hz，向左数第一个尖峰为1000Hz，是通带截止频率，向右数第一尖峰为1250Hz，是阻带截止频率。 信号频谱分析后的主要参数如表2所示。 声音信号的长度 fL fH 25000 200Hz 1700Hz 表2 声音信号的频谱分析的参数 2.2.3 2.2.3.1确定滤波器的参数 （1）模拟频率与数字频率之间为线性关系，，T为抽样周期，故 通带的截止频率； 阻带的截止频率； （2）滤波器的通带允许的最大衰减阻带必须达到的最小