第六章IIR滤波器的设计方法

合集下载

第六章IIR滤波器的设计方法

第六章IIR滤波器的设计方法

j4
数字滤波器的频率响应:
H
(e
j
)
1
0.318e j 0.4177e j 0.01831e
j
2
冲激响应不变法的优缺点
优点: h(n)完全模仿模拟滤波器的单位抽样响应ha(t)
时域逼近良好 保持线性关系:=T
线性相位模拟滤波器转变为线性相位数字滤波器
缺点: 频率响应混迭
只适用于限带的低通、带通滤波器
通带、阻带与过渡带:信号允许通过的频带 为通带,完全不允许通过的频带为阻带,通 带与阻带之间为过渡带。
滚降与滚降率:滤波器幅频特性在过渡 带的衰减和衰减速度称为滚降与滚降率。
阻带衰减:输入信号在阻带的衰减量
带内平坦度:通带和阻带内的平坦程度
4、数字滤波器的设计步骤
数字滤波器的设计三个步骤:
后各频率成分的衰减情况
( j) 为相频特性:反映各频率成分通过滤波
器后在时间上的延时情况
理想滤波器不可实现,只能以实际滤波器逼近
通带: c
1 1 H (e j ) 1
阻带: st H (e j ) 2
过渡带: c st
c :通带截止频率
st :阻带截止频率
:通带容限
Ha(s) Ha(-s)的零极点分布
由 Ha ( j) 2 确定Ha (s)的方法
由幅度平方函数得象限对称的s平面函数 对比 Ha ( j) 和 Ha (s) ,确定增益常数 由零极点及增益常数,得Ha (s)
低通

π
高通
H (e j ) π
2π ω

π
带通
H (e j ) π
2π ω

π
带阻
H (e j ) π

IIR滤波器设计报告

IIR滤波器设计报告

IIR滤波器设计报告IIR(Infinite Impulse Response)滤波器是数字信号处理领域中常用的一种滤波器类型。

与FIR(Finite Impulse Response)滤波器相比,IIR滤波器具有更高的滤波效率和更窄的频带宽度。

本文将介绍IIR滤波器的设计原理、设计步骤以及在实际应用中的一些注意事项。

一、IIR滤波器的设计原理IIR滤波器的设计基于递归差分方程的实现方法。

其基本结构包括反馈回路和前馈路径。

具体而言,IIR滤波器的输出值是输入值和过去输出值的加权和。

这种反馈结构使得IIR滤波器具有无限冲击响应的特性,即滤波器的输出值受到过去输出值的影响。

二、IIR滤波器的设计步骤1.确定滤波器的类型:根据实际需求确定滤波器是低通、高通、带通还是带阻类型。

2.确定滤波器的阶数:滤波器的阶数决定了滤波器对信号的响应速度和滤波器的复杂程度。

一般而言,阶数越高,滤波器的响应速度越快,但也会增加计算的复杂度。

3.确定滤波器的截止频率:根据实际需求确定滤波器的截止频率,即滤波器开始起作用的频率。

4. 计算滤波器的系数:根据滤波器的类型、阶数和截止频率,使用滤波器设计软件或公式来计算出滤波器的系数。

常用的设计方法包括巴特沃斯(Butterworth)滤波器设计、切比雪夫(Chebyshev)滤波器设计和椭圆(Elliptic)滤波器设计等。

5.实现滤波器:将滤波器的系数应用到差分方程或差分方程的转移函数中,从而实现滤波器。

三、IIR滤波器的应用注意事项1.阶数选择:较低的阶数可以实现基本的滤波效果,但可能无法满足更高的要求。

较高的阶数可以实现更精确的滤波效果,但同时也会增加计算的复杂度。

在实际应用中,需根据具体要求和系统的计算能力来选择适当的阶数。

2.频率响应:不同类型的IIR滤波器具有不同的频率响应特性。

在设计和选择滤波器的时候,需要根据实际应用需求来确定适合的滤波器类型。

3.稳定性:IIR滤波器可能会存在稳定性问题,即滤波器的输出会发散或产生震荡。

DSP 课件 IIR滤波器设计

DSP 课件 IIR滤波器设计

幅度响应
Type 2 Chebyshev Filter 1 N=3 N=5 N=7
Magnitude
0.8 0.6 0.4 0.2 0 0 1 2

