dsp:FIR滤波器的设计

设计一低通FIR滤波器,滤波器技术指标要求:通带:2KHZ,阻带:10KHZ,系统采样频率20KHZ,通带与阻带至少有38DB的衰减。回答下列问题:

(1)FIR差分方程的一般形式

(2)设计满足上述要求的FIR滤波器,确定各参数。

(3)画出程序流程图。

(4)用汇编语言或者C语言编写满足上述要求的FIR滤波器程序。(如用C语言则假设从数据键盘输入,结果显示在屏幕上;如用汇编语言则假设从PA0口输入,结果从PA1口输出。)

第1个回答  2014-06-13
题目:利用DSP的FIR滤波器设计

数字处理器(DSP)有很强的数据处理能力,它在高速数字信号处理领域有广泛的应用,例如数字滤波、音频处理、图像处理等。相对于模拟滤波器,数字滤波器没有漂移,能够处理低频信号,频率响应特性可做成非常接近于理想的特性,且精度可以达到很高,容易集成等。使用可编程的DSP芯片实现数字滤波可以通过修改滤波器的参数十分方便地改变滤波器的特性,下面主要说明利用TMS320VC54x DSP芯片设计实现FIR数字滤波器。
设计目的意义
一个实际的应用系统中,总存在各种干扰,所以在系统设计中,滤波器的好坏将直接影响系统的性能。使用DSP进行数字处理,可以对一个具有噪声和信号的混合信号源进行采样,再经过数字滤波,滤除噪声,就可以提取有用信号了。所以说,数字滤波器是DSP最基本的应用领域,熟悉基于DSP的数字滤波器能为DSP应用系统开发提供良好的基础。
技术指标
1、数字滤波器的频率参数主要有:①通带截频:为通带与过渡带的边界点,在该点信号增益下降到规定的下限。②阻带截频:为阻带与过渡带的边界点,在该点信号衰耗下降到规定的下限。③转折频率:为信号功率衰减到1/2(约3dB)时的频率,在很多情况下,也常以fc作为通带或阻带截频。④当电路没有损耗时,固有频率:就是其谐振频率,复杂电路往往有多个固有频率。
2、增益与衰耗
滤波器在通带内的增益并非常数。①对低通滤波器通带增益,一般指ω=0时的增益;高通指ω→∞时的增益;带通则指中心频率处的增益。②对带阻滤波器,应给出阻带衰耗,衰耗定义为增益的倒数。③通带增益变化量指通带内各点增益的最大变化量,如果通带增益变化量以dB为单位,则指增益dB值的变化量。
3、阻尼系数与品质因数
阻尼系数α是表征滤波器对角频率为ω0信号的阻尼作用,是滤波器中表示能量衰耗的一项指标,它是与传递函数的极点实部大小相关的一项系数。
4、灵敏度
滤波电路由许多元件构成,每个元件参数值的变化都会影响滤波器的性能。
5、群时延函数
在滤波器设计中,常用群时延函数评价信号经滤波后相位失真程度。
以上的几个技术指标是一般滤波器的特性,但在实际应用中,数字滤波器通常用来实现选频操作,因此在利用DSP实现数字滤波器设计中要求的技术指标主要为在频域中给出的幅频响应和相频响应。如下图所示

幅频响应和相频响应特性曲线
对于幅频响应,它的含义是信号通过系统之后的输出信号的幅度与它输入时的信号的幅度的比值,一般以分贝值表示。对于相频响应,含义是信号通过系统之后的输出信号的相位与它输入时的信号的相位之差,在运用线性相频响应指标进行滤波器设计具有如下优点:①只包含实数算法,不涉及复数运算;②不存在延迟失真,只有固定数量的延迟;③可以采用FFT算法,从而提高运行效率;④由于FIR滤波器的单位脉冲响应是有限长序列,故FIR滤波器没有不稳定的问题,且误差较小。
基本原理
利用DSP实现FIR滤波器的设计方法主要有窗函数法和频率抽样法,其中窗函数法是基本的设计方法,这里采用窗函数法设计FIR滤波器。设希望得到的滤波器理想响应为 ,那么FIR滤波器的设计就在于寻找一个传递函数

