《数字滤波网络演示教学》由会员分享,可在线阅读,更多相关《数字滤波网络演示教学(34页珍藏版)》请在金锄头文库上搜索。
1、通信教研室,1,2020/8/6,第4章 数字滤波网络,zhang,2020/8/6,通信教研室,2,数字滤波器基本结构,为了处理信号,必须要涉及和实现称之为滤波器的各种系统。首先关注一下在设计中这些滤波器是如何实现的。 基本单元: 加法器 乘法器 延迟单元(移位器或存储器),2020/8/6,通信教研室,3,基本单元:,加法器,乘法器,延迟单元(移位器或存储器),2020/8/6,通信教研室,5,由三种结构用于实现一个IIR滤波器,直接型-按差分方程直接实现 这类滤波器有两个部分:滑动平均部分和递归部分(等效为分子部分和分母部分)。这种实现导致两种形式:直接 I 型和直接 II 型 级联型:
2、H(z)分解为二阶环节乘积的级联实现。 并联型: H(z)分解为二阶环节之和的并联实现。,2020/8/6,通信教研室,6,直接 I 型结构,2020/8/6,通信教研室,7,直接 II 型,2020/8/6,通信教研室,8,MATLAB实现,在MATLAB中,直接型结构用两个行向量来描述:含有系数bn的向量 b 和含有系数an的向量 a ,其结构实现用filter函数。,2020/8/6,通信教研室,9,级联型:,H(z)分解为二阶环节乘积的级联实现。,2020/8/6,通信教研室,10,其中Hk(z),称为第k个二阶环节,,且有,2020/8/6,通信教研室,11,级联型MATLAB实现,
3、已知直接型滤波器系数bn 和an,必须求得系数b0,Bk,i和Ak ,i 。可由函数dir2cas来完成。,函数dir2cas 将向量b和a转换为K3的B和A矩阵,例(p178.m),2020/8/6,通信教研室,12,N=4的级联结构,2020/8/6,通信教研室,13,例题:(p178.m)的解,2020/8/6,通信教研室,14,并联型,H(z)分解为二阶环节和的级联实现。,2020/8/6,通信教研室,15,N=4的并联结构,C0,B1,0,B1,1,- A1,1,- A1,2,- A2,1,- A2,2,B2,0,B2,1,z-1,z-1,z-1,z-1,x(n),y(n),2020
4、/8/6,通信教研室,16,并联型MATLAB实现,-18,-10.05,-3.95,- 1,- 0.05,0.25,0.125,28.1125,-13.3625,z-1,z-1,z-1,z-1,x(n),y(n),例题:(p178.m)的并联解(p183.m),2020/8/6,通信教研室,17,FIR滤波器结构,设单位脉冲响应h(n)长度N,其系统函数H(z)和差分方程分别为,2020/8/6,通信教研室,18,FIR滤波器结构的4种结构,直接型:直接安插分方程给出 级联型:H(z)分解为二阶因式,以级联形式实现 线性相位型:利用具有线性相位FIR滤波器对称性的结构形式 频率采样型:基于频
5、率响应H(ejw)采样的并联结构形式,2020/8/6,通信教研室,19,直接型,y(n) =b0 x(n)+b1x(n-1) +b2x(n-2) +b3x(n-3) +b4x(n-4),y(n),x(n),b0,b1,b2,b3,b4,z-1,z-1,z-1,z-1,2020/8/6,通信教研室,20,级联型,H(z)分解为二阶因式,以级联形式实现,2020/8/6,通信教研室,21,级联型,FIR滤波器结构,y(n),x(n),B1,1,B1,2,z-1,b0,z-1,z-1,z-1,z-1,z-1,B2,1,B2,2,B3,1,B3,2,2020/8/6,通信教研室,22,线性相位型,H
6、(ejw)=const-aw,|w|pai,const=0 or pi/2,对称脉冲响应,反对称脉冲响应,2020/8/6,通信教研室,23,线性相位型FIR滤波器结构,对称脉冲响应(M=7),2020/8/6,通信教研室,24,线性相位型FIR滤波器结构,对称脉冲响应(M=6),2020/8/6,通信教研室,25,MATLAB实现,直接型:差分方程为,线性相位型:差分方程为,2020/8/6,通信教研室,26,级联型(p189.m),2020/8/6,通信教研室,27,频率采样型,重写频率采样公式,(5.4.2),2020/8/6,通信教研室,28,频率采样型FIR滤波器结构,2020/8/
7、6,通信教研室,29,实系数频率采样型FIR滤波器结构,利用DFT和旋转因子的对称性(5.4.2)可化为,N为奇,L=(N-1)/2;N为偶L=N/2-1;而,2020/8/6,通信教研室,30,实系数频率采样型FIR滤波器结构图,N=4实系数频率采样型FIR滤波器结构图,2020/8/6,通信教研室,31,实系数频率采样型FIR滤波器MATLAB实现,function dir2fs 例:设h(n)=1,2,3,2,1/9,画出频率采样结构 p194.m,2020/8/6,通信教研室,32,p194.m的解,y(n),z-5,x(n),0.2,0.5818,-1,z-1,1,-0.809,0.618,z-1,z-1,-0.809,-1,0.0848,-0.309,-1.618,z-1,z-1,-1,2020/8/6,通信教研室,33,状态变量法,2020/8/6,通信教研室,34,本章练习,教材第5章习题:1,2,5,13,