基于MATLAB的切比雪夫I型模拟低通滤波器设计.docVIP

课 程 设 计 课程设计名称：数字信号处理课程设计 专 业 班 级 ： 电信0604 学 生 姓 名 ： 学 号 ： 20064300430 指 导 教 师 ： 课程设计时间：2009.6.8-2009.6.14 数字信号处理 专业课程设计任务书 学生姓名 专业班级 电信0604 学号 20064300430 题 目 基于MATLAB的切比雪夫I型模拟低通滤波器设计 课题性质 其他 课题来源 自拟课题 指导教师 同组姓名 主要内容 根据已学的知识并结合MATLAB来设计一个切比雪夫I型模拟低通滤波器，技术指标如下： 通带截止频率：fp=1000Hz，通带最大衰减： 阻带截止频率：fs=1500Hz，阻带最小衰减： 画出滤波器的幅频、相频特性曲线。 任务要求 写出设计原理和设计思路，画出程序流程图 用MATLAB画出幅频特性图 用MATLAB画出相频特性图 用MATLAB画出零极点图 参考文献 1．程佩青著，《数字信号处理教程》，清华大学出版社，2001 2．Sanjit K. Mitra著，孙洪，余翔宇译，《数字信号处理实验指导书（MATLAB版）》，电子工业出版社，2005年1月 3．郭仕剑等，《MATLAB 7.x数字信号处理》，人民邮电出版社，2006年 审查意见 指导教师签字： 教研室主任签字： 年 月 日 说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页 内容包括： 一 设计内容与技术要求 设计一个切比雪夫I型模拟低通滤波器，满足指标如下：通带截止频率：fp=1000Hz，通带最大衰减：,阻带截止频率：fs=1500Hz，阻带最小衰减：,写出设计原理和设计思路，画出程序流程图，用MATLAB编写程序并画出幅频特性图，相频特性图和零极点图。 二 设计原理及设计思路 1 设计原理 Chebyshev滤波器是在通带或阻带上频率响应幅度等波纹波动的滤波器。在通带波动的为“ChebyshevI型滤波器”，在阻带波动的为“ChebyshevII滤波器”。滤波器在过渡带比Butterworth滤波器的衰减快，但频率响应的幅频特性不如后者平坦。Chebyshev滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波。 Chebyshev滤波器N阶多项式为 ： 当|x|1时 CN(x)= cos[Narcoses(x)] 当|x|1时 CN(x) = Cosh[Narccosh(x)] N不论为何值都有TN(1)=1; 当x1,TN(x)单调上升。此外，Chebyshey多项式满足下列递推公式 TN+1(x)=2x TN(x)- TN-1(x),N=1,2,… Chebyshey滤波器的滤波特性具有下列特点： （1）所有曲线在Ω=Ωc时通过点，因而把Ωc定义为Chebyshey滤波器的截至频率。 （2）在通带内≤1，|H(Ω)|在1和之间变化；在通带外，|Ω/Ωc|1,特性呈单调下降，下降速度为20NdB/dec。 （3）N为奇数，|H(Ω)|=。通带内误差分布是均匀的，实际上这种逼近称为最佳一致逼近。因此通带等波纹滤波器是在通带内以最大误差最小化对理想低通滤波器的最佳一致逼近。 （4）由于滤波器通带内有起伏，因而使通带内的相频特性也有相应的起伏波动，即相位是非线性的，这给信号传输时带来线性畸变，所以在要求群时延为常数时不宜采用这种滤波器。 切比雪夫I型滤波器的振幅平方函数为： 图（1） 切比雪夫滤波器的振幅平方特性 2 切比雪夫I型模拟低通滤波器设计步骤如下： （1）由通带最大衰减，阻带最小衰减、通带截止频率，阻带截止频率，确定波纹系数和系统的阶数； （2）计算出系统函数中的分子、分母； （3）画出幅频特性图，相频特性图和零极点图； 三 程序流程图 给定模拟低通滤波器技术指标 给定模拟低通滤波器技术指标 计算波纹系数和系统阶数 绘制幅频响应 调用切比雪夫I型滤波器的文件 绘制相频响应 绘制系统零极点 四 程序源代码 clc;clear all Rp = 1; % 通带最大衰减 As = 25; % 阻带最小衰减 fp=1000; % 通带截止频率 fs=1500; % 阻带截止频率 wp=2*pi*fp; ws=2*pi*fs; [N, Wn] =cheb1ord(wp,ws,Rp, As,s); [num,den] =cheby1(N,Rp,Wn,s); [h,w]=freqs(num,den); g=20*log10(abs(h)); g1=a