dsp数字滤波器c语言窗函数,数字滤波器的设计(DSP课程设计)资料.doc

大软

于 2021-05-21 23:42:05 发布

阅读量645

点赞数

文章标签： dsp数字滤波器c语言窗函数

数字滤波器的设计(DSP课程设计)资料.doc

摘要

在数字信号处理中,数字滤波器是一种被广泛使用的信号处理部件。针对有限长冲激响应(FIR)数字滤波器的原理,讨论了窗函数法设计线性相位FIR数字滤波器的基本思路,介绍了用MATLAool)来设计滤波器的方法及在DSP上的实现。开发环境——实现了达到目标要求的滤波器设计和仿真。

关键词：数字滤波器，

1.设计内容

1.1设计背景

数字滤波是数字信号处理的基本方法。数字滤波与模拟滤波相比有很多优点，它除了可避免模拟滤波器固有的电压漂移、温度漂移和噪声等问题外，还能满足滤波器对幅度和相位的严格要求。DSP(数字信号处理器)与一般的微处理器相比有很大的区别，它所特有的系统结构、指令集合、数据流程方式为解决复杂的数字信号处理问题提供了便利，本文选用TMS320C5作为DSP处理芯片，通过对其编程来实现滤波器。对数字滤波器而言,从实现方法上,有限长冲激响应(FIR滤波器和无限冲激响应(IIR)滤波器之分。由于FIR滤波器只有零点,因此这一类系统不像IIR系统那样易取得比较好的通带与阻带衰减特性。但是FIR系统有自己突出的优点:系统总是稳定的;易实现线性相位;允许设计多通带(阻带)滤波器。其中后两项是IIR系统不易实现的。

已知x1(n)=sin(2*pi*f11*n*T);x2(n)=0.7*sin(2*pi*f12*n*T)

x3(n)=0.5*sin(2*pi*f13*n*T);

x(n) = x1(n)+ x2(n)+ x3(n)

其中：f11=500Hz; f12=2000Hz; f13=4000Hz; fs=10000Hz

要求设计一个基于DSP的FIR高通滤波器，把f11和f12滤掉，保留f13。

FIR高通滤波器的系数由Matlab中的fdatool产生，信号也由Matlab产生。

1.2．2 技术指标：

根据设计要求可以确定如下技术指标：

(1)阻带衰减αp=40dB

(2)采样频率Fs=10000Hz

(3)过渡带宽度1000Hz

1.3设计思路

通过编写MATLAB程序得到信号x(n)，并导出数据文件input.dat作为信号处理的输入信号。根据设计任务及技术指标利用MATLAB中的Fdatool工具设计滤波器的系数。在ccs中创建滤波器工程，根据上述得到的滤波器系数及输入信号编写相应的c语言源程序进行编译仿真，得到滤波结果。

1.4相关介绍

1.4.1 MATLAB简介

MATLAB是“矩阵实验室”(MATrix LABoratoy)的缩写，是一种科学计算软件，主要适用于矩阵运算及控制和信息处理领域的分析设计，它使用方便，输入简捷，运算高效，内容丰富，因此很多专家在自己擅长的领域用它编写了许多专门的MATLAB工具包。其中FDATool是一个功能强大的数字滤波器设计与分析工具，它涵盖了信号处理工具箱中所有的滤波器设计方法。利用它可以方便地设计出满足各种性能指标的滤波器，并可查看该滤波器的各种分析波形。

1.4.2 CCS简介

CCS是一种针对TMS320系列DSP的集成开发环境,在Windows操作系统下，采用图形接口界面，提供有环境配置、源文件编辑、程序调试、跟踪和分析等工具。 CCS有两种工作模式，即软件仿真器模式：可以脱离DSP芯片，在PC机上模拟DSP的指令集和工作机制，主要用于前期算法实现和调试。硬件在线编程模式：可以实时运行在DSP芯片上,与硬件开发板相结合在线编程和调试应用程序。

2.FIR滤波器简介

设h ( n) ( n = 0， 1， 2 N - 1)为滤波器的冲激响应，输入信号为x ( n) ，则F IR 滤波器就是要实现下列差分方程:

式(1)就是FIR 滤波器的差分方程。F IR 滤波器的最主要的特点是没有反馈回路，因此它是无条件稳定系统。它的单位脉冲响应h ( n)是一个有限长序列。由上面的方程可见， F IR 滤波算法实际上是一种乘法累加运算，它不断地输入样本x ( n) ，经延时( Z 3 /1)做乘法累加，再输出滤波结果y(n)[1,2,3] 。对式(1)进行Z 变换，整理后可得FIR 滤波器的传递函数为:

由式(2)可以看出， FIR 滤波器的一般结构如图1 所示。

图2-1 FIR 滤波器的一般结构

FIR数字滤波器的设计方法主要有窗函数法和频率抽样设计法，其中窗函数法是基本而有效的设计方法。

3.信号产生程序的设计

3.1信号产生程序

%信号产生程序

%输入信号的三个频率分量

f11=500; %/Hz

f12=2000; %/Hz

f13=4000; %/Hz

fs=10000; %/采样Hz

N=1000 %数据个数