《数字信号处理》课程研究性学习报告_DSP3.docVIP

《数字信号处理》课程研究性学习报告 姓名 张卓 学号 同组成员 王欣然 王一明 马高飞 郭鹏飞 指导教师 胡健 时间 2011年5月29日 数字信号处理综合应用专题研讨 【目的】 (1) 能够灵活应用DFT（FFT）分析实际信号的频谱。 (2) 熟悉通过IIR和FIR数字滤波器进行实际系统设计的方法。 (3) 培养学生自主学习能力，以及发现问题、分析问题和解决问题的能力。 【研讨题目】 1．附件给出了一段含有噪声的音频信号。 （1）分析该信号的频谱特点。 （2）通过脉冲响应不变法设计一个滤波器对其进行处理，得到有用信息，自主确定各项指标。 （3）使用双线性变换法重做（2）。 （4）试定量比较上述两种滤波器的各项性能，画出能说明性能差异的相关图形，对比并解释。 （5）通过不同窗函数法设计FIR滤波器对其进行处理，比较不同窗函数的滤波效果，给出分析和结论。 (6) 请尝试采用其它的音频信号，混入不同的噪声，利用所学的滤波方法进行分析，会得到什么样的效果？ 【题目分析】 本题讨论用IIR和FIR数字滤波器进行实际系统设计的方法。 【仿真结果】 (1) (2) 脉冲响应不变法 双线性变换法 (3) Hamming窗 Blackman窗 Kaiser窗 (4) Butterworth (5) 【结果分析】 对几种滤波器得到的音频信号进行对比。 原始音频信号频谱在频率为10^4Hz以及4*10^4Hz处分别有两处毛刺，关于2.6*10^4Hz大致对称。 对于IIR滤波器，用Butterworth滤波器，指标：fp=6000Hz，fs=8000Hz，Ap=0.5dB，As=50dB。 对于FIR滤波器，用Hamming窗，指标：fp=6000;fs=8000。 用种滤波器滤波后均去处了毛刺，都达到了滤波的效果。 画出两种滤波器的幅度相应图形，可以看出，Butterworth滤波器的幅度相应较为平滑，没有等波纹振荡，通带比FIR滤波器的小，过渡带也比FIR滤波器要小。而FIR滤波器的幅度响应在阻带有所波动。 加Hamming窗的滤波器N=79大于Butterworth滤波器的N,IIR数字滤波器可以用较少的阶数获得较好的幅度相应，但由于其结构存在反馈，在实现时可能会造成系统不稳定，并且是非线性相位的。由于FIR数字滤波器的单位脉冲相应是有限长，因而系统是稳定的，并且很容易设计成线性相位。 【自主学习内容】 一、IIR 模拟模拟低通滤波器 （1）Butterworth型低通滤波器 [N,wc]=buttord(wp,ws,Ap,As,’s’) [num,den]=butter[N,wc,’s’] （2）Chebyshev I型低通滤波器 [N,wc]=cheb1ord(wp,ws,Ap,As,’s’) [num,den]=cheby1[N,Ap,wc,’s’] （3）Chebyshev II型低通滤波器 [N,wc]=cheb2ord(wp,ws,Ap,As,’s’) [num,den]=cheby2[N,As,wc,’s’] （4）椭圆型低通滤波器 [N,wc]=ellipord(wp,ws,Ap,As,’s’) [num,den]=ellip[N,Ap,As,wc,’s’] 模拟频率变换 （1）低通-高通 [numt,dent]=lp2hp(num,den,w0) （2）低通-带通 [numt,dent]=lp2bp(num,den,w0,B) （3）低通-带阻 [numt,dent]=lp2bs(num,den,w0,B) 3. 脉冲响应不变法 [numd,dend]=impinvar(num,den,Fs) 双线性变换法 [numd,dend]=bilinear(num,den,Fs) 二、FIR 窗函数法 w=hamming(N) w=blackman(N) w=kaiser(N,beta) [M,Wc,beta,ftype]=kaiserord(f,a,dev) 直接实现 h=firl(M,Wc,’ftype’,window) h=fi