《离散时间系统》PPT课件.ppt

,五、离散时间系统； 六、 LSI系统输入、输出关系； 七、 LSI系统的频率响应； 八、确定性信号的相关函数,Section B 离散时间系统基础,1.5离散时间系统,1. 系统定义 数字信号处理的任何处理都是依靠系统来完成的,所以系统是数字信号处理的核心，系统一般包括系统硬件和系统所完成的处理算法。 系统在数学上定义为将输入序列x(n)映射成输出序列y(n)的唯一性变换或运算。这种映射是广义的，实际上表示的是一种具体的处理，或是变换，或是滤波。,系统可以表示为 其中，符号T 表示系统的映射或处理，可以把T 简称为系统。 系统的图形表示如下图所示，输入x(n)称为系统的激励，输出y(n)称为系统的响应。由于它们均为离散时间信号，将系统T 称为离散时间系统或时域离散系统。,2 . 线形离散时间系统 满足叠加原理的系统，或满足齐次性和可加 性的系统称为线性系统。 设 y1(n)=Tx1(n)， y2(n)=Tx2(n) 对任意常数a,b，若 Tax1(n)+bx2(n)=aTx1(n)+bTx2(n) =a y1(n)+b y2(n) 则称T 为线性离散时间系统。,推广到一般情况，设 yk(n) = Txk(n) ， k=1，2，.N 线性系统满足 1kN 线性系统的特点是多个输入的线性组合的系统输出等于各输入单独作用的输出的线性组合。,例1 证明由线性方程表示的系统 是非线性系统。 证明 设 所以，该系统是非线性系统。,3 . 非时变离散时间系统 若满足下列条件，系统称为非时变（非移变）系统，或时不变（移不变）系统。 设 y(n) = Tx(n) 对任意整数k,有 y(n-k)=Tx(n-k) 即系统的映射T 不随时间变化，只要输入x(n)是相同的，无论何时进行激励，输出y（n）总是相同的，这正是系统非时变性的特征。,下图形象说明了系统非时变性的概念。,例2 设系统的映射 y = Tx(n) = nx(n) ，判断系统的线性和时不变性。 解 设 y1(n) = nx1(n), y2(n) =nx2(n) a1x1(n)+a2x2(n) = x(n) 则 Tx(n) = nx(n) = na1x1(n)+na2x2 = a1y1(n)+ a2y2(n) 所以，系统为线性系统。 设 y(n) = nx(n), x1(n) = x(n-k) y1(n) = nx1(n) = nx(n-k) 而 y(n-k) = (n-k)x(n-k)y1(n) 所以，系统为时变系统。,4 . 线性时不变离散系统 定义 同时具备线性和时不变性的系统称作线 性非时变系统或线性时不变系统。 它的重要意义在于，系统的处理过程可以统一采用这种系统的特征描述之一单位采样响应，以一种相同的运算方式卷积运算，进行统一的表示。 任何一个信号可以表示成单位采样序列的线性组合，即,系统对 的响应为 设系统对单位采样序列 的响应为 ， 即 称为系统的“单位采样响应”，它是描述系统的一个非常重要的信号。,根据时不变性，有 则系统输出y(n)可表示为 上式表明：当线性时不变系统的单位采样响应h(n)确定时，系统对任何一个输入x(n)的响应y(n)就确定了， y(n)可以表示成x(n)和h(n)之间的一种简单的运算形式。将上式的运算方式称作“离散卷积”，简称“卷积”，采用符号“*”表示，即,5 . 离散卷积的计算 卷积的计算一般采用两种方法：解析法和图解法，或是两种方法的结合。 例2.3 设线性时不变系统的单位脉冲响应和输入序列如下图所示，画出输出的波形。,图2.2 例2.2图解法,（2）采用解析法。 因为 所以 将x(n)的表达式代入上式，得到 两种方法结果一致。,6 . 离散卷积的运算规律 (1) 交换率 h(n)*x(n) = x(n)*h(n) 它的意义可以解释为，如果互换系统的单位采样响应h(n)和输入x(n)，系统的输出保持不变。 (2) 结合率 x(n)*h1(n)*h2(n)=x(n)*h2(n)*h1(n) =x(n)*h2(n)*h1(n) 它的意义可以解释为一种级联系统结构，级联顺序可以交换，或系统级联可以等效为一个系统，输出保持不变。,(3) 分配率 x(n)*h1(n)+x(n)*h2(n) = x(n)*h1(n)+ h2(n) 它的意义可以解释为一个并联系统结构，或并联系统可以等效为一个系统，输出保持不变。 (4) 与(n)卷积不变性 x(n)*(n) = x(n) 它的意义可以解释为输入通过一个零相位的全通系统。 (5) 与(n-k)卷积的移位性 x(n)*(n-k) = x(n-k) 它的意义可以解释为输入通过一个线性相位的全通系统。,2. 系统的稳定性和因果性,2.1 稳定性 稳定系统是有界输入产生有界输出的系统。 若 则 线性时不变离散系统是稳定系统的充要条件：,2.2 因果性 若系统 n时刻的输出，只取决于n时刻以及n时刻以前的输入序列，而与n时刻以后的输入无关，则称该系统为因果系统。 线性时不变离散系统是因果系统的充要条件：,例3 某线性时不变离散系统，其单位采样响 应为 试讨论其是否是因果的、稳定的。 解 讨论因果性: 该系统是非因果系统 讨论稳定性: 当 时系统稳定,当 时系统不稳定。,结论:因果稳定的线性时不变系统的单位取样响应是因果的,且是绝对可和的,即,1.6离散时间系统的 频率响应特性,离散系统频响特性的定义 频响特性的几何确定法,1离散系统频响特性的定义,正弦稳态（正弦序列作用下系统的稳态响应）,系统对不同频率的输入，产生不同的加权，这就是系统的频率响应特性。,由系统函数得到频响特性,输出对输入序列的相移,离散时间系统在单位圆上的z变换即为傅氏变换，即系统的频率响应特性:,输出与输入序列的幅度之比,：幅频特性,：相频特性,通过本征函数透视系统的频响特性,为输入序列的加权， 体现了系统对信号的处理功能。 是 在单位圆上的动态， 取决于系统的特性。,离散系统（数字滤波器）的分类,2频响特性的几何确定法,几点说明,小结,1. 系统的频响特性 ：幅频特性，输出与输入序列的幅度之比 ：相频特性，输出对输入序列的相移,3.因为 是周期为 的周期函数，所以系统的频响 特性 为周期为 的周期函数。,4. 是关于 的偶函数， 是关于 的奇函数。,2.系统的频率响应就是系统函数在单位圆上的动态， 因 而变化，影响输出的幅度与相位。,相关是研究两个信号之间，或一个信 号和其移位后的相关性，是信号分 析、检测与处理的重要工具；在随机 信号的理论中起到了中心的作用。,x(n), y(n),n = 0,1 2,xy =,x(n)y(n) n x2(n)y2(n) n n,| xy |1 相关系数,rxy(m) ryx(m),所以ryx(m) = rxy(m), x(n +i)y(n + j),相关函数的时间变量：, n=,rxy(m) =,3.,可正可负。,= rxy(n + j)(n +i)= rxy( j i) 1. 保持 x(n)不动，将 y(n) 往左， 或右移动 m 个抽样间隔，然后 将 x(n) 和y(n +m) 对应相乘与相 加，即得 rxy(m) ； 2. x(n)和 y(n) 的长度应一样长；,m, x(n)x(n +m), x (n)x(n + m), x (n)y(n + m),自相关函数：, , n= n= n=,rx(m) = rx(m) = rxy(m) =,实序列 复序列,rx(m) = rx(m), rx(m) = rx (m)； 性质： rx(0) rx(m)； Lim rx(m) = 0 m,卷积和相关的关系： = x(n m)y(n) = x(m) y(m) n=,定义: 两个序列的关系,描述LSI系统输入输出关系,计算: 任一序列都不需,要翻转,其中一个序列要翻转,rxy(m) = lim,rx(m) = lim,x(n)x(n +m),功率信号相关函数的定义：,1 N 2N +1,1 N 2N +1,互相关,自相关,对于能量信号 ：, n=,rx(m) =,自相关,功率信号自相关函数的性质： 1. 若 x(n)是周期的, 周期是N, 则,rx(m) = rx(m+ N),2. 若x(n)是实的, 则 rx(m) = rx(m),m,3. rx(0)取最大值, rx(0) = P x 为信号功率 4. rx(m) = 0,5. 若 x(n)是复信号, 则,sin(n)sin(n +m),sin (n), = n 0sin(n)cos(n),n = 0,1, N 1,2 N,例: x(n) = sin(n) =,2,1 2,1 N,1 N,cos(m),N1 n=0,N1 n=0,1 N1 N,rx(m) =,= cos(m),+sin(m),=,同频率余弦,例：信号的检测 x(n) = s(n)+u(n) （白噪声）,rx(m) =s(n)+u(n)s(n + m)+u(n + m) n,= rs(m)+ru(m),1 N,rx(m) =,实际计算相关函数时： x(n), n = 0,1,L,N 1,m = M M, M N 1 所以，rx(m) 的最大长度为 2N 1,to be continued,第二章 z变换和DTFT,主要内容：,1、z变换的定义及收敛域 2、z变换的反变换 3、z变换的基本性质和定理,2.1 z变换的定义及收敛域,信号和系统的分析方法有两种： 时域分析方法 变换域分析方法 连续时间信号与系统 LT FT 离散时间信号与系统 ZT FT,一、ZT的定义,z 是复变量，所在的复平面称为z平面,二、ZT的收敛域,对于任意给定序列x(n)，使其z变换X(z)收敛的所有z值的集合称为X(z)的收敛域。 级数收敛的充要条件是满足绝对可和,1）有限长序列,除0和两点是否收敛与n1和n2取值情况有关外，整个z 平面均收敛。,如果n20 ，则收敛域不包括点 如果n10 ，则收敛域不包括0点 如果n10n2，收敛域不包括0 、点,2）右边序列,因果序列的z变换必在处收敛 在处收敛的z变换， 其序列必为因果序列,3）左边序列,4）双边序列,例1,收敛域应是整个z的闭平面,例2：求x(n)=anu(n)的变换及其收敛域,给定z变换X(z)不能唯一地确定一个序列，只有同时给出收敛域才能唯一确定。 X(z)在收敛域内解析，不能有极点，故： 右边序列的z变换收敛域一定在模最大的有限极点所在圆之外 左边序列的z变换收敛域一定在模最小的有限极点所在圆之内,2.2 z反变换,实质：求X(z)幂级数展开式 z反变换的求解方法： 围线积分法（留数法） 部分分式法 长除法,z反变换: 从X(z)中还原出原序列x(n),1、围数积分法求解（留数法）,若函数X(z)zn-1在围数C上连续，在C以内有K个极点zk，而在C以外有M个极点zm，则有：,1、围数积分法求解（留数法）,根据复变函数理论，若函数X(z)在环状区域 内是解析的，则在此区域内X(z)可展开成罗朗级数，即 而 其中围线c是在X(z)的环状 收敛域内环绕原点的一条 反时针方向的闭合单围线。,若F(z)在c外M个极点zm，且分母多项式z的阶次比分子多项式高二阶或二阶以上，则：,利用留数定理求围线积分，令,若F(z)在围线c上连续，在c内有K个极点zk，则：,单阶极点的留数：,2、部分分式展开法求解IZT ：,常见序列的ZT参见书p.54页的表2-1,若函数X(z) 是z的有理分式，可表示为：,利用部分分式的z反变换和可以得到函数X(z) 的z反变换。,例2 设 利用部分分式法求z反变换。,解：,3、幂级数展开法求解（长除法）:,一般X(z)是有理分式，可利用分子多项式除分母多项式（长除法法）得到幂级数展开式，从而得到x(n)。,根据收敛域判断x(n)的性质，展开成相应的z的幂级数 将X(z) X(z)的 x(n) 展成z的 分子分母 按z的 因果序列 负幂级数 降幂排列 左边序列 正幂级数 升幂排列,例1,ROC1:,长除法示例,解：由Roc判定x(n)是因果序列，用长除法展成z的负幂级数,ROC2:,解：由Roc判定x(n)是左边序列，用长除法展成z的正幂级数,1、线性性,2.3 Z变换的基本性质和定理,R1R2,R,|a|R,R,2、序列的移位,3、z域尺度变换 （乘以指数序列）,4、 z域求导 （序列线性加权）,Z变换的基本性质（续）,5、翻褶序列,1/R,R,6、共轭序列,7、初值定理,8、终值定理,Z变换的基本性质（续）,9、有限项累加特性,ZT的主要性质参见书p.69页的表2-2,10、序列的卷积和,11、序列乘法,12、帕塞瓦定理,2.4 序列ZT、连