序列和离散系统教材课件

黄黄 海海哈尔滨理工大学哈尔滨理工大学2015年年3月月 DSP第一章第一章 离散时间信号与系统离散时间信号与系统Q数字信号处理概述数字信号处理概述Q离散时间信号与系统离散时间信号与系统QZ变换与离散时间傅里叶变换变换与离散时间傅里叶变换Q离散傅里叶变换离散傅里叶变换Q快速傅里叶变换（快速傅里叶变换（FFT）Q数字滤波器基本结构数字滤波器基本结构Q无限长单位冲激响应数字滤波器的设计方法无限长单位冲激响应数字滤波器的设计方法Q有限长单位冲激响应数字滤波器的设计方法有限长单位冲激响应数字滤波器的设计方法软件学院软件学院Slide 1 DSP本章主要内容本章主要内容F离散时间信号离散时间信号序列序列F离散时间系统离散时间系统F线性移不变系统线性移不变系统F常系数线性差分方程常系数线性差分方程F连续时间信号的抽样连续时间信号的抽样软件学院软件学院Slide 2 DSP离散时间信号是对模拟信号xa(t)进行等间隔采样获得的，采样间隔为T，得到：1.1 1.1 离散时间信号离散时间信号序列序列0txa(t)0 xa(nT)tT2T软件学院软件学院Slide 3 DSP这里n取整数。对于不同的n值，xa(nT)是一个有序的数字序列，该数字序列就是离散时间信号。注意，这里的n取整数，非整数时无定义，另外，在数值上它等于信号的采样值，即 离散时间信号的表示方法：公式表示法、图形表示法、公式表示法、图形表示法、集合符号表示集合符号表示法，如软件学院软件学院Slide 4 DSP二、常用序列二、常用序列1.单位抽样序列(n)01/t(t)0(1)t(t)1n0(n)软件学院软件学院Slide 5 DSP2.单位阶跃序列单位阶跃序列u(n)t0u(t)10nu(n)软件学院软件学院Slide 6 DSP(n)与与u(n)之间的关系之间的关系令n-k=m，有软件学院软件学院Slide 7 DSP3.矩形序列矩形序列RN(n)N为矩形序列的长度0nR4(n)123软件学院软件学院Slide 8 DSP4.实指数序列实指数序列，a为实数0n0a1a-1或-1a0,序列的幅值摆动0n-1a00na0 时，序列右移延迟延迟当 n00 时，序列左移超前超前x(n)n0n0 x(n-2)软件学院软件学院Slide 18 DSP4.序列的翻转序列的翻转n0 x(-n)x(-n)是x(n)的翻转序列。x(-n)是以纵轴（n=0）为对称轴将序列x(n)加以翻转。x(n)n0软件学院软件学院Slide 19 DSP5.尺度变换尺度变换x(n)n0n0 x(2n)是序列每隔m点取一点形成的，相当于时间轴n压缩了m倍。抽取序列抽取序列是序列相邻抽样点间补（m1)个零值点，表示零值插值。插值序列插值序列软件学院软件学院Slide 206.累加（等效积分）累加（等效积分）7.7.差分运算差分运算 前向差分前向差分 后向差分后向差分8.8.卷积和卷积和等效为等效为翻褶、移位、相乘和相加翻褶、移位、相乘和相加四个步骤。四个步骤。软件学院软件学院Slide 21 DSP1.2 1.2 线性移不变系统线性移不变系统离散时间系统Tx(n)y(n)在时域离散系统中，最重要、最常用的是线性时不变系统是线性时不变系统。系统可定义为将输入序列x(n)映射成输出序列y(n)的唯一变换或运算，并用T表示，即软件学院软件学院Slide 22 DSP1.2.1 线性系统线性系统若系统满足可加性与比例性,则称此系统为离散时间线性系统。其中a、b为任意常数。设软件学院软件学院Slide 23 DSP例例是否线性系统。证：所以，此系统是线性系统。软件学院软件学院Slide 24 DSP例例所代表的系统不是线性系统。证：但是所以，此系统不是线性系统。软件学院软件学院Slide 25 DSP1.2.2 时不变系统（移不变系统）时不变系统（移不变系统）时不变系统Tx(n)y(n)若则n0为任意整数。输入移动任意位（如n0位），其输出也移动这么多位，而幅值却保持不变。软件学院软件学院Slide 26 DSP例例证：所以，此系统是时不变系统。软件学院软件学院Slide 27 DSP例例证：所以，此系统不是时不变系统。同理，可证明所代表的系统不是时不变系统。软件学院软件学院Slide 28 DSP1.2.3 线性时不变系统输入与输出之间的关系线性时不变系统输入与输出之间的关系T(n)h(n)一一个个既既满满足足叠叠加加原原理理，又又满满足足时时不不变变条条件件的的系系统统，被被称称为为线线性性时时不不变变系系统统(linear shift invariant,LTI)。线性时不变系统可用它的单位抽样响应来表征。单位取样响应，也称单位冲激响应，是指输入为单位冲激序列时系统的输出，一般用h(n)来表示：软件学院软件学院Slide 29根据线性系统的叠加性质又根据时不变性质设系统的输入用x(n)表示，而因此，系统输出为通常把上式称为离散卷积或线性卷积。这一关系常用符号“*”表示：软件学院软件学院Slide 30 DSP线线性性时时不不变变系系统统的的一一个个重重要要特特性性是是它它的的输输入入与与输输出出序列之间存在着线性卷积关系：序列之间存在着线性卷积关系：v用单位取样响应h(n)来描述系统h(n)x(n)y(n)软件学院软件学院Slide 31 DSP线性卷积的计算线性卷积的计算计算它们的卷积的步骤如下：(1)折叠：先在变量坐标轴k上画出x(k)和h(k)，将h(k)以纵坐标为对称轴折叠成 h(-k)。(2)移位：将h(-k)移位n，得h(n-k)。当n为正数时，右移n；当n为负数时，左移n。(3)相乘：将h(n-k)和x(k)的对应取样值相乘。(4)相加：把所有的乘积累加起来，即得y(n)。软件学院软件学院Slide 32 DSP例已知x(n)和h(n)分别为：和a为常数，且1a，试求x(n)和h(n)的线性卷积。计算线性卷积时，一般要分几个区间分别加以考虑，下面举例说明。软件学院软件学院Slide 33 DSP解参看图，分段考虑如下：(1)对于n4，且n-60，即46，且n-64，即64，即n10。0nx(n)40nh(n)6n-6mh(n-m)n软件学院软件学院Slide 34图解说明图解说明0mx(m)40mh(m)6-6mh(0-m)06(1)n0n-6mh(n-m)n 0(2)0n4n-6mh(n-m)n04软件学院软件学院Slide 35(3)4n6n-6mh(n-m)n04 6n-6mh(n-m)n06(4)610n-6mh(n-m)n04(2)0n4n-6mh(n-m)n04图解说明图解说明软件学院软件学院Slide 36(2)在在0n4区间上区间上n-6mh(n-m)n040mx(m)4软件学院软件学院Slide 37(3)在在4n6区间上区间上n-6mh(n-m)n04 60mx(m)4软件学院软件学院Slide 38(4)在在6n10区间上区间上n-6mh(n-m)n06100mx(m)4软件学院软件学院Slide 39 DSP综合以上结果，y(n)可归纳如下：软件学院软件学院Slide 40 DSP卷积结果y(n)如图所示 6ny(n)1004软件学院软件学院Slide 41 DSP例例 设有一线性时不变系统，其单位取样响应为解：分段考虑如下：(1)对于n0；(2)对于0n N1；(3)对于nN。软件学院软件学院Slide 42 DSP(2)在0nN 区间上Slide 43 DSP(3)在nN区间上(1)(2)(3)y(n)软件学院软件学院Slide 44 DSP例例设有一线性时不变系统，其3142x(m)m0 1 2 3 4215h(m)m102 3 4解：m0-2-3-4-11h(-m)软件学院软件学院Slide 45-3-11 20mh(1-m)-23-11 20mh(2-m)-2ny(n)-11 20-23 4 5 665241322103142x(m)m0 1 2 3 4软件学院软件学院Slide 46 DSP对有限长序列相卷，可用竖乘法注：1.各点要分别乘、分别加且不跨点进位；2.卷和结果的起始序号等于两序列的其实序号之和。软件学院软件学院Slide 47 DSP由上面几个例子的讨论可见，h(n)x(n)y(n)设x(n)和h(n)两序列的长度分别是N 和M，线性卷积后的序列长度为(N+M-1)。软件学院软件学院Slide 48 DSP线性卷积满足以下运算规律：线性卷积满足以下运算规律：交换律h(n)x(n)y(n)x(n)h(n)y(n)软件学院软件学院Slide 49结合律分配律h1(n)x(n)y(n)h2(n)h1(n)*h2(n)x(n)y(n)h1(n)x(n)y(n)h2(n)+h1(n)+h2(n)x(n)y(n)软件学院软件学院Slide 50 DSP软件学院软件学院序列与单位取样序列的线性卷积等于序列本身：序列与单位取样序列的线性卷积等于序列本身：若序列与一个移位的单位取样序列若序列与一个移位的单位取样序列(n-nn-n0 0)进行线性卷进行线性卷积，就相当于将序列本身移位积，就相当于将序列本身移位n n0 0：Slide 51序列与单位阶跃序列的线性卷积等于序列累加：序列与单位阶跃序列的线性卷积等于序列累加：DSP1.2.4 1.2.4 系统的因果性和稳定性系统的因果性和稳定性在系统中，若输出y(n)只取决于只取决于n n时刻，以及时刻，以及n n时刻以前的时刻以前的输入输入，即称该系统是因果系统。对于线性时不变系统，具有因果性的充要条件是系统的单位取样响应满足：如因果系统是指输出的变化不领先于输入的变化的系统。软件学院软件学院Slide 52 DSP稳定系统稳定系统对一个线性时不变系统来说，系统稳定的充要条件是单位取充要条件是单位取样响应绝对可和样响应绝对可和，即稳定系统是指对于每个有有界界输输入入x(n)x(n)，都都产产生生有有界界输输出出y(n)y(n)的 系 统。即 如 果|x(n)|M(M为 正 常 数)，有|y(n)|+，则该系统被称为稳定系统。软件学院软件学院Slide 53 DSP例例设某线性时不变系统，其单位取样响应为式中a是实常数，试分析该系统的因果稳定性。解：由于n0时，h(n)=0，故此系统是因果系统。所以时，此系统是稳定系统。软件学院软件学院Slide 54 DSP例例设某线性时不变系统，其单位取样响应为式中a是实常数，试分析该系统的因果稳定性。解：(1)讨论因果性由于n0时，h(n)0，故此系统是非因果系统。(2)讨论稳定性所以时，此系统是稳定系统。软件学院软件学院Slide 551.3 1.3 线性常系数差分方程线性常系数差分方程一个N 阶线性常系数差分方程用下式表示：连续时间线性时不变系统 线性常系数微分方程离散时间线性时不变系统线性常系数差分方程求解差分方程的基本方法有三种：经典法求齐次解、特解、全解递推法求解时需用初始条件启动计算变换域法将差分方程变换到Z域进行求解软件学院软件学院Slide 56 DSP1.4 1.4 连续时间信号的抽样连续时间信号的抽样连续时间信号离散时间信号采样内插1.信号经过采样以后，将发生一些什么变化？例如，信号频谱将发生怎样变化；2.经过采样后信号内容会不会有丢失；3.如果信号没有被丢失，其反变换应该怎样进行，即由数字信号恢复成模拟信号应该具备那些条件等。软件学院软件学院Slide 57 DSP1.4.1 采样采样S0tT2T0tP(t)T0txa(t)最高频率为fc软件学院软件学院Slide 58 DSP理想采样一、理想采样xa(t)P(t)0txa(t)0t0tT1T软件学院软件学院Slide 59 DSP定义单位冲击函数单位冲击函数t0(t)(1)单位冲击函数有一个重要的性质：采样性若f(t)为连续函数，则有将上式推广，可得t0(t-t0)软件学院软件学院Slide 60 DSP软件学院软件学院理想采样过程理想采样过程Slide 61 DSP即即采样信号的频谱是原模拟信号频谱的周期延拓，采样信号的频谱是原模拟信号频谱的周期延拓，其延拓周期为其延拓周期为 s s 。软件学院软件学院Slide 68 DSP讨论：讨论：S S/2/2 C C S S2 2 S S3 3 S S 0 0-S S(c)c)-C C C C S S/2/2 0 0(a)a)最高最高截止截止频率频率 S S/2/2 0 0-S S2 2 S S S S(b)b)软件学院软件学院Slide 69 DSP称称NyquistNyquist采样率采样率称折叠频率称折叠频率 C C S S/2/2 S S 0 0-S S称称NyquistNyquist范围范围采样定理采样定理 ：要想采样后能够不失真地还原出原信号，则采样频率必须大于要想采样后能够不失真地还原出原信号，则采样频率必须大于两倍原信号频谱的最高截止频率（两倍原信号频谱的最高截止频率（s s 2 2 C C）。）。由上面的分析有，频谱发生混叠的原因有两个：由上面的分析有，频谱发生混叠的原因有两个：1.1.采样频率低采样频率低2.2.连续信号的频谱没有被限带连续信号的频谱没有被限带软件学院软件学院Slide 70 DSP0C 2C 3C 4C 可选可选 s=(3 4)C 低通采样软件学院软件学院Slide 71 DSP软件学院软件学院理想采样后信号频谱发生的变化理想采样后信号频谱发生的变化2信号在时域的特性和在频域的特性之间存在以下对应关系：2时域的周期性 频域的离散性2时域的离散性 频域的周期性2时域的非周期性 频域的连续性2时域的连续性 频域的非周期性Slide 72 DSP软件学院软件学院傅里叶变换的各种形式时域信号频域信号非周期连续连续非周期周期连续离散非周期非周期离散连续周期周期离散离散周期Slide 73 DSP频域分析频域分析且在 时，0 0 T T S S/2/2-S S/2/2 G(j)g(t)采样的恢复采样的恢复软件学院软件学院Slide 74 DSP时，时，000软件学院软件学院Slide 75 DSP时域分析时域分析g(t)时，0T软件学院软件学院Slide 76 DSP或或称为称为内插函数内插函数软件学院软件学院Slide 77 DSP软件学院软件学院Slide 78 DSP采样内插公式采样内插公式采采样样内内插插公公式式说说明明：只只要要满满足足采采样样频频率率高高于于两两倍倍信信号号最最高高截截止止频频率率，则则整整个个连连续续时时间间信信号号就就可可以以用用它它的的采采样样值值来来完完全代表，而不会丢失任何信息。全代表，而不会丢失任何信息。软件学院软件学院Slide 79 DSP软件学院软件学院KT(K-1)T(K-2)T(K+1)T(K+1)TSlide 80 DSP软件学院软件学院Matlab 脚本脚本Z%signalrecoverZt=-2:0.1:8;Zy=sinc(t);采样函数Zforn=1:8Zy1=n.*sinc(t-n);Zplot(t,y1),holdonZendZplot(t,t)Slide 81 DSP软件学院软件学院Slide 82