dsp带通滤波器设计.docx

Y Y z N 1 X z i 0 -- 一、 设计要求 试用 DSP设计 FＩＲ滤波器，分别实现带通的功能 , 具体要求如下： ① 滤波器的阶数≥ 5，截止频率自行选定，滤波系数用 MＡＴ LAB确定。 ② 编制 C54XDSP实现 FIR 滤波器的汇编源程序。 ③ 用软件仿真器完成上述程序的模拟调试。 ④ 以数据文件形式自行设定滤波器输入数据，以数据文件形式输出滤波结果， 并与输入数据进行比较分析。 用软件仿真器有关工具显示ＦＩ R 滤波器的输入输 出波形 , 以证明滤波器滤波性能。 二、 FI Ｒ滤波器的基本原理 数字滤波是将输入的信号序列 , 按规定的算法进行处理 , 从而得到所期望的 输出序列。 一个线性位移不变系统的输出序列 y n 和输入序列 x n 之间的关系 , 应满足常系数线性差分方程。 FIR 滤波器的差分方程为： N 1 y(n) h(m) x(n m) m0 m FIR 滤波器的传递函数为 : H z bi z i 直接由差分方程得出的实现结构如图 1 所示: 图 1 横截型 (直接型﹑卷积型 ) ＦＩ R 滤波器的结构图 由上面的公式和结构图可知， FIR 滤波算法实际上时一种乘法累加运算。它 不断地从输入端读入样本值 x n ，经过（ z 1 ）后做乘法累加 ,输出滤波结果 y n 三、 FIR 滤波器的设计 FIR 滤波器的设计方法主要有窗函数法和频率采样法。 其中 , 窗函数法是最基 本的方法。本次设计采用窗函数法。 一般是先给定所要求的理想滤波器频率响应 H d (ej ) , 由 H d(ej ) 导出 hd ( n) , 我们知道 , 理想滤波器的冲击响应 hd ( n)是无限长的非因果序列 , 而我们要设计的 是 hd ( n)是有限长的 hd (n) ，设： F Ｉ R 滤波器 , 所以要用有限长序列 hd (n)来逼近无限长序列 21hd (n) 2 1 H d (ej )ej d ?（３ -1） -- 2 2 -- 常用的方法是用有限长度的窗函数 ｗ （ ｎ )来截取 即： h(n) (n)hd (n) ?（3－２ ) 这里窗函数就是矩形序列 RN( n) ，加窗以后对理想低通滤波器的频率响应将 产生什么样的影响呢 ?根据在时域是相乘关系 , 在频域则是卷积关系 : ?H (ej ) 1 H d (ej )WR[ ej ( ) ]d ?（3-3） 其中， 为矩形窗谱 , H (ej ) 是 FIR 滤波器频率响应。 通过频域卷积过程看 的幅度函数 H( ω) 的起伏现象 , 可知 , 加窗处理 后 , 对理想矩形的频率响应产生以下几点影响 : (1)使理想频率特性不连续点处边沿加宽，形成一个过渡带，其宽度等于窗的频 率响应的主瓣宽度。 (2）在截止频率的两边的地方即过渡带的两边，出现最大的肩峰值 , 肩峰的两侧 形成起伏振荡，其振荡幅度取决于旁瓣的相对幅度 , 而振荡的多少，则取决于旁 瓣的多少。 (3）改变 N, 只能改变窗谱的主瓣宽度 , 改变 ω 的坐标比例以及改变的绝对值大 小 , 但不能改变主瓣与旁瓣的相对比例 (此比例由窗函数的形状决定 ) 。 (4)对窗函数的要求： ａ、窗谱主瓣尽可能窄，以获取较陡的过渡带； ｂ、尽量减小窗谱的最大旁瓣的相对幅度；即能量集中于主瓣 , 使肩峰和波纹减 小 , 增大阻带的衰减。 四、 带通滤波器设计指标 本次要设计一个Ｆ IＲ带通滤波器，其采样频率 fs＝3 ００ 0 ０Ｈ z，通带为 3000Hz～５ 000Ｈ z,即保留频率在３ 0 ００ Hz~ ５０ 00Hz 范围内的信号成分， 幅度 失真小于１ｄ B；阻带边界频率为 2000Ｈ z 和 60 ００Ｈ z,衰减大于 40dB。同时， F Ｉ R 滤波器的阶数 N=32。 五、 设计思路及实现步骤 1.使用ma ｔ l ａ b获得输入信号数据， 这里输入信号为 f １ =1500,f2 ＝４ 0 ００ ,f -- -- ３ =70 ００，采样频率为 fs= ３ 00