双线性变换法方向Elliptic型IIR数字低通滤波器.docVIP

课 程 设 计 课程设计名称：数字信号处理课程设计 专 业 班 级 ：电信1002 学 生 姓 名 ：周瑞峰 学号 ：2 指 导 教 师 ：李相国 课程设计时间：2018年6月 数字信号处理专业课程设计任务书 学生姓名 周瑞峰 专业班级 电信1002 学号 2 题 目 基于双线性变换法设计Elliptic型IIR数字低通滤波器。 课题性质 其它 课题来源 自拟课题 指导教师 李相国 同组姓名 主要内容 设计一个具有如下性能指标的Elliptic型IIR数字低通滤波器：=0.25，=0.4，=0.01，=0.001；分析最后设计结果性能。 任务要求 掌握双线性变换法设计数字滤波器的基本流程； 掌握Elliptic型模拟低通原型滤波器； 求出所设计滤波器的传输函数； 绘制所设计滤波器的幅频特性图。 参考文献 [1]胡广书. 数字信号处理—理论、算法与实现[M]. 北京: 清华大学出版社, 1997. [2]R. Lyons. Understanding Digital Signal Processing [M]. 2nd ed. Prentice Hall PTR., 2004. [3]A.V.奥本海姆, R.W.谢弗 and J.R.巴克. 离散时间信号处理[M]. 第二版. 西安交通大学出版社, 2001. [4]S. K. Mitra. Digital Signal Processing: A Computer-Based Approach[M]. 3rd ed. McGraw-Hill, 2005. 审查意见 指导教师签字： 教研室主任签字： 年 月 日 1 需求分析 此次实验是根据规定的要求=0.25，=0.4，=0.01，=0.001，利用双线性变换法设计一个椭圆形IIR 数字低通滤波器。根据题目可看出设计要求有几点： 1.双线性变换法2.椭圆型3.IIR 4.数字低通5.各项参数如上 2 概要设计 开始 开始 根据已给参数算出编程用的参数 根据已给参数算出编程用的参数 低通滤波器设计与实现 低通滤波器设计与实现用ellipord，ellip，filter函数实现） 低通滤波器设计与实现绘图 结束 3 运行环境 硬件：计算机一台 软件环境：WINDOWS XP MATLAB 4 开发工具和编程语言 开发工具：MATLAB 编程语言：c语言 5详细设计 参数配置程序： Wp=0.25*pi。 通带截止频率 Ws=0.4*pi。 阻带截止频率 Rp=-20*log10(1-0.01。 峰值通带波纹 Rs=-20*log10(0.001。 最小阻带衰减 椭圆型低通数字滤波器的生成程序： [N,Wc]=ellipord(Ws,Wp,Rp,Rs。 ellipord椭圆滤波器最小阶数和截止频率的确定函数。 语法：[N,Wc]=ellipord(wp,ws,Rp,Rs,‘s’ Rp、Rs分别为通带最大波纹和阻带最小衰减；wp、ws分别为为通带边界频率和阻带边界频率，单位为rad/s。这四个参数为滤波器的基本性能指标 返回值N为实现给定性能指标的最小阶数，Wc为数字滤波器的截止频率，单位为rad/s。s表示模拟滤波器，缺省时该函数适用于数字滤波器，本程序中‘s’缺省。此时wp及ws均为0~1之间的数值，即频率与抽样频率的比值。 [b,a]=ellip(N,Rp,Rs,Wc。 求传递函数的分子分母系数，b，a分别为传递函数的分子分母的系数，N是实现给定性能的最小阶数，Rp，Rs，分别为峰值通带波纹，最小阻带衰减，Wc，为数字滤波器的截止频率。 [hw,w]=freqz(b,a，N。 b和a分别为离散系统的系统函数分子、分母多项式的系数向量，返回量hw则包含了离散系统频响在 0~pi范围内N个频率等分点的值其中N为正整数），w则包含了范围内N个频率等分点。调用默认的N时，其值是512，本程序N用默认值。 subplot(3,2,1。 subplotm,n,p）或者subplotm n p）。 subplot是将多个图画到一个平面上的工具。其中，m表示是图排成m行，n表示图排成n列，也就是整个figure中有n个图是排成一列的，一共m行，如果m=2就是表示2行图。p表示图所在的位置，p=1表示从左到右从上到下的第一个位置。 plot(w/pi,20*log10(abs(hw。 plot(x,y?以x?元素为横坐标值，y?元素为纵坐标值绘制曲线。本程序中，x是由freqz函数得到的512个等分点除以pi的值，y中abs函数式求绝对值，即求hw）的绝对值，并依此为纵坐标作图。 grid。在画图中显示网格线 axis([0,1,-80,5]。设