数字信号处理第五章经典讲解
胡广书_数字信号处理题解及电子课件_第5章
π 3. 在实轴但不在圆上,无共轭,角度=0, π 4. 在实轴,但在圆上,无共轭,角度=0,
模<1; 模=1;
在单位圆内 四个零点同时存在, 构成四阶系统.
把该式展开,其系数也是对称的,是具 有线性相位的子系统。
无共轭零点, 有镜象零点
无镜象对称零点, 有共轭零点.
无镜象零点, 也无共轭零点.
− jω n
+
n = ( N +1) / 2
∑
N −1
h ( n )e
− jω n
N −1 + h e 2
N −1 − j ω 2
令:
m = N −1 − n
( N −3) 2 ( N −3) 2 m=0
H (e jω ) 的对称性,有 并利用
H (e jω ) =
则该系统具有线性相位。
上述对称有四种情况:
偶对称
N : even N : odd
第一类 FIR 系统
N : even N : odd
奇对称
第二类 FIR 系统
1. N 为奇数
H (e jω ) = ∑ h(n)e − jω n
n =0 N −1
第一类 FIR 系统
=
( N −3) 2
∑
n=0
h(n)e
幅频
一 阶 全 通 系 统
0 -0.5 -1 -1 0 -1 0 (a) 1
0
0.2 (b)
0.4
1
0.5 -2 0 -3 -4 -0.5
0
0.2 (c)
0.4
0
5
10 (d)
15
20
相频
抽样响应
1 1 0.5 0 -0.5 -1 -1 0 -2 0.5 -4 0 -6 -8 -0.5 0 (a) 1 0 0.5
《数字信号处理教程》程佩青课件 第五章 数字滤波器的基本结构
a12 Z-1 β12
y(n) a1MZ-1 β1M
a21 Z-1 β21
a22 Z-1 β22
…... a2MZ-1 β2M
例子
设IIR数字滤波器系统函数为:
H(z) 1 2z 1 2z 2 z 3 (1 z 1 )(1 z 1 z 2 )
1 2z 1 z 3
(1 z 1 )(1 z 1 z 2 )
则:
M1
M2
(1 gi z 1 ) (1 1i z 1 ,i z 2 )
H (z)
A
i 1 N1
i 1 N2
(1 pi z 1 ) (1 1i z 1 2i z 2 )
i 1
i 1
(3)基本二阶节的级联结构
M1
M2
(1 gi z 1 ) (1 1i z 1 2i z 2 )
Hi
(z)
1 1i 1 1i
z 1 z 1
2i z 2 2i z 2
一般用直接II型(典范型表示)
x(n)
y(n)
a1i Z-1 β1i a2i Z-1 β2i
(5)用二阶节级联表示的滤波器系统
整个滤波器则是多个二阶节级联
N 1 2
H (z) A Hi (z) i 1
x(n)
A
a11 Z-1 β11
a2 a N-1 aN
Z-1 Z-1 b2
Z-1 Z-1 b M+1
Z-1 Z-1
bM
a1 Z-1 b1 a2Z-1 b2
a N-Z1 -1 b M+1 aNZ-1 bM
合并
这就是直接II型的结构流图。
(4)直接II型特点
直接II型结构特点: (1)两个网络级联。 第一个有反馈的N节延时网络实现极点; 第二个横向结构M节延时网络实现零点。 (2)实现N阶滤波器(一般N>=M)只需N级延时 单元,所需延时单元最少。故称典范型。
数字信号处理-第五章
系 统 函 数
:
H (z) n M 0h (n )z n Y X ( (z z ) ) b 0 b 1 z 1 b 2 z 2 b M z M
单位脉冲响应的值等于差分方程系数:
h
h(n)=bn
n=0,1,·····,M
33
FIR数字滤波器的特点:
系统函数:
N1
H(z) h(n)zn n0
有N-1个零点分布于z平面 z=0处 是N-1阶极点
h
26
还可以如下式这样进行分解: H (z)1 1 0 0..6 4z z 1 11 1 0 0..5 3z z 1 1H 3(z)H 4(z)
h
27
级联型结构的特点:
调整某一路的分子系数能单独调整滤波器的一组 零点,而不影响其它零极点; 调整某一路的分母系数能单独调整滤波器的一组 极点,而不影响其它零极点;便于调整滤波器频 率响应性能
直接型
h
38
将H(z)进行因式分解,得到: H(z)=(0.6+0.5z-1)(1.6+2z-1+3z-2) 按照上式画出它的级联型结构如图所示。
级联型
h
39
5.5 线性相位网络结构
FIR滤波器单位抽样响应h(n)为实数,0 n N 1 且满足:
偶对称: h (n ) h (N 1 n ) 或奇对称 h (n ) h (N 1 n ) : 即对称中心在 (N-1) / 2处 则这种FIR滤波器具有严格线性相位。
bi zi
i0 N 1 ak zk
k 1
基本运算:加法,乘法(乘以常数),移位(时延)
h
3
信号流图由基本支路构成
1.基本支路箭头表示信号流向,两个圆点表示输入输出节点,箭头旁边的 符号表示增益(缺省为1)
第五章 时域离散系统的基本网络结构
本章的主要内容就是描述数字滤波器的基 本网络结构。(IIR、FIR)
引言
时域离散系统或网络可以用差分方程、单 位脉冲响应以及系统函数进行描述。
M
N
y(n) bi x(n i) ai y(n i)
i0
i 1
系统函数H(z)为
M
H (z)
(2) 流图环路中必须存在延时支路;
(3) 节点和支路的数目是有限的。
信号流图表达的系统含义
每个节点连接的有输入支路和输出支路,节点变 量等于所有输入支路的输出之和.
根据信号流图可以求出系统函数(节点法、梅逊 公式法)。
1(n) 2 (n 1) 2 (n) 2 (n 1) 2 (n) x(n) a12 (n) a21n y(n) b21(n) b12 (n) b02(n)
画出H(z)的直接型结构和级联型结构。
级联型
解: 将H(z)进行因式分解,得到: H(z)=(0.6+0.5z-1)(1.6+2z-1+3z-2)
其直接型结构和级联型结构如图所示。
x(n)
0.6
z- 1 0.5
1.6 z- 1
2 z- 1
3
y(n) x(n)
z- 1
z- 1
z- 1
0.96 2
2.8 1.5 y(n)
0 j
y(n)
1 j
z- 1 1j
1 j
z- 11 j
(a)
2 j
z-
1
2
j
(b)
一阶和二阶直接型网络结构 (a)直接型一阶网络结构;(b)直接型二阶网络结构
IIR的级联型例题
第五章 数字信号处理- 微弱信号处理
第五章微弱信号处理5.1 微弱信号检测技术中气体浓度检测仪中的应用微弱信号不仅意味着信号的幅度小,而且主要指被噪声淹没中的信号。
为了检测被背景噪音淹没的信号,就需要分析噪声产生的原因和规律,研究被测信号的特点、相关性以及噪声的统计特性,以寻找出从背景噪声中检测有用信号的方法。
因此,微弱信号检测技术的首要任务是提高信噪比。
它不同于一般的检测技术,它注重的不是传感器的物理模型和传感原理、相应的信号转换电路和仪表实现方法,而是如何抑制噪声和提高信噪比。
由于被测量的信号微弱,传感器的固有噪声、放大电路及测量仪器的固有噪声以及外界的干扰噪声往往比有用信号的幅度大得多,放大被测信号的过程同时也放大了噪声,而且必然还会附加一些额外的噪声,因此只靠放大是不能把微弱信号检测出来的。
只有在有效地抑制噪声的条件下增大微弱信号的幅度,才能提取出有用的信号。
为了表征噪声对信号的淹没程度,引入信噪比SNR来表示,它指的是信号的有效值S与噪音的有效值N之比。
而评价一种微弱信号检测方法的优劣,经常采用两种指标:一种是信噪改善比SNIR,它等于系统输出端的信噪比SNR和系统输入段oSNR之比。
SNIR越大,表明系统抑制噪声的能力越强。
i另一个指标是检测分辨率,指的是检测仪器指示值可以响应与分辨的最小输入值的变化值。
检测分辨率不同于检测灵敏度,后者表示的是检测系统标定曲线的斜率,定义为输出变化量y∆之比。
一般情况下,∆的输入变化量x∆与引起y灵敏度越高,分辨率越好。
但提高系统的放大倍数虽可提高灵敏度,但却不一定能提高分辨率,因为分辨率要受噪声和误差额制约。
5.1.1 本检测系统的噪声源广义的噪声是扣除被测信号真实值以后的各种测量值,可以分为两类:一是干扰;另一被称为电子噪声(狭义)。
干扰是指被非被测信号或非测量系统所引起的噪声。
从理论上讲,干扰是属于理想上可排除的噪声。
不少干扰源发出的干扰是有规律的,有些具有周期性,有些只是瞬时值。
《信号、系统与数字信号处理》第五章 Z变换与离散系统的频域分析
同理
sinh0nun
1 2
e0n
e0n
un
1 z
2
z
e0
z z e0
z2
z sinh0 2z cosh0
1
z max e0 , e0
2、双边z变换的移位 n0 0
若 xn X z
RX
z
R X
则 x n n0 z n0 X z
RX
z
R X
证明: Z x n n0
n
xT t nT estdt
n
xnT esnT
n
令 z esT 引入新的复变量, 将上式写为
X s s xnT zn
n
此式是复变量 z 的函数(T 是常数),记为
X z xnzn
n
x 2z2 x 1z x0 x1z1 x2z2
Z xn 2un z2 X z z1x1 x 2
3) 若 xn 为因果序列 xnun X z
则 xn mun zm X z
m0
xn
mun
zm
X
z
m1 k 0
xk
z
k
例5-9 求周期序列的单边z变换
解: 周期序列 xn xn rN
m0
令 n 0 ~ N 1 的主值区序列为 x1 n ,
( z 1)
4、指数序列加权
若 xn X z RX z RX
则 an xn X a1z
RX a 1z RX
证:Z an xn an xnzn
n
xn a1z n X z / a
n
RX a 1z RX
a
R X
z
a
R X
利用
数字信号处理第五章-IIR数字滤波器的设计
2、由模平方函数确定系统函数
模拟滤波器幅度响应常用幅度平方函数表示:
| H ( j) |2 H ( j)H *( j)
由于冲击响应h(t)为实函数,H ( j) H *( j)
| H ( j) |2 H ( j)H ( j) H (s)H (s) |s j
H (s)是模拟滤波器的系统函数,是s的有理分式;
分别对应:通带波纹和阻带衰减(阻带波纹)
(4种函数)
只介绍前两种
31
32
33
无论N多大,所 有特性曲线均通 过该点
特性曲线单调减小,N越大,减小越慢 阻
特性曲线单调减小,N越大,减小越快
34
20Nlog2:频率增加一倍,衰减6NdB
35
另外:
36
无论N多大,所 有特性曲线均通 过Ωc点: 衰减3dB, Ωc 为 3dB带宽
8
根据
(线性相位滤波器)
非线性相位滤波器
9
问题:
理想滤波器的幅度特性中,频带之间存 在突变,单位冲击响应是非因果的;
只能用逼近的方法来尽量接近实际的要 求。
滤波器的性能要求以频率响应的幅度特 性的允许误差来表征,如下图:
10
p
11
低通滤波器的频率响应包括:
通带:在通带内,以幅度响应的误差δp逼近 于1;
20
3、数字滤波器设计的基本方法
利用模拟理论进行设计 先按照给定的技术指标设计出模拟滤波 器的系统函数H(s),然后经过一定的变 换得到数字滤波器的系统函数H(z),这实 际上是S平面到Z平面的映射过程: 从时域出发,脉冲响应不变法 从频域出发,双线性变换法 适合于设计幅度特性较规则的滤波器, 如低通、高通等。
由于系统稳定, H(s)的极点一定落在s的左半 平面,所以左半平面的极点一定属于H(s),右 半平面的极点一定属于H(-s)。
精品文档-数字信号处理(吴瑛)-第5章
第5章 数字滤波器概论
5.3 实际滤波器的设计指标
5.3.1 图5.3.1是理想低通滤波器的幅频响应,该理想低通滤波
器具有截止频率ωd。可以看出,理想滤波器在通带内幅度为常 数(非零),在阻带内幅度为零。另外,一般理想滤波器 要求具有线性相位(在第8章讨论),这里假设相频响应 θ(ω)=0
h(n) sin(nd )
第5章 数字滤波器概论
1. 根据H(ejω) 一般数字滤波器从滤波功能上分类,和模拟滤波器一样, 可以分成低通、高通、带通和带阻等滤波器。它们的理想幅频 响应如图5.2.2
第5章 数字滤波器概论
图5.2.2 (a) 低通; (b) 高通; (c) 带通; (d) 带阻
第5章 数字滤波器概论
需要注意的是,数字滤波器的频率响应H(ejω)都是以2π 为周期的,滤波器的低通频带处于2π的整数倍处,而高通频 带处于π
5.3.2 当滤波器形状为非理想时,要用一些参数指标来描述其关
键特性。图5.3.5表示低通滤波器的幅频响应。滤波器的通带 定义了滤波器允许通过的频率范围。在阻带内,滤波器对 信号严重衰减。ωp和ωs分别称为通带截止频率(或通带上限频 率)和阻带截止频率(或阻带下限频率)。参数δ1定义了通带波 纹(Pass Band Ripple),即滤波器通带内偏离单位增 益的最大值。参数δ2定义了阻带波纹(Stop Band Ripple),即 滤波器阻带内偏离零增益的最大值。
截短脉冲响应自然会对频率响应产生影响。截短后,滤波 器幅频响应曲线不再是理想矩形,通带不再平坦,有过渡带, 同时阻带衰减不再为零。图5.3.4给出了因果脉冲响应 的幅频响应。当然,脉冲响应保留的采样点越多,即滤波器阶
第5章 数字滤波器概论 图5.3.4 非理想低通滤波器因果脉冲响应的幅频响应
数字信号处理 第五章07
an example
2
2
1.5
h(n)
1.5
a(n)
1
1
0.5
0.5
0
0
-0.5
-0.5
-1
-1
-1.5
-1.5
-2
0
2
4
6
-2
0
2
4
6
4 2 0 -2 -4
H()
0 -5 -10 -15 -20
()
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
0
0.2
0.4
0.6
0.8
h(n)
2 2
d(n)
1
1
0
0
-1
-1
-2
-2
-3
0
2
4
6
-3
0
2
4
6
5
0
H()
0
-5 -10 -15
()
-5
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
-20
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
5
4 2
arg[H(ej)]
0
0
|H(ej)|
-5 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
h(n ) e
j n
e
j (N n 1)
数字信号处理 第五章
第五章离散时间信号的数字处理Q5.1运行程序P5.1,产生连续时间序号及其抽样形式,并显示它们。
clf;t = 0:0.0005:1;f = 13;xa = cos(2*pi*f*t);subplot(2,1,1)plot(t,xa);gridxlabel('时间, msec');ylabel('振幅');title('连续时间序号 x_{a}(t)');axis([0 1 -1.2 1.2])subplot(2,1,2);T = 0.1;n = 0:T:1;xs = cos(2*pi*f*n);k = 0:length(n)-1;stem(k,xs);grid;xlabel('时间 n');ylabel('振幅');title('离散事件序号 x[n]');axis([0 (length(n)-1) -1.2 1.2])Q5.2 正弦信号的频率是多少赫兹?抽样周期是多少秒?正弦信号的频率f=13Hz,抽样周期T=0.1s。
Q5.3 解释两个axis命令的效果。
给x,y轴标刻度。
Q5.4 以比在程序P5.1中列出的抽样周期低的两个抽样周期和高的两个抽样周期的四个其他值,运行程序P5.1.评论你的结果。
T=0.04s T=0.08sT=0.15s T=0.3s由上图可以发现:当取的T越小时,得到的图形越接近原图形。
Q5.5 通过将正弦信号的频率分别变为3HZ和7HZ,重做习题Q5.1。
相应的等效离散时间信号与习题Q5.1中产生的离散时间信号之间有差别么?若没有,为什么没有?f=3Hz f=7Hz由图可以看出,变换频率得到的两个图没有区别,因为他们的抽样周期一样。
Q5.6 运行程序P5.2,产生离散时间信号x[n]及其连续时间等效ya[t],并显示它们。
clf;T = 0.1;f = 13;n = (0:T:1)';xs = cos(2*pi*f*n);t = linspace(-0.5,1.5,500)';ya = sinc((1/T)*t(:,ones(size(n))) - (1/T)*n(:,ones(size(t)))')*xs;plot(n,xs,'o',t,ya);grid;xlabel('时间, msec');ylabel('振幅');title('重构的连续时间序号 y_{a}(t)');axis([0 1 -1.2 1.2]);图1 图2Q5.7 在程序P5.2中,t的范围和时间增量的值是什么?在图中,t的范围是什么?改变t的范围,显示上述程序所计算的全范围ya[t]并再次运行程序P5.2,。
中国石油大学《数字信号处理》第五章 有限离散傅立叶变换与快速傅立叶变换
X ( f ) x(n)e
n 0
N 1
j2nf
第一节 离散傅立叶变换
对于一个周期信号,只需研究其在一个周期内的 变化过程和规律。
因此,对XΔ(f )只要分析[0,1/Δ]范围内的频谱即可。
频谱XΔ(f )在一个周期内[0,1/Δ]是连续的,即所 对应的采样点数是无穷多个。实际上在计算有限离 散序列的频谱时,只能计算有限个点上的值,因此 1 把XΔ(f )在区间[0,1/Δ]分成N等份,每份的间隔为 N 则得到有限离散频谱:
对于非周期连续时间信号x(t),其频谱X(f )为一 个无限长的非周期连续频谱函数 。
x(t ) X ( f )
F
若用采样间隔Δ对信号x(t)进行离散采样,则得到离 散序列x(nΔ) ,其频谱XΔ(f )是周期连续,可以表示 为:以1/Δ为周期将X(f )向两边延拓而得到的。则有
m X( f ) X ( f ) m
第一节 离散傅立叶变换 二、离散傅立叶变换(DFT)
上面讨论的离散傅立叶级数中xp(n)和Xp(k)都是周期性 的序列,是无限长的。
1、有限长离散序列和周期离散序列的关系
有限长序列x(t)
周期延拓 T T N
周期性的序列xp(n)
离散采样
x(t ) x p (t ) x p (n)
第一节 离散傅立叶变换
数字信号处理
Digital Signal Processing(D S P )
中国石油大学
数字信号处理
Digital Signal Processing(D S P )
第五章 有限离散傅立叶变换与快速傅立叶变换
第五章 有限离散傅立叶变换与快速傅立叶变换
[工学]第五章 数字信号处理1AD_DA_采样定律、能量泄漏
西安工业大学机电学院
(2)采样定理
上述情况表明,如果fs>2fm,就不发生频混现象,因此对 采样脉冲序列的间隔Ts须加以限制,即采样频率fs(1/Ts) 必须大于信号x(t)中的最高频率fm的两倍,即
f s > 2 fm 或 ω s > 2 ω m
为了保证采样后的信号能真实地保留原始模拟信号 的信息,采样信号的频率必须大于原信号中最高频 率成分的2倍。这是采样的基本法则,称为采样定 理。
西安工业大学机电学院
3.3 采样定理
采样是将采样脉冲序列p(t)与连续时间信 号x(t)相乘,取离散点x(nt)的值的过程。
西安工业大学机电学院
X(0), X(1), X(2), ……, X(n)
西安工业大学机电学院
1)采样信号的频谱 采样:等间隔的脉冲序列p(t)与连续时间信号x(t)相乘。 理想脉冲采样过程如下:
当采样信号的频率低于被采样信号的最高频率时,采样所得的信 号中混入了虚假的低频分量,这种现象叫做频率混叠。
西安工业大学机电学院
(1) 频混现象 频混现象又称频谱混叠效应:因采样频率过低造成信号失 真的现象。
f
频带范围为-fm~fm
f s f
其周期为fs=1/Ts
f s f
f
西安工业大学机电学院
•当采样周期Ts较小时,fs>2fm,周期谱图相互分离, 如图b; •当采样周期Ts较大时,fs<2fm,周期谱图相互重叠, 即谱图之间高频与低频部分发生重叠,如图 c,此即 频混现象,这将使信号复原时丢失原始信号中的高频 信息。
X ( f kf )
s
•一个连续信号经过理想采样以后,它的频谱将沿着 频率轴每隔一个采样频率fs ,重复出现一次,即其 频谱产生了周期延拓。
同济大学数字信号处理课件第五章1数字滤波器的基本结构3
B1 z z1B1 z1
B2 z B1 z k2z1B1 z 1 k1z1 k1k2z1 k2z2
B2
z
k2 B1
z
z 1B1
z
k2
k1k2 z1
k2 z 1
z 2
2020/6/12
B2 z z2B2 z1
课件
8
L Bm z zmBm z1
a(3) 1
1.8313708,
a(3) 2
1.4319595,
a(3) 3
0.448
二阶系统:A2 (z) 11.4886262z1 0.7650549z2
a(2) 1
1.4886262,
a(2) 2
0.7650549
一阶系统: A1(z) 1 0.8433879z1
a(1) 1
0.8433879
z
2
B 2020/6/12 m1 z zBm z z课k件mBm1 z
3
6
Bm
z
Bm1
z
km z
B 1 m1
z
Bm
z
km
Bm1
z
z
B 1 m1
z
Bm1 z zBm z zkmBm1 z
1 2
3
(3)代入(1)得(4)
Bm1
z
1 1 km2
Bm
z
km Bm
z
4
(4)代入(3) 得:
b12
k2b12
0.8433879
k1
得一阶系统:
B1(z) 1 0.8433879z1
2020/6/12
课件
13
2、全极点系统(IIR系统)的格型结构
数字信号处理教程(程佩青)第五章习题讲解
考虑分子分母的组合及级联的次序,共有以下 四种级联型网络:
2023/11/3
信号处理
3、给出以下系统函数的并联型实现:
H z= 5.2 1.58z1 1.41z2 1.6z3 1 0.5z1 1 0.9z1 0.8z2
解:对此函数进行因式分解并展成部分分式,得
H z 5.2 1.58z1 1.41z2 1.6z3 1 0.5z1 1 0.9z1 0.8z2
1 0.5z1 1 0.9z1 0.8z2
则 A4
11 1 21 0 0.5 2023/11/311 21 0
12 1.4 22 1 信号处理12 0.9 22 0.8
11 1 21 0 11 0.5 21 0
12 1.4 22 1 12 0.9 22 0.8
H z
1 z1
1
1
6
10
z 1
1
1
6
10
z 1
1 7 10
1 7 10
2 20
2 20
1 1 10 z1 1 1 10 z1
6
6
1
1.61 0.69 z 1
1
0.61 0.36 z 1
2023/11/3
信号处理
H
z
1
1.61 0.69 z 1
1
0.61 0.36 z 1
H z=
1 z1
1 z1 1 z1 z2
5 3z3 1 z1 z2
2023/11/3
信号处理
由 H z 5 3z3 1 z1 z2
得
Hk
H
z
5 3e jk
z2 k N
j k
j 2 k
1 e 3 e 3
程佩青_数字信号处理_经典版(第四版)_第5章_5.1-5.3
N阶FIR 数字滤波器
N 1
N 1
H (z) h(k)zk bk zk
k 0
k 0
M
yn bk xn k k 0
h(n)* x(n) (5.3.1)
M
N
特点: yn bk xn k ak yn k
k 0
k 1
(1)系统单位抽样响应h(n)是有限长的。
(5.1.1)
(2)必须全部的ak = 0(k=1,2,…,N),即结构上没有输出 到输入的反馈(没有递归结构,只有非递归结构)。
5.2 IIR滤波器的基本结构
5.2.1 IIR滤波器的特点
h(n)是无限长的, h(n)、 ak 、bk 均为实数 有某一个ak 0(k=1,2,3,…,N),即一定是递归结构 至少有一个ak 0(k=1,2,3,…,N),故H(z)在有限z平 面上一定有极点存在,即不一定稳定
若只有b0 0,其它bk = 0(k=1,2,3,…,N),称为全极 点型IIR滤波器,或自回归(AR)系统
试画出其直接型、
H(z)
33
级联型和并联型结构。
(1 1 z 1)(1 1 z 1 1 z 2 )
3
22
解: 并联型
将系统函数H(z)表达为 H (z) 2
1 z 1
部分分式之和的形式
1 1 z 1 1 1 z 1 1 z 2
1
3
22
2
一阶节
1/ 3
z 1
yy([nk)]
x(xn[k) ]
x[k] A
11 21
z 1 11 z 1 21
1L
z
1
1L
2L
z
1
2L
数字信号处理-第五章数字滤波器的基本结构(new)
H ( z) A
将两个一阶因子组合成二阶因子,则
数字信号处理-第五章 数字滤波器络结构及 FIR数字滤波器的基本网络结构
数字信号处理-第五章 数字滤波器的基本结构
滤波器表示方式
(1)系统函数
k b z k M
Y ( z) H ( z) X ( z)
1 ak z k
k 1
k 0 N
1 ak z k
k 1
k 0 N
N2 M N Ak Bk (1 g k z 1 ) k G z k 1 1 * 1 1 c z ( 1 d z )( 1 d z ) k 1 k 1 k 0 k k k N1
一般IIR滤波器满足
N1
数字信号处理-第五章 数字滤波器的基本结构
5.2 无限长单位冲激响应(IIR)滤波器的基本结构)
IIR滤波器有以下几个特点: (1)系统的单位冲激响应 (2)系统函数
h( n)
是无限长的
H ( z)
在有限z平面(
0 z
)上有极点存在
(3)结构上存在输出到输入的反馈,也就是结构是递归的 1、直接Ⅰ型 一个IIR滤波器的有理系统函数为:
x n
3 1.5 -1.5 0.5
z 1 z 1 z 1
-3.5 2.5
y n
数字信号处理-第五章 数字滤波器的基本结构 级联型:
3z 3 3.5z 2 2.5z 3 3.5z 1 2.5z 2 1 H ( z) 2 2 z z 1 z 0.5 1 z z 1 0.5z 1
数字信号处理第5章
第5章 数字滤波器的基本结构5.1 学习要求1 掌握IIR 数字滤波器的基本网络结构,包括直接型、级联型和并联型;2 掌握FIR 数字滤波器的基本网络结构,包括直接型、级联型和频率抽样型;3 了解数字信号处理中的量化效应和数字信号处理的实现。
5.2 学习要点5.2.1 数字滤波器的结构特点与表示方法一个数字滤波器可以用系数函数表示为:01()()()1Mkk k N kk k b zY z H z X z a z -=-===-∑∑ (5-1) 直接由此式可得出表示输入输出关系的常系数线性差分方程为:1()()()N Mk k k k y n a y n k b x n k ===-+-∑∑ (5-2)由式(5-2)看出,实现一个数字滤波器需要几种基本的运算单元—加法器、单位延时和常数乘法器。
这些基本的单元可以有两种表示法:方框图法和信号流图法,如图5-1所示。
用方框图表示较明显直观,用流图表示则更加简单方便。
z ⊕aa单位延时乘常数相加方框图表示法信号流图表示法图5-1 基本运算过程的表示5.2.2 无限长单位脉冲响应(IIR)滤波器的基本结构无限长单位脉冲响应(IIR)滤波器有以下几个特点:(1) 系统的单位脉冲响应()h n 是无限长的;(2) 系统函数()H z 在有限z 平面(0z <<∞)上有极点存在; (3) 结构上存在着输出到输入的反馈,也就是结构上是递归型。
同一种系统函数()H z 的基本网络结构有直接I 型、直接Ⅱ型、级联型和并联型四种。
1直接I 型直接型按式(5-2)差分方程式将输入采样值(序列))(n x 延迟并乘以系数k b ,将输出采样(序列))(n y 延迟并乘以系数k a ,再把它们加起来,这种结构称为直接I 型,结构流图如图5-2所示。
由图可看出,总的网络)(z H 由Mkk k b z-=∑和11Nkk k a z-=-∑两部分网络级联组成,第一个网络实现零点,第二个网络实现极点,从图中又可看出,直接I 型结构需要N M +级延时单元。
数字信号处理第五章chhy
1 1 H1 ( z ) 1 H1 ( z ) 1 0.8 z 1 0.15z 2 1 2 1 1 0.8 z 0.15z 2 1.5 2.5 1.5 2.5 H 2 ( z) H 2 ( z ) 1 0.3z 1 1 0.5 z 1 1 1 2 1 0.3z 1 1 0.5 z 1 1 1 1 1 H3 ( z) H 3 ( z ) 1 0.3z 1 1 0.5 z 1 1 1 3 1 0.3z 1 1 0.5 z 1
5.3 IIR系统的基本网络结构
例:已知IIR数字滤波器的系统函数,画出该滤波器的直接型结构。
8 4 z 1 11z 2 2 z 3 H ( z) 5 3 1 1 z 1 z 2 z 3 4 4 8
5 3 1 y (n 1) y (n 2) y (n 3) 4 4 8 5 3 1 y (n) y (n 1) y (n 2) y (n 3) 8 x(n) 4 x(n 1) 11x(n 2) 2 x(n 3) y (n)
i 0 i 1
M
N
对输入信号的直接算法,
其系统函数H(z)为: M
刻以前的y(n-i),可以 y (n ) bi x (n i ) ai y (n
i 0 递推出y(n) i 1
已知x(n)、ai、bi和n时 M
N
N i i 0 YY z z Y i Y( 1) ii aNz NaaYY ( ()z )zjiz,Hj((,zH z )( z( ))) ( z)0 i bi z i j i i i0 i ) ) X X z i ( z )z Y ( z )Y ( zY (bi bi( z()z (z z z ja j Y z , H z ) ( z bi X ) i 0 H ( z ) X ( )) X ( z )NN j jN i 1 j 1 1 X (z z N i 0 i 0 i 0 j j i=11 a j j z i X ( z ) 11 a1z a 1 j1 ai z j 1 j 1
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
3、模拟滤波器和数字滤波器
经典滤波器从功能上分又可分为:
1、低通滤波器( LPAF/LPDF)
(Low pass analog filter / Low pass digital filter)
2、高通滤波器( HPAF/HPDF)
(High pass analog filter / High pass digital filter)
x(n) b0
第一部分是一 个对输入x(n) 的 M节延时链结构 。即每个延时 抽头后加权相 加,即是一个 横向网络
第五章 数字滤波器结构 Digital Filter Structures
第一节 引言
一、什么是数字滤波器
顾名思义: 其作用是对输入信号起到滤波的作用; 即DF是由差分方程描述的一类特殊的 离散时间系统。
功能: 把输入序列通过一定的运算变换成输出 序列。不同的运算处理方法决定了滤波 器的实现结构的不同。
三、数字滤波器表示方法
? 表示方法: 方框图表示法、流图表示法 ? 三种运算: 相加、乘以常数、延时 ? 基本运算单元: 加法器、单位延时、乘常数的
乘法器。
1、方框图、流图表示法
相加
方框图表示法
乘常数
a
延时
z-1
信号流图表示法
a
z-1
例:二阶数字滤波器:
y(n) ? a1 y(n ? 1)? a2 y(n ? 2)? b0x(n)
IIR DF类型有: 直接型、级联型、并联型。
直接型结构: 直接I型、直接II型(正准型、典范型)
1、 IIR DF系统函数及差分方程
一个N阶IIR DF有理的系统函数可能表示为:
M
?? H (z) ?
bi Z ? i
i?0 N
1 ? ai Z ? i
?
Y (z) X (z)
i?1
注:以下我们讨论 M<=N情况。
|H(ej? )|
|Y(ej? )|
有用
?c ?
?c ?
?c ?
? 如果信号和噪声的频谱相互重叠,那么经典滤波器将 无能为力,此时可以设计现代滤波器来解决。
2、现代滤波器
? 它主要研究内容是从含有噪声的数据记录(又称时间序列) 中估计出信号的某些特征或信号本身。一旦信号被估计出, 那么估计出的信号将比原信号会有 高的信噪比。
第二节 IIR DF的基本结构
一、IIR DF特点
1、单位冲激响应 h(n) 是无限长的: n→∞ 2、系统函数 H(z)在有限长 z平面(0<|z|< ∞)
有极点存在。 3、结构上存在输出到输入的反馈,也即结构上
是递归型的。 4、因果稳定的 IIR滤波器其全部极点一定在单位圆内。
二、IIR DF基本结构
则这一系统差分方程为:
N
M
? ? y(n ) ? ai y(n ? i ) ? bi x (n ? i )
i?1
i? 0
2、直接I型
? 直接I型流图
N
M
? ? y( n ) ? ai y(n ? i ) ? bi x (n ? i )
i?1
i? 0
IIR DF的差分方程就代表了一种最直接的计算公式,用流图表 现出来的实现结构即为直接I型结构(即由差分方程直接实现)。
其方框图及流图结构如下:
x(n) b0 a2
y(n) a1 z-1
z-1
x(n) b0
a1
z-1
a2
z-1
y(n)
说明:可通过流图或方框图看出系统的运算步骤和运算结构。 以后我们用流图来分析数字滤波器结构。
四、数字滤波器的分类
滤波器的种类很多,分类方法也不同。
1、从功能上分;低通、带通、高通、带阻。 2、从实现方法上分: FIR、IIR 3、从设计方法上来分: Butterworth (巴特沃斯)、
3、带通滤波器( BPAF/BPDF)
(Bandpass analog filter / Bandpass digital filter)
4、带阻滤波器( BSAF/BSDF)
(Bandstop analog filter / Bandstop digital filter)
4、模拟滤波器的理想幅频特性
二、数字滤波器的工作原理
设:x(n) 是系统的输入, X(ej? )是其傅立叶变换; y(n) 是系统的输出, Y(ej? )是其傅立叶变换;
则:
x(n)
y(n)
h(n)
LTI 系统的输出为:
?
? y(n) ? x(m )h(n ? m ) ? F ?1[ X (e j? )H (e j? )] m ? ??
Chebyshev( 切比雪夫)、 Ellips (椭圆)等。 4、从处理信号分:经典滤波器、现代滤波器
1、经典滤波器
? 假定输入信号x(n) 中的有用成分和希望去除的成分,各自 占有不同的频带。
? 当x(n) 经过一个线性系统(即滤波器)后即可将欲去除的 成分有效地去除。
|X(ej? )| 无用
? 滤波器的基本特性 (如有限长冲激响应 FIR与无限 长冲激响应 IIR)决定了结构上有不同的特点。
? 不同结构所需的 存储单元及乘法次数不同 ,前者影 响复杂性,后者影响 运算速度 。
? 有限精度(有限字长)实现情况下,不同运算结构 的误差及稳定性 不同。
? 好的滤波器结构应该易于 控制滤波器性能 ,适合于 模块化实现,便于时分复用。
? 现代滤波器把信号和噪声都视为随机信号,利用它们的统计 特征(如自相关函数、功率谱等)导出一套 最佳估值算法, 然后用硬件或软件予以实现。
? 现代滤波器理论源于维纳在40年代及其以后的工作,这一类 滤波器的代表为维பைடு நூலகம்滤波器,此外,还有卡尔曼滤波器、线 性预测器、自适应滤波器。
注:本课程主要讲经典滤波器
………
H(e j? )
? 3? ? 2? ? ? ? ? c ? c ? 2?
BPAF
………
? 3? ? 2? ? ?
BPAF
………
? 3? ? 2? ? ?
H(e j? )
? c1? c2 ? 2?
H(e j? )
? c1? c2 ? 2?
………
?
………
3? ?
………
3? ?
………
3? ?
五、研究数字滤波器结构意义
LPAF
H ( j? )
HPAF
?? c
?c
?
H ( j? )
BPAF
?? c
?c
?
H ( j? )
BSAF
? c1 ? c2
?
H ( j? )
? ? c2 ? ? c1 ? c1 ? c2
?
5、数字滤波器的理想幅频特性
LPAF
………
H(e j? )
? 2? ? ? ? ? c ? c ? 2?
HPAF