数字信号处理第6章概论
数字信号处理-第六章

第六章数字滤波器结构6.1:级联的实现num = input('分子系数向量 = ');den = input('分母系数向量 = ');[z,p,k] = tf2zp(num,den);sos = zp2sos(z,p,k)Q6.1使用程序P6.1,生成如下有限冲激响应传输函数的一个级联实现:H1(z)=2+10z^(-1)+23z^(-2)+34z^(-3)+31z^(-4)+16 z^(-5)+4z^(-6)画出级联实现的框图。
H1(z)是一个线性相位传输函数吗?答:运行结果:sos = zp2sos(z,p,k)Numerator coefficient vector = [2,10,23,34,31,16,4]Denominator coefficient vector = [1]sos =2.0000 6.0000 4.0000 1.0000 0 01.0000 1.00002.0000 1.0000 0 01.0000 1.0000 0.5000 1.0000 0 0级联框图:H1(z)不是一个线性相位传输函数,因为系数不对称。
Q6.2使用程序P6.1,生成如下有限冲激响应传输函数的一个级联实现:H2(z)=6+31z^(-1)+74z^(-2)+102z^(-3)+74z^(-4)+31 z^(-5)+6z^(-6)画出级联实现的框图。
H2(z)是一个线性相位传输函数吗?只用4个乘法器生成H2(z)的一级联实现。
显示新的级联结构的框图。
Numerator coefficient vector = [6,31,74,102,74,31,6]Denominator coefficient vector = [1]sos =6.0000 15.0000 6.0000 1.0000 0 01.00002.00003.0000 1.0000 0 01.0000 0.6667 0.3333 1.0000 0 0级联框图:H2(z)是一个线性相位传输函数。
《数字信号处理导论_第6章》

每一种又有模拟(AF)、数字(DF)两种滤波器.
对数字滤波器, 从实现方法上, 有IIR滤波器和
FIR滤波器之分, 转移函数分别为: FIR DF:
直接设计
IIR DF:
利用模拟滤 波器设计
3. 滤波器的技术要求
低通:
:通带允许的最大衰减;
:阻带内应达到的最小衰减
单位 (dB)
若幅度下降到 0.707, 则幅平方下降 0.5 (半 功率点):
k 1,2,..., N
5)滤波器的设计步骤:
确定技术指标: p p s s
注: p c
根据技术指标求出滤波器阶数N及 :
2 10
1
0.1 p
1
1 ch 100.1 s 1 N 1 s ch c
k 1, 2,..., N
c cr 1 rad / s
为归一化系统的系统函数 H an ( s) 去归一化,得 H ( s ) H ( s ) a an
cr s H an cr s s c c
4)滤波器的设计步骤:
确定技术指标: p p s s
由H a ( s)
s 0
H a ( j)
0
,得K0 4
4( s 2 25) 4s 2 100 H a (s) 2 ( s 7)( s 6) s 13s 42
2、Butterworth 低通逼近
幅度平方函数: H a ( j)
2
1 1 c
• 极点在s平面呈象限对称,分布在Buttterworth圆上,共2N点 • 极点间的角度间隔为 / N rad • 极点不落在虚轴上 • N为奇数,实轴上有极点,N为偶数,实轴上无极点
数字信号处理_第六章

