数字信号处理课程实验报告.docVIP

数字信号处理课程实验报告 课题名称： IIR滤波器相位校正实验 学院：计算机与电子信息学院 专业：电子信息工程 班级：创新121 姓名：程小宁 一、实验内容与分析 1、实验目的和内容 1）利用MATLAB设计一个IIR滤波器； 2）结合课本关于全通滤波器特性知识（课本p128），在IIR滤波器后级联一个全通相位滤波器进行相位校正，使此滤波器最终实现线性相位特性； 3）分别使用相位校正前后两滤波器实现对某一信号的处理； 4）画出IIR滤波器、全通滤波器、相位校正后滤波器的幅度频率特性曲线、相位频率特性曲线，信号时域波形、信号的幅度频率特性曲线、相位频率特性曲线； 5）详述实验设计原理，分析相位校正前后两类滤波器对信号处理后的区别。 2、实验的分析 1）、IIR滤波器的设计 通过对实验内容的理解，我们首先需要设计一个IIR滤波器，对课本第六章的学习我们知道IIR数字滤波器有两种设计方法：间接设计法和直接设计法。间接设计法中有巴特沃斯滤波器，切比雪夫I型、II型滤波器，椭圆滤波器和贝塞尔滤波器五种。我们选择设计切比雪夫II型低通滤波器，其中的技术指标为：通带边界频率fp=1000Hz，阻带边界频率fs=2000 Hz，阻带最小衰减As=40 dB，通带最大衰减Ap=1 dB。 全通滤波器的设计 全通滤波器的幅度特性是在整个频带上均等于常数，或者等于1.信号通过全通滤波器后，其输出的幅度特性保持不变，仅相位发生变化。由于IIR滤波器后需要级联一个全通相位滤波器，使整个系统实现线性相位特性，为了求解全通滤波器的参数，我们先假设整个系统具有线性相位特性，再根据已经设计好了的切比雪夫II型滤波器的系统参数，求解全通滤波器的参数。 二、实验的过程 切比雪夫II型滤波器的设计过程 在确定了滤波器的参数之后，我们运用cheb2ord函数计算模拟低通滤波器的最小阶数；然后用cheby2计算滤波器传输函数的系数。然后运用脉冲响应不变法将模拟低通滤波器转换成数字滤波器。这样我们就设计出了满足给定参数的切比雪夫II型滤波器。此部分的程序为： [N2,ws2]=cheb2ord(wp,ws,Rp,As,s); [B2,A2]=cheby2(N2,As,ws2,s); [bz,az]=impinvar(B2,A2,Fs); 程序中，图1为切比雪夫II型滤波器的幅频特性图和相频特性图。 全通滤波器的设计过程 因为题目要求是设计一个IIR滤波器级联一个全通滤波器，使其最终能够满足线性相位，所以我们先假设滤波器满足线性相位，反求全通滤波器的参数。可知全通滤波器不改变信号的幅度，只改变信号的相位。假设全通滤波器的幅度为1，则由级联公式，得：,H1为切比雪夫II型滤波器的系统函数，其幅度为，相位角为angle(H1)，H2为所求全通滤波器的系统函数，幅度为1，H3为最终系统的系统函数，幅度为切比雪夫II型滤波器的幅度与全通滤波器幅度的乘积，即H3的幅度为，相位角为2*pi*f/Fs。所以可得全通滤波器的系统函数为：H2=exp(j*(2*pi*f/Fs-unwrap(angle(H))))。 用切比雪夫II型滤波器对某一信号进行处理 假设输入信号为：x=sin(2*pi*f1*t)+0.5*sin(2*pi*f2*t)。其中输入信号含有的频率f1=1000Hz，f2=3000Hz，利用filtfilt函数对输入信号进行滤波，即可以得到通过切比雪夫II型滤波器的输出信号。对于输入信号本身，我们利用快速傅里叶算法（FFT）可以得到输入信号的幅频特性图额相频特性图。所以程序中图2和图3分别为输入信号通过切比雪夫II型滤波器之后和输入信号本身的幅频、相频特性图。 将切比雪夫II型滤波器处理过的信号再通过全通滤波器 当我们得到通过切比雪夫II型滤波器之后的输出信号，再次用filtfilt函数对其进行滤波，同上面做法不同的是，此时filtfilt函数的参数是全通滤波器的参数，所以最终的系统的输出函数可以得到线性相位。 最终输出信号 三、实验中遇到的困难和解决办法 本实验的思路并不难，但是要做出来还是需要一定的时间，首先第一步设计IIR滤波器时很顺利，但是后来的全通滤波器，我们就遇到了障碍。我们以为任意一个全通滤波器都可以，然而后来做出来之后整个系统并不能实现线性相位。在去图书馆查阅了相关资料以及询问老师和学长之后，我们想到了可以利用结论反推条件，也就是先假设整个系统是可以实现线性相位的，再反求全通滤波器的参数，这样，我们就设计出了最终的总系统。 四、实验总结和存在的不足 通过本次实验，我们对数字信号处理这门课有了更深一层的理解，虽然平时上课的时候我们也做过一些课后的MATLAB仿真作业，但是远比这次的实验要简