第6章 滤波器
现代信号处理--清华胡广书讲义-第6章滤波器组基础
150第6章 滤波器组基础6.1 滤波器组的基本概念一个滤波器组是指一组滤波器,它们有着共同的输入,或有着共同的输出,如图6.1.1所示。
图6.1.1 滤波器组示意图,(a )分析滤波器组,(b )综合滤波器组。
假定滤波器)(0z H ,)(1z H ,…,)(1z H M -的频率特性如图6.1.2(a )所示,)(n x 通过这些滤波器后,得到的)(0n x ,)(1n x ,…,)(1n x M -将是)(n x 的一个个子带信号,它们的频谱相互之间没有交叠。
若)(0z H ,)(1z H ,…,)(1z H M -的频率特性如图6.1.2(b )所示,那么,)(0n x ,)(1n x ,…,)(1n x M -的频谱相互之间将有少许的混迭。
由于)(0z H ,)(1z H ,…,)(1z H M -的作用是将)(n x 作子带分解,因此我们称它们为分析滤波器组。
将一个信号分解成许多子信号是信号处理中常用的方法。
例如,若图6.1.1中的2=M ,那么,在图6.1.2中,)(0z H 的频率特性将分别占据2~0π和ππ~2两个频段,前者对应低频段,后者对应高频段。
这样得到的)(0n x 将是)(n x 的低频成份,而)(1n x 将是其高频)(0n x )(1n x )(1n x M -)(n x(ˆ0x (ˆ1x)(ˆ1n xM -)(ˆn x151成份。
我们可依据实际工作的需要对)(0n x 和)(1n x 作出不同的处理。
例如,若我们希望对)(n x 编码,设)(n x 的抽样频率为20KHz ,若每个数据点用16bit ,那么每秒钟需要的码图6.1.2 分析滤波器组的频率响应,(a )无混迭,(b )稍有混迭流为320Kbit 。
若)(n x 是一低频信号,也即)(n x 的有效成份(或有用成份)大都集中在)(0n x 内,)(1n x 内含有很少的信号能量。
这样,我们可对)(0n x 仍用16bit ,对)(1n x 则用8bit ,甚至是4bit ,由于)(0n x 和)(1n x 的带宽分别比)(n x 减少了一倍,所以,)(0n x 和)(1n x 的抽样频率可降低一倍。
数字信号处理 第六章
各种数字滤波器的理想幅度频率响应 数字滤波器的设计步骤 理想滤波器的逼近 数字滤波器的系统函数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 数字滤波器的设计方法
第6章 匹配滤波器(第18讲111109).
(2)对于τ的变化
MF是用t0的改变来适应 采样峰值出现的时刻改变,采样时刻t0随之改变即可
信息传输与处理实验室 蒋铃鸽
5. t0 的选择—系统的因果性
作为线性时不变系统的MF要满足因果性
hm
(t
)
=
⎧ ⎨
si
⎩
(t0
− 0
t
)
t≥0 t<0
设在t=0注入MF,持续期为 [0,T ]
hm
(
t
)
=
⎧⎨si ⎩
(
t0 − 0
t
)
0 ≤ t < t0 其余
为把 si (t) 全部纳入 [0,t0 ] t0 ≥ T
t0 = T
即选择在输入信号码元周期T的末尾
信息传输与处理实验室 蒋铃鸽
si (t )
t
0
T
si (−t )
不可实现部分
0
hm (t )
信息传输蒋与铃处鸽6理.实1验最室 大信噪比准则—MF
6.1.1 瞬时信噪比表示
接收信号: x (t ) = si (t ) + ni (t )
∫ A( f ) =
( ) ∞
−∞ si
t
e− jωt dt
白噪声 Sni ( f ) = N0 / 2
经滤波后: y (t ) = so (t ) + no (t )
=
T
∫0
⎣⎡si
(τ
)
+
ni
(τ
)⎤⎦
si
(T
−
t
+
τ
现代信号处理基础及应用6章-白化滤波器
内积空间:
设有 M 个两两正交的随机矢量 ε1, ε2, , εM ,满足
εi , εj 0, i j
令 Y=ε1, ε2,
εM , 是由这 M 个随机矢量张成的线
性子空间,那么随机矢量就是该内积空间的正交基底。
根据正交分解定理,对于任何随机矢量 x , 相对于线性子空间 Y ,可唯一分解为两个互 相正交的部分,即
D(
z)
可能不是因果的,
D(z) G(z)
就不是因果的;
D(z) (3) G(z) 对应的是一个因果稳定的 IIR 滤波器,而所
设计的 H(z) 是一个 n 阶的 FIR 滤波器。
上述因素都会使滤波器的实际输出 y = g* h 不一定等
于期望输出 d 。
设 d l2 g l2 ,且 g 是因果的。则滤波器实际输出 y 与期望
lim n
δ-
gh
2
1
1 2
例 6.2 假设信道的传递函数为G(z) 1 z2 ,它是非最小相 4
设 x = x1 + x2 ,其中 x1 与 y 相关, x2 与 y 不相关,由 于 Rxy R[ xyT ] E[( x1 x2 ) yT ] Rx1y Rx2 y Rx1y
所以, xˆ Rxy Ry1 y Rx1y Ry1 y xˆ1,因此, xˆ 实际上就 是对 xˆ1的估计,即对 x 中与 y 相关部分的估计。所 以相关抵消器的输出中与 y 相关的部分 x1 得到了 尽可能大的抵消。
Yn ε1, ε2 , , εn y1, y2,
, yn
用符号 yˆn n1 来表示 yn 在子空间Yn1 上的正交投影即
n 1
1
yˆn n1 E yni E ii i
数字信号处理第六章
1)幅度函数特点:
H a ( j)
2
1 1 c
2
2N
0
c
H a ( j) 1 H a ( j) 1/ 2 1 3dB 3dB不变性
2
c 通带内有最大平坦的幅度特性,单调减小
c 过渡带及阻带内快速单调减小
3、逼近情况
1)
s平面虚轴
2)
z平面单位圆
s平面
左半平面
z平面 单位圆内 单位圆外 单位圆上
右半平面
虚轴
例7.4
已知模拟滤波器的传输函数为
1 H a ( s) 2 2s 3s 1
采用双线性变换法将其转换为数字滤波 器的系统函数,设T=2s 解 将s代入Ha(s)可得
H ( z ) H a ( s ) s 2 1 z 1 ,T 2
i 1,2,..., m
例6.4.1试分别用脉冲响应不变法和双 线性不变法将图6.4.4所示的RC低通滤波器 转换成数字滤波器。 解 首先按照图6.4.4写出该滤波器的传 输函数Ha(s)为 1
H a ( s)
s
,
RC
利用脉冲响应不变法转换,数字滤波器的系统函 数H1(z)为
低通
0 高通
0 带通 0
带阻
0
全通 0
通带
阻带 过渡带 平滑过渡
三、DF频响的三个参量 1、幅度平方响应
2、相位响应
3、群延迟
它是表示每个频率分量的延迟情况;当其为常 数时, 就是表示每个频率分量的延迟相同。 四、DF设计内容 1、按任务要求确定Filter的性能指标; 2、用因果稳定LSI的系统函数去逼近这一性 能要求; 3、选择适当的运算结构实现这个系统函数; 4、用软件还是用硬件实现。
数字信号处理 第6章
H ( z ) h( n) z n
n 0
N 1
(6.1.2)
(6.1.1)式中的H(z)称为N阶IIR数字滤波器系统函数; (6.1.2) 式中的H(z)称为N-1阶FIR数字滤波器系统函数。这两种 数字滤波器的设计方法有很大区别,因此下面分成两章分 别进行学习。
第6章 无限脉冲响应数字滤波器的设计
s 20 lg
| H (e j0 ) |
j s
dB
(6.1.4b)
p 20 lg | H (e
j p
) | dB
(6.1.5)
s 20 lg | H (e js ) | dB
(6.1.6)
第6章 无限脉冲响应数字滤波器的设计
当幅度下降到 2 / 2 时,标记ω=ωc,此时 p 3dB,称 ωc为3 dB通带截止频率。ωp、ωc和ωs统称为边界频率, 它们是滤波器设计中所涉及到的很重要的参数。对其他 类型的滤波器,(6.1.3b)式和(6.1.4b)式中的H(ej0)应改 成
拟滤波器得到系统函数Ha (s),然后将Ha(s)按某种方法转
换成数字滤波器的系统函数H(z)。这是因为模拟滤波器的 设计方法已经很成熟,不仅有完整的设计公式,还有完善
的图表和曲线供查阅; 另外,还有一些典型的优良滤波
器类型可供我们使用。直接法直接在频域或者时域中设计 数字滤波器,由于要解联立方程,设计时需要计算机辅助 设计。FIR滤波器不能采用间接法,常用的设计方法有窗 函数法、频率采样法和切比雪夫等波纹逼近法。
第6章 无限脉冲响应数字滤波器的设计
图6.1.3所示的单调下降幅频特性,p和s别可以表
示为
p 20 lg
| H (e j0 ) | | H (e
第6章 数字滤波器的基本结构
由于滤波器的系数应为实数 ; 将分子、 将分子、分母中的共轭复根因子合并为二阶实系 数因子,得到: 数因子,得到:
H ( z) = K ⋅
返回到本节向导
(1 − pm z ) ∏ (1 +β m z −1 +β2 m z −2 ) ∏ 1
−1 m =1 N1
M1
M2
∏ (1 − ck z −1 ) ∏ (1 +α1k z −1 +α2k z −2 )
返回到本节向导
6.2.2
直接型结构
5 2 3 + z −1 + z −2 3 3 H ( z) = 1 1 1 1 + z −1 + z −2 − z −3 6 3 6
例6.2-2 已知 3 阶IIR数字滤 6.2IIR数字滤 波器的系统函数; 波器的系统函数;
求:直接Ⅰ型、直接Ⅱ型和转置直接Ⅱ型结构; 直接Ⅰ 直接Ⅱ型和转置直接Ⅱ型结构;
表明:滤波器可以由若干一阶和二阶子系统级联 表明: 组成, 组成,从而构成滤波器的级联型结构 ; 将分子、分母中一阶因子(即实零、极点因子) 将分子、分母中一阶因子(即实零、极点因子) 两两合并为实系数二阶因子,得到: 两两合并为实系数二阶因子,得到:
H ( z) = K
返回到本节向导
∏
k =1
N0
返回到本节向导
6.2.2 例 6.2-2 6.2-
直接型结构
返回到本节向导
6.2.2
直接型结构
直接型结构的优点:简单、直观 直接型结构的优点:简单、
缺点: 缺点: 系数 bm 和 ak 对滤波器性能的控制关系不
直接,调整困难; 直接,调整困难; 极点分布对系数变化的灵敏度高, 零、极点分布对系数变化的灵敏度高, 对有限字长效应敏感, 对有限字长效应敏感,易引起不稳定现 象和较大的误差; 象和较大的误差; 滤波器的阶次越高,上述问题就越明显; 滤波器的阶次越高,上述问题就越明显; 产生上述缺点的原因: 产生上述缺点的原因: 不明显; 不明显; 的改变会影响所有零点或极点的分布; 且 bm、ak 的改变会影响所有零点或极点的分布;
现代信号课件第6章自适应滤波课件
自适应滤波器的计算复杂度较高 ,尤其是在处理大规模数据时,
计算量会变得非常大。
计算复杂度问题可能导致滤波器 实时性差、功耗大等问题,限制
了其在某些领域的应用。
解决计算复杂度问题的方法包括 优化算法、采用并行计算等技术
。
自适应滤波器的未来发展方向
未来自适应滤波器的发展方向主要包 括提高性能、降低计算复杂度、拓展 应用领域等方面。
自适应滤波器的发展历程
20世纪50年代
20世纪60年代
自适应滤波器的概念开始出现,最早的应 用是在通信领域。
线性自适应滤波器的研究取得突破性进展 ,如最小均方误差(LMS)算法和递归最 小二乘(RLS)算法等。
20世纪70年代
21世纪初
非线性自适应滤波器开始受到关注,如神 经网络和模糊逻辑等。
随着数字信号处理技术的发展,自适应滤 波器的应用领域不断扩展,涉及通信、雷 达、图像处理、医学成像等多个领域。
稳态误差
自适应滤波器的稳态误差 越小,说明其跟踪期望信 号的能力越强。
鲁棒性
自适应滤波器的鲁棒性越 好,说明其对输入信号的 异常变化和噪声干扰的抵 抗能力越强。
03
自适应滤波器的实现方 法
最小均方误差算法
最小均方误差算法(LMS)是一种常用的自适应滤波算法,其基本思想是使滤波器 的输出信号与期望信号之间的均方误差最小。
。
05
自适应滤波器的挑战与 展望
自适应滤波器的稳定性问题
稳定性是自适应滤波器的核心 问题之一,它关系到滤波器的 性能和可靠性。
稳定性问题主要表现在系Байду номын сангаас参 数的变化和噪声的影响,可能 导致滤波器性能下降甚至失稳 。
解决稳定性问题的方法包括改 进算法、增加系统稳定性约束 条件等。
模电第六章 基于集成运算放大器的有源 滤波器分析与设计
相频响应
arctg
1 0 /
0 / Q
2
第六章
基于集成运算放大器的有源 滤波器分析与设计
三、二阶Sallen-Key带通滤波器
高通
反馈
设 Y 1 1/ R 1
Y2 1 R2 Y3 sC3 Y4 sC4 Y5 1 R5
得到二阶有源带通滤波电路
5、设计有源滤波器比设计LC滤波器更具灵活性,也可得到电 压增益。
第六章
基于集成运算放大器的有源 滤波器分析与设计
4.滤波器的用途 滤波器主要用来滤除信号中无用的频率成分,例 如,有一个较低频率的信号,其中包含一些较高频率 成分的干扰。滤波过程如图所示。
第六章
基于集成运算放大器的有源 滤波器分析与设计
第六章
基于集成运算放大器的有源 滤波器分析与设计
稳态响应
H ( j ) H (0 ) 1 jQ 0 0
幅频响应
H ( j ) H (0 ) 1 Q2 0 0
2
相频响应
arctgQ
低通
第六章
基于集成运算放大器的有源 滤波器分析与设计
A1 A0 通带 O 测评 通带 阻带 阻带
有源带通滤波电路可理解为
由低通和高通串联得到
1
1 低通特征角频率 1 R1C 1 1 高通特征角频率 2 R2 C 2
必须满足
A2 A0
阻 碍 阴
通带 阻 碍 测评 O 2 阴 阻 碍 A A0 阴 通带 阻带 O 阻 碍
低通(LPF) 高通(HPF) 带通(BPF) 带阻(BEF) 全通(APF)
第六章
基于集成运算放大器的有源 滤波器分析与设计
第6章FIR数字滤波器的设计
表6-1a 四种线性相位FIR滤波器的特性 类型 h(n) h(n)=h(N-1-n) N为奇数 h(n)=h(N-1-n) N为偶数
H ( )
( )
1型
关于 0, ,2 偶对称
( )
2型
关于 0,2 偶对称 关于 奇对称
N 1 2
第一类线性相位
H()
1 H ( ) d ( n) sin ( n ) 2 n 1 o N N 其中:d ( n) 2h( n) n 1,2,3, , 2 2 由此看出:
N /2
2
1 ()由于sin ( n ) 在 0,处为0, 1 2 2 即H ( )在 0,2处为零。即H ( z )在z 1处有一零点。 H ( )对 0,处呈奇对称,对 呈偶对称。 2 (2 )此类型不能用于设计 低通、带阻滤波器。
0
N 1 2
N 1 π
N/2 1 H () b(n) cos n 2 n1
N-1 n H() o
2
2型
情 况 2
b(n)
0
N 2
n
19
奇对称单位冲激响应
相位响应
h(n)=-h( N-1-n)
3型
情 况 3
7
H (e j ) sin 4e j 3 | sin 4 | e j ( )
1 0.5 0 -0.5 -1 1 0.5 0 -0.5 -1
0
0.5
1
0
0.5
1
0
1 0.5 0
-1
-2
-0.5 -1
-3
数字信号处理程佩青第三版课件第六章IIR滤波器的设计方法
6.2 最小与最大相位延时系统、最小 与最大相位超前系统
LSI系统的系统函数:
M
M
(1cmz1)
(zcm)
H(z)KmN1
Kz(NM)
m1 N
(1dkz1)
(zdk)
k1
ቤተ መጻሕፍቲ ባይዱk1
频率响应:
M
(ejcm)
H(ej)Kej(NM)m N 1
H(ej)ejarg[H(ej)]
(ejdk)
可k 整1理ppt
H (ej)H (ej)ej(j)
H ( e j ) 为幅频特性:表示信号通过该滤波器
后各频率成分的衰减情况
( j) 为相频特性:反映各频率成分通过滤波
器后在时间上的延时情况
可整理ppt
理想滤波器不可实现,只能以实际滤波器逼近
通带: c
11H(ej)1
阻带: st H(ej) 2
过渡带: cst
c :通带截止频率 s t :阻带截止频率
1 :通带容限
2 :阻带容限
可整理ppt
通带最大衰减: 1
12 0lgH H ((e ejj 0 c)) 2 0lgH (ej c) 2 0lg (11)
阻带最小衰减: 2
220lgH H ((e ej j0 st)) 20lgH (ejst) 20lg2
第六章 IIR滤波器的设计
可整理ppt
1
主要内容
理解数字滤波器的基本概念 了解最小相位延时系统 理解全通系统的特点及应用 掌握冲激响应不变法 掌握双线性变换法 掌握Butterworth、Chebyshev低通滤波器的特点 了解利用模拟滤波器设计IIR数字滤波器的设计过程 了解利用频带变换法设计各种类型数字滤波器的方法
第6章---锁相环路的应用
cos
2ct
mA 2
sin[(2c
)t
]
mA 2
sin[(2c
)t
]
(6-10) (6-11)
《锁相技术》
第6章 锁相环路的应用
《锁相技术》
图 6-7 AM信号的PLL同步解调
第6章 锁相环路的应用
《锁相技术》
图 6-7 AM信号的PLL同步解调
第6章 锁相环路的应用
二、模拟调频和调相信号的调制与解调
f p
1( 2
f3
f2)
1 (120.8 72.5) 2
24.15kHz
《锁相技术》
第6章 锁相环路的应用
二、频率特性 锁相环路对输入高频信号的带通特性是由环路传 递函数的低通特性所决定的。设输入信号被正弦音频 信号调频,则输入瞬时频率为
i (t) c sin t
(6-1)
式中ωc是载频;
第6章 锁相环路的应用
《锁相技术》
图 6-2 锁相环路跟踪特性的测量
第6章 锁相环路的应用
当输入频率下降时得到图中实线,在
fi=f3=1208kHz处环路捕获,在fi=f1=41kHz处失锁。由 此可算得环路的同步带
捕获带
f H
1( 2
f4
f1 )
1 (161.5 41)0 2
60.25kHz
U c {sin ct
mA 2
cos[(c
)t
]
mA 2
cos[(c
)t
]
(6-8)
《锁相技术》
第6章 锁相环路的应用
2调制器 用集成锁相环路很容易构成一个性能良 好的AM调制器。这时,环中的相乘器不再作鉴相器应 用,而是直接用它的相乘功能;压控振荡器也不再作 被控振荡器,而是直接产生载波信号。 由此构成如图 6-4框图。
数字信号处理--第6章无限脉冲响应数字滤波器的设计
1 0 0.1a p 1 k sp 1 0 0.1as 1 0 .0 2 4 2
sp
2 2
fs fp
2.4
N lg 0.0242 4.25, N 5 lg 2.4
2019/10/17
数字信号处理
(2) 按照(6.2.12)式,其极点为
j3
s0 e 5 ,
滤波器幅频特性。其幅度平方函数用A2(Ω)表示:
A2()Ha(j)2 12C1N 2( p)
(6.2.19)
2019/10/17
数字信号处理
2019/10/17
图6.2.5 切比雪夫Ⅰ型滤波器幅频特性
数字信号处理
式中,ε为小于1的正数,表示通带内幅度波动的 程度,ε愈大,波动幅度也愈大。Ωp称为通带截止频率。 令λ=Ω/Ωp,称为对Ωp的归一化频率。CN(x)称为N阶切 比雪夫多项式,定义为
数字信号处理
2019/10/17
数字信号处理
2019/10/17
数字信号处理
例6.2.1 已知通带截止频率fp=5kHz,通带最大衰减 αp=2dB , 阻 带 截 止 频 率 fs=12kHz , 阻 带 最 小 衰 减 αs=30dB,按照以上技术指标设计巴特沃斯低通滤波器。
解 (1) 确定阶数N。
(6.2.3) (6.2.4)
以上技术指标用图6.2.2表示。图中Ωc称为3dB截止 频率,因 H a (j c ) 1 /2 , 2 0 lg H a (j c ) 3 d B
2019/10/17
数字信号处理
2019/10/17
图6.2.2 低通滤波器的幅度特性
数字信号处理
滤波器的技术指标给定后,需要设计一个传输函
第6章无限冲激响应(IIR)滤波器设计
单位 (dB)
若在 ωp 处幅度下降到 0.707, 则幅平方下降 0.5 (半功率点): 若在 ωs 处幅度下 降到 0.01: | H(e j0 )| 1 αs = 20lg = 20lg = 40dB jω
| H(e
p
| H(e j 0 )|2 1 αp =10lg =10lg = 3dB jωp 2 0.5 | H(e )|
(b) | H(e jω ) |;
(c) | Y (e jω ) |
x(n)通 过系统 (n)后使 h 输出y(n)中 不再含有| ω |> ωc的 频率成 ,而使| ω |< ωc的成分 失真的 分 "不 通过".
2
2. 滤波器的分类
x(n) = s(n) + u(n)
若
加法性噪声
x(n) 中的有用成分 s(n) 和希望去除的成分
5
3. 滤波器的技术要求
低通:
δ2 :阻带容限;
决 定 于
δ1 : 通带容限;
ωp : 通带截止频率(又称通带上限频率) ωs : 阻带下限截止频率
αp :通带允许的最大衰减;
αs :阻带内应达到的最小衰减
6
3. 滤波器的技术要求
低通:
δ2 :阻带容限;
决 定 于
δ1 : 通带容限;
αp :通带
11
数字滤波器设计的一般步骤:
1. 给定所设计的滤波器的技术指标:
ωp , ωs , αp , αs , fs
LP, HP
ωsl , ωsh , ω1, ω3, αp , αs , fs 阻带起始频率
(下 带 止 (上 带 止 (通 下 截 )(通 上 截 ) 阻 截 ) 阻 截 ) 带 限 止 带 限 止
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
X
13
6.2.2 巴特沃思低通滤波器 Butterwoth
设计滤波器的一般工程方法: 设计滤波器的一般工程方法:
利用逼近理论寻找可实现的( s) =
∑a s
i =0 m i
bk s k ∑
k =0
∏ (s + z ) = ∏ (s + p )
i i k k
( ω)
α p = −20 lg H (ωc )
α s = −20 lg H (ωs )
11
X
6.2.1 概述
当系统的傅里叶变换存在时: 当系统的傅里叶变换存在时: H ∗ (ω ) = H (−ω ) 2 H (ω ) = H (ω ) H ∗ (ω ) H (ω ) = H ( s ) s = jω = H (ω ) H (−ω ) H (−ω ) = H (− s ) s = jω 2 ∴ H (ω ) = H ( s ) H (− s ) s = jω 意义:
2/2
|H(ω)| 过 渡 带 ωc ωs
通带
阻带 ω
一般设|H(ω)|的峰值等于 的峰值等于1 一般设 的峰值等于
7
X
6.1.3 滤波器的技术要求
截止频率ω 截止频率 c
1
1 2的频率
2/2
|H(ω)|
|H(ω)|下降 下降 -3dB频率 频率 一个或多个截止频率
过 渡 带 通带 ωc
20 lg(
所有极点以jω轴对称分布在巴特沃思圆上 所有极点以 轴对称分布在巴特沃思圆上 虚轴jω上无极点 n为偶数时,jω轴上无极点 为偶数时, 轴上无极点 为偶数时 n为奇数时,在s=±ωc的实轴上有两个极点 为奇数时, 为奇数时
18
极点6个 极点 个,间距 π/3
X
(3) 巴特沃思低通滤波器的传递函数
因滤波器必须是稳定的 取左s平面的全部极点作为 平面的全部极点作为H(ω)的极点 取左 平面的全部极点作为 的极点 平面是的极点属于H(-ω) 右s平面是的极点属于 平面是的极点属于 ωcn 故: H ( s ) = n ∏ (s − sk ) n为偶数时 为偶数时 H ( s) = n / 2
设计模拟滤波器的依据: 设计模拟滤波器的依据:
给定的工作损耗α 给定的工作损耗 p、αs~|H(ω)|2——幅度平方函数 频率选择性取决于传递函数|H(ω)|2 H(s) 频率选择性取决于传递函数 n H(s)必须是稳定的时不变系统 必须是稳定的时不变系统 i
实系数有理函数 分子阶数n≤分母阶数 分母阶数m 分子阶数 分母阶数 极点分布在左半s平面 极点分布在左半 平面 H(ω)有共轭对称性 有共轭对称性
保证S 保证 k=0时, 时 H(0)=1
ω
k =1 n c
∗ ( s − sk )( s − sk ) ∏ k =1
ω ωc
最大平坦幅值滤波器
15
X
(1) 巴特沃思低通滤波器的幅频特性
性能参数
衰减函数 通带衰减函数
1 α = −10 lg 2n 1 + ω ωc
ω 2 n = 10 lg 1 + ωc
y (t ) = x(t ) ∗ h(t )
ωc
y(t)=s(t) y(m)=s(m)
4
滤波器的截至频率
X
6.1.2 滤波器的分类
按照选频特性
低通滤波器 高通滤波器 带通滤波器 带阻滤波器 全通滤波器
|H(ω)| |H(ω)|
LP
ωc
|H(ω)|
HP
ω
|H(ω)|
ωc
ω
BP
ωc1 ωc2 ω
|H(ω)|
9
X
6.1.3 滤波器的技术要求
相移φ(ω) 相移
信号通过滤波器后 的相位滞后 jφ (ω) H(ω) = H(ω) e
1
2/2
|H(ω)|
过 渡 带 通带 ωc ωs 阻带 ω
相位滞后是频率的函数
群延迟τ 群延迟 g
d 相移对频率的导数(变化率 τ 变化率): 相移对频率的导数 变化率 :g = − φ(ω) dω
1 ) = −3dB 2
阻带 ωs ω
带宽(通带 带宽 通带)B 通带
|H(ω)|从1(0dB)下降到 从 下降到 带通滤波器上下截止频率之间的区域: = ωc2 −ωc1 带通滤波器上下截止频率之间的区域: B
1 的通频带宽度 2 (-3dB)的通频带宽度
中心频率ω 中心频率 0
带通滤波器上下截止频率的几何平均值 ω0 = ωc1 ⋅ωc2
1+ ( jωc )
1+ (
ω 2n ) ωc
1
令分母多项式=0, 得到极点: 令分母多项式= 得到极点:
1+ ( s jωc )2n = 0
j( 2 k −1 π π+ ) 2n 2
2n
2n
−1 = e
π
j(
2 k −1 π) 2n
⇒ s = jωc
sk = ωc e
2n个极点均匀分布 2n个极点均匀分布 间隔 n ,第1个极点 位置 π + π
2
滤波器
定 义: 滤除噪声、 滤除噪声、提取特征信息的系统 物理形式:
模拟滤波器(R、 模拟滤波器(R、L、C) (R 数字滤波器(软件算法) 数字滤波器(软件算法)
X
6.1.1 滤波器的基本原理
信号与噪声通常占据不同的频带 滤波器实质上是一种选频器件
使一种频率的信号分量(噪声 大幅度衰减 使一种频率的信号分量 噪声)大幅度衰减 噪声 使另一种频率的信号分量(信号 信号)顺利通过 使另一种频率的信号分量 信号 顺利通过 滤除噪声、 滤除噪声、获得有用信息
2n 2
17
−1
sk = ω
2k − 1 π ϕs = π+ c 2n 2 k = 1, 2,⋯ , 2n
X
(2) 巴特沃思低通滤波器的极点分布
巴特沃思圆
极点分布图
sk = ωc
2k − 1 π ϕk = π+ 2n 2 k = 1, 2,⋯ , 2n
2阶 阶
3阶 阶
极点4个 间距 极点 个,间距π/2
LP
ωc1 ωc2 ω
按滤波器元件性质
无源滤波器(R、 、 无源滤波器 、L、C) 有源滤波器(含运放 含运放) 有源滤波器 含运放
ω
5
X
6.1.3 滤波器的技术要求
理想滤波器与实际滤波器的区别
|H(ω)| |H(ω)|
过渡带 波纹度
物理上 不可实现
ωc
ω
ωc
ω
截至频率处突变 通带内幅频特性常数 阻带内幅频特性为0 阻带内幅频特性为
12
X
例题5-1 求给定滤波器的平方幅度函数的 (1 − ω 2 ) 2 最小相位滤波器的传递函数: (ω )2 = H 2
代替ω 解:以-s2代替 2,得
(1 + s ) H ( s) H (− s ) = (4 − s 2 )(9 − s 2 )
2 2
(4 + ω )(9 + ω 2 )
jω
◎j × -3 × -2 × × 0 2 3 ◎ -j
(s + j) (s − j) = ( s + 2)( s − 2)( s + 3)( s − 3)
2 2
σ
一对2重共轭零点: 一对 重共轭零点:±j 重共轭零点 两对实极点: 两对实极点:±2、 ±3 、 取左半s平面的极点 、-3 一对共轭零点±j 确定 取左半 平面的极点-2、 一对共轭零点± 平面的极点 滤波器传递函数H(s) 滤波器传递函数
通带波动∆ 通带波动
通带内最大值与最小值之差
8
X
6.1.3 滤波器的技术要求
衰减函数α 衰减函数
又称工作损耗
α = 20lg
H(0) H(ω)
2
1
2/2
|H(ω)|
过 渡 带 通带 ωc ωs 阻带 ω
= −20lg H(ω) = −10lg H(ω)
描述幅频特性的衰减程度 理想滤波器通带衰减= ,阻带衰减= 理想滤波器通带衰减=0,阻带衰减=∞ 实际滤波器衰减在0~ 之间 实际滤波器衰减在 ~∞之间 通带最大衰减α 阻带最小衰减α 通带最大衰减 p、阻带最小衰减 s αp = −20lg H(ωc ) αs = −20lg H(ωs )
H (ω ) = H (ω ) e jϕh (ω )
Y(ω) X(ω) S(ω) N(ω)
离散系统
y ( n) = x ( n) ∗ h( n)
Y (Ω) = H (Ω ) X (Ω)
H (Ω) = H (Ω) e jϕh ( Ω )
设计滤波器即设计合适的H(ω),满足滤波效果
x(t)= s(t) + n(t) 或 x(n)=s(m)+n(m) h(t)或h(n)
ω 2 n s 阻带衰减函数 α |ω =ωs = 10 lg 1 + ≥ αs ωc 确定滤波器的阶数n 确定滤波器的阶数
n≥ lg 10
0.1α s
ω 2 n α p = α |ω =ωc = 10 lg 1 + c = 3dB ω c
信号分析与处理 Signals analysis & processing
第6章 滤波器
华侨大学机电及自动化学院
6.1 滤波器概述