j2
j2
s p 3 c e 3,s p 4 c ,s p 5 c e3,k 0 ,1 ,2
系统函数为: H a(s)(ssp3)(s sc 3p4)(ssp5)
或
1 H a (s) (s/ c)3 2 (s/ c)2 2 (s/ c) 1
令
p
s c
,则
pk sk /c
有,归一化的三阶滤波器的系统函数
22 0 lgH H ((e e j j0 st)) 2 0 lgH (ej st) 2 0 lg2
其中: H(ej0) 1
当 H (ej c)2/20 .7 0 7时,1 3dB
称 c 为3dB通带截止频率
3.滤波器的设计方法
直接设计法
在时域或频域直接设计数字滤波器
间接设计法
先根据指标要求设计对应的模拟滤波器 再将模拟滤波器转换为数字滤波器
通带: c
11H(ej)1
阻带: st H(ej) 2
过渡带: cst
c :通带截止频率
s t :阻带截止频率
:通带容限
1
2 :阻带容限
通带最大衰减: 1
1 2 0 lgH H ( (e e jj 0 c ) ) 2 0 lgH (ej c) 2 0 lg (1 1 )
阻带最小衰减: 2
数字滤波器分为: 无限长单位脉冲响应滤波器(IIR DF) [递归系统] 有限长单位脉冲响应滤波器(FIR DF)[非递归系统]
M
bj z r
H (z) j0 N 1 ak z k k 1
N1
H(z) h(n)zn n0
(6.1.1) (6.1.2)
2. 数字滤波器的技术指标
数字滤波器的频率响应函数H(ejω)用下式表示: H(ejω)=|H(ejω)|ejθ(ω)
数字信号处理第三版 教材第六章习题解答

6.2 教材第六章习题解答1. 设计一个巴特沃斯低通滤波器,要求通带截止频率6p f kHz =,通带最大衰减3p a dB =,阻带截止频率12s f kHz =,阻带最小衰减3s a dB =。
求出滤波器归一化传输函数()a H p 以及实际的()a H s 。
解:(1)求阶数N 。
lg lg sp spk N λ=-0.10.30.1 2.51011010.0562101101p s asp a k --==≈--332121022610s sp p πλπΩ⨯⨯===Ω⨯⨯将sp k 和sp λ值代入N 的计算公式得lg 0.05624.15lg 2N =-=所以取N=5(实际应用中,根据具体要求,也可能取N=4,指标稍微差一点,但阶数低一阶,使系统实现电路得到简化。
) (2)求归一化系统函数()a H p ,由阶数N=5直接查表得到5阶巴特沃斯归一化低通滤波器系统函数()a H p 为54321() 3.2361 5.2361 5.2361 3.23611a H p p p p p p =+++++或 221()(0.6181)( 1.6181)(1)a H p p p p p p =+++++ 当然,也可以按(6.12)式计算出极点:121()22,0,1,2,3,4k j Nk p ek π++==按(6.11)式写出()a H p 表达式41()()a k k H p p p ==-代入k p 值并进行分母展开得到与查表相同的结果。
(3)去归一化(即LP-LP 频率变换),由归一化系统函数()a H p 得到实际滤波器系统函数()a H s 。
由于本题中3p a dB =,即32610/c p rad s πΩ=Ω=⨯⨯,因此()()a a cH s H p s p ==Ω5542332453.2361 5.2361 5.2361 3.2361c c c cc cs s ss s Ω=+Ω+Ω+Ω+Ω+Ω对分母因式形式,则有()()a a cH s H p s p ==Ω52222(0.6180)( 1.6180)()c c c c cc s s s s s Ω=+Ω-Ω+Ω-Ω+Ω如上结果中,c Ω的值未代入相乘,这样使读者能清楚地看到去归一化后,3dB 截止频率对归一化系统函数的改变作用。
现代数字信号处理概论_6