3
特点:1、通带内单调下降;2、阻带内具有等波纹
阶数估计:
H a ( j s )
2
1 2 1 2TN ( s / p )
1 1 2 2 1 1 cosh[N cosh ( s / p ) A
Butterworth Filter 1 N=2 N=4 N = 102源自3、滚降的陡峭度(参数N)
2 1,1 c时, 2 G ( 2 ) 10 log10 1 G (1 ) 6 NdB
2N
Magnitude
0.8 0.6 0.4 0.2 0 0 1 2
模拟低通滤波器设计
滤波器特性
理想低通滤波器:
T
H r ( j )
c
c

典型的模拟低通滤波器幅度响应指标
滤波器特性:
通带 0 p 中,
1 p H a ( j) 1 p , for p
阻带 s 中,
H a ( j) s , for s
归一化形式(c=1)
1
s N a1 s n 1 a 2 s n 2 a N 1 s a N
巴特沃斯滤波器的特性:
1、 =0处的最大平坦幅度特性(前2N-1阶导数为0) 2、-3dB截止频率(参数c )
G ( c ) 10 log10 H a ( j c ) 10 log10 (1 / 2) 3dB c : 3dB截止频率
椭圆滤波器的幅度响应

IIR低通滤波器设计

IIR低通滤波器设计

一.数字滤波器1.1 数字滤波器的概念滤波器是指用来对输入信号进行滤波的硬件和软件。

数字滤波器是对数字信号实现滤波的线性时不变系统。

数字滤波器可以理解为是一个计算程序或算法,将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列,并在转化过程中,使信号按预定的形式变化。

数字滤波实质上是一种运算过程,实现对信号的运算处理。

数字滤波器和模拟滤波器相比,因为信号的形式和实现滤波的方法不同,数字滤波器具有比模拟滤波器精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配等优点。

输入数字信号(数字序列)通过特定的运算转变为输出的数字序列,因此,数字滤波器本质上是一个完成特定运算的数字计算过程,也可以理解为是一台计算机。

描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则,使其按照这个规则完成对输入数据的处理。

时域离散系统的频域特性:)()()(ωωωj j j e H e X e Y = (1-1)其中)(ωj e Y 、)(ωj e X 分别是数字滤波器的输出序列和输入序列的频域特性(或称为频谱特性), )(ωj e H 是数字滤波器的单位取样响应的频谱,又称为数字滤波器的频域响应。

输入序列的频谱)(ωj e X 经过滤波后)()(ωωj j e H e X ,因此,只要按照输入信号频谱的特点和处理信号的目的, 适当选择)(ωj e H ,使得滤波后的)()(ωωj j e H e X 满足设计的要求,这就是数字滤波器的滤波原理。

1.2数字滤波器的分类按照不同的分类方法,数字滤波器有许多种类,但总起来可以分成两大类:经典滤波器和现代滤波器。

经典滤波器的特点是其输入信号中有用的频率成分和希望滤除的频率成分占有不同的频带,通过一个合适的选频滤波器滤除干扰,得到纯净信号,达到滤波的目的。

但是,如果信号和干扰的频谱相互重叠,则经典滤波器不能有效地滤除干扰,最大限度地恢复信号,这时就需要现代滤波器,例如维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器。

数字信号处理 第六章

数字信号处理 第六章

各种数字滤波器的理想幅度频率响应 数字滤波器的设计步骤 理想滤波器的逼近 数字滤波器的系统函数H(z) IIR滤波器设计方法
6.1 引言

数字滤波器的设计步骤:


按任务要求,确定滤波器性能要求。 用一个因果稳定的离散线性移不变的系统函数去逼 近这一性能要求。逼近所用系统函数有无限冲激响 应(IIR)系统函数与有限长单位冲激响应(FIR) 系统函数两种。 利用有限精度算法来实现这个系统函数。 实际的技术实现。

零极点分布对系统相角的影响

相位“延时”(或相位“滞后”)系统

最小相位延时系统 最大相位延时系统 最大相位超前系统 最小相位超前系统

相位“超前”(或相位“领先”)系统

当全部零点在单位圆外时,相位变化最大,又是负数, 当全部零点在单位圆外时,相位变化最小, 当全部零点在单位圆内时,相位变化最大, 当全部零点在单位圆内时,相位变化最小, 故称为最小相位超前系统。 故称为最大相位超前系统。 故称为最大相位延时系统。 故称为最小相位延时系统。
2、可实现Ha(s)Ha(-s)零极点分布
j
σ
1、零极点中一半属Ha(s),另一 半属Ha(-s)。如要求系统稳定, 则左半平面极点属于Ha(s)。 2、挑选零点时,不加任何限制, 则Ha(s)的解不唯一。 3、如限定Ha(s)是最小相位的, 则只能取所有左半平面的零极 点作为Ha(s)的零极点,Ha(s) 的解唯一。 4、虚轴上的零点阶数减半分配给 Ha(s)。 5、稳定系统虚轴上无极点,临界 稳定时虚轴上才会有极点。
第6章 无限冲激响应IIR 数字滤波器的设计方法
刘笑楠
第6章 无限冲激响应IIR 数字滤波器的设计方法

IIR数字滤波器的设计方法

IIR数字滤波器的设计方法
6.2 IIR滤波器设计的特点
将IIR滤波器的系统函数用极、零点表示:
M
M
bk zk
(1 ck z1)
H(z)
k 0 N

A
k 1 N
1 ak zk
(1 dk z1)
k 1
k 1
M≤N
对系统函数的设计就是确定各系数ak, bk或零极点ck,dk和A, 使滤波器满足给定的性能要求
14
第12讲 无限长单位脉冲响应(IIR)数字滤波器的设计方法 6.3.2 巴特沃思低通逼近 (最平幅度逼近)
巴特沃思低通滤波器在通带内有最大平坦的幅度特 性,因而又称为最平幅度特性滤波器
巴特沃思低通滤波器幅度平方函数定义
|
H
a
(
j)
|2

1

(
1 / c
)2
N
式中,N为正整数,代表滤波器的阶数。Ωc为 3dB截止频率。当Ω=Ωc时,衰减为 3 dB
器• Ha(s)Ha(-s)的极点为
sk
1
(1)2N ( jc )
ej

1 2

22kN1
c
k=1, 2, …, 2N
• Ha(s)Ha(-s)的2N个极点等间隔分布在半径为Ωc的圆(巴特沃 思圆)上,极点间的角度间隔为π/N rad
16
第12讲 无限长单位脉冲响应(IIR)数字滤波器的设计方法
|Ha(jΩ)|2单调减小,N越大,通带内特性越平坦,过渡带越窄
15
第12讲 无限长单位脉冲响应(IIR)数字滤波器的设计方法
巴特沃思滤波器的极(零)点分布 (公式法求解低通Ha(s))
|
H

脉冲响应不变法设计iir数字滤波器

脉冲响应不变法设计iir数字滤波器

脉冲响应不变法设计iir数字滤波器以脉冲响应不变法设计IIR数字滤波器引言:数字滤波器是数字信号处理中的重要组成部分,用于处理和改变数字信号的频率特性。

脉冲响应不变法(Impulse Invariance Method)是一种常用的IIR数字滤波器设计方法,其基本原理是通过将连续时间域中的模拟滤波器的脉冲响应与采样脉冲响应进行匹配,从而实现滤波器的设计。

一、脉冲响应不变法基本原理脉冲响应不变法的基本原理是将连续时间域中的模拟滤波器的脉冲响应与离散时间域中的数字滤波器的脉冲响应进行匹配。

在这种方法中,首先需要确定模拟滤波器的脉冲响应,然后通过采样得到数字滤波器的脉冲响应,最后将其离散化得到数字滤波器的差分方程。

二、脉冲响应不变法的设计步骤1. 确定模拟滤波器的脉冲响应:选择适当的模拟滤波器类型,并设计其频率响应。

根据滤波器的阶数和截止频率,确定模拟滤波器的差分方程。

2. 采样得到数字滤波器的脉冲响应:通过将连续时间域中的模拟滤波器的脉冲响应与采样脉冲进行卷积,得到数字滤波器的脉冲响应。

3. 离散化得到数字滤波器的差分方程:将数字滤波器的脉冲响应离散化,得到数字滤波器的差分方程。

根据差分方程,可以计算数字滤波器的各个系数。

三、脉冲响应不变法的优缺点脉冲响应不变法具有以下优点:1. 设计方法简单:通过匹配模拟滤波器和数字滤波器的脉冲响应,可以直接得到数字滤波器的差分方程,设计方法相对简单。

2. 精度较高:脉冲响应不变法可以保持模拟滤波器的频率响应特性,因此可以实现较高的滤波器精度。

3. 适用范围广:脉冲响应不变法适用于各种模拟滤波器类型和滤波器规格的设计。

然而,脉冲响应不变法也存在一些缺点:1. 频率响应失真:由于采样过程中的截断和抽样误差,脉冲响应不变法可能导致数字滤波器的频率响应失真。

2. 高阶滤波器设计困难:对于高阶滤波器的设计,脉冲响应不变法可能会导致数字滤波器的稳定性问题和数值计算问题。

四、脉冲响应不变法的应用领域脉冲响应不变法广泛应用于数字信号处理领域,特别是在音频信号处理、图像处理和通信系统中的滤波器设计中。

IIR滤波器的设计方法

IIR滤波器的设计方法

把H(z)单位圆外的零点:z 1/ z0, 1/z0*, z0 1 映射到单位圆内的镜像位置:z z0, z0*
构成Hmin(z)的零点。
而幅度响应不变:
P231 图6-6
H (e j ) Hmin (e j ) Hap (e j ) Hmin (e j )
2)级联一个全通系统可以使非稳定滤波器变成 一个稳定滤波器
相位响应:
(e
j
)
arctan
Im[H Re[H
(e (e
j j
)] )]
H *(e j ) H (e j ) e j (e j )
H (e j ) H *(e j )
e2 j (e j )
(e j )
1 2j
H (e j )
ln
H
*
(e
j
)
1 2j
ln
H(z)
H
(
z
1
)
z
阻带最小衰减: 2
H (e j0 )
2 20lg H (e jst ) 20lg H (e jst ) 20lg2
其中: H (e j0 ) 1
当 H (e jc ) 2 / 2 0.707 时,1 3dB
称 c 为3dB通带截止频率
6、表征滤波器频率响应的特征参量
幅度平方响应
H (e j ) 2 H (e j )H *(e j )
2π ω

