(数字滤波器的设计.docVIP

实验报告 课程名称 数字信号处理实验 实验名称 数字滤波片的设计 姓名 田红亮 学号 201005020216 专业班级 电信1002班 实验日期 成绩 指导老师 （①实验目的②实验原理⑶主要仪器设备④实验内容与步骤⑤实验数据记录与处理⑥实验结果与分析⑦问题建议） 一、实验目的 1. 了解IIR数字滤波器的特点； 2． 掌握IIR数字滤波器的设计方法； 3. 掌握脉冲响应不变法和双线性变换法设计数字滤波器的方法。 二、实验所采用的功能函数 巴特沃斯滤波器阶数选择函数 （1）[N,wc]=buttord(wp，ws，αp，αs) 作用： 计算巴特沃斯数字滤波器的阶数N和3dB截止频率wc, wc为数字频率，单位rad。 说明： 调用参数wp，ws分别为数字滤波器的通带、阻带截止频率的归一化值，要求：0≤wp≤1，0≤ws≤1。αp，αs分别为通带最大衰减和组带最小衰减（dB）。当ws≤wp时，为高通滤波器；当wp和ws为二元矢量时，为带通或带阻滤波器，这时wc也是二元向量。 （2）[N,Ωc]=buttord(Ωp，Ωs，αp，αs，‘s’) 作用： 计算巴特沃斯模拟滤波器的阶数N和3dB截止频率Ωc。 说明：Ωp，Ωs，Ωc均为实际模拟角频率。 模拟频率f：每秒经历多少个周期，单位Hz，即1/s，信号的真实频率，可用于模拟信号和数字信号；模拟角频率Ω：每秒经历多少弧度，单位rad/s，通常只于模拟信号；数字频率w：每个采样点间隔之间的弧度，单位rad，通常只用于数字信号。关系：Ω=2pi*f；w = Ω*T=2pi*f/F。(F=1/Ts为采样频率，Ts为采样间隔) 实验报告 课程名称 数字信号处理实验 实验名称 数字滤波片的设计 姓名 田红亮 学号 201005020216 专业班级 电信1002班 实验日期 成绩 指导老师 （①实验目的②实验原理⑶主要仪器设备④实验内容与步骤⑤实验数据记录与处理⑥实验结果与分析⑦问题建议） 完整巴特沃斯滤波器设计函数 （1）格式： [b，a]=butter（N，wc，‘ftype’） 作用： 计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数向量b、a。 说明： 调用参数N和wc分别为巴特沃斯数字滤波器的阶数和3dB截止频率的归一化值，一般是调用buttord格式（1）计算N和wc。系数b、a是按照z-1的升幂排列。 （2）格式：[B，A]=butter（N，Ωc，‘ftype’，‘s’） 作用：计算巴特沃斯模拟滤波器系统函数的分子、分母多项式系数向量。 说明：调用参数N和Ωc分别为巴特沃斯模拟滤波器的阶数和3dB截止频率(实际角频率)，可调用buttord（2）格式计算N和Ωc。系数B、A按s的正降幂排列。 tfype为滤波器的类型： ◇ftype=high时，高通；Ωc只有1个值。 ◇ftype=stop时，带阻；Ωc=[Ωcl,Ωcu]，分别为带阻滤波器的通带3dB下截止频率和上截止频率。 ◇ ftype缺省时：若Ωc只有1个值，则默认为低通；若Ωc有2个值，则默认为带通；其通带频率区间ΩclΩ Ωcu。 3. 求离散系统频响特性的函数freqz()[H,w]=freqz(b,a,N) 说明：b和分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H则包含了离散系统频响在 0pi范围内N个频率等分点的值（其中N为正整数），w则包含了范围内N个频率等分点。调用默认的N时，其值是512。可以先调用freqz()函数计算系统的频率响应，然后利用abs()和angle()函数及plot()函数，绘制出系统的频响曲线。 实验报告 课程名称 数字信号处理实验 实验名称 数字滤波片的设计