《现代数字信号处理》课件第1章.ppt_163文库

资源描述

1、第1章基础理论第1章基础理论1.1 概述1.2 离散时间信号与系统1.3 信号抽样、量化和编码1.4 基本的信号变换方法第1章基础理论1.1 概概述述1.1.1 信号信号通信的目的主要是传递信息，信号则是信息的载体。信号定义为随着时间、空间或者其他自变量变化的物理量。虽然信号可以用很多方法表示，例如光信号、声信号、电信号、磁信号等，但是在所有情况下，信息都可以包含在以某种方式变化的一个图形中，各种信号可表示为一个或几个独立变量的函数。第1章基础理论根据时间变量的取值方式不同，信号可以分为连续时间信号和离散时间信号。时间取值连续的为连续时间信号；不连续的为离散时

2、间信号。根据信号的周期性，可以分为周期信号和非周期信号。对于信号x(n)，若有x(n)=x(n+kN),k和N均为正整数，则称x(n)为周期信号，周期为N；否则，称x(n)为非周期信号。当然，一个非周期信号也可以视为周期信号，其周期认为是无穷大的。根据信号在任意时刻的取值是否能够精确确定，可以分为确定性信号和随机信号。顾名思义，随机信号在任意时刻的取值是随机的，不能给出精确的预测。第1章基础理论根据信号的能量和功率大小，可以将信号分为能量信号和功率信号。如果信号能量有限，则称为能量信号；如果信号功率有限，则称为功率信号。周期信号、准周期信号及随机信号，由于其时间是无限的，所以它们总是功

3、率信号。一般认为在有限区间内存在的确定性信号有可能是能量信号。根据信号变量的个数，可以将信号分为一维信号、二维信号和多通道信号。第1章基础理论1.2 离散时间信号与系统离散时间信号与系统1.2.1 离散时间信号离散时间信号一个信号x(t)，可以代表一个实际的物理信号，也可以是一个数学函数。例如，x(t)=A sin(2ft)既是正弦信号，又是正弦函数。因此，在信号处理中，信号与函数往往是通用的。随机信号与随机过程也是通用的。在函数中，如果t是时间轴上的连续变量，那么称x(t)为连续时间信号，又称模拟信号。如果t仅在时间轴的离散点上取值，那么称其为离散时间信号，记为x(n)。第1章基

4、础理论x(n)在时间上是离散的，其幅度可以在某一个范围内连续取值。但是目前的信号处理装置是以计算机或专用信号处理芯片来实现的，它们都以有限的位数来表示幅度，因此，其幅度也要量化，即取离散值。在时间和幅度上都取离散值的信号称为数字信号。第1章基础理论1.2.2 典型的离散信号典型的离散信号在离散时间信号中，我们经常遇到一些典型的离散信号，下面对它们进行简要介绍。(1)单位抽样序列:0 ,0 0 ,1)(nnn(1.2.1)(n)为单位抽样序列，该序列在离散信号和离散系统的分析与综合中有重要的作用，其地位犹如单位冲激信号(t)在连续时间系统中的地位。(n)在数字信号处理中也起着重要的作

5、用，它的许多性质与单位冲激信号相似，只是定义不同。(t)又称Kronecker函数，是以积分为基础定义的，即第1章基础理论(1.2.2)1d)(tt且t0时，(t)=0。(t)表示在极短的时间内产生的巨大冲激，而(n)在n=0时定义为1。图 1.1(a)和(b)分别给出了(n)和(t)的波形。第1章基础理论图 1.1 单位抽样信号和单位冲激信号第1章基础理论(2)单位阶跃序列:(1.2.3)0 ,0 0 ,1)(nnnu该序列只在时间轴的正半轴上取值，且恒为数值1。若序列y(n)=x(n)u(n)，那么y(n)的自变量n的取值就限定在n0的右半轴上。x(n)=A sin

6、(2fnTs+j)(1.2.4)式中:f是频率，单位是Hz。令=2f，则的单位是rad/s,是相对连续信号x(t)的连续角频率变量，当f由至时，也由至，令(1.2.5)ssss1 ,22TffffT第1章基础理论fs称为抽样频率。显然，当f从0增至fs时，从0增至2；当 f 从0减至fs时，从0减至2；当f再增加或减少fs的整数倍时，重复02。的单位是rad，我们称为圆频率或圆周频率，它是相对离散信号x(n)的频率分量。于是式(1.2.4)的正弦序列可以表示为x(n)=A sin(n+j)(1.2.6)(4)复正弦序列:x(n)=ejn=cos(n)+j sin(n)(1.2.7)上式

7、称为欧拉恒等式。复正弦ejn在数字信号处理中有着重要的应用，它不但是离散信号做傅立叶变换的基函数，同时也是离散系统的特征函数。第1章基础理论(5)指数序列:x(n)=a|n|(1.2.8)式中，a为常数且|a|0，00，这样，x(n)变成复值信号，即x(n)=r|n|ej0|n|。若r0)到n内的输入样本确定，那么该系统称为具有长为N的存储。如果N=0，那么系统是静态的；如果0N，那么系统称为有限存储系统；如果N=，那么系统称为无限存储系统。第1章基础理论由以下输入输出方程y(n)=ax(n)(1.2.12)y(n)=nx(n)+bx3(n)(1.2.13)所描述的系统都是静态

