第九讲matlab在信号处理中的应用二.pptVIP

Z变换 信号与系统的变换域分析 连续时间信号与系统：拉普拉斯变换和傅里叶变换 离散时间信号与系统：z变换和离散时间傅里叶变换 Z变换的作用和数字滤波器的功能类似：将输入信号通过一定的运算变化成输出序列，它把描述离散系统的差分方程转化为简单的代数方程。 Z变换定义 若序列为x(n),则幂级数 称为序列x(n)的z变换，其中z为变量,它表示z平面，用极 坐标表示为： 。 收敛域： 即只要满足绝对可和的条件该级数收敛。 本章我们只讲右边序列即无限长序列左边为0. 滤波器：所谓的滤波器是指对输入信号起滤波的作用的装置。为了处理信号，必须设计和实现称之为滤波器（在某些领域或称为频谱分析仪）的各种系统。滤波器设计受到它的类型（IIR或FIR）和实现结构等因素的影响。所以有必要了解滤波器是如何实现的。当输入、输出是离散信号，滤波器的冲激响应是单位抽样响应 h（n）时，这样的滤波器称作数字滤波器（ Digital Filter，DF ）。DF是由差分方程来描述的一种特殊的离散时间系统。 常系数线性差分方程与系统函数 假设常系数线性差分方程的一般形式为： 若系统起始状态为零，直接两边取z变换，利用移位特性 得： 9.1 数字滤波器结构的表示方法 数字滤波器实现的方法： （1）利用计算机编程，即软件实现; （2）数字信号处理器（DSP）,即专用硬件实现。 数字滤波器的系统函数可以表示为： 其中ak和bk都是常量，当ak不等于0时，有反馈环路是 递归型结构。若ak=0则没有反馈，称为非递归结构。 补充知识 实现数字滤波器的基本运算单元 单位延时 常数乘法器 加法器 这些单元有两种表示法： 1）方框图 2）信号流图 方框图与信号流图 若二阶数字滤波器表示为： 方框图 信号流图 其中1,2,3,4,5称为网络节点，x(n)为输入节点（源节点），y(n) 处为输出节点（阱节点）节点之间用有向支路相连接，任一节点 的节点值等于它的所有输入支路的信号之和。而输入支路的信号 值等于这一支路起点处节点信号值乘以支路上的传输系数。如果 支路上不标注传输系数值，则认为其传输系数为1。 分支节点：只有一个输入，一个输出或多个输出。如2,3,4。 相加器：有两个或两个以上的输入。如1,5。 9.2无限长单位冲激响应滤波器 无限长单位冲激响应（IIR）滤波器的特点： 1）系统的单位冲激响应h(n)无限长； 2）系统函数H(z)在有限z平面内有极点存在； 3）从输出到输入存在反馈环路（回路），即存 在递归型结构。 1.直接型 直接型结构是按照给出的系统差分方程或者系统函数直接实现的。 例9-1二阶数字滤波器的系统差分方程为： 试构造它的直接型I型和直接II型方框图，并观测它们 的单位冲激响应。 转置定理 如果将线性时（移）不变网络中所有支路方向倒转，并将输入x(n)和输出y(n)相互交换，则其系统函数H(z)不改变。 1）将全部路径箭头方向颠倒过来。 2）将全部支路节点用加法器替换，将全部加法器节点替换成支路节点。 3）将输入和输出节点交换。 所得的结构称为转置直接型结构，即直接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.级联型结构 级联型结构是把系统函数H（z）进行因式分解为较低的二阶节，然后将系统函数表示成二阶节的乘积，其中每个二阶节用直接II型结构实现。 用多个二阶节级联表示滤波器系