■ 现代滤波器 ,则不是建立在频率领域 ,而是通过随机 过程的数学手段 ,通过对噪声和信号的统计特性做一 定的假定 ,然后通过合适的数学方式 ,来提高信噪比。 譬如Kalman滤波器中 , 总会假定状态噪声和测量噪声 是不相关的; 而在Weiner滤波器中还必须假定信号是 平稳的 ,等 。总之各有所用,要针对不同的问题采用 不同滤波器 。譬如,要滤除工频50HZ的影响,显然不 宜采用Kalman滤波器 ,可以采用限波器就可以了。
■ 金年文 , 韦岗 ,现代数字信号处 理简明教程 ,清华大学出版社 , 2004 年1月
教材选择
■ 个人可以选择所列参考书目中任意一本 作为教材 , 都是经典教材 , 本人所选为 何子述的教材 , 利于讨论
■ 课件为主 , 教材为辅 ■ 考核中涉及到的均从相应配套习题集中
选择(姚天任 , 何子述 , 丁玉美等)
■ 20世纪40年代 , 维纳确定了关于最佳滤波器研究的基础 。 即假定 线性滤波器的输入为有用信号和噪声之和 , 两者均为广义平稳过 程 ,且知它们的二阶统计特性 。维纳根据最小均方误差准则(滤 波器的输出信号与需要信号之差的均方值最小),求得了最佳线 性滤波器参数 , 这种滤波器被称为维纳滤波器 。在维纳研究的基 础上 , 人们还根据最大输出信噪比准则 、统计检测准则以及其他 最佳准则求得的最佳线性滤波器 。实际上 ,在一定条件下 , 这些 最佳滤波器与维纳滤波器是等价的 。 因而 ,讨论线性滤波器时, 一般均以维纳滤波器作为参考。
些序列 ■ 目的
■ 估计信号的特征参数(脑电图和心电图分析 ,或语 音传输分析和语音识别系统中 )
数字信号处理ppt第六章

