IIR数字滤波器设计课程设计.docVIP

理工大学计算机科学与技术学院 课程设计 课程名称 课题名称专 业 通信工程 班 级 学 号 姓 名 指导教师 年 月 日一 设计目的 4 二设计原理及要求 4 三 功能模块 5 3.1语音信号分析 5 3.2噪语音信号合成 5 3.3滤波器的设计 5 四 程序及结果 6 4.1程序 6 4.2运行结果 9 五 流程图 12 六 心得体会 13 七 参考文献 14 一 、设计目的 综合运用数字信号处理的理论知识进行频谱分析和滤波器设计，通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。1．熟悉用双线性变换法设计IIR数字滤波器的原理和方法； 2．了解用脉冲响应不变法设计IIR数字滤波器的原理和方法； 3．掌握双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点； 4．掌握数字滤波器的计算机仿真方法要求IIR数字滤波器的系统函数为的有理分式： 设计IIR滤波器的系统函数，就是要确定的阶数N及分子分母多项式的系数和，使其满足指定的频率特性。，因此IIR滤波器设计的方法之一是：先设计一个合适的模拟滤波器，然后将模拟滤波器通过适当的 由于模拟滤波器的设计有许多简单而严谨的设计公式和大量的图表可以利用变换转换成满足给定指标的数字滤波器。 1、 MATLAB 的使用，掌握 MATLAB 的程序设计方法。 2、 Windows 环境下语音信号采集的方法。 3、数字信号处理的基本概念、基本理论和基本方法。 4、TLAB 设计 FIR 和 nR 数字滤波器的方法。 5、用 MATLAB 对信号进行分析和处理出采样后的语音信号的时域波形和频谱图MATLAB软件平台下，给原始的语音信号叠加上噪声，噪声类型分为如下几种：（1）白噪声；（2）单频噪色（正弦干扰）；（3）多频噪声（多正弦干扰）；（4）其它干扰，可设置为低频、高频、带限噪声，或Chirp干扰、冲激干扰。绘出加噪后的语音信号时域和频谱图，原始信号对比， [N,Wn]=buttord(wp,ws,Rp,Rs); （4）设计巴特沃兹滤波器，调用系统函数： [B,A]=butter(N,Wn,…); （5）数字滤波器的频率响应：[H,W]=freqz(B,A,…); （6）显示频率响应。 （7）滤波后的的语音信号时域波形、信号频谱。 四、程序及结果： 4.1程序： [x1,fs,bits]=wavread(e:\hello酷狗.wav);%打开酷狗时的音频文件 [x1,fs,bits]=wavread(e:\SAH.wav); t=0:(size(x1)-1); figure(1); subplot(4,1,1);plot(t,x1); title(原始语音信号时域波形); subplot(4,1,2);y1=fft(x1,40000);plot(abs(y1)); title(原始语音信号频谱); N=length(x1)-1;t=0:1/fs:N/fs; d=[0.3*cos(2*pi*5000*t)];x2=x1+d; sound(x2,fs,bits); subplot(4,1,3);plot(x2); title(加单频余弦信号的语音信号时域波形); subplot(4,1,4);y2=fft(x2,40000);plot(abs(y2)); title(加单频余弦信号的语音信号频谱); Fs=500;t=0:1/Fs:1; fp=30;fs=60; wp=2*fp/Fs;ws=2*fs/Fs;Rp=1;Rs=40; figure(2); [N,wc]=buttord(wp,ws,Rp,Rs); [B,A]=butter(N,wc); [H1,W1]=freqz(B,A);y1=filter(B,A,x1); subplot(311);plot(W1*Fs/(2*pi),abs(H1)); title(低通滤波器频率响应); xlabel(频率);ylabel(幅度); subplot(312);plot(y1