FIR有限冲击滤波系统设计

合集下载

fir滤波器的设计方法

fir滤波器的设计方法

fir滤波器的设计方法一、引言二、基本概念1.数字信号2.离散时间信号3.FIR滤波器三、FIR滤波器的设计方法1.窗函数法(1)矩形窗函数法(2)汉宁窗函数法(3)汉明窗函数法(4)布莱克曼窗函数法2.最小二乘法3.频率抽样法四、FIR滤波器设计实例五、总结一、引言数字信号处理在现代通信技术中得到了广泛的应用,其中滤波器是数字信号处理的重要组成部分。

FIR滤波器是一种常用的数字滤波器,具有无限冲击响应和线性相位特性。

本文将介绍FIR滤波器的基本概念和设计方法,并给出一个实例。

二、基本概念1.数字信号数字信号是在时间轴上取样后离散化的模拟信号。

在计算机中,数字信号由一系列离散的数值表示。

2.离散时间信号离散时间信号是以时间为自变量且取值为离散值的函数。

通常使用序列表示,如x(n)。

3.FIR滤波器FIR滤波器是一种数字滤波器,其系统函数是有限长冲击响应的线性时不变系统。

FIR滤波器的输出只与当前和过去的输入有关,与未来的输入无关。

FIR滤波器具有无限冲击响应和线性相位特性。

三、FIR滤波器的设计方法1.窗函数法窗函数法是一种常用的FIR滤波器设计方法。

它通过在频域上对理想低通滤波器进行截止频率处理得到所需的频率响应,并使用窗函数将其转换为时域上的序列。

(1)矩形窗函数法矩形窗函数法是最简单的FIR滤波器设计方法。

它将理想低通滤波器在频域上乘以一个矩形窗函数,得到所需频率响应后再进行反变换得到时域上的系数序列。

(2)汉宁窗函数法汉宁窗函数法是一种常用的FIR滤波器设计方法。

它将理想低通滤波器在频域上乘以一个汉宁窗函数,得到所需频率响应后再进行反变换得到时域上的系数序列。

(3)汉明窗函数法汉明窗函数法是一种常用的FIR滤波器设计方法。

它将理想低通滤波器在频域上乘以一个汉明窗函数,得到所需频率响应后再进行反变换得到时域上的系数序列。

(4)布莱克曼窗函数法布莱克曼窗函数法是一种常用的FIR滤波器设计方法。

FIR滤波器的设计及特点

FIR滤波器的设计及特点

FIR滤波器的设计及特点FIR滤波器(Finite Impulse Response Filter)是一种数字滤波器,它的特点是其冲激响应是有限长度的。

FIR滤波器通过对输入序列做线性加权的运算来实现滤波的效果。

FIR滤波器的设计需要确定滤波器的系数以及长度,其设计方法有很多种,其中比较常用的有窗函数法、频率采样法以及最小二乘法。

FIR滤波器的设计方法之一是窗函数法,它是根据所设定的频率响应曲线来进行设计的。

具体的步骤是:首先,在频率域上设定所需的频率响应曲线;然后,将该曲线转换到时域上,得到滤波器的单位冲激响应;最后,对单位冲激响应进行加窗处理,得到最终的滤波器系数。

在窗函数法中,常用的窗函数有矩形窗、汉宁窗、哈宁窗等,不同的窗函数会导致滤波器具有不同的性能,如频域主瓣宽度、滤波器的过渡带宽度等。

另一种常用的FIR滤波器设计方法是频率采样法,它是通过在频率域上进行采样来确定滤波器的系数。

在频域上,滤波器的频率响应可以表示为幅度特性和相位特性。

通过选取一组频率,在这些频率上等幅响应,并且在其余的频率上衰减至零,然后对这些采样点进行IFFT运算,即可得到滤波器的系数。

频率采样法的特点是可以直观地设计滤波器,但是在采样点之间的频率响应无法得到保证,会产生幅度插值误差。