8、的或者无记忆的。对于这种系统，在计算当前输出时，无须存储任何过去的输入或输出值。另一方面，由以下输入输出关系y(n)=x(n)+3x(n1)(1.2.14)(1.2.15)nkknxny0)()(0)()(kknxny(1.2.16)第1章基础理论所描述的系统是动态系统或者是具有存储器的系统。式(1.2.14)和式(1.2.15)描述的系统属于有限存储系统，而式(1.2.16)描述的系统属于无限存储系统。我们注意到，静态或者无记忆系统通常由形式为y(n)=x(n)，n(1.2.17)的输入输出方程来描述，这里并不包括延迟元件(存储器)。第1章基础理论2.时不变系统与时变系统时不

9、变系统与时变系统可以把一般的系统类别再分为两大类:时不变系统和时变系统。如果系统的输入输出特性不随着时间变化，那么该系统称为时不变系统。为了详细说明这一点，假设有一个弛豫系统，当它被输入信号x(n)激励时，产生输出信号y(n)。因此，可以写为y(n)=x(n)(1.2.18)第1章基础理论定义定义1.1 一个弛豫系统是时不变或者平移不变系统，当且仅当这意味着，对任何输入信号x(n)和任意时间平移量k，有为了确定一个给定系统是否为时不变系统，需要进行前面定义所指定的测试。基本方法是，用一个任意的输入序列x(n)激励系统，产生一个输出序列,记为y(n)。下一步，将输入序列延迟k个单位并重新

10、计算输出。通常，输出可以写为y(n，k)=x(nk)()(nynx)()(knyknx第1章基础理论如果对于所有可能的k值，输出y(n，k)=x(nk)，那么该系统就是时不变的。另一方面，如果即使只对一个k值，输出y(n，k)x(nk)，那么该系统就是时变的。第1章基础理论例例 1.3 确定图 1.10所示的系统是时变还是时不变的。图 1.10 系统判断实例第1章基础理论解解:(a)这个系统的输入-输出方程为y(n)=x(n)=x(n)x(n1)现在，如果输入在时间上延迟了k个单位，并将它作用到系统，那么从结构图中显然可以看出该系统的输出为y(n，k)=x(nk)x(n

11、k1)另一方面，如果将y(n)从时间上延迟k个单位，那么就有y(nk)=x(nk)x(nk1)因此该系统是时不变的。同理可以判断(b)、(c)、(d)是时变系统。第1章基础理论3.线性系统与非线性系统线性系统与非线性系统一般的系统类别还可再分为线性系统和非线性系统。满足叠加性的系统就是线性系统。简单说明一下，叠加性要求系统对信号加权和的响应等于系统对每个独立输入信号对应的响应(输出)的加权和。因此，得到以下的线性定义。定义定义1.2 一个系统是线性的，当且仅当对任意输入序列x1(n)和x2(n)，以及任意常数a1和a2，满足式(1.2.20)。a1x1(n)+a2x2(n)=a1x1(

12、n)+a2x2(n)(1.2.20)图 1.11给出了叠加性的图形表示。第1章基础理论图 1.11 叠加性的图形表示第1章基础理论关系式(1.2.20)表现的叠加性可以分成两部分。第一，假设a2=0，那么式(1.2.20)简化为a1x1(n)=a1x1(n)=a1y1(n)(1.2.21)其中，y1(n)=x1(n)。关系式(1.2.21)展示了线性系统的乘法或者缩放特性。即如果系统对输入x1(n)的响应为y1(n)，那么系统对a1x1(n)的响应就仅仅是a1y1(n)。因此，对输入的任意缩放会导致相应输出的等同缩放。第二，假设在式(1.2.20)中有a1=a2=1，那么 x1

13、(n)+x2(n)=x1(n)+x2(n)=y1(n)+y2(n)(1.2.22)第1章基础理论这个关系式展示了线性系统的可加性。可加性和乘法性组成了适用于线性系统的叠加性。式(1.2.20)所包含的线性条件，通过归纳，可以扩展到信号的任意加权组合。通常，可以得到(1.2.23)()()()(1111nyanynxanxMkkkMkkk其中 yk(n)=xk(n)，k=1，2，M1(1.2.24)从式(1.2.21)可以看到，如果a1=0，那么y(n)=0。换言之，一个弛豫线性系统，零输入产生零输出。如果零输入产生非零输出，那么该系统可能是非弛豫的或者是非线性的。如果一个弛豫系统不满足

14、定义所给出的叠加性，那么就称它为非线性的。第1章基础理论例例 1.4 判断下列输入-输出方程所描述的系统是线性的还是非线性的。(a)y(n)=nx(n);(b)y(n)=x(n2);(c)y(n)=x2(n);(d)y(n)=Ax(n)+B;(e)y(n)=ex(n)。解解:(a)对于两个输入序列x1(n)和x2(n)，对应的输出为 y1(n)=nx1(n)y2(n)=nx2(n)第1章基础理论两个输入序列的线性组合产生的输出为y3(n)=na1x1(n)+a2x2(n)=a1nx1(n)+a2nx2(n)另一方面，两个输出的线性组合产生的输出为a1nx1(n)+a2nx2(n

15、)所以该系统是线性的。同理可以判断，(b)是线性的，(c)、(e)是非线性的，(d)不满足判断条件，即系统不是弛豫的。第1章基础理论4.因果系统与非因果系统因果系统与非因果系统我们先从因果离散时间系统的定义开始。定义定义1.3 如果一个系统在任意时刻n的输出(即y(n)仅依赖于当前和过去的输入(即x(n)，x(n1),x(n2)，)，而与将来的输入(即x(n+1)，x(n+2)，)无关，那么这个系统就称为是因果的。以数学形式表示，因果系统的输出应该满足如下方程:y(n)=Fx(n)，x(n1)，x(n2)，(1.2.25)其中，F是某个任意函数。第1章基础理论例例 1.5 判断

16、由下列输入-输出方程描述的系统是因果的还是非因果的。(a)y(n)=x(n)x(n1);(b);(c)y(n)=ax(n);(d)y(n)=x(n)+3x(n+4);(e)y(n)=x(n2);(f)y(n)=x(2n);(g)y(n)=x(n)。nkkxny)()(第1章基础理论解解:(a)、(b)、(c)所描述的系统很明显是因果的，因为它们的输出只依赖于当前和过去的输入。(d)、(e)、(f)所描述的系统很明显是非因果的，因为它们的输出依赖于将来的输入。(g)所描述的系统也是非因果的，我们注意到，当选择n=1时，得到的输出为y(1)=x(1)。因此，在n=1时的输出依赖于n=1时的

17、输入，在时间上这是将来两个单位后的值。第1章基础理论5.稳定系统与不稳定系统稳定系统与不稳定系统定义定义1.4 一个任意的弛豫系统称为有界输入有界输出(BIBO)稳定，当且仅当每个有界输入产生有界的输出。输入序列x(n)和输出序列y(n)有界的条件转换到数学上，就意味着存在某些有限数，如Mx和My，并且对所有的n，使得|x(n)|Mx，|y(n)|My(1.2.26)如果对某些有界输入序列x(n)，输出是无界(无限)的，那么该系统归为不稳定系统。第1章基础理论1.2.7 离散时间系统的互连离散时间系统的互连离散时间系统可以相互连接以构建更大的系统。有两种基本的系统互连方法:级联

18、(串联)或者并联。这些互连如图 1.12所示。注意，这两个互连的系统是不同的。图 1.12 系统的级联和并联第1章基础理论下面证明c是时不变的。为了证明时不变性，假设1和2是时不变的，那么)()(11knyknx并且)()(21knykny)()(12knyknxc因此因而c是时不变的。第1章基础理论1.2.8 离散时间线性时不变系统的分析离散时间线性时不变系统的分析在1.2.6节中，我们依据许多特性或者范畴对系统进行了分类，即线性、因果、稳定和时不变。有了这样的分类，就可集中分析一类重要的线性时不变(LTI)系统。特别地，我们将会说明这类系统的时域特征是通过它们对单位采样序列

19、的响应来进行描述的。我们还会说明，任意输入信号都可以分解并表示为单位采样序列的加权和。因此，由于系统的线性时不变性，系统对任意输入信号的响应可以表示为单位采样响应的形式。将系统的单位采样响应、任意的输入信号与输出信号关联起来的一般表达形式，称为卷积和或者卷积公式，我们也将推导这个公式。因此，我们就可以计算线性时不变系统对任意输入信号的输出。第1章基础理论1.2.9 线性系统的分析方法线性系统的分析方法有两种分析线性系统对给定输入信号反应或者响应的方法。一种是基于系统输入输出方程的直接求解方法，通常输入输出方程的形式为y(n)=Fy(n1)，y(n2)，y(nN)；x(n)，x(n1)，

20、x(nM)其中，F在方括号内表示各个量的某些函数。特别地，对LTI系统，在后面我们将看到，它的输入输出关系的一般形式为MkkNkkknxbknyany01)()()(第1章基础理论第二种分析线性系统对给定输入信号响应的方法是:首先将输入信号分解成基本信号的和，基本信号的选择要使得系统对每个信号分量的响应易于求解；然后利用系统的线性特性，将系统对基本信号的响应相加，就得到了系统对给定输入信号的总响应。第二种方法即是本节将讲述的方法。为了详细说明，假设输入信号x(n)分解为基本信号分量xk(n)的加权和，使得(1.2.31)kkknxcnx)()(其中，ck是信号x(n)分解后的一组幅度(

21、权系数)。现在，假设系统对基本信号分量xk(n)的响应为yk(n)，从而yk(n)xk(n)(1.2.32)第1章基础理论假设系统是弛豫的，由于线性系统的缩放特性，于是系统对ckxk(n)的响应为ckyk(n)。系统对输入x(n)的总响应为(1.2.33)kkknxcnxny)()()(kkknxc)(kkknyc)(=在式(1.2.33)中，我们利用了线性系统的可加性。第1章基础理论虽然基本信号的选择看似是任意的，但依据是对想要求解的输入信号的分类。如果对输入信号的特性不做限制，那么将输入信号分解为单位采样(冲激)序列的加权和可以证明是运算方便且十分普遍的。另一方面，如果我们

22、的讨论限制在输入信号的一个子集上，那么或许存在其他的基本信号集，使得求解输出的运算更加方便。例如，如果输入信号x(n)是周期为N的周期信号，那么运算方便的基本信号集是指数信号xk(n)=ejkn，k=0，1，N1(1.2.34)其中，频率k是谐波相关的，即(1.2.35),2kNwk1,1,0Nk第1章基础理论1.2.10 离散时间信号分解为冲激信号离散时间信号分解为冲激信号假设对于任意的输入信号x(n)，我们希望将它分解成单位采样序列的和。为了能够利用前一节建立的符号，选择基本信号xk(n)为xk(n)=(nk)(1.2.36)其中，k表示单位采样序列的延迟。为了处理具有无限长非零值

23、的任意信号x(n)，该组的单位冲激也必须是无限的，以构成无限数量的延迟。第1章基础理论现在假设将两个序列x(n)和(nk)相乘，因为(nk)除了在n=k处的值为单位1外，其他处的值均为零，所以乘法的结果是另外一个序列，它除了在n=k处的值为x(k)外，其他处的值均为零，正如图 1.13所示。因此x(n)(nk)=x(k)(nk)(1.2.37)是一个序列，它除了在n=k处的值为x(k)外，其他处的值均为零。如果重复计算x(n)和(nm)的乘积，其中m是另一个延迟量(mk)，那么这个结果也是一个序列，它除了在n=m处的值为x(m)外，其他处的值均为零。因此x(n)(nm)=x(m)(nm

24、)(1.2.38)第1章基础理论图 1.13 信号x(n)和移位单位采样序列的乘积第1章基础理论换言之，信号x(n)和单位冲激在某个延迟k(即(nk)处的每个乘积，其实就是取信号x(n)在单位冲激的非零处的延迟值x(k)。因此，如果要重复计算所有可能延迟k的乘积，kM。这个结构是递归系统(1.2.48)()()(1nxknwanwNkk及其后的非递归系统(1.2.49)Mkkknwbny0)()(的级联。第1章基础理论图 1.16 式(1.2.45)描述的系统的直接型结构第1章基础理论我们看到，如果NM，那么该结构需要的延迟器的数量等于系统的阶N。然而，如果MN

25、，那么需要的存储器由M指定(可以修改图 1.16以解决这种情况)。因此，直接型结构需要M+N1个乘法器和max(M，N)(M，N的最大值)个延迟器。因为要实现式(1.2.45)描述的系统需要最少数量的延迟器，所以该结构有时也称为规范形式。如果令系统参数ak=0，k=1，N，就会出现式(1.2.45)的特殊情况，那么系统的输入输出关系就简化为(1.2.50)Mkkknxbny0)()(第1章基础理论这是一个非递归线性时不变系统，该系统可以视为只将最近的M+1个输入信号样本通过bk集中的适当的系数bk加权后再相加。换言之，系统的输出基本上是输入信号的加权滑动平均。因此，该系统有时也称为滑动

26、平均(MA)系统。这样的系统是冲激响应hk等于系数bk的FIR系统，即(1.2.51)其他 ,00 ,)(Mkbkhk如果再回到式(1.2.45)并令M=0，那么一般线性时不变系统就简化为完全递归系统，它的差分方程为(1.2.52)()()(01nxbknyanyNkk第1章基础理论由二阶差分方程描述的线性时不变系统是由式(1.2.45)、式(1.2.49)或式(1.2.52)描述的更具一般系统特征的一个重要子类。在后面讨论量化效果时，我们会解释这些系统重要性的原因。可以肯定地说，迄今为止，二阶系统通常是用来实现高阶系统的基本运算单元。描述最一般二阶系统的差分方程为 y(n)=a1y(

27、n1)a2y(n2)+b0 x(n)+b1x(n1)+b2x(n2)(1.2.53)这是由式(1.2.45)通过令N=2以及M=2得到的。实现这种系统的直接型结构如图 1.17(a)所示。如果令a1=a2=0，那么式(1.2.53)就简化为第1章基础理论y(n)=b0 x(n)+b1x(n1)+b2x(n2)(1.2.54)这是式(1.2.50)描述的FIR系统的一个特殊情况，实现该系统的结构如图 1.17(b)所示。如果令式(1.2.53)中b1=b2=0，那么就得到完全二阶递归系统，描述它的差分方程为y(n)=a1y(n1)a2y(n2)+b0 x(n)(1.2.55)这是式(1.

28、2.52)的一种特殊情况。实现该系统的结构如图 1.17(c)所示。第1章基础理论图 1.17 二阶系统的实现结构第1章基础理论1.2.13 FIR系统的递归和非递归实现系统的递归和非递归实现基于系统的冲激响应h(n)是有限长的还是无限长的，我们已经能够区分FIR和IIR系统。我们还能够区分递归和非递归系统。基本上，因果递归系统是由形式为y(n)=Fy(n1)，y(nN)；x(n)，x(nM)(1.2.56)的输入输出方程来描述的；而特殊的线性时不变系统是由差分方程(1.2.57)MkkNkkknxbknyany01)()()(第1章基础理论来描述的。另一方面，因果非递

29、归系统并不依赖于输出的过去值，所以它是由形式为y(n)=Fx(n)，x(n1)，x(nM)(1.2.58)的输入输出方程来描述的；而对于特殊的线性时不变系统，它是由式(1.2.57)中的差分方程来描述的，其中ak=0，k=1，2，N。在FIR系统的情况下，我们已经看到，通常可以非递归地实现这些系统。实际上，令式(1.2.57)中的ak=0，k=1，2，N，就得到这样一个非递归FIR系统，它的输入输出方程为(1.2.59)Mkkknxbny0)()(第1章基础理论正如式(1.2.51)所指出的，该系统的冲激响应只等于系数bk。因此，每个FIR系统可以非递归地实现。另一方面，任何FIR系统

30、也可以递归地实现。虽然我们在后面会给出这个陈述的一般证明，但现在我们将会以一个简单的例子来说明这一点。假设有一个形式为(1.2.60)MkknxMny0)(11)(的FIR系统，现在要计算输入信号x(n)的滑动平均。很明显，该系统是一个冲激响应为(1.2.61)MnMnh0,11)(第1章基础理论的FIR系统。图 1.18画出了该系统的非递归实现结构。现在，假设将式(1.2.60)表示为(1.2.62)MkknxMny0)1(11)()1()(11MnxnxM)1()(11)1(MnxnxMny此时，式(1.2.61)表示了该FIR系统的递归实现。这个滑动平均系统的递归实现结构如图1.

31、19所示。第1章基础理论图 1.18 FIR滑动平均系统的非递归实现第1章基础理论图 1.19 FIR滑动平均系统的递归实现第1章基础理论1.3 信号抽样、量化和编码信号抽样、量化和编码1.3.1 抽样抽样抽样是连续时间信号到离散时间信号的转换过程，通过对连续时间信号在离散时间点处取样本值获得。因此，如果x(t)是抽样器的输入，那么输出是x(nT)=x(n)，其中T称为抽样间隔。抽样器可以看成是一个电子开关，开关每隔T秒闭合一次(对于理想抽样，闭合时间应该无穷短，对于实际抽样，闭合时间是秒，T)，使输入信号得以抽样，得到连续信号的抽样输出信号。第1章基础理论1.

32、3.2 理想抽样的抽样定理理想抽样的抽样定理在0的极限情况(当T时，就可认为是理想抽样)下，抽样脉冲序列p(t)变成冲激函数序列T(t)，而各冲激函数准确地出现在抽样瞬间上，面积为1，抽样后输出理想抽样信号的面积(即积分幅度)则准确地等于抽样信号xa(t)在抽样瞬间的幅度。理想抽样信号过程如图 1.20所示。第1章基础理论图 1.20 连续时间信号抽样第1章基础理论冲激函数序列T(t)为(1.3.1)nmTmTtt)(理想抽样输出为()ax t(1.3.2)()()()aaTx tx tT将式(1.3.1)代入式(1.3.2)得(1.3.3)()()()aaTmx tx ttm

33、T根据(t)函数的性质，可以得到(1.3.4)()()()aaTmx tx mttmT第1章基础理论下面讨论理想抽样信号后信号频谱发生的变化。由于在连续时间信号与系统中已经学过，时域相乘则频域为卷积运算，因此若各个信号的傅立叶变换分别表示为Xa(j)=DTFTxa(t)T(j)=DTFTT(t)(j)DTFT()aaXxt 其中，DTFT表示离散时间信号的傅立叶变换，则有(1.3.5)1(j)(j)(j)2aTaXX Xa(j)可表示为(1.3.6)j(j)DTFT()()edtaaaXxtxtt 第1章基础理论现在来求T(j)=DTFTT(t)。由于T(t)是周期函数，可表示

34、成傅立叶级数，即 ktkkTsAtje)(称此级数的基频为抽样频率，即,1TfsTs2而系数可表示成jj2222j2211()ed()ed11 ()edsssTTktktTTkTmTktTAtttmTtTTttTT第1章基础理论以上结果考虑到在|t|T/2的区间内，只有一个冲激(t)，而当m0时，(tmT)都在积分区间之外，且利用关系式(0)()()dff ttt因而有ktkTsTtje1)(1.3.7)由此得出jj11(j)DTFT()DTFTeDTFTessktktTTkktTT 由于DTFTejkst=2(ks)(1.3.8)第1章基础理论所以图 1.21 周期冲激序列T

35、(t)和它的傅立叶变换T(j)(1.3.9)sss2(j)DTFT()()()TTkktkkT T(t)和T(j)的图形如图 1.21所示。第1章基础理论根据时域乘积等于频域卷积可以算出(1.3.10)12(j)jjaakXXkTT 由此看出，一个连续时间信号经过理想抽样后，其频谱将以抽样频率s=2/T为间隔重复，即频谱产生了周期延拓，如图 1.22所示，因为频谱是复数，故只画出了其幅度。也就是说，理想抽样信号的频谱是频率的周期函数，其周期为s，而频谱的幅度则受1/T加权。由于T是常数，所以除了一个常数因子区别外，每一个延拓的谱分量都和原频谱分量相同。因此只要各延拓分量与原频谱分量不发

36、生频率上的交叠，就有可能恢复出原信号。也就是说，如果xa(t)是带限信号，其频谱如图 1.22(a)所示，且最高分量h不超过s/2，那么原信号的频谱和各次延拓分量的频谱彼此不重叠，如图1.22(b)所示。第1章基础理论这时采用一个截止频率为s/2的理想低通滤波器，就可以得到不失真的原信号频谱，也就是说，可以不失真地恢复出原来的连续信号，即有(1.3.11)ss1j,2(j)0,2aaXTX 如果信号的最高频谱h超过s/2，则各周期延拓分量产生频谱的交叠，称为频谱的混叠现象，如图 1.22(c)所示。由于Xa(j)一般是复数，所以混叠现象也是复数相加。我们将抽样频率之半称为折叠频率，即(

37、1.3.12),212sTfT2s第1章基础理论图 1.22 频谱的周期延拓第1章基础理论由此可以得出结论:若xa(t)的频带宽度是有限的，要想抽样后x(n)=xa(nT)能够不失真地还原出原信号xa(t)，则抽样频率必须大于或等于两倍信号频谱的最高频率，这就是奈奎斯特抽样定理，即fs2fh(1.3.13)若xa(t)的频带宽度不是有限的，为了避免混叠，一般在抽样器前加一个保护性的前置低通滤波器，称为防混叠滤波器，其截止频率为fs/2，可以滤除高于fs/2的频谱分量。用同样的方法可以证明，理想抽样后，信号的拉普拉斯变换在s平面上沿虚轴周期延拓，也就是说在s平面虚轴上是周期函数，

38、即()aXs第1章基础理论(1.3.14)s1()(j)aakXsXskT其中()()edstaaXsxtt()()edstaaXsxtt第1章基础理论1.3.3 实际抽样实际抽样实际情况中，抽样脉冲不是冲激函数，而是有一定宽度的矩形脉冲p(t)，这时奈奎斯特抽样定理是否有效？由于p(t)是周期函数，故仍可以展成傅立叶级数(1.3.15)sj()ektkkp tC同样可以求出p(t)的傅立叶系数Ck,假定p(t)的幅度为1,则有(1.3.16)2jss0j22je22sinde1de1sssktkTTtkkkkTtTttpTC第1章基础理论如果、T一定，则随着k的变化，

39、Ck的幅值|Ck|将按xxkksssin22sin而变化，其中x=ks/2。进而可以求出实际抽样时的信号的频谱为(1.3.17)s(j)(jj)akakXC Xk 由此看出，和理想抽样一样，抽样数据信号的频谱是连续信号频谱的周期延拓，因此，如果满足奈奎斯特抽样定理，则不会产生频谱的混叠失真。和理想抽样不同的是，这里频谱分量的幅度有变化，其包络是随着频率增加而逐渐下降的，如图 1.23所示。第1章基础理论图 1.23 实际抽样时的包络变化第1章基础理论由于s22siness2jkkkkTC又由于包络的第一个零点出现在022sinsskk这就要求2sk第1章基础理论所以Tk

40、第1章基础理论1.3.4 正弦信号的抽样正弦信号的抽样连续时间正弦信号是一种很重要的信号，不管是理论上还是信号处理的实际应用中，它都有着广泛的应用。例如，常用正弦信号加白噪声作为输入信号来研究某一实际系统或某一算法的性能。因此，正弦信号的抽样显得十分重要。设连续时间正弦信号为x(t)=A sin(0t+j)=A sin(2f0t+j)(1.3.18)第1章基础理论对于两个不同频率的正弦信号x1(t)、x2(t)，如果用同一抽样频率对其抽样，抽样出的序列可能是一样的，则我们无法判断它来源于x1(t)还是x2(t)。例如设x1(t)、x2(t)分别为x1(t)=cos(240t)

41、x2(t)=cos(2140t)如果用fs=100 Hz对两个信号进行抽样，可以看出，x1(t)的抽样满足抽样定理，x2(t)的抽样不满足抽样定理，抽样后的序列分别为54cos402cos)(s1nfntx54cos514cos1402cos)(s2nnfntx第1章基础理论它们都是5点的周期序列，其基本周期内的序列值分别为1，0.809，0.309，0.399，0.809,这时我们无法判断这个序列来自x1(t)还是x2(t)。第1章基础理论1.4 基本的信号变换方法基本的信号变换方法1.4.1 Z变换的定义与收敛域变换的定义与收敛域1.Z变换的定义变换的定义若序列为x(n)，

42、则幂级数(1.4.1)nnznxzX)()(称为序列x(n)的Z变换，其中z为变量。亦可将x(n)的Z变换表示为Zx(n)=X(z)(1.4.2)第1章基础理论2.Z变换的收敛域变换的收敛域显然，只有当式(1.4.1)的幂级数收敛时，Z变换才有意义。对任意给定序列x(n)，使其Z变换(即式(1.4.1)收敛的所有z值的集合称为X(z)的收敛域。按照级数理论，式(1.4.1)的级数收敛的充要条件是满足绝对可和的条件，即要求(1.4.3)Mznxnn|)(|要满足此不等式，|z|值必须在一定范围之内才行，这个范围就是收敛域。不同形式的序列，其收敛域形式不同。下面以有限长序列为例进行说明。第

44、，表示为x(n)=Z1X(z)(1.4.5)由式(1.4.1)可看出，式(1.4.5)实质上是求X(z)的幂级数的展开式。求Z反变换的方法通常有三种：围线积分法(留数法)、部分分式展开法和长除法。这里以围线积分法为例进行介绍。围线积分法是求Z反变换的一种有用的分析方法。根据复变函数理论，若函数X(z)在环状区Rx|z|Rx+(Rx0,Rx+)内是解析的，则在此区域内X(z)可以展开成罗朗级数，即(1.4.6)xxnnnRzRzCzX|,)(第1章基础理论而(1.4.7),2,1,0,d)(j211nzzzXCncn其中围线c是在X(z)的环状解析域(即收敛域)内环绕原点的一条逆时针方向

45、的闭合单围线，如图1.25所示。将式(1.4.7)与式(1.4.1)的Z变换定义相比较可知，x(n)就是罗朗级数的系数Cn，故式(1.4.7)可写成(1.4.8),(,d)(j21)(1xxncRRczzzXnx式(1.4.8)就是用围线积分的Z反变换公式。第1章基础理论图 1.25 围线积分路径第1章基础理论式(1.4.8)的正确性可以用以下的柯西积分定理来证明:为整数kkkRRzzkkkckck,0,00,1de2dReej21dj21jj)1j(11其中，z=Rej，RxRRx+。将式(1.4.8)右端写成zzmxzzzmxzzzXcmnmncmmcndj21)(d)(j

46、21d)(j211)(11 若将此式中的nm看成式(1.4.9)中的k，则有)(dj21)(1)(nxzzmxcmnm第1章基础理论即),(),(d)(j211xxcnRRcnxzzzX直接计算围线积分比较麻烦，一般都采用留数定理来求解。按留数定理，若函数F(z)=X(z)zn1 在围线c上连续，在c以内有K个极点zk，而在c以外有M个极点zm，(M、K为有限值)，则有(1.4.10)kzzncnkzzXzzzX11)(Resd)(j21或(1.4.11)mzzncnmzzXzzzX11)(Resd)(j21第1章基础理论式(1.4.11)应用的条件是X(z)zn1在z=有二阶

47、或二阶以上零点，即要分母多项式z的阶次比分子多项式z的阶次高二阶或二阶以上。其中符号ResX(z)zn1z=z表示函数F(z)=X(z)zn1在点z=zk(c以内极点)的留数。式(1.4.10)说明，函数F(z)沿围线c逆时针方向的积分等于F(z)在围线c内部各极点的留数之和；式(1.4.11)说明，函数F(z)沿围线c顺时针方向的积分等于F(z)在围线c外部各极点的留数之和。由于cczzFzzFd)(d)(k(1.4.12)所以由式(1.4.10)及式(1.4.11)可得(1.4.13)mkzznmzznkzzXzzX)(Res)(Res11第1章基础理论将式(1.4.10)及式(1

48、.4.13)分别代入式(1.4.8)，可得(1.4.14a)kzznknczzXzzzXnx)(Resd)(j21)(11(1.4.14b)mzznmnczzXzzzXnx)(Resd)(j21)(11同样，应用式(1.4.14b)，必须满足X(z)zn1的分母多项式z的阶次比分子多项式z的阶次高二阶或二阶以上。第1章基础理论根据具体情况，既可以采用式(1.4.14a)，也可以采用式(1.4.14b)。例如，如果当n大于某一值时，函数X(z)zn1在z=处，也就是在围线外的外部可能有多重极点，这时选c的外部极点计算留数就比较麻烦，而选c的内部极点求留数则较简单。如果当n小于某值时，X(

49、z)zn1在z=0处，也就是在围线的内部可能有多重极点，这时选用c外部的极点求留数就方便得多。第1章基础理论现在来讨论如何求X(z)zn1在任一极点zr处的留数。如果zr是X(z)zn1的单(一阶)极点，则有ResX(z)zn1z=zr=(zzr)X(z)zn1z=z(1.4.15)如果zr是X(z)zn1的多重(l阶)极点，则有r(1.4.16)rrzznlrlzznzzXzzllzzX)()()!1(d)!1(1)(Res111第1章基础理论1.4.3 傅立叶变换的几种可能形式傅立叶变换的几种可能形式1.连续时间、连续频率连续时间、连续频率傅立叶变换傅立叶变换这就是连续时间

50、非周期信号x(t)的傅立叶变换关系，所得到的是连续的非周期的频谱密度函数X(j)。在“信号与系统”课的内容中，已知这一变换对为(1.4.17)ttxXtde)()(jj(1.4.18)Xtxtde)(j21)(j这一变换对的示意图(只说明关系，并不代表实际的变换对)见图1.26。由图可以看出，时域连续函数造成频域是非周期的谱，而时域非周期函数造成频域是连续的谱密度函数。第1章基础理论图 1.26 连续时间连续频率的谱线第1章基础理论2.连续时间、离散频率连续时间、离散频率傅立叶级数傅立叶级数设x(t)代表一个周期为T0的周期性连续时间函数，x(t)可展开成傅立叶级数，其傅立叶级

展开阅读全文