π
π
2π ω
2、LP到其他滤波器的变换 由LP实现的HP
LP实现的BP
LP实现的BRF
3、 滤波器的性能指标
带宽:当幅度降低到0.707时的宽度称为
滤波器的带宽(3dB带宽)

IIR低通滤波器的设计

IIR低通滤波器的设计

IIR低通滤波器的设计一.设计题目IIR低通滤波器的设计二.设计目的通过本次课程设计,结合所学的理论知识,进一步掌握冲激响应不变法IIR 低通滤波器的设计原理,并利用MATLAB进行设计。

三.设计要求1) 采用冲激响应不变法设计一巴特沃斯数字滤波器;2)数字滤波器的技术指标ωp=0.15π,ωst =0.35π,R p=1dB, R s=20dB,采样=1000HZ;频率fs3) 设计中需要有设计的理论内容;4)设计中必须有每一步的MATLAB处理的结果。

四.设计原理1.数字滤波器设计的基本思路1) 利用ω=ΩT,将数字低通的技术指标转化为模拟低通的技术指标,将ω、stω转换成pΩ、,stΩ而R P ,R S不变;p2)利用巴特沃斯逼近法,求出模拟滤波器的系统函数H a(s);3)利用冲激响应不变法,将模拟滤波器数字化,得到数字滤波器的系统函数H(z);4)利用MATLAB软件实现所设计的H(z)。

数字滤波器的设计步骤如下图所示:图1 数字滤波器设计步骤2.滤波器的技术指标无论是数字滤波器还是模拟滤波器,他们技术指标的建立都是以所谓的“固有衰减”参数为参照。

以数字滤波器为例,固有衰减参数定义为:()()20lgjjH edBH eωδω=-IIR滤波器的设计包括三个步骤:①给出所需要的滤波器的技术指标;②设计一个H(z)使其逼近所需要的技术指标;③实现所设计的H (z)。

目前IIR数字滤波器设计的最通用的方法是借助于模拟滤波器的设计方法。

模拟滤波器设计已经有了一套相当成熟的方法,它不但有完整的设计公式,而且还有较为完整的图表供查询,因此,充分利用这些已有的资源将会给数字滤波器的设计带来很大方便。

IIR数字滤波器的设计步骤是:①按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标;②根据转换后的技术指标设计模拟低通滤波器H a(s);③再按一定规则将H a(s)转换成H(z)。

3.模拟低通滤波器的概述用模拟—数字变换法设计IIR数字滤波器,首先必须设计一个模拟滤波器,它有许多不同的类型,本设计中采用的是巴特沃斯滤波器。

IIR滤波器的设计方法

IIR滤波器的设计方法

