fir数字滤波器设计与软件实现_FIR数字滤波器设计(中)

1f1377708b3f872203e6fc15c15b2bfa.gif大侠好,欢迎来到FPGA技术江湖,江湖偌大,相见即是缘分。大侠可以关注FPGA技术江湖,在“闯荡江湖”、"行侠仗义"栏里获取其他感兴趣的资源,或者一起煮酒言欢。

今天给大侠带来FIR数字滤波器设计,由于篇幅较长,分三篇。今天带来第二篇,FIR数字滤波器设计基础,包括FIR数字滤波器的特点、线性相位条件以及基本结构。话不多说,上货。

83918a22840a4658befcb48154eb827b.gif

数字滤波器的输入输出均为数字信号,信号通过数字滤波器后,可以改变频率成分的相对比例或滤除某些频率成分。数字滤波器可以分为IIR数字滤波器和FIR数字滤波器。

本篇只介绍FIR数字滤波器的设计,可以根据所给定的频率特性直接设计FIR数字滤波器。FIR数字滤波器在保证幅度特性满足要求的同时,能够做到严格的线性特性。

本篇采用了窗函数法、频率采样法以及基于firls函数和remez函数的最优化方法设计FIR滤波器。对FIR滤波器进行了详细的理论分析,并且对应于每种方法都给出了设计实例。通过编写MATLAB语言程序,运行程序,得到幅频和相频特性图。

对于窗函数和firls函数设计的滤波器,还通过建立Simulink系统模块进行仿真,观察滤波器滤波情况。

FIR数字滤波器设计基础

一、FIR数字滤波器的特点

FIR滤波器在保证幅度特性的同时,很容易做到严格的线性相位特性。

在数字滤波器中,FIR滤波器的最主要特点是没有反馈回路,故不存在不稳定的问题。

同时,在幅度特性可以任意设置的同时,保证了精确的线性相位。

稳定和线性相位是FIR滤波器的突出优点。

另外还有以下特点:

  • 设计方式是线性的;
  • 硬件容易实现;
  • 滤波器过渡过程具有有限区间;
  • 相对IIR滤波器而言,阶次较高,其延迟也要比同样性能的IIR滤波器大得多。

二、FIR数字滤波器的线性相位条件

设滤波器单位脉冲响应的长度为N,系统函数为关系式(2-1),如下:

f656c252f05271ed3966d577e0e44112.png

由此式可见,H(z)是81b7edc2fa250568c123de1f9fae7309.png的(N-1)次多项式,它在Z平面上有(N-1)个零点,原点z=0是(N-1)阶重极点,位于r =1的单位圆内,系统永远稳定。稳定性和线性相位特性是FIR滤波器的突出优点。

 FIR滤波器的设计任务是选择有线长度的h(n),使传输函数402e9e115cde86fde5fcb2fa5e76025f.png满足要求。

线性相位条件为关系式(2-2),如下:

8aa9508b5466699e7ef0da42de280ef0.png

对于长度为N的h(n),传输函数为关系式(2-3),如下:

17d2510f87195491bd28668ec7b4350e.png

上式中,8d976baae5e8bbd8fc16697ab641f00c.png称为幅度特性,ad4440ae54bfc77c3e030bc2c4a17260.png称为相位特性。线性相位是指相位函数ad4440ae54bfc77c3e030bc2c4a17260.png满足如下特性:

a33e37cda32c191dca77afa3868ee257.png

44233985013541cb1df4c865aaaba14d.png

8b0c89848912cb664b3ab1141880ac29.png 是起始相位, ad116aa631698acac834b6c92c8a29fa.png 为常数,一般称第一种情况为第一类线性相位,称第二种情况为第二类线性相位。

满足第一类线性相位的充要条件是:h(n)为实序列,并且对(N-1)/2偶对称,

即:

7704a11bddd6625983bbb211e545bc50.png

满足第二类线性相位的充要条件是:h(n)为实序列,并且对(N-1)/2奇对称。

即:

9c62818de6ec0169782521f5b77d5dc4.png

三、FIR数字滤波器的基本结构

FIR滤波器的基本结构有以下几种:直接型、级联型、线性相位型、频率采样型。

1、直接型

设FIR滤波器的单位冲击响应h(n)为一个长度为N的序列,则滤波器系统函数为关系式(2-4),如下所示:

6d93b303b6690a7771a298cef99634f5.png

表示这一系统输入输出关系的差分方程为关系式(2-5),如下所示:

fb9c9aafde1a0a433de1188e90012cc4.png

直接由差分方程可得出对应的网络结构如图2-1所示:

678aa26e50621d7a3332da709bb3832d.png

图2-1  FIR滤波器的直接型结构

直接型结构的优点:简单直观,乘法运算量较少。

缺点:调整零点较难。

2、级联型

当需要控制滤波器的传输零点时,可将H(z)分解为实系数二阶因子的乘积形式,见关系式(2-6),如下所示:

c85fd04ef1cb9a566ac5af33a3b141d2.png

