第六章 数字滤波器的基本结构
《数字信号处理教程》程佩青(第三版)清华大学出版社课后答案
9
c) 设 x(n) = δ (n) + δ (n − 1)
i)向n > 0 处递推
y3 (1) = ay3 (0) + x3 (1) = 1 y3 (2) = ay3 (1) + x3 (2) = a y3 (3) = ay3 (2) + x3 (3) = a2
┇ y3 (n) = ay3 (n − 1) + x3 (n) = a n−1 ∴ y3 (n) = a n−1 , n ≥ 1 ii)向 n < 0 处递推
解:(1 )
n
y(n) = ∑ x(m ) m = −∞
n
y1 (n ) = T [x1 (n )] = ∑ x1 (m ) m = −∞
y2 (n ) = T [x2 (n )] =
n
∑ x2 (m )
m = −∞
n
ay1(n)+ by2 (n) = ∑[ax1(m) + bx2 (n)] m = −∞
10
T [ax1(n)+ bx2 (n)] =
n
∑
[ax1
(n
)
+
bx2
(n
)]
m = −∞
T[ax1(n) + bx2(n)] = ay1(n) + by2(n)
∴ 系统是线性系统
解:(2) y(n) =
[x(n )] 2
y1(n)
= T [x1(n)] = [x1(n)] 2
y2 (n) = T [x2 (n)] = [x2 (n)] 2
试判断系统是否是线性的?是否是移不变的?
分析:已知边界条件,如果没有限定序列类型(例如因果序列、反因果序列等), 则递推求解必须向两个方向进行(n ≥ 0 及 n < 0)。
第六章 无限冲击响应数字滤波器设计
分子分母同除Ω , 得 Ha (s) = :
N c
1 s − sk ) ∏(Ωc Ωc k =0
N −1
因为
令λ = ΩΩc 称为归一化频率;p = jλ 称为归一化复变量
归一化巴特沃斯的传输函数为
Ha ( p) = 1
s = jΩ Ωc Ωc
∏( p − p )
k =0 k
N −1
Ha ( p) =
sk = (−1)
( jΩc ) = Ωce
1 2k −1 jπ [ + ] 2 2N
, k =1 ⋯,2N ,2
jΩ
例如,N=3时, 例如, 时
s0 =Ωce
s3 =Ωce
j 2π 3
s0
s5
s1 = −Ωc s2 =Ωce
− j 2π 3
s1
s2
s4
s3
− j 1π 3
s4 = Ωc
s5 =Ωce
Ha ( jΩ)
2
由于一般滤波器的单位冲响应为实数, 由于一般滤波器的单位冲响应为实数,其传递函数是对称 的,有:
* Ha ( jΩ) = Ha (s)Ha (−s) s= jΩ = Ha ( jΩ)Ha ( jΩ) 2
确定 Ha ( jΩ) 极、零点,并将左半S平面极点分配给 Ha (s) , 零点,并将左半 平面极点分配给 此系统是因果稳定的。 得到滤波器的传递函数 Ha (s) ,此系统是因果稳定的。
Ha ( jΩ)
1
1 2
Ωc
巴特沃斯幅度特性和N的关系 巴特沃斯幅度特性和 的关系
Ω
3、巴特沃斯滤波器的极、零点分布 、巴特沃斯滤波器的极、 由于 H (s)H (−s) = a a 在
数字信号处理 第六章
各种数字滤波器的理想幅度频率响应 数字滤波器的设计步骤 理想滤波器的逼近 数字滤波器的系统函数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滤波器具有较好的幅频特性,但相位特性较差,且存 在稳定性问题。
在实际应用中,应根据具体需求选择合适的滤波器类型 和设计方法。
04
数字滤波器的实现
数字滤波器的实现步骤
确定滤波器参数
设计滤波器系数
根据实际需求,确定滤波器的阶数、截止 频率等参数。
根据滤波器类型和参数,计算滤波器系数 。
实现滤波器算法
描述滤波器实现的难易程度,包括运算量和 存储需求。
数字滤波器的基本结构
直接实现型
将输入信号直接与滤波器系数进行运算,得到输 出信号。
级联实现型
将滤波器分解为若干个简单滤波器的级联,以降 低计算复杂度。
并行实现型
将滤波器分解为若干个简单滤波器的并行运算, 以提高处理速度。
03
数字滤波器的设计方法
窗函数法
验证滤波器效果
根据滤波器系数,编写滤波器算法,实现 信号的滤波处理。
对滤波后的信号进行验证,确保满足设计 要求。
数字滤波器的编程实现
选择编程语言
根据实际需求,选择适合的编程语言,如C、 Python等。
设计滤波器函数
根据滤波器算法,编写滤波器函数,实现信 号的滤波处理。
测试滤波器函数
对滤波器函数进行测试,确保其正确性和稳 定性。
数字信号处理教程答案
数字信号处理教程 课后习题及答案目录第一章 离散时间信号与系统 第二章 Z 变换第三章 离散傅立叶变换 第四章 快速傅立叶变换 第五章 数字滤波器的基本结构第六章 无限长单位冲激响应(IIR )数字滤波器的设计方法 第七章 有限长单位冲激响应(FIR)数字滤波器的设计方法 第八章 数字信号处理中有限字长效应第一章 离散时间信号与系统1 。
直接计算下面两个序列的卷积和)n (h *)n (x )n (y =请用公式表示。
分析:①注意卷积和公式中求和式中是哑变量m ( n 看作参量),结果)(n y 中变量是 n ,; )()()()()(∑∑∞-∞=∞-∞=-=-=m m m n x m h m n h m x n y ②分为四步 (1)翻褶( -m ),(2)移位( n ),(3)相乘,00 , 01()0 , ,()0,n n n a n N h n n n n x n n n β-⎧≤≤-=⎨⎩⎧≤⎪=⎨<⎪⎩其他; )( )( 4n y n n y n 值的,如此可求得所有值的)相加,求得一个( ③ 围的不同的不同时间段上求和范一定要注意某些题中在 n 如此题所示,因而要分段求解。
2 。
已知线性移不变系统的输入为)n (x ,系统的单位抽样响应 为)n (h ,试求系统的输出)n (y ,并画图)(5.0)(,)1(2 )()4()(5.0)(,)2( )()3()()(,)( )()2()()(,)( )()1(3435n u n h n u n x n R n h n n x n R n h n R n x n R n h n n x n n n =--==-=====δδ分析:①如果是因果序列)(n y 可表示成)(n y ={)0(y ,)1(y ,)2(y ……},例如小题(2)为()∑∑∑+-=+-=--+===-=-+≥nN n m mn n nN n m mn n m nn m m n h m x n y N n n 111N -00)()()( , 1)3(αββααβ全重叠时当()()()()βααβαβαβαββααβαβαβ==≠--=--=---+++--,)(,100111n n N N n N n n N n n nN n y ∑∞-∞=-==m m n h m x n h n x n y )()()(*)()(:解0)()1(0=<n y n n 时当, 1)2(00部分重叠时当-+≤≤N n n n ()∑∑∑==--===-=nn m mnn n n m mn n m nn m m n h m x n y 0)()()(αββααβ()()βαβαβαβααβαβαβ≠--=--=-+-++-,100111nn n n n n n n())(,1)(00βαα=-+=-n n n y n n)(n y ={1,2,3,3,2,1} ;②)()(*)( , )()(*)(m n x n x m n n x n x n -=-=δδ ;③卷积和求解时,n 的分段处理。
第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 的改变会影响所有零点或极点的分布;
DSP第6章数字滤波器基本网络结构
4.Ideal Amplitude-frequency character (模拟滤波器的理想幅频特性)
H ( j)
LPAF
c
s c H ( j) 2
HPAF
c c
H ( j) c
s 2
BPAF
s H ( j) c 2
BSAF
c 2 c1 c1
is described by a differential function. (顾名思义:其作用是对输入信号起到滤波 的作用;即DF是由差分方程描述的一类特 殊的离散时间系统) • Function: transfer the input sequence to the output sequence. • (它的功能:把输入序列通过一定的运算变 换成输出序列。不同的运算处理方法决定 了滤波器的实现结构的不同。)
第六章 DF (Digital Filter) 数字滤波器结构
第一节 Introduction: Represent of the structure on Digital Filter 引 言:数字滤波器结构的 表示方法
一、What is Digital Filter ( 什么是数字滤波器 ) Digital filter is used to filter the input signal. It
第二节 Basic Structure of IIR (IIR DF的基本结构)
1.h(n) is infinite. (单位冲激响应h(n)是无限长的n→∞)
一、Definition (IIR DF定义)
2.H(z) has pole(i.e. H(z) has denominator). (系统函数H(z)在有限长Z平面(0<|Z|<∞)有 极点存在(即H(z)有分母))
第六章数字滤波器的基本结构
可以表示为
M
bi zi
H(z)
i0 N
ai zi
i 1
一个线性时不变离散时间系统只包含三种基本运算单元: 加法器、标量乘法器和单位延时器。通常采用如图6.1 所示的符号表示基本运算单元。通过这些基本单元的组 合可以构成离散时间系统的方框图。
单位延迟器 Z-1
常数乘法器
加法器
图6.1 基本运算单元的方框图
4
4
8
11x(n 2) 2x(n 3)
x(n)
8
y(n)
Z-1
5/4
-4
-3/4 Z-1 11
1/8 Z-1 -2
直接型结构的优点是简单直观,所用的延 时 器数量少。缺点是系数ak,bk对滤波器的性能 控制作用不明显,更严重的是极点对系数的变 化过于灵敏,易出现不稳定或较大误差,以至 于积累的误差很大,因此对于三阶以上的IIR 滤波器。几乎都不使用直接型结构。而是采用 级联型、并联型等其他形式的结构。
k2
z-1
kN-1
z-1
kN-1 z-1
r1
r2
rN-1
N阶全零点格形网络
y(n)
kN kN
2 全极点系统的格形网络
eN x(n)
rN
eN-1
eN-2
e2
e1
kN
kN-1
-k-N z-1
-kN-1
rN-1
k2
z-1
-k2 z-1
r2
r1
N阶全极点格形网络
e0
y(n)
k1 -k1 z-1
r0
第六节 数字滤波器的MATLAB实现
16 20z1 1 z1 0.5z2
x(n)
第六章 无限长单位冲激响应(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 )
数字信号处理课件——06数字滤波器的结构
(6-3)
由于H z 的系数 ai , bi都是实系数,因此零极点 ci , di 只有两种情 况:或者是实根,或者是共轭复根。即
H( z )
( 1 gi z ) ( 1 hi z 1 )( 1 hi* z 1 )
1
M1
M2
( 1 pi z 1 ) ( 1 qi z 1 )( 1 qi* z 1 )
2 1
其结构如图6-6所示。
x n
y2 n
a0
z 1
y n
x n x n
H1 z H2 z
y1 n
y2 n
H2 z H1 z
y n
b1
z 1 z 1
y n 1
a1
b2
y n
y n 2
z 1
N
(6-1)
aN 1
bN 1
用差分方程可以表达为
y (n) ai x(n i ) bi y (n i )
i 0 i 1 N N
xn N
z 1
aN
bN
z 1
y n N
图6-5 N 阶数字滤波器的信号流图表达
(6-2)
y (n) ai x(n i ) bi y (n i )
a2
bN 1
aN 1
y n N
bN
z 1
z 1
aN
图6-6 直接型的变形
即信号先经过反馈网络 H2 z ,其输出为中间变量
H2 1 1 bi z i
i 1 N
y2 n
y2 n bi y2 n i x n
数字滤波器的基本结构
数字网络的信号流图表示
① 通路:沿同一方向传输的连通支路 ② 环路:闭合的通路 ③ 环路增益 : 环路中所有支路增益之积 ④ 前向通路 :从输入节点到输出节点通过 任何节点仅一次的通路 ⑤ 前向通路增益:前向通路中所有支路增 益之积
29
二阶数字滤波器的例子: y(n) a1y(n 1) a2 y(n 2) b0x(n)
级联型 I I R 数字滤波器
并联型
直接Ⅰ型 直接Ⅱ型
转置型
34
N
M
y(n) ak y(n k) bm x(n m)
k 1
m0
x(n)
b0
y(n)
Z 1
b1 x(n 1)
Z 1
x(n 2)
b2
Z 1 a1
y(n 1)
Z 1
a2
y(n 2)
Z 1 bM
x(n M )
Z 1
aN 1
y(n N 1)
M2
(1 pm z1) (1 qm z1)(1 qm z1)
A
m1 N1
m1 N2
(1 ck z1) (1 dk z1)(1 dkz1)
k 1
k 1
44
将共轭因子组合成实系数的二阶因子,两 个一阶构成一个二阶有:
M1
M2
(1 pm z1) (1 1m z1 2m z2 )
H (z)
A
M
bm zm
H(z)
m0 N
1 ak zk
k 1
式中 N N1 2N2
N1
Ak
k 1 1 ck z1
N2 k 1
Bk (1 gk z1)
(1
d
k
z
1
数字信号处理教程 程佩青 课后题答案
第一章 离散时间信号与系统2.任意序列x(n)与δ(n)线性卷积都等于序列本身x(n),与δ(n-n 0)卷积x(n- n 0),所以(1)结果为h(n) (3)结果h(n-2) (2(4)3 .已知 10,)1()(<<--=-a n u a n h n,通过直接计算卷积和的办法,试确定单位抽样响应为 )(n h 的线性移不变系统的阶跃响应。
4. 判断下列每个序列是否是周期性的,若是周期性的,试确定其周期:)6()( )( )n 313si n()( )()873cos()( )(ππππ-==-=n j e n x c A n x b n A n x a分析:序列为)cos()(0ψω+=n A n x 或)sin()(0ψω+=n A n x 时,不一定是周期序列,nmm m n n y n - - -∞ = - ⋅ = = ≥ ∑ 2 31 2 5 . 0 ) ( 01当 3 4n m nm m n n y n 2 2 5 . 0 ) ( 1⋅ = = - ≤ ∑ -∞ = - 当 aa a n y n a a an y n n h n x n y a n u a n h n u n x m m nnm mn -==->-==-≤=<<--==∑∑--∞=---∞=--1)(11)(1)(*)()(10,)1()()()(:1时当时当解①当=0/2ωπ整数,则周期为0/2ωπ;②;为为互素的整数)则周期、(有理数当 , 2 0Q Q P QP =ωπ ③当=0/2ωπ无理数 ,则)(n x 不是周期序列。
解:(1)0142/3πω=,周期为14 (2)062/13πω=,周期为6 (2)02/12πωπ=,不是周期的 7.(1)[][]12121212()()()()()()[()()]()()()()[()][()]T x n g n x n T ax n bx n g n ax n bx n g n ax n g n bx n aT x n bT x n =+=+=⨯+⨯=+所以是线性的T[x(n-m)]=g(n)x(n-m) y(n-m)=g(n-m)x(n-m) 两者不相等,所以是移变的y(n)=g(n)x(n) y 和x 括号内相等,所以是因果的。
六章节FIR数字滤波器设计-资料
2
推导方法与前面类似,可得:
h( N 1) 0 2
(N1)/2
H() c(n)sinn()
n1
c(n)2h(N1n),n1,2,, N1
2
2
上式表明,当 0,,2时,H()0,相当于 H (z)在 z
=1和z = -1有两个零点,并且由于sinn()对 0,,2 呈奇对称,因而 H() 对 0,,2也呈奇对称。
4
FIR数字滤波器的性质
H(ej)H(z)|zej
其求和项全为实数
j(N1)N1
e 2 h(n)c
o s[(nN1)]
n0
2
将 H (ej)表示成 ( 相 )和位 幅函 H 度 ()的 数 函 形 ,即 数 式
H (ej)H ( )ej()
则:H()
N1
h(n)cos[(n
N1)]
n0
2
() (N1)
版权所有 违者必究 第六章第1讲
1
§1 FIR数字滤波器的性质
【问题的引入】
优点:能借助模拟滤波器已有成果设计;
IIR数字滤波器
简单方便
缺点:相位非线性
优点:严格线性相位;系统十分稳定(传输 FIR数字滤波器 函数为全零点型);进行滤波时可采用FFT
缺点:幅度特性较差;设计繁琐
版权所有 违者必究 第六章第1讲
2
2
特即点H (:z)当在z =-1时为,零c点o,s[(且m由于12)c]os0[,(m故H12()])对0,
呈奇对称,因而 H() 对 也呈奇对称。
因此这种情况不适合做在 处不等于零的滤波
器,如高通滤波器。
版权所有 违者必究 第六章第1讲
数字信号处理第六章8设计IIR滤波器的频率变换法
频率变换法可以根据不同的设计需求,灵活地调 整滤波器的频率响应特性,实现多种多样的滤波 效果。
数学基础坚实
频率变换法基于严格的数学推导,具有坚实的理 论基础,保证了设计结果的准确性和可靠性。
频率变换法的局限性
对初值敏感
频率变换法对初值的选择较为敏感,初值选取不当可能导致设计 结果偏离预期目标。
模糊、锐化、边缘检测等效果。
控制系统
03
在控制系统中,频率变换法可用于设计IIR滤波器,以改善系统
的动态性能和稳定性。
04
CATALOGUE
结论
频率变换法的总结
频率变换法是一种有效的设计IIR滤波器的方法 ,通过将原始滤波器的频率响应转换为易于实 现的形式,从而简化了滤波器的设计过程。
该方法基于频率域的变换,通过一系列的数学 推导,将原始滤波器的频率响应转换为易于实 现的形式,从而得到滤波器的系数。
选择合适的频率变换方法
根据滤波器类型和性能指标,选择合 适的频率变换方法,如巴特沃斯、切 比雪夫等。
计算滤波器系数
根据选择的频率变换方法,计算滤波 器系数。
验证滤波器性能
通过仿真或实验验证滤波器的性能是 否满足设计要求。
设计实例
选择巴特沃斯频率变换方法。 通过仿真验证滤波器性能,确保满足设计要求。
03
在应用中,需要注意滤波器的 实时性要求,选择合适的算法 和实现方式以满足实时处理的 需求。
THANKS
感谢观看
数字信号处理第六 章8设计iir滤波器的 频率变换法
目录
• 引言 • 设计IIR滤波器 • 频率变换法在IIR滤波器设计中的应用 • 结论
01
CATALOGUE
引言
背景介绍
第6章IIR FILTER-1
4 j 5
s2 e , s4 e
7 j 5
j
6 j 5
37
(2)查表求极点: P1,2=-0.3090j0.9511; p3,4=-0.8090 j0.5857; p0=-1.000 所以
Ha ( p )
1
( p p )
k k 0
38
4
使用zp2tf命令
X(ejω )H(ejω )符合人们的要求,这就是数字滤波器的滤
波原理。
3
数字滤波器分类 同模拟滤波器一样,线性数字滤 波器按照频率响应的通带特性可划分为 低通、高通、带通和带阻几种形式。它 们的理想模式如图6-1所示。(系统的 频率响应H(ejω )是以2π 为周期的。)
4
H (e j )
(a)
y(n)是y(n-m)以及x(n-m)的函数。
IIR 数字滤波器一定包含某种形式的反馈。
7
m 0 M
m 0
N
m 1
IIR系统函数为:
b0 b1z b2 z ... bMz H ( z) 1 2 N 1 a1z a 2 z ... aNz
1
2
M
j
14
6.2 模拟滤波器设计
模拟滤波器的理论和设计方法已发展得 相当成熟,且有若干典型的模拟滤波器可供 选择。
这些滤波器都有严格的设计公式、现成的 曲线和图表供设计人员使用。
在matlab的help 中输入Classic IIR Filter Design
15
常用模拟原型滤波器及特点page153
由于滤波器冲激响应ha(t)是实函数,因而Ha(jΩ )满足
所以
数字滤波器的基本结构
四 并联型型结构
1.结构
将H(Z)展成部分分式形式:
H (Z ) A0 H1(Z ) H2 (Z ) H N (Z )
A0
N i 1
1
Ai pi z i
H( Z ) A0
E
i1
1
Ai pi z
i
F
i1
0i 1iZ 1 1 1iZ 1 2iZ 2
N E 2F
A0
x(n) A1
p1
极点:zk ej2πk/N , k 0, 1,,N 1
频率取样型结构流图
x(n)
z N
H[0] 1/ N y(n)
WN0
z 1
WN1
H[1]
z 1
H [N 1]
WN( N 1)
z 1
2.频率取样型结构优缺点
优x点(n)
H[0] 1/ N y(n)
H(k)是W z N=2πk/N 在处WN0的频z1 率响应值,
y(n N )
2 .特点
第一个网络实现零点,即实现x(n)加权延时:
M
bi x(n i)
i0
第二个网络实现极点,即实现y(n)加权延时:
N
ai y(n i)
k 1
可见,第二网络是输出延时,即反馈网络。
*共需(M+N)个存储延时单元。
二 直接型II结构
1.结构
H2(z)
H
(
z)
1
M
b
i0
z-1
x(n)
x(n 1)
2)信号流图法
x( n ) Z 1 x( n 1 )
2. 乘常数 1) 框图表示
y(n) a
2)信号流图法
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1
第四节 有限长脉冲响应基本网络结构
有限冲激响应(FIR)滤波器的单位抽样 有限冲激响应(FIR)滤波器的单位抽样 响 应为有限长度,一般采用非递归形式实现。 通常的FIR数字滤波器有横截型、级联型和频 通常的FIR数字滤波器有横截型、级联型和频 率采样型。我们假设其单位脉冲响应h(n)长 率采样型。我们假设其单位脉冲响应h(n)长 度 为N,因此它的传递函数一般具有如下形式。
H(Z) = ∑h(n)Z =∏α0i + 1i z +α2i z ) ( α
n 1 2 n=0 i=1
N1
M
其中M=N/2取整数值,这样就可以用二阶节级 其中M=N/2取整数值,这样就可以用二阶节级 联起来构成,如图所示: 联起来构成,如图所示:
第五节 数字滤波器的格型结构
1 全零点系统(FIR系统)的格型结构 全零点系统(FIR系统)
H ( z ) = ∑ h( n) z n
n =0
N 1
1. 横截型 我们从式可以看出,是的N 我们从式可以看出,是的N-1阶多项式。在 处有个零点可以位于z 处有个零点可以位于z平面上任意位置。系 统的差分方程为:
y (n) = ∑ h(i ) x(n i )
i =0 N 1
2. 级联型 当需要控制滤波器的传输零点时,可将传递函 数分解为二阶实系数因子的形式:
1. 直接型 IIR数字滤波器的系统函数一般表示为 IIR数字滤波器的系统函数一般表示为
bi Z i ∑ 1 ∑ a k Z k
k =1 i =0 N M ∞
H ( z) =
= ∑ h( n) Z n
n =0
例题1 例题1
IIR数字滤波器的系统函数H(z)为 IIR数字滤波器的系统函数H(z)为
84z +11z 2z H(z) = 5 1 3 2 1 3 1 z + z z 4 4 8
1 2 3
画出该滤波器的直接型结构。
解 由H(z)写出差分方程如下: H(z)写出差分方程如下:
5 3 1 y ( n ) = y ( n 1) y ( n 2) + y ( n 3) + 8 x ( n ) 4 x ( n 1) 4 4 8 + 11 x ( n 2) 2 x ( n 3)
IIR滤波器的差分方程为: IIR滤波器的差分方程为:
y (n ) = ∑ bi x(n i ) ∑ ai y (n i )
i =0 i =1 M N
FIR滤波器的差分方程为:
y (n ) = ∑ bi x(n i )
i =0 M
IIR滤波器的常用结构形式有:直接型、级 IIR滤波器的常用结构形式有:直接型、级 联型和并联型 FIR滤波器的常用结构形式有:直接型、 FIR滤波器的常用结构形式有:直接型、 级联型和频率采样型。
第六章 数字滤波器的基本结构
第一节 第二节 第三节 第四节 第五节 第六节 小结 引言 离散时间系统的表示方法 无限长脉冲响应基本网络结构 有限长脉冲响应基本网络结构 数字滤波器的格型结构 数字滤波器的MATLAB实现 数字滤波器的MATLAB实现
第一节 引言
数字滤波器是一种常用的离散时间系统, 数字滤波器是一种常用的离散时间系统 , 其 功能就是采取不同的软件或硬件把输入序列通 过一定的运算变换成输出序列。 过一定的运算变换成输出序列。 本章研究的重点是根据系统函数( 本章研究的重点是根据系统函数 ( 或差分方 程),设计不同运算结构实现数字滤波器,以 设计不同运算结构实现数字滤波器, 及数字滤波器的两大类型---无限长脉冲响应滤 及数字滤波器的两大类型---无限长脉冲响应滤 波器和有限长脉冲响应滤波器。 波器和有限长脉冲响应滤波器。
M
H (Z ) =
∑ bk z
k =0 N k =1
k
=A
(1 p k z )∏ (1 q k z )(1 q k z 1 ) ∏
1 1
M1
M
2
1 ∑ ak z k
(1 c k z 1 )∏ (1 d k z 1 )(1 d k z 1 ) ∏
k =1 k =1
k =1 N1
k 数 1 2 3 8 4 z + 11z 2 z H ( z) = 1 2 3 1 1.25z + 0.75z 0.125z 试画出其级联型网络结构。
解 将H(z)分子分母进行因式分解,得到 H(z)分子分母进行因式分解,得到
(2 0.379 z 1 )(4 1.24 z 1 + 5.264 z 2 ) H ( z) = (1 0.25z 1 )(1 z 1 + 0.5z 2 )
3. 并联型 如果将级联形式的H(z)展成部分分式形式, 如果将级联形式的H(z)展成部分分式形式,则得到 IIR并联型结构。 IIR并联型结构。
H (z ) = a i z i ∑
i =0 N
1 ∑ bi z i
i =1
N
Ai = A0 + ∑ 1 d i z i i =1
N
例题3 例题3
2.无限脉冲响应(IIR)滤波器的基本结构 .无限脉冲响应(IIR)滤波器的基本结构 3.有限脉冲响应(FIR)滤波器的基本结构 .有限脉冲响应(FIR)滤波器的基本结构 4.格型结构的数字滤波器的分类和特点
画出例题2中的H(z)的并联型结构 画出例题2中的H(z)的并联型结构
8 4 z + 11z 2 z H ( z) = 1 2 3 1 1.25z + 0.75z 0.125z
试画出其级联型网络结构。
1
2
3
解: 将H(z)分子分母进行因式分解,得到 H(z)分子分母进行因式分解,得到
8 16 + 20 z H ( z ) = 16 + + 1 1 2 1 0.5 z 1 z + 0.5 z
第二节 离散时间系统的表示方法
1 离散时间系统结构的方框图表示法
数字滤波器是根据预定要求对序列进行变换,实 现某种“算法”的离散时间系统,因此其系统函数 M 可以表示为 i
H ( z) =
∑b z
i =0 N i =1 i
ai z i ∑
一个线性时不变离散时间系统只包含三种基本运算单元: 加法器、标量乘法器和单位延时器。通常采用如图6 加法器、标量乘法器和单位延时器。通常采用如图6.1 所示的符号表示基本运算单元。 所示的符号表示基本运算单元。通过这些基本单元的组 合可以构成离散时间系统的方框图。 合可以构成离散时间系统的方框图。
第三节 无限长脉冲响应基本网络结构
IIR数字滤波器的单位脉冲响应无限长,其系统 IIR数字滤波器的单位脉冲响应无限长, 函数H(z) 在有限z 平面上存在极点, 函数 H(z)在有限 z 平面上存在极点 , 因而结构上存 在 反馈支路,也就是说结构上是递归的。IIR基本网络 反馈支路,也就是说结构上是递归的。IIR基本网络 结构可分为三种,即直联型、级联型和并联型。 结构可分为三种,即直联型、级联型和并联型。
N阶全零点格形网络
2 全极点系统的格形网络
N阶全极点格形网络
第六节 数字滤波器的MATLAB实现 数字滤波器的MATLAB实现
数字滤波器的类型分为IIR和FIR两种。IIR 数字滤波器的类型分为IIR和FIR两种。IIR 滤波器的特征是具有无限持续时间脉冲响应, 一般情况下是递归型滤波器结构。FIR滤波器 一般情况下是递归型滤波器结构。FIR滤波器 的特征是有限持续时间脉冲响应,一般情况下 是非递归型滤波器结构。
直接型结构的优点是简单直观,所用的延 时 器数量少。缺点是系数a 器数量少。缺点是系数ak,bk对滤波器的性能 控制作用不明显,更严重的是极点对系数的变 化过于灵敏,易出现不稳定或较大误差,以至 于积累的误差很大,因此对于三阶以上的IIR 于积累的误差很大,因此对于三阶以上的IIR 滤波器。几乎都不使用直接型结构。而是采用 级联型、并联型等其他形式的结构。
2. 级联型 因为IIR滤波器是递归运算, 因为IIR滤波器是递归运算,所以对高阶系统有 误差积累和容易不稳定的缺点。 误差积累和容易不稳定的缺点。可以采取使用低阶 系统来级联方法实现,即将系统函数H(z)的分子、 系统来级联方法实现,即将系统函数H(z)的分子、 分母按零、极点进行因式分解。 分母按零、极点进行因式分解。
小结
1.数字滤波器结构的表示方法 任何一个数字滤波器的输出等于系统输人各 延 迟信号与输出各延迟信号的线性组合,在数字网 络中只有三种基本运算单元——乘法器、加法器 络中只有三种基本运算单元——乘法器、加法器 和延迟器。因此可以用三种基本运算——乘法、 和延迟器。因此可以用三种基本运算——乘法、 加法和单位延迟——流图组合来表示数字滤波器 加法和单位延迟——流图组合来表示数字滤波器 的算法结构。用信号流图表示滤波器的结构比其 他方法简单、方便。
图6.1 基本运算单元的方框图
2 离散时间系统结构的信号流图表示法
由节点与有向支路构成的能表征系统功能与信号流动方向的图,称为 系统的信号流图。离散时间系统结构的信号流图表示法基本单元的符号 如图6.2所示,其运算结构和方框图表示是相同的。 如图6.2所示,其运算结构和方框图表示是相同的。
图6.2 基本单元的信号流图表示