实验六 IIR 滤波器的设计方法⏹ 实验学时:2学时 ⏹ 实验类型:设计、验证一、实验目的IR 数字滤波器的常用指标和设计过程的理解;响应不变法和双线性变换法设计IIR 滤波器的原理,设计方法,步骤; 现象,比较脉冲响应不变法和双线性变换法设计IIR 滤波器的特点熟悉Matlab 信号处理工具箱中的常用函数二、实验原理与方法IIR 数字滤波器设计步骤:(一)先设计模拟低通原型滤波器(二)AF 数字化为DF1. 脉冲响应不变法用数字滤波器的单位脉冲响应序列)(n h 模仿模拟滤波器的冲激响应)(t h a ,让)(n h 正好等于)(t h a 的采样值,即)()(nT h n h a =,其中T 为采样间隔,如果以)(s H a 及)(z H 分别表示)(t h a 的拉式变换及)(n h 的Z 变换,则)2(1)(m Tj s H T z H m a e z sT ∑∞-∞==+=π2.双线性变换法S 平面与z 平面之间满足以下映射关系:11112--+-⋅=zz T s (***) s 平面的虚轴单值地映射于z 平面的单位圆上,s 平面的左半平面完全映射到z 平面的单位圆内。

双线性变换不存在混叠问题。

双线性变换是一种非线性变换,这种非线性引起的幅频特性畸变可通过预畸而得到校正。

双线性变换法特别适合用于设计常见的选频性滤波器,下面以双线性变换法设计低通数字滤波器为例,将设计步骤归纳如下:(1)根据需要确定数字滤波器的性能指标:通带截止频率p f 、阻带截止频率p f 、通带最大衰减αp 、阻带内的最小衰减αs ,采样频率s f ; (2) 确定相应的数字角频率,T f p p πω2=;T f s s πω2=; (3) 计算经过预畸的相应模拟低通原型的截止频率,)2(2pp tg T ω=Ω,)2(2s s tg T ω=Ω;(4)根据Ωp ,αp 和Ωs ,αs 计算模拟低通原型滤波器的传递函数)(s H a ; (5)用上述(***)表示的双线性变换公式代入)(s H a ,求出)(z H ;分析滤波器特性,检查其指标是否满足要求。

第六章 无限长单位冲激响应(IIR)数字滤波器的设计方法

第六章 无限长单位冲激响应(IIR)数字滤波器的设计方法

6.1 引 言
数字滤波器的性能要求
通带截止频率 误差容限 阻带截止频率
具有误差容限的的滤波器具有三个特征范围: 通带 过渡带
阻带
1 1 H (e j ) 1,
c
c st
H (e j ) 2 , st
第六章 IIR数字滤波器的设计方法
z e j
H ( z)H ( z 1) 的零极点特征:
1 若 z re ji 是H(z)的极点,则 z e ji 是H(z-1)的极点。 r 即 H ( z)H ( z 1) 的极点是以单位圆镜像对称的,同时也是共轭的。
jIm[z]
满足上述条件的极点可能有几种情况。 对于可实现的系统,系统函数的极点都在单位圆内。 对于零点的分析类似极点,只是系统函数的零点没
0 1
Re[z]
4
有只在单位圆内的限制。
第六章 IIR数字滤波器的设计方法
6.1 引 言
M
设计IIR数字滤波器的方法 数字滤波器可用系统函数表示
H ( z)
1 ak z k
k 0
k 0 N
bk z k
对IIR系统,N>0,且一般有MN。
设计的目的就是要求出ak和bk,使对应的传输函数逼近所要求的特性。
对于因果稳定的LSI系统,其单位冲击响应 h(n)为实函数,因而满足共轭 对称条件,即 H (e j ) H (e j )
第六章 IIR数字滤波器的设计方法
6.1 引 言
幅度平方响应:
H (e
j 2
) H (e j ) H (e j ) H (e j ) H (e j ) H ( z ) H ( z 1 )

IIR滤波器设计

IIR滤波器设计

IIR滤波器的DSP实现一、IIR滤波器的基本结构IIR滤波器差分方程的一般表达式为:式中x(n)为输入序列;y(n)为输出序列;和为滤波器系数.若所有系数等于0,则为FIR滤波器.IIR滤波器具有无限长的单位脉冲响应,在结构上存在反馈回路,具有递归性,即IIR滤波器的输出不仅与输入有关,而且与过去的输出有关.将上式展开得出y(n)表达式为:在零初始条件下,对上式进行z变换,得到:设N=M,则传递函数为:上式可写成:上式具有N个零点和N个极点.若有极点位于单位圆外将导致系统不稳定.由于FIR滤波器所有的系数均为0,不存在极点,不会造成系数的不稳定.对于IIR滤波器,系统稳定的条件如下: 若||<1,当n→时,h(n)→0,系统稳定;若||>1,当n→时,h(n)→,系统不稳定.IIR滤波器具有多种形式,主要有:直接型(也称直接I型)、标准型(也称直接II型)、变换型、级联型和并联型.二阶IIR滤波器,又称为二阶基本节,分为直接型、标准型和变换型.对于一个二阶IIR滤波器,其输出可以写成:1.直接型(直接I型)根据上式可以得到直接二型IIR滤波器的结构图.如图1所示.共使用了4个延迟单元().直接型二阶IIR滤波器还可以用图2的结构实现.图2 直接I型二阶IIR滤波器此时,延时变量变成了w(n).可以证明上图的结构仍满足二阶IIR滤波器输出方程.前向通道:反馈通道:将1.2式代入1.1式可得:2.标准型(直接II型)从图2可以看出,左右两组延迟单元可以重叠,从而得到标准二阶IIR滤波器的结构图,如图3所示.由于这种结构所使用的延迟单元最少(只有2个),得到了广泛地应用,因此称之为标准型IIR滤波器.图3 标准型二阶IIR滤波器二、IIR滤波器的设计IIR滤波器的设计可以利用模拟滤波器原型,借鉴成熟的模拟滤波器的设计结果进行双线性变换,将模拟滤波器变换成满足预定指标的数字滤波器,即根据模拟设计理论设计出满足要求的传递函数H(s),然后将H(s)变换成数字滤波器的传递函数H(z).设计IIR滤波器的基础是设计模拟滤波器的原型,这些原型滤波器主要有:①巴特沃兹(Butterworth)滤波器,其幅度响应在通带内具有最平特性;②切比雪夫(Chebyshev)滤波器,在通带内具有等波纹特性,且阶数小于巴特沃兹滤波器.③椭圆(Elliptic)滤波器,在通带内具有等波纹特性,且阶数最小.将模拟滤波器转换为数字滤波器常用的方法是双线性变换,其作用是完成从s平面到z平面的一个映射.其关系为:双线性变换的基本性质如下:①s平面上的轴映射到z平面的单位圆上;②s平面的左半平面映射到z平面的单位圆内;③s平面的右半平面映射到z平面的单位圆外.考虑到s平面上的虚轴映射为z平面的单位圆,令,它代表一个可变的模拟频率.其z平面上相应的数字频率为,即,将以上两式代入2.1式得:对上式求解得:模拟频率和相应的数字频率之间的映射关系如图4所示:图4 和之间的映射关系当在0~1之间变化时,在0~之间变化,为采样频率:当时,其对应的在之间.双线性变换会造成频率失真,通常采用预畸变来补偿频率失真.双线性变换设计的步骤如下:①选择一个合适的模拟传递函数H(s);②对截止频率或预定的数字频率进行畸变,并根据式2.3求得相应的模拟频率;③用对H(s)中的频率进行换算,即:④用式2.1计算H(z):三、二阶IIR滤波器的DSP实现1.标准型二阶IIR滤波器的实现在二阶IIR滤波器结构中,标准型结构是最常见的滤波器结构,其结构如图5所示:图5 标准型二阶IIR滤波器由结构图可以写出反馈通道和前向通道的差分方程:反馈通道:前向通道:由以上两式对二阶IIR滤波器进行编程,其中乘法-累加运算可采用单操作数指令或双操作数指令,数据和系数可存放在DARAM中,如图6所示:图6 双操作数数据存放和系数表2.直接型二阶IIR滤波器的实现二阶IIR滤波器可以用直接型结构来实现.在迭代运算中,先衰减后增益,系统的动态范围和鲁棒性要好些.直接型二阶IIR滤波器的结构如图7所示:图7 直接型二阶IIR滤波器直接型二阶IIR滤波器的脉冲传递函数为:差分方程为:为了实现直接型滤波,可在DARAM中开辟4个循环缓冲区,用来存放变量和系数,并采用循环缓冲区方式寻址.这4个循环缓冲区的结构如图8所示:图8 循环缓冲区结构四、总体设计方案1.利用MATLAB来确定IIR滤波器的参数;2.启动CCS,在CCS中建立一个汇编源文件、建立一个C源文件和一个命令文件,并将这三个文件添加到工程,再编译并装载程序;3.设置波形时域观察窗口,得到其滤波前后波形变化图;4.设置频域观察窗口,得到其滤波前后频谱变化图.五、主要参数f1=500Hz f2=1000Hz f s=2000六、源程序1.汇编源文件程序(.asm).global _iir,_init,_outdata_iirbset frctsub #1,t0mov t0,mmap(csr)add t0,ar0mov xar2,xdpmov ar2,cdpmov #0,ac0rpt csrmacmz *ar0-,*cdp+,ac0add t0,ar1mov xar3,xdpmov ar3,cdprpt csrmacmz *ar1-,*cdp+,ac0mar *ar1+mov hi(ac0),*ar1mov hi(ac0),t0bclr frctret_initsub #1,t0mov t0,mmap(csr)mov #0,ac0rpt csrmov ac0,*ar0+rpt csrmov ac0,*ar1+ret_outdatamov t1,ac0sub #2,ac0mov ac0,mmap(csr)add ac0,ar0rpt csrdelay *ar0-mar *ar0+mov t0,*ar0ret2.C源文件程序(.c)#includ e "math.h"#define signal_1_f 500#define signal_2_f 10000#define signal_sample_f 25000#define pi 3.1415926#define IIRNUMBER_L 2#define bufer_L 256int N_L=IIRNUMBER_L;int data_in[bufer_L];int out[bufer_L] ;int x[IIRNUMBER_L+1];int y[IIRNUMBER_L+1];int k=0;int bufer=bufer_L;int fBn[IIRNUMBER_L]={0,0x634a};int fAn[IIRNUMBER_L]={0xe5c,0xe5c};extern int iir(int *x,int *y,int *fAn,int *fBn,int N_L); extern int init(int *,int *,int);extern int outdata(int *,int,int);void inputwave();void main(){int iirout;inputwave();init(x,y,N_L);while(1){x[0]=data_in[k];iirout=iir(x,y,fAn,fBn,N_L);outdata(out,iirout,bufer);k++;if(k>=bufer_L){k=0;}}}void inputwave(){fl oat wt1;fl oat wt2;int i;for(i=0;i<=bufer_L;i++){wt1=2*pi*i*signal_1_f;wt1=wt1/signal_sample_f;wt2=2*pi*i*signal_2_f;wt2=wt2/signal_sample_f;data_in[i]=(cos(wt1)+cos(wt2))/2*32768;}}3.命令文件程序(.cmd)-stack 0x0500-sysstack 0x0500-heap 0x1000-c-u _Reset-l rts55.libMEMORY{PAGE 0:RAM(RWIX): origin=0x000100, length=0x01ff00ROM(RIX): origin=0x020100, length=0x01ff00VECS(RIX): origin=0xffff00, length=0x000200 PAGE 2:IOPORT(RWI):origin=0x000000, length=0x020000 }SECTIONS{.text >ROM PAGE 0.data >ROM PAGE 0.bss >RAM PAGE 0.const >RAM PAGE 0.sysmem >RAM PAGE 0.stack >RAM PAGE 0.cio >RAM PAGE 0.sysstack >RAM PAGE 0.switch >RAM PAGE 0.cinit >RAM PAGE 0.pinit >RAM PAGE 0.vectors >VECS PAGE 0.ioport >IOPORT PAGE 2}七、实验结果及分析1.输入波形:①时域波形:②频域波形:2.输出波形:①时域波形:②频域波形:八、心得体会课程设计是培养学生综合运用所学知识,发现,提出,分析和解决实际问题,锻炼实践能力的重要环节,是对学生实际工作能力的具体训练和考察过程.随着科学技术发展的日新日异,数字技术已经成为当今计算机应用中空前活跃的领域, 在生活中可以说得是无处不在.通过这次课设,我掌握了数字滤波器的设计过程,了解了IIR的原理和特性,学习IIR滤波器的DSP 实现原理.熟悉设计IIR数字滤波器的原理和方法,学习使用CCS的波形观察窗口观察输入输出信号波形和频谱变化情况.通过这次课程设计使我懂得了理论与实际相结合是很重要的,只有理论知识是远远不够的,只有把所学的理论知识与实践相结合起来,从理论中得出结论,才能真正为社会服务,从而提高自己的实际动手能力和独立思考的能力.DSP原理及应用课程设计报告设计课题: IIR滤波器的DSP实现专业班级: 信息07-2班姓名: 王晓萱学号: 07100637同组人: 王鑫、蔚金昌、吴杰、徐扬、李峰。

第六章IIR滤波器的设计方法

第六章IIR滤波器的设计方法
一. IIR滤波器的一般设计方法:
1.累试法
滤波器的幅度特性和相位特性完全由其零点和极点位置所决定。
如果在单位圆内处
设置一对共轭极点
j0 Z r e ,频响 在 0 将有一峰值。

极点 0 离单位圆愈远,频响在 0 处的峰值比较平缓。 极点 0 越接近单位圆,频响在 0 处的峰值就越尖锐。 如果通带太窄,可以把极点向原点平移,如果通带太宽,可以将极点

2 m Ha j T m


表达式表明:数字滤波器的频响并不是简单的重现模拟滤波器的频响,
H e j 为模拟滤波器频响

的周期重复 ,是以 j Ha
为周期的周期延。 s
需要强调指出,在周期重复模拟滤波器频响的过程中所存在的问题:
首先分析一下
变函数的映射变换,这种映射变换应该满足两个基本的要求: ① H z 的频响应该模仿 H a s 的频响
即要求
②是因果稳定的映射 指 H a s 的因果稳定性通过映射后, H z 仍保持因果稳定。
§6.5 脉冲响应不变法
根据容限设计好一个模拟滤波器后,就可对此模拟系统进行模仿。 数字滤波器从什么角度去模仿模拟滤波器呢?第一种方法是脉冲响应 不变法。 1. 脉冲响应不变法
1
N j 1 1 j C

2N
N —— 为整数,表示滤波器的阶次;
c——
为截止频率。
当 0 时, 当 时, c
H a j 0 1
2
H a j C
2
1 2
又称为滤波器的3bB带宽
巴特沃斯低通滤波器的特点:
①在 0 ,所以巴特沃斯滤波器通带内具 处,即靠近零频处,衰减为 0 有最大平坦的振幅特性,故得名为最平坦响应滤波器。

IIR滤波器的设计和实现

IIR滤波器的设计和实现

IIR滤波器的设计和实现摘要目前IIR滤波器的设计方法主要有两种,数模转换法和计算机辅助设计法。

数模转化法就是借助模拟滤波器进行设计。

先按照数字滤波器的设计指标得到模拟滤波器的系统函数,再通过一定的转换方式将其转换为数字滤波器。

而计算机辅助设计法就是将滤波器系数作为要求解的未知数,根据设计指标设联立方程,最终求解得出滤波器系数。

在IIR滤波器的设计当中,使用PSO算法能简化步骤。

种群中每个粒子的位置都是每一组滤波器系数解,而适应函数则是该系数所代表的滤波器幅频响应和理想滤波器之间的均方差之和。

这样适应函数越小,滤波器越接近设计目标。

经过多次迭代,我们求得的全局最佳值即是最好的一组滤波器系数。

关键词:IIR滤波器,双线性转换法,脉冲响应不变法,巴特沃兹滤波器,粒子群优化算法Design and implementation of IIR filterAbstractAt present, there are mainly two methods of IIR filter design, DAC and CAD. The analog to digital conversion method is designed by means of analog filter. First, according to the design index of the digital filter, the system function of the analog filter is obtained, and then converted to digital filter through a certain conversion mode. The computer aided design method is to use the filter coefficient as the unknowns of the required solution, and set up a joint equation according to the design index. Finally, the filter coefficient is obtained.In the design of IIR filter, the PSO algorithm can simplify the steps. The position of each particle in the population is the filter coefficient solution for each group, while the adaptive function is the sum of the amplitude frequency response of the filter and the mean square variance between the ideal filter. The smaller the adaptive function, the closer the filter is to the design target. After many iterations, the best global value is the best set of filter coefficients.Keywords:IIR filter,Bilinear transformation method,impulse invariance,butterworth,PSO目录摘要 (1)Abstract (2)第一章绪论 (4)1.1 研究数字滤波器的背景 (4)1.1.1 数字滤波器 (4)1.1.2 IIR滤波器 (4)1.2 国内外研究现状 (5)1.3 本文的研究内容 (5)第二章采纳数模转换法设计IIR滤波器 (6)2.1 IIR滤波器的简介 (6)2.1.1 IIR滤波器的特性 (6)2.1.2 IIR滤波器的结构 (7)2.1.3 模拟滤波器的特性 (10)2.2 采纳双线性变换法设计IIR滤波器 (11)2.2.1 双线性变换法的基本原理 (11)2.2.2 基于双线性变换法的IIR滤波器设计实例 (13)2.3 采纳脉冲响应不变法设计IIR滤波器 (13)2.3.1 脉冲响应不变法的基本原理 (13)2.3.2 基于脉冲响应不变法的IIR滤波器设计实例 (15)2.4 本章小结 (15)第三章基于PSO算法设计IIR滤波器 (16)3.1 粒子群优化算法 (16)3.1.1 PSO算法的进展 (16)3.1.2 PSO算法的基本原理 (16)3.2 IIR滤波器设计的PSO算法 (18)3.2.1设计原理 (18)3.2.2设计实例 (19)3.3 本章小结 (20)结论 (21)参考文献 (22)致谢 (23)第一章绪论1.1 研究数字滤波器的背景1.1.1 数字滤波器现代社会中,人们的日常生活越来越数字化,越来越多的数字信号需要处理。

手把手教系列之IIR滤波器设计

手把手教系列之IIR滤波器设计

⼿把⼿教系列之IIR滤波器设计[导读]:在嵌⼊式系统中经常需要采集模拟信号,采集模拟信号的信号链中难免引⼊⼲扰,那么如何滤除⼲扰呢?今天就来个⼀步⼀步描述如何设计部署⼀个IIR滤波器到你的系统。

何为IIR滤波器?⽆限冲激响应(IIR: Infinite Impulse Response)是⼀种适⽤于许多线性时不变系统的属性,这些系统的特征是具有⼀个冲激响应h(t),该冲激响应h(t)不会在特定点上完全变为零,⽽是⽆限期地持续。

这与有限冲激响应(FIR: Finite Impulse Response)系统形成对⽐,在有限冲激响应(FIR)系统中,对于某个有限T,在时间t> T时,冲激响应确实恰好变为零。

线性时不变系统的常见⽰例是⼤多数电⼦和数字滤波器。

具有此属性的系统称为IIR系统或IIR滤波器。

这是常见的教科书式数学严谨定义,很多⼈看到这⼀下就蒙了,能说⼈话吗?线性时不变系统理论俗称LTI系统理论,源⾃应⽤数学,直接在核磁共振频谱学、地震学、电路、信号处理和控制理论等技术领域运⽤。

它研究的是线性、⾮时变系统对任意输⼊信号的响应。

虽然这些系统的轨迹通常会随时间变化(例如声学波形)来测量和跟踪,但是应⽤到图像处理和场论时,LTI系统在空间维度上也有轨迹。

因此,这些系统也被称为线性⾮时变平移,在最⼀般的范围理论给出此理论。

在离散(即采样)系统中对应的术语是线性⾮时变平移系统。

由电阻、电容、电感组成的电路是LTI系统的⼀个很好的例⼦。

⽐如⼀个运放系统在⼀定频带范围内满⾜信号的时域叠加,输⼊⼀个100Hz和200Hz正弦信号,输出频率是这两种信号的线性叠加。

⽤数学对LTI系统描述:线性:输⼊x_1(t),产⽣响应y_1(t),⽽输⼊x_2(t),产⽣响应y_2(t),那么缩放和加和输⼊a_1x_1(t)+a_2x_2(t),产⽣缩放、加和的响应a_1y_1(t)+a_2y_2(t),其中a_1和a_2是标量,对于任意的有:输⼊\sum_{k=0}^{N}a_kx_k(t),产⽣响应\sum_{k=0}^{N}a_ky_k(t)时不变性:指如果将系统的输⼊信号延迟\tau秒,那么得到的输出响应也相应延时\tau秒。

IIR数字滤波器的设计及软件实现

IIR数字滤波器的设计及软件实现

IIR 数字滤波器的设计及软件实现一.实验目的(1)熟悉用双线性变换法设计IIR 数字滤波器的原理与方法; (2)学会用MATLAB 信号处理工具箱中的滤波器设计函数(或滤波器设计分析工具FDAtool )设计各种滤波器,学会根据滤波需求确定滤波器指标参数; (3)掌握IIR 数字滤波器的MATLAB 实现方法;(4)通过观察滤波器输入、输出信号的时域波形及其频谱,建立数字滤波的概念。

二.实验原理设计IIR 数字滤波器一般采用间接法(脉冲响应不变法和双线性不变法),应用最广泛的是双线性变换法。

基本的设计过程是:①将给定的数字滤波器指标转换成模拟滤波器的指标; ②涉及模拟滤波器;③将模拟滤波器的系统函数转换成数字滤波器的系统函数。

MATLAB 信号处理工具箱中的各种IIR 数字滤波器设计函数都是采用双线性变换法。

本实验的数字滤波器的MATLAB 实验是调用MATLAB 信号处理工具箱的函数filter 对给定的输入信号x(n)进行滤波,得到滤波后的输出信号y(n). 三.实验容及步骤1.信号处产生函数mstg 产生由三路抑制载波调幅信号相加构成的复合信号st ,该函数还会自动回图显示st 的时域波形和幅频特性曲线,由后图可见,三路信号时域混叠无法在时域分离。

但频域是分离的,所以可通过滤波的方法在频域分离。

2.将st 中三路调幅信号分离,通过观察st 的幅频特性曲线,分别确定可以分离st 中三路抑制载波单频调幅信号的三个滤波器(低通滤波器、带通滤波器、高通滤波器)的通带截止频率和阻带截止频率。

且滤波器的通带最大衰减为0.1dB ,阻带最小衰减为60bB 。

提示:抑制载波单频调幅信号的数学表示式为()()()()()()[]tt t t t s f f f f f f c c c 0002cos 2cos 212cos 2cos )(++-==ππππ 其中,()t fc2cos π称为载波,fc为载波频率,()t f2cos π称为单频调制信号,f 0为调制正弦波信号频率,且满足fc>f 0。

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

H (e ) K
e
m 1 N k 1
M
j
cm
j e dk
各零矢量模的连乘积 各极矢量模的连乘积

相角:
N H (e j ) M j j arg arg[ e c ] arg[ e d k ] ( N M ) m k 1 K m1
N H (e j ) M j j arg arg[e cm ] arg[e d k ] ( N M ) k 1 K m1
令: 单位圆内零点数为mi
单位圆外的零点数为mo
单位圆内的极点数为pi
mi mo M
单位圆外的极点数为po
2)全部零点在单位圆外: mi 0, mo M
arg[] 2 ( N M ) 为最小相位超前系统
最小相位延时系统的性质
1)在
H (e )
j
相同的系统中,具有最小的相位滞后
2)最小相位延时系统的能量集中在n=0附近,而总 能量相同 m m 2 2 h(n ) hmin (n ) m N 1
H (e ) H (e ) e
*
j
j
j (e
j
)
H (e j ) 2 j ( e j ) e * j H (e )
j 1 H (e ) 1 H ( z ) j (e ) ln * j ln 1 2 j H (e ) 2 j H ( z ) z e j H (e j )
6.3 全通系统
对所有,满足: H ap (e j ) 1
称该系统为全通系统