去逼进 ,设

这里 就是傅立叶级数的系数。在这种逼近中,最直接的一种方法就是从单位脉冲响应 入手,使 逼近理想的单位脉冲响应 。由于 是一个无限长序列,因此,最简单的方法就是对 做截尾处理,即得到一个近似的传递函数

上式中,Q就是最终确定FIR滤波器的阶数,Q越大,近似程度就越高。对 截尾,实际上就是对 乘上一个矩形窗口 ,即

令z= ,则

其脉冲响应系数为 , ,…, , , ,…, , 。为使 具有因果性,延时Q个样值,可得:

令n+Q=k,上式成为

令 ,N=2Q,得

式中, 是脉冲响应系数,这里 …, ,…, 。
一般来说,FIR数字滤波器输出 的Z变换形式 与输入 的Z变换形式之间的关系如下:

实现结构如下图所示:

Z变换结构图
从上面的Z变换和结构图可以很容易得出FIR滤波器的差分方程表示形式,即对上式进行反Z变换得:

上式为FIR数字滤波器的时域表示方法,其中x(n)是在时间n的滤波器的输入抽样值,根据上式即可对滤波器进行设计。
硬件设计
1、DSP芯片
根据设计原理,实现的核心器件采用美国德州仪器公司生产的低功耗定点数字信号处理器芯片TMS320C5402。选择该芯片主要是因为它是目前最常用的低成本DSP芯片,而且包括以下主要特点:
⑴运算速度快,最快可达532MIPS;
⑵多总线结构,片内共有8 条总线(1条程序存储器总线、3条数据存储总线和4条地址总线);
⑶CPU采用冯? 诺依曼并行结构设计,使其能在一条指令周期内,高速地完成多项算术运算;
⑷片内集成了4K×16bitROM和16K×16bit的双存取RAM;
⑸丰富的片上外围电路(通用I/O 引脚,定时器,时钟发生器, HPI 接口,多通道缓冲串行口McBSP)使其与外部接口方便;
⑹3.3V I/O电压,1.8V核点压,工作电流平均值为75mA,其中核45mA,I/O约30mA;
⑺144脚BGA封装,使体积减少,功耗降低。
2、AD和DA电路
在本数字滤波器系统中选择了TI公司的TLV1570芯片作为模数转换器件,8通道10位2.7到5.5 V低电压模数转换芯片。TLVl570在3V电压下的采样频率为625KSPS,输入信号最高频率不能超过300K。
由于模数转换选择了10位器件,为了简化程序代码,减少DSP 的运算工作量,在本数字滤波器系统中选择了TI公司的TLV5608芯片,它是一款8通道10位2.7到5.5V低电压数模转换芯片。
3、电源电路
根据DSP芯片工作的电压电流需求,及芯片采用双电源供电对加电顺序的要求,考虑使用TI公司的电源转换芯片TPS73HD318,其输出电压为一路3.3V、一路1.8V,每路电源的最大输出电流为750mA,能满足本系统的供电需求。而且TPS73xx具有非常低的静态电流,能使稳压器输出稳定。
4、时钟电路
C54xx系列的时钟端子为X1和X2/CLKIN,采用无源晶振提供时钟信号,由于DSP有一组端子可以用来调整其工作频率的高低,故对晶振频率大小的选定没有特别的要求,这里选用10Mhz的晶振。
5、复位电路
为了克服DSP系统因时钟频率较高导致在运行时可能发生的干扰和被干扰的现象,最好是使用具有监视(Watchdog)功能的自动复位电路,于是采用专门的自动复位芯片MAX706。MAX706的电源为3.1V~5.0V,低电平复位输出,复位门限为3.08V。
6、未用端子处理
根据使用DSP芯片的相关原则,以及芯片手册具体决定未用端子是接上拉电阻还是悬空。
7、基于上述的各部分电路组成,可以得出DSP数字滤波器的整体硬件电路连线图,如下所示

