离散时间信号与系统教程(ppt 103页).ppt

第一章离散时间信号与系统,2,本章目录,离散时间信号序列,离散时间系统,线性常系数差分方程,连续时间信号的取样,Matlab实现,3,1.1引言,信号信号与信息信号的表示信号的分类,系统系统的作用系统的分类系统的描述与分析,4,信号与信息,信号是信息的表现形式,信息则是信号的具体内容交通灯信号传递的信息：红灯停而绿灯行。,信号是传递信息的函数数学上表示成一个或多个独立变量的函数一维变量：时间或其它参量语音信号表示为一个时间变量的函数静止图像信号表示为两个空间变量的亮度函数,5,信号的分类,连续时间信号：连续时间域内的信号幅度可以是连续数值，或是离散数值,离散时间信号：离散时间点上的信号幅度同样可以是连续数值，或是离散数值,特殊形式：模拟信号和数字信号模拟信号：时间和幅度都是连续数值的信号，实际中与连续时间信号常常通用。数字信号：时间和幅度都离散化的信号。,6,本章主要内容,离散时间信号的基本概念,离散时间系统的定义及其性质,线性常系数差分方程及其求解方法,理想取样：连续时间信号数字处理的概念和基本方法,Matlab实现,7,1.2离散时间信号序列,序列的定义及表示,序列的基本运算,几种常用序列,序列的周期性,用单位脉冲序列表示任意序列,8,1.2.1序列的定义及表示,序列的定义数字序列：离散时间信号一般只在均匀间隔的离散时间nT上给出数值,序列的表示x=x(n)，-n+(1.1)图1.1图形表示用单位脉冲序列表示,9,序列表示,x=x(n)，-n+,n代表nT,nT指均匀间隔的离散时间点,T采样时间间隔,n为非整数时没有定义，不能认为此时x(n)的值是零,10,图1.1序列的图形表示,11,1.2.2序列的基本运算,和积移位标乘翻转,累加差分时间尺度变换序列的能量卷积和,12,基本运算序列的和,设序列为x(n)和y(n)，则序列z(n)=x(n)+y(n)(1.2)表示两个序列的和，定义为同序号的序列值逐项对应相加。,13,例：序列的和,例1.1设序列,计算序列的和x(n)+y(n)。,解：,14,例：序列求和图示,15,基本运算序列的积,设序列为x(n)和y(n)，则序列z(n)=x(n)y(n)(1.3)表示两个序列的积，定义为同序号的序列值逐项对应相乘。,16,例：序列的积,例1.1设序列,计算序列的和x(n)y(n)。,解：,17,例：序列求积图示,x(n),18,基本运算序列的移位,设序列为x(n)，则序列y(n)=x(n-m)(1.4)表示将序列x(n)进行移位。,m为正时x(n-m)：x(n)逐项依次延时(右移)m位x(n+m)：x(n)逐项依次超前(左移)m位m为负时，则相反。,19,例：序列的移位,例1.1设序列,计算序列的和x(n+1)。,解：,20,例：序列移位图示,x(n),21,基本运算序列的标乘,设序列为x(n)，a为常数(a0)，则序列y(n)=ax(n)(1.5)表示将序列x(n)的标乘，定义为各序列值均乘以a，使新序列的幅度为原序列的a倍。,22,例：序列的标乘,例1.1设序列,计算序列的和4x(n)。,解：,23,基本运算序列的翻转,设序列为x(n)，则序列y(n)=x(-n)(1.6)表示以n=0的纵轴为对称轴将序列x(n)加以翻转。,24,例：序列的翻转,例1.2设序列,计算序列的和4x(n)。,解：,25,基本运算序列的累加,设序列为x(n)，则序列(1.7)定义为对x(n)的累加，表示将n以前的所有x(n)值求和。,26,基本运算序列的差分,前向差分：将序列先进行左移，再相减x(n)=x(n+1)-x(n)(1.8),后向差分：将序列先进行右移，再相减x(n)=x(n)-x(n-1)(1.9),由此，容易得出x(n)=x(n-1),27,多阶差分运算,二阶前向差分,二阶后向差分,单位延迟算子D，有Dy(n)=y(n-1),y(n)=y(n)-y(n-1)=y(n)-Dy(n)=(1-D)y(n),=1-D,k阶后向差分,(按二项式定理展开),二阶后向差分,28,基本运算时间尺度(比例)变换,设序列为x(n)，m为正整数，则序列抽取序列y(n)=x(mn)(1.10),x(mn)和x(n/m)定义为对x(n)的时间尺度变换。,插值序列(1.11),29,抽取序列,x(mn)：对x(n)进行抽取运算不是简单在时间轴上按比例增加到m倍以1/m倍的取样频率每隔m-1个点抽取1点。保留x(0),30,插值序列,x(n/m)：对x(n)进行插值运算表示在原序列x(n)相邻两点之间插入m-1个零值点保留x(0),31,基本运算序列的能量,设序列为x(n)，则序列(1.12)定义为序列的能量，表示序列各取样值的平方之和；若为复序列，取模值后再求平方和。,32,基本运算序列的卷积和,设序列为x(n)和z(n)，则序列(1.13)定义为x(n)和z(n)的卷积和。卷积和又称为离散卷积或线性卷积，是很重要的公式。,33,卷积和计算的四个步骤,翻转：x(m)，z(m)z(-m)移位：z(-m)z(n-m)n为正数时，右移n位n为负数时，左移n位相乘：z(n-m)x(m)（m值相同）相加：y(n)=z(n-m)x(m),34,对应点相乘！,例：卷积和计算,例1.3设序列,求y(n)=x(n)*z(n)。,解：,n0时，x(m)与z(n-m)没有重叠，得y(n)=0。,0n4时，,对应点相乘！,35,例：卷积和计算,4n6时，,4n6时，,n10时，x(m)与z(n-m)没有重叠，得y(n)=0。,36,1.2.3几种常用序列,单位脉冲序列单位阶跃序列矩形序列实指数序列正弦序列复指数序列,37,单位脉冲序列,(n)只在n=0时取确定值1，其它均为零(n)类似于(t),(n-m)只有在n=m时取确定值1，而其余点取值均为零,38,单位阶跃序列,u(n)类似于u(t)u(t)在t=0时常不定义，u(n)在n=0时为u(0)=1,(n)和u(n)的关系：,(n)=u(n)-u(n-1),39,单位矩形序列,N为矩形序列的长度,和u(n)、(n)的关系：,40,实指数序列,a为实数,当|a|1时序列收敛,当|a|1时序列发散,41,正弦序列,A为幅度为数字域角频率为起始相位,x(n)由x(t)=sint取样得到,x(n)=Asin(n+),归一化:=T=/fs(与线性关系),42,复指数序列,为数字域角频率,用实部与虚部表示,用极坐标表示,=0时，序列具有以2为周期的周期性,43,1.2.4序列的周期性,对于序列x(n)，如果对所有n存在一个最小的正整数N，满足x(n)=x(n+N)则序列x(n)是周期序列，最小周期为N。以正弦序列为例讨论周期性,设x(n)=Asin(n+),则有x(n+N)=Asin(n+N)+=Asin(N+n+),若满足条件N=2k，则,x(n+N)=Asin(n+N)+=Asin(n+)=x(n),44,周期性讨论,N、k为整数，k的取值满足条件，且保证N最小正整数。其周期为,2/为整数时，取k=1，保证为最小正整数。此时为周期序列，周期为2/。,例1.4序列，因为2/=8，所以是一个周期序列，其周期N=8。,45,周期性讨论,2/为有理数而非整数时，仍然是周期序列，周期大于2/。,例1.5序列，2/=8/3是有理数，所以是周期序列，取k=3，得到周期N=8。,2/为无理数时，任何k都不能使N为正整数，这时正弦序列不是周期序列。,例序列,指数为纯虚数的复指数序列的周期性与正弦序列的情况相同。,46,1.2.5用单位脉冲序列表示任意序列,任何序列都可以用单位脉冲序列的移位加权和来表示，即,x(n)可看成是x(n)和(n)的卷积和，式中,例1.6,47,1.3离散时间系统,离散时间系统的定义及表示,线性时不变系统,单位脉冲响应与卷积和,线性时不变系统的性质,因果系统和稳定系统,48,1.3.1离散时间系统的定义及表示,离散时间系统定义为将输入序列x(n)映射成输出序列y(n)的惟一变换或运算。以T表示这种运算y(n)=Tx(n),对变换T加以不同的约束条件，所定义的系统就具有不同的特性和功能。线性时不变系统:最重要、最常用，可表征许多物理过程。,49,1.3.2线性时不变系统,线性系统满足叠加原理叠加原理包含可加性和齐次性两方面性质,时不变系统系统的响应与输入信号施加于系统的时刻无关运算关系在整个运算过程中不随时间而变化,线性时不变系统既满足叠加原理，又满足时不变性的系统,50,线性系统,设系统的输入序列与输出分别为,可加性:如果系统的输入之和与输出之和满足,齐次性(或比例性):设a为常数，系统的输入增大a倍，输出也增大a倍,线性系统与非线性系统,51,例：证明一个线性系统,注意：必须证明系统同时满足可加性和齐次性，且信号及比例常数都可以是复数。,例1.7试分析下列系统的线性(1)y(n)=2x(n)-3，(2)y(n)=x(Mn)，其中M为正整数。,不满足叠加原理，非线性系统,满足叠加原理，线性系统,52,时不变系统,输入序列x(n)移动任意m位后，输出序列y(n)也移动m位，数值却保持不变。,m为任意常整数,时不变系统也称为移不变系统,53,例：证明一个时不变系统,例1.7试分析下列系统的时不变性(1)y(n)=2x(n)-3，(2)y(n)=x(Mn)，其中M为正整数。,二者相等，具有时不变性,时变系统,54,1.3.3单位脉冲响应与卷积和,单位取样响应(单位脉冲响应)h(n)=T(n)线性时不变系统输入为(n)时对应的输出,线性时不变系统都可以用它的单位脉冲响应h(n)来表征已知h(n)可得到线性时不变系统对任意输入的输出,55,推导卷积和表达式,(n)表示x(n),系统输出,叠加原理,时不变性,卷积和表达式:表示线性时不变系统的输出等于输入序列和单位脉冲响应的卷积。,56,1.3.4线性时不变系统的性质,交换律结合律分配律,可以推广到多个系统的情况，由卷积和的定义可以很容易加以证明。,57,1.3.5因果系统和稳定系统,因果系统系统某时刻的输出y(n)只取决于此时刻x(n)和以前的输入x(n-1)，x(n-2)，而和此时刻以后的输入x(n+1)，x(n+2)，无关。,先因后果因果系统的响应不会出现于外加输入之前。,非因果系统当前的输出还取决于未来的输入，不符合因果关系。,58,因果性的充分必要条件,线性时不变系统具有因果性的充要条件h(n)=0，n0,证明充分条件若n0时，h(n)=0，则,因而n0时刻的输出,可见，y(n0)只与mn0时的x(m)有关，因而是因果系统。,59,因果条件证明,证明利用反证法证明必要条件,假设因果系统，n0时h(n)0，则,在所设条件下，第二个求和式中至少有一项不为零，y(n)将至少和mn时的某一个x(n)值有关，这不符合因果性，假设不成立。,60,例：判断因果系统,例1.8判断差分系统的因果性。(1)前向差分系统:y(n)=x(n+1)-x(n);(2)后向差分系统:y(n)=x(n)-x(n-1)。,解因为前向差分系统的y(n)决定于x(n+1)，故系统为非因果的。,而后向差分系统定义为y(n)=x(n)-x(n-1)，显然是因果的。,61,讨论因果系统可实现性,因果系统是物理可实现的系统；非因果系统是不可实现的系统。,在具有较大延时的情况下，可以用因果系统去逼近非因果系统。例如语音处理、气象、地球物理学等。,非因果系统在理论上是存在的。例如，理想低通滤波器以及理想微分器都是非因果系统，但它们是不可实现的。,62,稳定系统,稳定系统系统的每个有界输入，对应产生的输出都有界。如果输入满足|x(n)|M+(M为正常数)，有输出|y(n)|P+(P为正常数)。,判断系统不稳定只要找出一个特别的有界输入，对应的输出是无界的，则该系统就是不稳定的。,判断系统稳定必须证明所有有界输入，其输出都是有界的。,63,稳定性的充分必要条件,线性时不变系统具有稳定性的充要条件是其单位脉冲响应绝对可和，即,证明充分条件若式成立，对于所有n都有|x(n)|M，得,即输出y(n)有界，系统不稳定。,64,稳定条件证明,证明利用反证法证明必要条件,假设系统稳定，但单位脉冲响应不绝对可和,定义一个有界输入,计算输出，有,即y(0)无界，系统不稳定，因此假设不成立。,65,例：判断稳定系统,例1.9判断累加器系统的稳定性,解考虑有界输入x(n)=u(n)，累加器的输出为,虽然n为有限值时，系统输出也为有限值，但对于所有n值(包括+)不存在有限值P，使得(n+1)P+，故系统输出无界。,66,例：判断因果稳定系统,例1.10已知线性时不变系统的单位脉冲响应,解因为n0时，u(-n-1)=1，所以h(n)0，故系统是非因果系统。,所以|a|1时系统稳定，|a|1时不稳定。,式中a为实常数，讨论其因果性和稳定性。,收敛序列：模值随n加大而减小，如|a|1时h(n)；发散序列：模值随n加大而加大，如|a|1时h(n)。,因为,67,1.4线性常系数差分方程,离散时间系统的数学模型差分方程,线性常系数线性差分方程求解,68,1.4.1离散时间系统的数学模型差分方程,差分方程是描述函数序列差分之间关系的方程，由序列及其各阶差分进行线性叠加组成。例如，对于一个二阶差分方程,将=1-D代入方程，得到,展开得到二阶线性常系数差分方程,69,线性常系数差分方程的一般形式,线性时不变系统的数学模型,式(1.44)不必是因果。假设是因果系统，变换得到,线性：x(n-r)和y(n-k)项都只有一次幂且不存在它们的相乘项，也没有相互交叉项,常系数：决定系统特征的系数均为常数,阶数：y(n-k)项变量k的最大值与最小值之差。,70,差分方程的方框图表示,理论上表示系统，也能在计算机上实现系统。例如，一阶差分方程,b0 x(n)表示将输入x(n)乘上常数b0-a1y(n-1)表示将序列y(n)延时一位后乘以常数-a1两个结果相加就得到y(n)序列图中代表相加器，代表乘法器，z-1代表延时一位的延时单元。,71,1.4.2线性常系数差分方程求解,差分方程的确定解不仅与差分方程的形式有关，而且还与其初始条件有关。差分方程求解实际上求系统的全响应零输入响应：y1(n)零状态响应：y2(n)全响应：全解y(n)=y1(n)+y2(n)求解差分方程时域求解法变换域求解法,72,例如递推法求特解,例1.13已知一个因果线性时不变系统的差分方程y(n)=ay(n-1)+x(n)，设初始条件y(n-1)=0，求系统的单位脉冲响应。,解令x(n)=(n)，于是有,由于系统具有因果性，递推如下,由此求出,h(n)=ah(n-1)+(n),73,1.5连续时间信号的取样,74,1.5.1理想取样,实际取样:=0;%生成阶跃序列stem(n,x);axis(-3,3,0,1.5);,89,矩形序列生成函数,functionx,n=rectseq(n0,n1,n2,N)%单位矩形序列生成函数%调用方式x,n=rectseq(n0,n1,n2,N)n=n0:n2;%生成位置向量x=(n-n1)=0%生成矩形脉冲序列,90,矩形序列,x,n=rectseq(-3,-1,4,5);stem(n,x);axis(-3,5,0,1.5);,91,实指数序列,n=0:10;%生成位置向量x=(0.6).n;%生成实指数序列stem(n,x);axis(0,10,0,1.5);,92,正弦序列3sin(0.1n+/3),n=0:1:20;%生成位置向量x=3*sin(0.1*pi*n+pi/3);%生成正弦序列stem(n,x);axis(0,20,-4,4);,93,复指数序列,n=-2:10;x=exp(0.2-0.5j)*n);%复指数序列subplot(1,2,1),stem(n,real(x);%用空心圆画点line(-5,10,0,0);%画横坐标subplot(1,2,2),stem(n,imag(x),filled);%用实心圆画点%line(-5,10,0,0),94,翻转:调用fliplr,n=-3:3;%生成一个序列x=0,0,1,0.5,0.25,0.125,0;stem(n,x);x=fliplr(x);%x排列次序左右翻转n=-fliplr(n);%向量n对n=0翻转stem(n,x);,95,序列的能量,conj求共轭复数sum求总和,E=sum(x.*conj(x);,abs求幅值sum求总和,E=sum(abs(x).2);,96,卷积和：调用conv,x=3,-3,7,0,-1,5,2;%序列x的非零区间-4n2h=2,3,0,-5,2,1;%序列x的非零区间-1n4,%调用conv计算卷积和y=conv(x,h);,运行结果：无位置信息,y=635619-313018-27-192,97,卷积和函数：convextd.m,functiony,ny=convextd(x,nx,h,nh)%序列y为序列x和序列h的卷积%ny，n