一、DF按频率特性分类 可分为低通、高通、带通、带阻和全通,
其特点为:
(1)频率变量以数字频率 ω 表示,ω = ΩT ,
Ω 为模拟角频率,T为抽样时间间隔; (2)以数字抽样频率 ωs = 2πfs ⋅T = 2π 为周期; (3)频率特性只限于 ω ≤ ω s / 2 = π 范围,这
3、由 A2 (Ω) = H a ( jΩ) 2 确定 H a (s)的方法
(1)求 H a (s)H a (−s) = A2 (Ω) Ω2 =−S 2
(2)分解 Ha (S)Ha (−S),得到各零极点,将左半面的 极点 归于 Ha (S),对称的零点任一半归 Ha (S)。若要求 最小相位延时,左半面的零点归 Ha (S)(全部零极点 位于单位圆内)。
将2、技Q∴计术算2H指0所a标l(g需j,ΩH的代)a阶2入( j=数Ω上1及式)/[3=1,d+B−可截(1得Ω0Ω止lC频g)[21率N+]Ω(CΩΩC )2N ]
{−10lg[1+ ( 2π×103 )2N ] ≥ −1 −10lg[1+ (3π×Ω1C03 )2N ] ≤ −15 ΩC
解上述两式得:
它是表示每个频率分量的延迟情况;当其为常数时, 就是表示每个频率分量的延迟相同。
四、DF设计内容 1、按任务要求确定Filter的性能指标; 2、用IIR或FIR系统函数去逼近这一性能要求; 3、选择适当的运算结构实现这个系统函数; 4、用软件还是用硬件实现。
五、IIR数字filter的设计方法
1、借助模拟filter的设计方法 (1)将DF的技术指标转换成AF的技术指标; (2)按转换后技术指标、设计模拟低通filter的 Ha (s); (3)将 H a (s) → H (z) (4)如果不是低通,则必须先将其转换成低通
精品课件-数字信号处理—理论与实践-第6章

第 6 章 数字滤波器的结构
因而在设计具体的实现算法时要分析和考虑选择什么样的网 络结构才合适。
一般来说, 设计好数字滤波器的结构后, 我们就可以通过 两种方法来具体实现数字滤波器:
(1) 将数字滤波器所要完成的运算编成程序, 利用计算 机进行软件实现;
(2) 设计专用的数字硬件、 专用的数字信号处理器或采 用通用数字信号处理器(DSP)进行硬件实现。
y(n)=a1y(n-1)+a2y(n-2)+b0x(n) 它对应的方框图结构如图6-3所示。
第 6 章 数字滤波器的结构
图6-2 基本运算的方框图表示法
第 6 章 数字滤波器的结构
图6-3 二阶数字滤波器的方框图结构
第 6 章 数字滤波器的结构
2. 信号流图法的特点是简单、 方便。 和方框图法相对应,
三种基本运算的信号流图表示如图 6-4 所示。
图6-4 基本运算的信号流图表示法
第 6 章 数字滤波器的结构
信号流图在本质上与方框图表示法等效, 只是符号上有差 异。 比如, 图6-3的二阶数字滤波器用信号流图表示的结 构如图6-5所示。 图中, 1, 2, 3, 4, 5称为网络节点, x(n)处为输入节点或称源节点, y(n)处为输出节点或称阱节点。 节点之间用有向支路相连接, 支路上的传输系数如果为常数, 则表示乘法运算; 如果没有标注传输系数, 则表示其传输系数 为1; 如果是延时算子z-1, 则表示单位延时。
第 6 章 数字滤波器的结构
图6-5 图6-3的二阶数字滤波器的信号流图结构
第 6 章 数字滤波器的结构
源节点没有输入支路, 阱节点没有输出支路, 其余网络节 点均可以有多条输入支路和多条输出支路。 每一个节点的节点 值都等于它的所有输入支路的信号之和。 这样, 通过分析各节 点的值, 就可以清楚地得到该网络的传输特性。 比如图6-5所 表示的二阶数字滤波器的各节点的值为
数字信号处理第六章

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、用软件还是用硬件实现。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2c
Sa
c
(t
t0
) cos
0t
这是以等效低通滤波器的冲激响应为包络的正弦调幅信号。
16
6.2 模拟滤波器的基本概念与设计方法
6.2.1 系统的物理可实现性
从时域上看,一个物理可实现系统的冲激响应h(t)应满足:
传输,而将频率高于c 的信号完全抑制。
H ( j)
H ( j)
e j ()
e
jt0
c
0 c
或
e jt0 [u( c ) u( c )]
H( j)
1
H ( j)
1
c
0 c
c 0
c
()
u( c ) u( c )
() t0
c --- 3dB 截止频率
0
t0
5
6.1 信号的传输与滤波
Y ( j ) Ke jt0 X ( j ) H ( j ) X ( j )
H ( j ) Ke jt0 (6.1 2)
H( j)
K
即
H ( j ) K ( ) t0
(6.1 3)
( )
t0
3
6.1 信号的传输与滤波
无失真传输系统应满足如下两个条件: (1)系统的幅频特性在整个频率范围内为常数;
(2)系统的相频特性在整个频率范围内应与 成正比变化。
6.1.2 理想滤波器
理想滤波器:在通带(pass-band)内,滤波器的幅频 特性为常数,相频特性呈线性;而在阻带(stop-band)内, 滤波器的幅频特性立即降1、理想低滤波器
理想低通滤波器是将频率低于 c 的所有信号予以无失真地
g(t) c c (t-t0 ) sin x 1 dx
-
x c
1 0 sin xdx 1 c (tt0 ) sin xdx
- x
0
x
7
6.1 信号的传输与滤波
g(t) 1 0 sin xdx 1 c (tt0 ) sin xdx
- x
0
x
1 1 c (tt0 ) sin xdx
0
t0
t
h(t ) 的特点:响应超前于激励(非因果系统)
6
6.1 信号的传输与滤波
(2) 阶跃响应
h(t)
c
Sa[c
(t
t0
)]
g(t) t h( )d c t sin[c ( t0 )]d
- c ( t0 )
令 x c ( t0 ), 则
d
1 c
dx,
当
时,x ,
当
t
时,
x c (t t0 )
x1(t)
1
0
t
Y1( j)
X1( j)
0
2 / 2 / 4 /
1 H( j)
0 2 / 2 / 4 /
c
0 c 4 /
14
6.1 信号的传输与滤波
2、理想带通滤波器
H LP
(
j
)
e 0
j t0
c c
HBP ( j ) HLP ( j ) ( 0 ) ( 0 )
15
c (t t0 )
t t0 / c
1 1 Si() 1.0895 2
10
6.1 信号的传输与滤波
(3) 矩形脉冲响应
x1(t) u(t) u(t ) 1
t
0
y1(t)
1
1
2
c c
0
t0
tr
c 6 t
t0
y1(t)
1
{Si[c
(t
t0
)]
Si[c
(t
t0
)]}
矩形脉冲响应除了比矩形脉冲输入延迟一段时间t0外, 矩形脉冲响应的波形也不再是矩形脉冲,即产生了失真。
2 0
x
sin x x
1
y sinxdx Si( y)
0x
--------- 正弦积分
2
2
0
x
Si( y)
/2
g(t)
1 2
1
Si[c
(t
t0
)]
0 2
y
/ 2
8
6.1 信号的传输与滤波
g(t)
1 2
1
Si[c
(t
t0
)]
h(t)
c
t0 c
Si( y)
/2
0 2
y
/ 2
0
t0
1. 时域条件
y(t) Kx(t t0 ) (6.11)
其中: K ------- 常数, t0 ------- 滞后时间
x(t) 1
y(t) x(t) 线性系统 y(t) K
t
t0
t
2
6.1 信号的传输与滤波
2. 频域条件
y(t) Kx(t t0 ) (6.11)
对式(6.1.1)两边取傅氏变换,得:
g (t )
t
g(t)的特点:
1
a) 响应波形的前沿是倾斜的,
1 2
0
c c
t0 tr
响应信号的建立需要一段时间。
t
tr
2
c
1 fc
1 B
-------阶跃响应的建立
(上升)时间
9
6.1 信号的传输与滤波
g (t ) 1
1 2
0
c c
t0
tr
tr
2
c
1 fc
1 B
b)理想低通滤波器阶跃响应 的建立(上升)时间与滤波
(1) 冲激响应
H ( j) H (
j)
e j ()
e jt0
c
0 c
h(t) F 1[H ( j)]
x(t)
1 H ( j)e jtd
2
(1) 0 h(t)
t
1 c e jt0 e jt d
2 c
c
t0 c
c
Sa[c
(t
t0
)]
0h(t)
t0
t
(1)
当 c 时,h(t) (t t0 )
6.1 信号的传输与滤波
h(t) F 1 HBP ( j ) 2F 1 HLP ( j )F 1 ( 0 ) ( 0 )
F
1
H LP
(
j
)
c
Sa
c
(t
t0
)
F 1
( 0 ) ( 0 )
1 2
ej 0t
e j 0t
1
cos
0t
h(t)
2
c
Sa
c
(t
t0
)
1
cos
0t
器的截止频率成反比。
c)响应与激励相比有波纹。
最大波峰的高度约为跳 变值的8.95%左右(波峰值 t 为1.0895),它与 c 无关。
-------- 吉伯斯现象
阶跃响应g(t)的第一个极大值发生
在 t t0 / c 处,将它代入到式
(6.1-8)中,得到阶跃响应的极大值:
gmax
1 1 Si 2
11
6.1 信号的传输与滤波
a)失真的程度既与理想低通滤波器的频带宽度有关, 也与矩形脉冲的频带宽度或脉冲宽度有关。
为矩形脉冲的频带宽度
12
6.1 信号的传输与滤波
b)不同 c的理想低通滤波器对矩形脉冲的响应
13
6.1 信号的传输与滤波
在下述情况下,会产生吉伯斯现象: 1)激励有跳变; 2)系统的带宽为有限值。
第6章 傅里叶变换的应用
6.1 信号的传输与滤波 6.2 模拟滤波器的基本概念与设计方法 6.3 信号的采样 6.4 调制与解调 6.5 MATLAB在信息处理与通信中的应用
1
6.1 信号的传输与滤波
6.1.1 无失真传输
信号无失真传输是指响应信号与激励信号相比,只有幅度 大小和出现时间的不同,而没有波形上的变化。