最小二乘法是一种通过最小二乘准则来设计滤波器的方法。

它在时域上通过对输入序列和输出序列之间的误差进行最小化,得到最优的滤波器系数。

最小二乘法可以看作是一种优化问题的求解方法,需要解决一个线性规划问题,因此需要求解线性方程组来确定滤波器的系数。

1.稳定性:FIR滤波器是一种无反馈结构的滤波器,其零点可以完全控制在单位圆内,因此具有稳定性保证。

2.线性相位特性:FIR滤波器的冲激响应通常是对称的,因此它不会引入相位失真,可以保持输入信号的相位。

3.精确控制频率响应:FIR滤波器的频率响应可以通过设计滤波器系数来精确控制,具有很高的灵活性。

4.零相移滤波:由于线性相位特性,FIR滤波器可以实现零相移的滤波效果,适用于对输入信号相位要求较高的应用。

数字信号处理:第7章下 有限冲激响应(FIR) 滤波器的设计方法

数字信号处理:第7章下 有限冲激响应(FIR) 滤波器的设计方法

N sin k
N
取N=64;
当L=1时,sinc函数值恒 为1,说明该滤波器对所 有频率的信号的增益为1, 即不衰减也不放大;
显然,L=1时离散矩形函 数退化为单位脉冲函数, 而单位脉冲函数的频率 响应为常数1;
另外,单位脉冲函数使 得输出即为输入,信号 不但幅度不改变,也不 会发生延迟。
不同L时sinc函数(N=64)
幅频响应的最常见方式是 显示振幅的绝对值,因此, 在振幅发生变号的地方, 其对应的相频响应发生了π 的奇数倍的相移。
R E CT(k )
A(k )e-j (k )
sin kL N
sin k
jk ( L1)
eN
N
RECT(k ) e-j(k)
因此,一个因果的线性相 位系统的相频响应φ(k)可 能不是严格线性的,而是 分段线性的,在断点处一 般有π的奇数倍的跳变。
n0
n0
令 2k q 有: N
RECT(k)
L1
e- jqn
n0
1 e-jqL 1- e-jq
e e e jqL / 2 jqL / 2 jqL / 2 e jq / 2 e jq / 2 e jq / 2
jq( L1)
e2
sin qL 2
sin q
2
R E CT(k )
jk ( L1)
频域Gibbs现象依然存在
Gibbs现象
1898年,波兰物理学家Albert Michelson发明了一个能够计算傅 里叶级数之分析及综合方程的机器 。当输入波形为方波时,综合出的 图形在方波跳变处发生振荡。 Michelson以为这是一个bug,修 复机器后便能消除这个现象。
1899年,美国物理学家J. Willard Gibbs指出,这些振荡是一种数学 现象,并且将在对非连续函数应 用傅里叶级数时必然产生。

实验四FIR数字滤波器的设计

实验四FIR数字滤波器的设计

实验四FIR数字滤波器的设计
FIR(有限冲击响应)数字滤波器是一种常见的数字信号处理器件,
可以用于滤波、降噪等应用。

下面是一种FIR数字滤波器的设计流程:
1.确定滤波器的需求:首先确定需要滤除的频率范围和滤波的类型,
例如低通、高通、带通、带阻等等。

2.设计滤波器的频率响应:根据滤波器的需求,设计其理想的频率响应。

可以使用窗函数、最小二乘法等方法获得一个理想的滤波器响应。

3.确定滤波器的阶数:根据设计的频率响应,确定滤波器的阶数。


数越高,滤波器的响应越陡峭,但计算复杂度也会增加。

4.确定滤波器的系数:根据滤波器的阶数和频率响应,计算滤波器的
系数。

可以使用频域窗函数或时域设计方法。

5.实现滤波器:根据计算得到的滤波器系数,实现滤波器的计算算法。

可以使用直接形式、级联形式、传输函数形式等。

