第九讲MATLAB在信号处理中的应用(二)题材.ppt

第九讲 MATLAB在信号处理中的应用（二） 　　　;主要内容;Z变换;Z变换定义;滤波器：所谓的滤波器是指对输入信号起滤波的作用的装置。为了处理信号，必须设计和实现称之为滤波器（在某些领域或称为频谱分析仪）的各种系统。滤波器设计受到它的类型（IIR或FIR）和实现结构等因素的影响。所以有必要了解滤波器是如何实现的。当输入、输出是离散信号，滤波器的冲激响应是单位抽样响应 h（n）时，这样的滤波器称作数字滤波器（ Digital Filter，DF ）。DF是由差分方程来描述的一种特殊的离散时间系统。 ;常系数线性差分方程与系统函数;9.1 数字滤波器结构的表示方法; 实现数字滤波器的基本运算单元 单位延时 常数乘法器 加法器 这些单元有两种表示法： 1）方框图 2）信号流图;方框图与信号流图;2;9.2无限长单位冲激响??滤波器;例9-1二阶数字滤波器的系统差分方程为： 试构造它的直接型I型和直接II型方框图，并观测它们 的单位冲激响应。;转置定理;对调;直接II型; MATLAB提供了专门求离散系统单位响应，并绘 制其时域波形的函数impz()（数字滤波器冲激响应），该函数有以下几种调用格式： impz(b,a) %以默认格式绘制向量a和b定义的离 散系统的单位响应 impz(b,a,n) %绘制0~n时间范围内的单位响应波形 impz(b,a,n1:n2) %绘制n1~n2时间范围内的单位 响应波形 y=impz(b,a,[n1:n2]) %并不绘制波形，而是求出单 位响应的数值解;将差分方程化简为： a=[1 -1/3 -1/4]; b=[1 1 1]; impz(b,a,[0:10]) title(单位脉冲响应) xlabel(Time(s)) axis([0,10,-3,3]) ;2.级联型结构;用多个二阶节级联表示滤波器系统;例9-2用级联型结构实现以下系统函数 ;级联型结构的仿真建模;3.并联型结构;例9-3用并联型结构实现以下用系统函数表示为IIR滤波器。 法一：分解为： ;并联型结构仿真模型;例9-4用并联型结构实现以下用系统函数表示为IIR滤波器。 法二：解：自定义函数文件dir2par(b,a).m b=[6 1.6 4.8 -1.3];a1=[1 -0.5];a2=[1 0.9 0.8]; a=conv(a1,a2) [C,B,A]=dir2par(b,a) C = 3.2500 B = -0.2500 -2.5250 3.0000 0 A = 1.0000 0.9000 0.8000 1.0000 -0.5000 0;并联型结构仿真模型;9.3有限长单位冲激响应（FIR）滤波器的结构;FIR的系统函数及差分方程;1.直接型（卷积型、横截型结构）;举例;直接型 FIR滤波器 仿真模型;2.级联型结构;FIR 级联型仿真模型;3.FIR滤波器 线性相位型;举例：例3-12;IIR和FIR数字滤波器的比较;习题：;令n’=N-1-n 代入;其中h(0)=h(N-1),h(2)=h(N-2)……;当N=奇数时，有一中间项h((N-1)/2)无法合并，需提出来：;其中h(0)=h(N-1),h(2)=h(N-2)……,h((N-3)/2)=h((N-1)/2;当h(n)=偶对称时，即h(n)=h(N-1-n),可求出:;当h(n)=奇对称时，即h(n)=-h(N-1-n),可求出:;Z-1;h(N-1);举例;解1） 2） 3） 4）其线性相位型结构如右图所示