上式中,752b8fafb6b74a3eb463eec17e5104d5.png2c77b93dd67ce5642d744c21c7eda2c5.png的 2d7532bbd7747fd9ef7d942cec0991ca.png变换,bc7b7a960471a27dc044543279a3be92.pngf3e0626df274f672031ee7d2e90c4595.png31ca57eb8e34678a19e9a9e7c839ce85.png为实数。级联型结构如图2-2所示:

059cd9d1f5c17aea54aba5c7d3c5eaf8.png

图2-2  FIR滤波器的级联型结构

该结构的优点:调整零点比直接型方便。

缺点:cb08d7ecd0cb08fcf6ff5debfd4cef2b.png中的系数比直接型多,因而需要的乘法器多。当cb08d7ecd0cb08fcf6ff5debfd4cef2b.png的阶次高时,也不易分解。

3、线性相位型结构

FIR滤波器的线性相位结构有偶对称和奇对称,不论6006e787b8c9a614041ce59ec71837aa.png为偶对称还是奇对称都有:

当N为偶数时,系统函数为关系式(2-7),如下所示:

420d34b63eb05a01aba8944ee24b9909.png

当N为奇数时,系统函数为关系式(2-8),如下所示:

de2991162a4776dc7d5f4b04344b51ca.png

对这两种情况,都可以用FIR直接型实现,其信号流图如图2-3所示。

76c61b2008aa7fa4424688ac0eba474c.png

(a)N为偶数

9ce76a7804a6d84afa3899e6420fbe80.png

(b)N为奇数

图2-3  线性相位型结构

    这种结构在本质上是直接型,但乘法次数比直接型省了一半。

4、频率采样型

频率采样型结构是一种用系数将滤波器参数化的一种实现结构。一个有限长序列可以由相同长度频域采样值惟一确定。

系统函数在单位圆上作N等分取样就是单位取样相应h(n)的离散傅里叶变换。188576437b7e2822f8bdb01a62dc27cc.png与系统函数之间的关系可用内插公式表示,为关系式(2-9),如下所示:

c69e78fe6a13c902bab520656ea6131e.png

上式中a456e0953d644f127921e00c7174a066.png18d5b60ce03beef693ebfd67f0f6d4ca.png

这样,abce6c731d75fd75ed5ca8f245755ffc.png是由梳状滤波器和N个一阶网络的并联结构进行级联而成的,其网络结构(信号流图)如图2-3所示。是一个梳状网络,其零点为:

faabe3fd4c8a79a6b8b48d46e2b15f7d.png,    k= 0, 1,2…,N-1

刚好和极点一样,等间隔地分布在单位圆上。理论上,极点和零点相互抵消,保证了网络的稳定性。

84f6666bbedf5760b6b7c4dab553b2f5.png

图2-5  FIR滤波器的频率采样结构

频率采样结构的优点:

1)在频率采样点c2e5edad19852c87f8a4afdab2c7da11.png175710c99fadc38764d5f274051c192b.png,只要调整b39c90b65d60ad06d983a24108c619b0.png就可以有效地调整频响特性。

2)只要47d85eec58095471513f8c843e0d390c.png长度N相同,对于任何频响,其梳状滤波器部分和N个一阶网络部分完全相同,只是各支路增益b39c90b65d60ad06d983a24108c619b0.png不同。相同部分便于标准化、模块化。

缺点:

1)寄存器长度都是有限的,零、级点可能不能正好抵消,造成系统不稳定。

2)当N很大时,其结构很复杂,需要的乘法器和延时单元很多。

83918a22840a4658befcb48154eb827b.gif

第二篇就到这里,下一篇带来第三篇,FIR数字滤波器设计,包括窗函数法设计FIR滤波器、频率采样法设计FIR滤波器以及基于firls函数和remez函数的最优化方法设计FIR滤波器。

END后续会持续更新,带来Vivado、 ISE、Quartus II 、candence等安装相关设计教程,学习资源、项目资源、好文推荐等,希望大侠持续关注。大侠们,江湖偌大,继续闯荡,愿一切安好,有缘再见!

往期推荐

  • 通信的数学原理-香农-英文 电子版

  • 数字电子技术基础(第六版)教材 闫石 电子版

  • FPGA 之 SOPC 系列 汇总篇

54a0d9d752ceb75f00706dd6ce94cb1e.gif

6331f844e58211d650f121bd9ab3cc32.png

FPGA技术江湖广发江湖帖

无广告纯净模式,给技术交流一片净土,从初学小白到行业精英业界大佬等,从军工领域到民用企业等,从通信、图像处理到人工智能等各个方向应有尽有,QQ微信双选,FPGA技术江湖打造最纯净最专业的技术交流学习平台。

FPGA技术江湖微信交流群

9dbe4bff2cebb5718a306ccfbe204395.png

加群主微信,备注职业+方向+名字进群

FPGA技术江湖QQ交流群

78334433e4d19687e736c774208342c2.png

备注地区+职业+方向+名字进群

5423e4e3c9819ff9890decbc792ae00b.png

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值