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

fir数字滤波器设计与软件实现_FIR数字滤波器设计(上)_第1张图片大侠好,欢迎来到FPGA技术江湖,江湖偌大,相见即是缘分。大侠可以关注FPGA技术江湖,在“闯荡江湖”、"行侠仗义"栏里获取其他感兴趣的资源,或者一起煮酒言欢。

今天给大侠带来FIR数字滤波器设计,由于篇幅较长,分三篇。今天带来第一篇,数字滤波器介绍,包括数字滤波器概述、分类以及设计指标。话不多说,上货。

fir数字滤波器设计与软件实现_FIR数字滤波器设计(上)_第2张图片

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

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

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

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

数字滤波器

一、数字滤波器的概述

所谓数字滤波器,是指输入输出均为数字信号,通过一定的运算关系,改变输入信号中所含频率成分的相对比例,或滤除某些频率成分的器件。数字滤波器具有稳定性高、精度高、灵活性大等突出优点。对于数字滤波器而言,若系统函数为H(z),其脉冲响应为h(n),输入时间序列为x(n),则它们在时域内的关系式(1-1)如下:

1f1ca0f75019f08f31c1b748619c6315.png

在Z域内,输入和输出存在如下关系式(1-2):

d50d2c4800c29b0376e8cccb959c2917.png

上式中, X(z)、Y(z)分别为x(n)和y(n)的Z变换。

在频域内,输入和输出则存在如下关系(1-3):     

b1e53537f503bcb4898ab58835128e1c.png

上式中,是数字滤波器的频率特性;

41b7b053375d8d1abe99f0f1b3fb03c2.png99e36bf0c2803155783ab66c815bee39.png分别为x(n)和y(n)的频谱,而ded4de9a50733359174d9ebac21baf59.png为数字角频率。

二、数字滤波器的分类

数字滤波器可以有很多种分类方法,但总体上可分为两大类。一类称为经典滤波器,即一般的滤波器,其特点是输入信号中的有用成分和希望滤除的成分占用不同的频带,通过合适的选频滤波器可以实现滤波。

例如,若输入信号中有干扰,信号和干扰的频带互不重叠,则可滤出信号中的干扰得到纯信号。但是,如果输入信号中信号和干扰的频带相互重叠,则干扰就不能被有效的滤除。

另一类称为现代滤波器,如维纳滤波器、卡尔曼滤波器等,其输入信号中有用信号和希望滤除的频带成分重叠。对于经典滤波器,从频域上也可以分为低通、高通、带通和带阻滤波器。

从时域特性上看,数字滤波器还可以分为有限脉冲响应(FIR,finite impulse response)数字滤波器和无限脉冲响应(IIR, infinite impulse response)数字滤波器[5]。

对于有限脉冲响应(FIR)数字滤波器,其输出y(n)只取决于有限个过去和现在的输入,x(n),x(n-1),…,x(n-m),滤波器的输入输出关系可表示为表达式(1-4),如下:

4844469c883a6eb800b28d238d7ab962.png

对于无限脉冲响应(IIR)数字滤波器,它的输出不仅取决于过去和现在的输入,而且还取决于过去的输出,其差分方程为表达式(1-5),如下:

552614534d1cf677728f676bd965338b.png

该差分方程的单位冲激响应是无限延续的。

三、数字滤波器设计指标

设数字滤波器的传输函数用下式(1-6)表示:

a48f9e730be13146fc46917e55630c21.png

式中,f363925447e5a1ffc5927b3d60531366.png为幅频特性,81dd8fea5e0e29b41ad35892674c1634.png为相频特性。幅频特性表示信号通过滤波器后各频率成分的衰减情况,相频特性则反映各频率成分通过滤波器后在时间上的延时情况。通常,选频滤波器的指标要求都以幅频特性给出,对相频特性不作要求,如果需要对输出波形有严格要求,如语音合成、波形传输等,则要求设计线性相位数字滤波器。

数字滤波器的参数指标是04f568fd2016ad0663b95c09777bcfa3.png1033b07f254d4a2dd5888dce442a30cf.pngcdd64b5ffb24d41d57369f590730dd13.png36e46714e813a22d5db55a338b79b5c4.png04f568fd2016ad0663b95c09777bcfa3.png1033b07f254d4a2dd5888dce442a30cf.png分别称为通带截止频率和阻带截止频率。通带和阻带内允许的衰减一般用分贝数表示,通带内允许的最大衰减用表示,阻带内允许的最小衰减用表示,和分别定义为关系式(1-7)和关系式(1-8):

80e4d8116d6928d6b5585659b600268f.png

ca76501bd8ecf34b1dd433a77c243b38.png

式中均假定 a891ec39527cb1109970f3bd84ea64a1.png 已被归一化为1。

fir数字滤波器设计与软件实现_FIR数字滤波器设计(上)_第3张图片

第一篇就到这里,下一篇带来第二篇,FIR数字滤波器设计基础,包括FIR数字滤波器的特点、线性相位条件以及基本结构。

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

往期推荐

  • 在FPGA中何时用组合逻辑或时序逻辑

  • 嵌入式中位宽的计算

  • FPGA Vivado设计流程

fir数字滤波器设计与软件实现_FIR数字滤波器设计(上)_第4张图片

c08e32dc280042099d08dedc3f4e6064.png

FPGA技术江湖广发江湖帖

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

FPGA技术江湖微信交流群

fir数字滤波器设计与软件实现_FIR数字滤波器设计(上)_第5张图片

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

FPGA技术江湖QQ交流群

fir数字滤波器设计与软件实现_FIR数字滤波器设计(上)_第6张图片

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

9ce0edf5a6ac452aae9fa8d446945683.png

你可能感兴趣的:(fir数字滤波器设计与软件实现,fir滤波器,fir滤波器c++程序,iir数字滤波器,iir数字滤波器设计及软件实现,matlab,带通滤波器)