6.评估滤波器的性能:使用所设计的FIR滤波器对输入信号进行滤波,评估其滤波效果。

可以使用频率响应曲线、幅频响应、群延时等指标进行
评估。

7.调整滤波器设计:根据实际的滤波效果,如果不满足需求,可以调
整滤波器的频率响应和阶数,重新计算滤波器系数,重新实现滤波器。

以上是FIR数字滤波器的基本设计流程,设计过程中需要考虑滤波器
的性能、计算复杂度、实际应用需求等因素。

数字信号处理实验FIR数字滤波器的设计

数字信号处理实验FIR数字滤波器的设计

数字信号处理实验:FIR数字滤波器的设计1. 引言数字滤波器是数字信号处理的关键技术之一,用于对数字信号进行滤波、降噪、调频等操作。

FIR (Finite Impulse Response) 数字滤波器是一种常见的数字滤波器,具有线性相应和有限的脉冲响应特性。

本实验旨在通过设计一个FIR数字滤波器来了解其基本原理和设计过程。

2. FIR数字滤波器的基本原理FIR数字滤波器通过对输入信号的每一个样本值与滤波器的冲激响应(滤波器的系数)进行线性加权累加,来实现对信号的滤波。

其数学表达式可以表示为:y(n) = b0 * x(n) + b1 * x(n-1) + b2 * x(n-2) + ... + bN * x(n-N)其中,y(n)表示滤波器的输出,x(n)表示滤波器的输入信号,b0~bN表示滤波器的系数。

FIR数字滤波器的脉冲响应为有限长度的序列,故称为有限冲激响应滤波器。

3. FIR数字滤波器的设计步骤FIR数字滤波器的设计主要包括以下几个步骤:步骤1: 确定滤波器的阶数和截止频率滤波器的阶数决定了滤波器的复杂度和性能,而截止频率决定了滤波器的通带和阻带特性。

根据实际需求,确定滤波器的阶数和截止频率。

步骤2: 选择滤波器的窗函数窗函数是FIR滤波器设计中常用的一种方法,可以通过选择不同的窗函数来实现不同的滤波器特性。

常用的窗函数有矩形窗、汉宁窗、汉明窗等。

根据实际需求,选择合适的窗函数。

步骤3: 计算滤波器的系数根据选择的窗函数和滤波器的阶数,使用相应的公式或算法计算滤波器的系数。

常见的计算方法有频率采样法、窗函数法、最小二乘法等。

步骤4: 实现滤波器根据计算得到的滤波器系数,可以使用编程语言或专用软件来实现滤波器。

步骤5: 评估滤波器性能通过输入测试信号,观察滤波器的输出结果,评估滤波器的性能和滤波效果。

常见评估指标有滤波器的幅频响应、相频响应、群延迟等。

4. 实验步骤本实验将以Matlab软件为例,演示FIR数字滤波器的设计步骤。

信号与系统复习资料第7章有限单位冲激响应FIR数字滤波器的设计方法

信号与系统复习资料第7章有限单位冲激响应FIR数字滤波器的设计方法

x1 y1 A[a sin(0 t 1 )] x2 y2 A[a sin(20 t 2 )] y(t ) y1 (t ) y2 (t )
不产生相移失真的条件:x1, x2 经过系统后产生的延 迟应相等。 ( ) 0 0 1 0 1 由 d ( ) 2 2 0 2 0 0 2
h ( N 1 n) z
n 0 ( N 1 m )
N 1
n
H(z-1)
m0
h ( m) z
m0
N 1
z
1
( N 1)
h ( m) z
N 1
m
即:
H ( z) z
( N 1)
H (z )
-13-
FIR 滤波器
H (z) z ( N 1) H (z 1 )
对低通、高通、带通及带阻均适用。
-24-
线性相位滤波器的幅度特点 2. h(n)偶对称,N为偶数
若N为奇数,则
N 1 由于N为偶数 h( ) 不存在,因此 2
前N/2项之 和的2倍
FIR 滤波器
N 1 ( N 3)/ 2 N 1 H ( ) h( ) 2h(n) cos ( n) 2 2 n 0
( N 1) h(n) cos (n ) 2 n 0 cos x 为耦函数
N 1
H (e ) e
显然
j
j
( N 1) 2
e j ( )
( N 1) h(n) cos ( n) 2 n 0
N 1 H ( )
H ( )