一阶全通系统:
z 1 a H ap ( z ) 1 az 1 a为实数
0 a 1
极点:z
a
z 1/ a 零点:
a为复数
z 1 a* H ap ( z ) 1 1 az
第六章 IIR滤波器的设计
主要内容


理解数字滤波器的基本概念
了解最小相位延时系统
理解全通系统的特点及应用
掌握冲激响应不变法 掌握双线性变换法 掌握Butterworth、Chebyshev低通滤波器的特点 了解利用模拟滤波器设计IIR数字滤波器的设计过程

了解利用频带变换法设计各种类型数字滤波器的方法

N 阶数字全通滤波器
N 1 *
1 ( N 1) N d d z ... d z z z ak N N 1 1 H ( z ) 1 ( N 1) N 1 1 d z ... d z d z 1 a z k 1 1 N 1 N k
0 a 1 * 极点:z a 零点: z 1/ a
零极点以单位圆为镜像对称

实系数二阶全通系统
z a z a H ap ( z ) 1 1 1 az 1 a * z
*
1
1
a 1
极点: z a,a*
* z 1/ a ,1/a 零点:
两个零点(极点)共轭对称 零点与极点以单位圆为镜像对称
相位延时系统
1)全部零点在单位圆内: mi M , mo 0
arg[] 0 为最小相位延时系统
2)全部零点在单位圆外: mi 0, mo M arg[] 2 M 为最大相位延时系统