程序设计
1、设计思路
⑴在DSP进行数字滤波运算前首先要进行初始化,只有正确设置了DSP的初始状态才能保证芯片能正常运行。本系统主要进行以下两方面的初始化:
①寄存器初始化:状态寄存器ST0、状态寄存器ST1、处理器模式控制寄存器PMST、软件等待状态寄存器SWWSR、组交换控制寄存器BSCR和时钟模式寄存器等。
②中断矢量表初始化:根据DSP芯片对各中断矢量的设置位置编写一个子程序;设置PMST控制寄存器;连接时将矢量表重定位到IPTR指定的地址。
⑵其次就是FIR 数字滤波的子程序设计,主要步骤如下:
①查询SPCR11寄存器的第二位,当为1时说明read ready,将DRR11的值读入AR3所指向的地址,该值为最新的采样值。
②将最新的采样值减去200h,然后AR3的值减1。
③执行MAC指令。
④将累加器的值送给变量Y,并将Y加上200h。
⑤查询SPCR20寄存器的第二位,当为1时说明writeready,将Y值赋给DXR10,该值为滤波器输出值。
⑥循环执行上面步骤。
2、程序流程图
依据上述程序设计思路可以得到利用DSP实现FIR滤波器设计的程序流程图,如下

3、程序代码
由于初始化程序部分过于庞大繁杂,这里只给出用MAC指令编程实现FIR低通滤波器的程序片断:
FILT_task1
LD Store_SICX,A
STLM A,ar4
STM #1,ar0 ;间址
STM #28,bk
LD DEM_Out,A
STL A,*ar4+% ;输入信号:实部
STM #Coef_Tab1,ar5 ;滤波器实部系数地址
LD #0,A
STM #27,brc
RPTB SICXU-1
MAC *AR4+0%,*AR5+,A
SICXU LD A,-16,A ;低通滤波结果
LD C7FFF,B
MIN A
NEG B
MAX A
STL A,DEM_Out
LDM AR4,A
STL A,Store_SICX
RET
Coef_Tab1
.word 100 ;h(0)
.word 7 ;h(1)

… ;脉冲响应系数

.word -248
.word -71 ;h(N-1)
.end
总结
通过利用DSP的FIR滤波器设计,对DSP芯片的使用,以及利用DSP芯片组成的基本系统的相关电路有了比较深的认识。熟悉DSP芯片的系统设计和应用开发流程,并利用图书馆、网络、询问同学等方式查找资料和解决相关的难题,这是最基础的工作,也是最关键的步骤。这样做可以培养自己的动手解决问题的能力和独立思考的处事方法,使自己具有技术人员的气质和工作态度,为将来的就业增加优势。
数字滤波器是DSP的典型应用,学会了有助于触类旁通,利于进一步的学习研究,能做到理解其他基于DSP的系统的功能和工作原理。掌握了基于DSP的应用开发,开阔了视野,增长了知识,是进入现代数字信号处理领域重要技能,乃至大规模集成电路的开发也是会用到的基础,今后要予以重视并积极努力去学习。
求采纳为满意回答。追问

你把原网页地址发给我看下吧

FIR滤波器的设计方法有哪些?
FPGA有着规整的内部逻辑阵列和丰富的连线资源,特别适合于数字信号处理任务。相对于串行运算为主导的通用DSP芯片来说,其并行性和可扩展性更好,利用FPGA乘累加的快速算法,可以设计出高速的FIR数字滤波器。

基于dsp的 f.i.r低通滤波器设计
利用DSP实现FIR滤波器的设计方法主要有窗函数法和频率抽样法,其中窗函数法是基本的设计方法,这里采用窗函数法设计FIR滤波器。设希望得到的滤波器理想响应为 ,那么FIR滤波器的设计就在于寻找一个传递函数 去逼进 ,设这里 就是傅立叶级数的系数。在这种逼近中,最直接的一种方法就是从单位脉冲响应 入手,使 逼近理想...

