DSP课程设计 MATLAB实现切比雪夫滤波器分析与设计.docVIP

目录 1 课题综述 1 1.1 课题来源 1 1.2 预期目标 1 1.3 面对的问题 1 1.4 需要解决的关键技术 1 2 系统分析 2 2.1 涉及的基础知识 2 2.2 解决的基本思路 2 2.3 总体方案 2 2.4 功能模块框图 2 3 详细设计 3 3.1 巴特沃斯低通滤波特性(MATLAB) 3 3.2 巴特沃斯高通滤波特性(MATLAB) 4 3.3 切比雪夫低通滤波特性(MATLAB) 4 3.4 切比雪夫高通滤波特性(MATLAB) 4 4 程序调试 5 4.1 巴特沃斯低通滤波特性 5 4.2 巴特沃斯高通滤波特性 6 4.3 切比雪夫低通滤波特性 8 4.4 切比雪夫高通滤波特性 9 5 运行与测试 10 5.1 选择音频文件（WAV） 10 5.2 滤波后音频特点 10 6 全文代码设计 10 总 结 14 致 谢 15 参 考 文 献 16 1 课题综述 1.1 课题来源 随着数字集成电路，设备和系统技术的快速进步，通过数字方法进行信号处理已变得越来越有吸引力。大规模一般用途的计算机和特殊用途硬盘的高效性，已使得实时滤波既实用又经济。目前主要有两类滤波器，模拟滤波器和数字滤波器，它们在物理组成和工作方式上完全不同，而模拟滤波器的技术发展已相当成熟，所以研究的重点基本上放在了数字滤波器上。滤波器的功能是用来移除信号中不需要的部分，比如随机噪声；或取出信号中的有用部分，如位于某段频率范围内的成分。目前随着计算机技术和数字信号处理器芯片的发展，使我们更为便利地识别和提取各种各样的信号。因此研究不同数字滤波器的设计原理和稳定性分析对于满足军事、航空、民营等等各个领域的信号处理要求具有十分重要的意义。 1.2 预期目标 能够理解并掌握无限脉冲响应数字滤波器（IIR）的机理，分析IIR滤波器的结构特性，观察IIR滤波器的频域特性。能够通过MATLAB实现巴特沃斯，切比雪夫的高通低通滤波器的幅频响应，相频响应，以及零极点的图像的勾画。 1.3 面对的问题 通过MATLAB实现切比雪夫滤波器的设计，以及通过MATLAB对高通，带通，带阻滤波器的设计。 1.4 需要解决的关键技术 了解巴特沃斯以及切比雪夫滤波器理论上的设计过程，从理论推至代码程序的实现。掌握高通低通带通带阻滤波器的幅频响应，相频响应，以及零极点的图像的意义，并且能够寻找到实现他们的MATLAB函Wp,Ws,Rp,Rs技术指标，求出滤波器阶数N。根据阶数N查表，得到归一化传输函数Ha(p)将Ha(p)去归一化。将p=s/ Ωc代入Ha(p)，得到实际的滤波器传输函数H(S)。 2.3 总体方案 2.3.1 巴特沃斯滤波器MATLAB实现 其调用格式为： [n W]=buttord(OP,OS,Rp,Rs,s); [B A]=butter(n,W,s); [num,den]=impinvar(B,A,f1); 2.3.2 切比雪夫滤波器MATLAB实现 其调用格式为： [n W]=cheb1ord(wp1,ws1,Rp,Rs,s); [B A]=cheby1(n,1,W,type,s); [num,den]=bilinear(B,A,f1); 功能模块框图 开始 ↓ 读入数字滤波器技术指标 ↓ 将指标转换成归一化模拟低通滤波器的指标 ↓ 设计归一化的模拟低通滤波器阶数N和3db截止频率 ↓ 模拟域频率变换，将G(P)变换成模拟带通滤波器H(s) ↓ 用脉冲响应不变法将H(s)转换成数字带通滤波H(z) ↓ 输入信号后显示相关结果 图 2.1 3 详细设计 3.1 巴特沃斯低通滤波特性(MATLAB) Wp=2*pi*f2/f1; Ws=2*pi*f3/f1; OmigaP=Wp*f1; OmigaS=Ws*f1; [n W]=buttord(OmigaP,OmigaS,Rp,Rs,s); [B A]=butter(n,W,s); [num,den]=impinvar(B,A,f1); 3.2 巴特沃斯高通滤波特性(MATLAB) Wp=2*pi*f2/f1; Ws=2*pi*f3/f1; OmigaP=Wp*f1; OmigaS=Ws*f1; [n W]=buttord(OmigaP,OmigaS,Rp,Rs,s); [B A]=butter(n,W,high,s); [num,den]=impinvar(B,A,f1); 3.3 切比雪夫低通滤