逆因果稳定系统 z r , r 1 n > 0时,h(n) = 0
全部极点在单位圆外:po = N,pi = 0

滚降与滚降率:滤波器幅频特性在过渡 带的衰减和衰减速度称为滚降与滚降率。

阻带衰减:输入信号在阻带的衰减量

带内平坦度:通带和阻带内的平坦程度
4、数字滤波器的设计步骤
数字滤波器的设计三个步骤:
(1) 按要求确定滤波器的性能参数; (2) 用一个因果稳定的离散线性移不变系统的系 统函数去逼近去逼近这一性能要求;
6.1 引言
数字滤波器: 是指输入输出均为数字信号,通过一定运算 关系改变输入信号所含频率成分的相对比例或者 滤除某些频率成分的器件。
优点:
高精度、稳定、体积小、重量轻、灵活,不要 求阻抗匹配,可实现特殊滤波功能
1、滤波器的基本概念
(1) 滤波器的功能
滤波器的功能是对输入信号进行滤波 以增强所需信号部分,抑制不要的部分。
器后在时间上的延时情况
理想滤波器不可实现,只能以实际滤波器逼近

通带:
c
st
1 1 H (e ) 1
j
阻带:
过渡带:
H (e j ) 2
c st
1 :通带容限
2 :阻带容限
c :通带截止频率 st :阻带截止频率
N H (e j ) M j j arg arg[e cm ] arg[e d k ] ( N M ) k 1 K m1
当 0 2 , 2
j Im[ z ]
0
Re[ z ]
位于单位圆内的零/极矢量角度变化为2 位于单位圆外的零/极矢量角度变化为 0
H ( e j )
π
2π ω
带阻
2π π π 2π ω
2、LP到其他滤波器的变换
由LP实现的HP
LP实现的BP
LP实现的BRF
3、 滤波器的性能指标