关于数字信号处理中FIR滤波器的设计的问题
DSP设计主要讲究的是在线算法,用时域的方式解决频域的问题。这个框图本身的效果并不是为了实现一个DSP的设计,而是为了解决两个离线的数据如何快速得到其卷积。例如多项式f(x),g(x)求乘法,平常需要O(n^2)的时间复杂度,n为多项式项数。而采用此结构(其实多项式乘法也是卷积),则需要:FFT算法O(nl...

fir数字滤波器设计原理是什么
滤波器输出的数据都是一串序列,要使它能直观地反应出来,还需经过数模转换,因此由FPGA构成的FIR滤波器的输出须外接D\/A模块。FPGA有着规整的内部逻辑阵列和丰富的连线资源,特别适合于数字信号处理任务,相对于串行运算为主导的通用DSP芯片来说,其并行性和可扩展性更好,利用FPGA乘累加的快速算法,可以...

FIR滤波器工作原理
FPGA,即现场可编程门阵列,是一种理想的FIR滤波器实现平台。它拥有整齐的内部逻辑阵列和丰富的连接资源,特别适合于数字信号处理任务。相比以串行运算为主的通用DSP芯片,FPGA在并行性和可扩展性上具有显著优势。利用FPGA的快速乘累加算法,能够设计出运行速度极快的FIR数字滤波器,满足高效处理的需求。

急!!!用窗函数法设计FIR滤波器的主要特点是什么?
没有输出到输入的反馈,但有些结构中(例如频率抽样结构)也包含有反馈的递归部分。设FIR滤波器的单位冲激响应h (n)为一个N点序列,0 ≤ n ≤N —1,则滤波器的系统函数为H(z)=∑h(n)*z^-k。就是说,它有(N—1)阶极点在z = 0处,有(N—1)个零点位于有限z平面的任何位置。

脉冲响应函数与频率响应函数的关系怎样?
DSP需要更多的计算时间,对DSP的实时性有影响。以下都是低通滤波器的设计。FIR的设计:FIR滤波器的设计比较简单,就是要设计一个数字滤波器去逼近一个理想的低通滤波器。通常这个理想的低通滤波器在频域上是一个矩形窗。根据傅里叶变换我们可以知道,此函数在时域上是一个采样函数。通常此函数的表达式为...

dsp中的FFT是如何实现信号的倒序输入
DSP处理器与通用处理器的比较 考虑一个数字信号处理的实例,比如有限冲击响应滤波器(FIR)。用数学语言来说,FIR滤波器是做一系列的点积。取一个输入量和一个序数向量,在系数和输入样本的滑动窗口间作乘法,然后将所有的乘积加起来,形成一个输出样本。类似的运算在数字信号处理过程中大量地重复发生,...

MATLAB下利用FDATool进行FIR滤波器设计
Matlab的FDATool在FIR滤波器设计中展现出强大功能,能便捷生成滤波器系数,适用于嵌入式、FPGA及DSP系统中的FIR滤波器应用。以下步骤演示在Matlab环境下使用FDATool进行FIR滤波器设计的过程。启动FDATool可直接在命令窗口输入“fdatool”或通过界面打开。设计流程如下:选择一个低通滤波器案例,配置参数如下:采样...

哪些是fir数字滤波器的基本网络结构类型
要使它能直观地反应出来。还需经过数模转换,因此由FPGA构成的FIR滤波器的输出须外接D\/A模块。FPGA有着规整的内部逻辑阵列和丰富的连线资源,特别适合于数字信号处理任务,相对于串行运算为主导的通用DSP芯片来说,其并行性和可扩展性更好,利用好FPGA乘累加的快速算法,可以设计出高速的FIR数字滤波器。

相似回答