有限冲击响应FIR数字滤波器设计zhjConcise资料重点

有限冲击响应FIR数字滤波器设计zhjConcise资料重点

2)h ] sin[(n M (n M 2)
2)l ]
相当于用一个截止频率在h 处的低通滤波器 减去一个截止频率在 l 处的低通滤波器。
带阻:令:
H d (e j )
e jM 0
2
l, h
其它
hd
(n)
sin[(n
M 2
)l
]
sin[(n M 2
(n M
)
2)
]
sin[(n
则:
hd (n)
1
2
1
2
H d (e j )en(cn) n
特点: 无限长 非因果 偶对称
解决方法: 截短, 移位
保留
即: h(n) hd (n M 2)
, n 0,1,..., M
隐含着使用 了窗函数
M
于是: H ( z) h(n)z n
n0
注意:H (z) 是因果的,且是线性相位的,即
2
c 0 c
hd (n) sin[(n M
2) ] sin[(n M (n M 2)
2) c ]
相当于用一个截止频率在 处的低通滤波器
(实际上是全通滤波器)减去一个截止频率
在 c 处的低通滤波器。
带通:
令:
H d (e j )
e jM 0
2
l h
其它
hd (n) sin[(n M
二、 FIR DF 设计的窗函数法的特点:
优点:1. 无稳定性问题; 2. 容易做到线性相位; 3. 可以设计各种特殊类型的滤波器; 4. 方法特别简单。
缺点:1. 不易控制边缘频率; 2. 幅频性能不理想; 3. h(n) 较长;
改进:1. 使用其它类型的窗函数; 2. 改进设计方法。

有限冲击响应滤波器 (FIR)

有限冲击响应滤波器 (FIR)

§4.5 [实验4.5] 有限冲击响应滤波器 (FIR) 实现一、实验目的1. 掌握FIR 滤波器的原理和窗函数设计法;2. 掌握用C 语言编写DSP 程序的方法。

二、实验设备1. 一台装有CCS 软件的计算机;2. DSP 实验箱的TMS320C5416主控板;3. DSP 硬件仿真器。

三、实验原理数字滤波是DSP 的最基本的应用领域之一。

对于许多应用来说,数字滤波一般具有如下的差分方程形式:()()()∑-+∑-=-=-=110N k k N k k k n y b k n x a n y式中,X(n) 为输入序列,Y(n)为输出序列,A k 和B k 为滤波器系数,N 是滤波器的阶数。

若式中所有的B k 均为零,且通常把系数A k 记为h k , 则有:()()∑-=-=1N k k k n x h n y上式就是FIR 滤波器的差分方程了。

FIR 滤波器的最主要的特点是没有反馈回路,因此它是无条件稳定系统。

它的单位脉冲响应h(n)是一个有限长序列。

由上面的方程可见,FIR 滤波算法实际上是一种乘法累加运算,它不断地输入样本x(n),经延时 (z –1),做乘法累加,再输出滤波结果 y(n)。

要设计一个FIR 滤波器就是要求出它的冲击响应系数h(n),设计方法主要有窗函数法和频率抽样法,本实验要求掌握窗函数法,这也是最基本的方法。

理想的低通滤波器的频率响应Hd (w)是一个矩形,这意味着它在时域上是无限长的序列,这在实际上是不可能实现的。

因此我们要采取某种方法截断 Hd(n),可以用一个有限长度的窗函数序列w(n)与之相乘。

这个窗函数序列的形状和长度都会对最后系统的频率响应特性产生影响,因此对窗函数的分析和选择是设计FIR 滤波器的关键问题所在。

本实验举了五种常用的窗函数为例,通过设置参数可以得到加上不同窗后的冲击响应序列h(n),并且可以观察到其幅频响应图。

关于根据给定频率要求进行FIR 滤波器设计的详细原理,以及在求得符合要求的h(n)后如何对输入信号序列进行滤波,请读者参考数字信号处理的有关资料。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

0 《现代信号处理课程设计》 课程设计报告

设计题目: FIR有限冲击滤波系统设计 1

目 录 第1章 绪论..............................................错误!未定义书签。 1.1设计背景.............................................错误!未定义书签。 1.2设计内容及目的.... ................................错误!未定义书签。 1.3设计思路简介.........................................2 第2章 系统开发平台与环境.................................2 1.1 CCS开发环境.........................................错误!未定义书签。 2.2 ICETEK-F2812 –A 评估板..............................错误!未定义书签。 第3章 FIR滤波器设计过程 ...............................4 3.1 FIR滤波器设计总框图 .................................4 3.2 FIR滤波器设计的原理及其方法..........................5 第4章 系统软件设计.......................................8 4.1 程序流程图...........................................8 4.2程序源代码 ...........................................9 第5章 系统仿真..........................................13 2

5.1仿真设置 ............................................13 5.2启动CCS ...........................................19 5.3 仿真图..............................................19 第6章 总结..............................................错误!未定义书签。 参考文献.................................................20 1

第1章 绪论 1.1设计背景 从广义上讲,任何对某些频率(相对于其他频率来说)进行修正的系统称为滤波器。严格地讲,对输入信号通过一定的处理得到输出信号,这个处理通常是提取信号中某频率范围内的信号成分,把这种处理的过程称为滤波。实现滤波处理的运算电路或设备称为滤波器。在许多科学技术领域中,广泛应用线性滤波和频谱分析对信号进行加工处理,模拟滤波是处理连续信号,数字滤波则是处理离散信号,而后者是在前者的基础上发展起来的。我们知道,无源或有源模拟滤波器是分立元件构成的线性网络,他们的性能可以用线性微分方程来描述,而数字滤波器是个离散线性系统,要用差分方程来描述,并以离散变换方法来分析。这些方程组可以用专用的或通用的数字计算机进行数字运算来实现。因此,数字滤波器的滤波过程是一个计算过程,它将输入信号的序列数字按照预定的要求转换成输出数列。

数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类:无限冲激响应(IIR)滤波器和有限冲激响应(FIR)滤波器。由于FIR滤波器的最大优点是可以实现线性相位滤波,信号无失真与滤波处理的条件是在信号的有效频谱范围内系统幅频响应为常数,相频响应为频率的线性函数(即具有线性相位)。在数字通信和图像传输处理等应用场合都要求滤波器具有线性相位特性。另外FIR滤波器是全零点滤波器,硬件和软件实现结构简单,不用考虑稳定性问题,在数字信号处理领域得到广泛应用。FIR滤波器的设计方法有许多种,如窗函数设计法、频率采样设计法和最优化设计法等。本次设计主要利用窗函数法设计FIR滤波器,并在DSP上实现滤波。 1.2设计内容及目的 设计内容主要是在给定滤波器的性能指标下,采用窗函数法和双线性变换设计滤波器,并画出滤波器的频率响应;然后用自己设计的滤波器对采集的信号进行滤波,画出滤波后信号的时域波形和频谱,并对滤波前后的信号进行对比,分析信号的变化。设计的目的为: 2

1、 掌握用窗函数法设计FIR数字滤波器的原理和方法 2、 熟悉线性相位FIR数字滤波器特性 3、 了解各种窗函数对滤波器特性的影响 1.3设计思路简介

在TMS320C54x系统开发环境CCS(Code Composer Studio)下对FIR滤波器的DSP实现原理进行讨论。利用C语言设计相应的滤波器,通过

实验仿真,从输入信号和输出信号的时域和频域曲线可看出在DSP上实现的FIR滤波器能完成预定的滤波任务。

第2章 系统开发平台与环境

1.1 CCS开发环境 CCS提供了配置、建立、调试、跟踪和分析程序的工具,它便于实时、嵌

入式信号处理程序的编制和测试,它能够加速开发进程,提高工作效率。

CCS提供了基本的代码生成工具,它们具有一系列的调试、分析能力。CCS支持如下图1.1所示的开发周期的所有阶段。

图2.1 CSS开发周期阶段 2.2 ICETEK-F2812 –A 评估板 ICETEK–F2812-A 评估板的内存映射 ICETEK–F2812-A 评估板(以下简称评估板)的内存映射与TMS320F2812A 芯片不同的是,所有的寄存器和存储器全部映射在XINTF ZONE 2 译码的空间内,占TMS320F2812A 芯片的80000h-0FFFFFh 地址单元中。因此,如果您需要使用DSP 芯片的管脚外扩硬件设备的话,请避免使用这部分地址资源。对于这部分地址单

设计概念性规划

编程和编译创建工程文件、编写源代码和配置文件调试语法检查、探测点设置和日志保存等分析

实时调试、统计和跟踪 3

元,基本上分为两个部分,从 80000h 到0BFFFFh 共256K 的地址单元分配给片外扩展存储器。这里有一个声明,一般来说,在出厂时,板上所连接的存储器一般为64K 的存储器,因此这部分的实际有效内存是80000h-8FFFFh。为了特殊的需求,这里最多可以扩展到256K 的实际存储器。从 C0000h 到FFFFFh 是外围其他设备的寄存器接口,这里包含一个12 位的模数转换器DAC7528,以及四位的状态显示数码管和四位可读入数字量的开关。ICETEK–F2812-A 评估板实物图和原理框图分别如图1.2、1.3所示。

图2.2 ICETEK–F2812-A 评估板实物图

图2.3 ICETEK–F2812-A 评估板原理框图 4

第3章 FIR滤波器设计过程 3.1 FIR滤波器设计总框图 FIR滤波器设计总框图如图3.1所示

开始

初始化DSP产生输入数据

读入新数据并存放到栈顶

滤波处理保存滤波结果到输出单元

待滤波的样本点数处理完成

结束否

图3.1 FIR滤波器设计总框图 5

3.2 FIR滤波器设计的原理及其方法 FIR数字滤波器设计原理是利用数字的方法,按预定的要求对信号进行变换,把输入的信号变成一定的输出信号,从而达到改变信号频谱的目的。 FIR滤波器(有限冲激响应滤波器)的差分方程是:

10)()()(Nkknxkhny

在数字滤波中,我们主要讨论离散时间序列。如图3.2所示。设输入序列为nx,离散或数字滤波器对单位抽样序列n的响应为nh。因n在时域离散信号和系统中所起的作用相当于单位冲激函数在时域连续信号和系统中所起的作用。

图3.2 数字滤波器原理 数字滤波器的序列ny将是这两个序列的离散卷积,即



kknxkhny (1.4)

同样,两个序列卷积的z变换等于个自z变换的乘积,即 zXzHzY (1.5)

用Tjez代入上式,其中T为抽样周期,则得到 TjTjTjeXeHeY (1.6)

式中TjeX和 TjeY 分别为数字滤波器输入序列和输出序列的频谱,而TjeH为单位抽样序列响应nh的频谱。由此可见,输入序列的频谱TjeX

过滤波后,变为TjTjeXeH ,按照TjeX的特点和我们处理信号的目的,选取适当的TjeH使的滤波后的TjTjeXeH符合我们的要求。 其直接型结构如图3.3所示:

数字滤波器nh,H(z)

zx

nx



ny

zy

相关文档
最新文档