带宽:当幅度降低到0.707时的宽度称为
滤波器的带宽(3dB带宽)

通带、阻带与过渡带:信号允许通过的频带 为通带,完全不允许通过的频带为阻带,通 带与阻带之间为过渡带。
可以证明,一阶全通节在任何频率上,其频率响应的模 都为1:
z a 1 az 1
1
* z e j
e re 1 re j e j j 1 re e 1 re j ( )
j
j
j ( )
1 r cos( ) jr sin( ) 1 1 r cos( ) jr sin( )
6、表征滤波器频率响应的特征参量

幅度平方响应
H (e ) H (e ) H (e )
*
j
2
j
j
H (e ) H (e
1
j
j
) H ( z) H ( z )
j Im[ z ]
1
z e j
H ( z ) H ( z ) 的极点既是共轭的,又是以单
位圆成镜像对称的
1/ a*

群延迟响应
相位对角频率的导数的负值
d (e ) (e ) d
j
j
dH ( z ) 1 Re z j dz H ( z ) z e
若滤波器通带内 (e ) = 常数,
j
则为线性相位滤波器
7、IIR数字滤波器的设计方法
用一因果稳定的离散LSI系统逼近给定的性能要求:
(3) 用有限精度的运算实现;实现可以采用通用 计算机,也可以采用DSP。
5、数字滤波器的技术要求

选频滤波器的频率响应:
H (e ) H (e ) e
j
j
j
j ( j )
H (e ) 为幅频特性:表示信号通过该滤波器
后各频率成分的衰减情况
( j ) 为相频特性:反映各频率成分通过滤波
Kz ( N M )
(z c
m 1 N k 1
M
m
)
(z d )
k
j j arg[ H ( e j )]
频率响应:
H (e ) Ke
j
j ( N M ) m 1 N k 1
j ( e cm ) j
(e
H (e ) e
dk )

模:
j
H (e j ) arg 2 mi 2 pi 2 ( N M ) K 2 2 mi 2 ( N M ) 0
相位超前系统 1)全部零点在单位圆内: mi M , mo 0
arg[] 2 N
为最大相位超前系统

n 0 N 1
n 0

h(n )
2
hmin (n )
n 0
n 0 N 1
2
3)最小相位序列的 4)在
j
hmin (0) 最大:hmin (0) h(0)
hmin (n) 唯一 H (e ) 相同的系统中,
5)级联一个全通系统,可以将一最小相位系统转变成一 相同幅度响应的非最小相位延时系统
a) 时域说明
b) 频域说明
(2) 四种基本的滤波器
四种基本滤波器为低通(LP)、高通 (HP)、带通(BP)和带阻滤波器(BRF):
(3) 四种基本滤波器的数字表示
H ( e j )
低通
2π π
H ( e j )
π
2πBiblioteka ω高通 2π π
H ( e j )
π
2π ω
带通
2π π
则:
pi po N
H ( e j ) arg 2 ( N M ) 2 mi 2 pi K 2

因果稳定系统
j
z r , r 1 n < 0时,h(n) = 0
全部极点在单位圆内:po = 0,pi = N
H (e ) arg 2 mi 2 pi 2 ( N M ) K 2 2 mi 2 M 2 mo 0
二、一般转换方法
1、 DLF ALF 2、 DHF AHF ALF 3、 DBF ABF ALF 4、 DBSF ABSF ALF
相关文档
最新文档