数字信号处理 第1章时域离散信号和时域离散系统

1、第1章时域离散信号和时域离散系统 第第1章章 时域离散信号和时域离散系统时域离散信号和时域离散系统 1.1引言引言 1.2时域离散信号1.3时域离散系统1.4时域离散系统的输入输出描述法线性常系数差分方程1.5模拟信号数字处理方法习题与上机题第1章时域离散信号和时域离散系统 1.1引言信号通常是一个自变量或几个自变量的函数。如果仅有一个自变量，则称为一维信号；如果有两个以上的自变量，则称为多维信号。本书仅研究一维数字信号处理的理论与技术。物理信号的自变量有多种，可以是时间、距离、温度、位置等，本书一般把信号看做时间的函数。针对信号的自变量和函数值的取值情况，信号可分为以下三种。第1章时域离散信

2、号和时域离散系统 如果信号的自变量和函数值都取连续值，则称这种信号为模拟信号或者称为时域连续信号，例如语言信号、温度信号等；如果自变量取离散值，而函数值取连续值，则称这种信号称为时域离散信号，这种信号通常来源于对模拟信号的采样; 如果信号的自变量和函数值均取离散值，则称为数字信号。我们知道，计算机或者专用数字信号处理芯片的位数是有限的，用它们分析与处理信号，信号的函数值必须用有限位的二进制编码表示，这样信号本身的取值不再是连续的，而是离散值。这种用有限位二进制编码表示的时域离散信号就是数字信号，因此，数字信号是幅度量化了的时域离散信号。第1章时域离散信号和时域离散系统 例如：，这是一个模拟信号

3、，如果对它按照时间采样间隔T=0.005s进行等间隔采样，便得到时域离散信号x(n)，即 =， 0.0， 0.6364， 0.9， 0.6364， 0.0， -0.6364， 0.9， -0.6364， 显然, 时域离散信号是时间离散化的模拟信号。 ) 50sin(9 . 0)()(anTtxnxnTt) 50sin(9 . 0)(attx第1章时域离散信号和时域离散系统 如果用四位二进制数表示该时域离散信号，便得到相应的数字信号xn，即xn=，0.000，0.101，0.111，0.101，0.000，1.101，1.111，1.101，显然，数字信号是幅度、时间均离散化的模拟信号，或者说是

4、幅度离散化的时域离散信号。 信号有模拟信号、时域离散信号和数字信号之分，按照系统的输入输出信号的类型，系统也分为模拟系统、时域离散系统和数字系统。当然，也存在模拟网络和数字网络构成的混合系统。 第1章时域离散信号和时域离散系统 数字信号处理最终要处理的是数字信号，但为简单，在理论研究中一般研究时域离散信号和系统。时域离散信号和数字信号之间的差别，仅在于数字信号存在量化误差，本书将在第9章中专门分析实现中的量化误差问题。 本章作为全书的基础，主要学习时域离散信号的表示方法和典型信号、时域离散线性时不变系统的时域分析方法，最后介绍模拟信号数字处理方法。第1章时域离散信号和时域离散系统 1.2时域离

5、散信号时域离散信号 实际中遇到的信号一般是模拟信号，对它进行等间隔采样便可以得到时域离散信号。假设模拟信号为xa (t)，以采样间隔T对它进行等间隔采样，得到：（1.2.1） nnTxtxnxnTt- )()()(aa第1章时域离散信号和时域离散系统 这里, x(n)称为时域离散信号，式中的n取整数，将 代入上式, 得到： 显然, x(n)是一个有序的数字，因此时域离散信号也可以称为序列。注意这里n取整数，非整数时无定义。时域离散信号有三种表示方法：),2(),(),0(),(,)(aaaaTxTxxTxnx, 3 , 2 , 1 , 0 ,n第1章时域离散信号和时域离散系统 1） 用集合符号

6、表示序列数的集合用集合符号表示。时域离散信号是一个有序的数的集合，可表示成集合: x(n)= xn, n=, 2, 1, 0, 1, 2, 例如，一个有限长序列可表示为x(n)=1, 2, 3, 4, 3, 2, 1; n=0, 1, 2, 3, 4, 5, 6也可简单地表示为x(n)1, 2, 3, 4, 3, 2, 1集合中有下划线的元素表示n=0时刻的采样值。 第1章时域离散信号和时域离散系统 2） 用公式表示序列例如：x(n)=a|n|0a1, n3） 用图形表示序列例如, 时域离散信号x(n)=sin(n/5)，n=5, 4, , 0, , 4, 5, 图1.2.1就是它的图形表示。

7、 这是一种很直观的表示方法。为了醒目，常常在每一条竖线的顶端加一个小黑点。 第1章时域离散信号和时域离散系统 图1.2.1x(n)=sin(n/5)的波形图第1章时域离散信号和时域离散系统 实际中要根据具体情况灵活运用三种表示方法，对于一般序列，包括由实际信号采样得下面介绍用MATLAB语言表示序列。MATLAB用两个参数向量x和n表示有限长序列x(n)，x是x(n)的样值向量，n是位置向量（相当于图形表示方法中的横坐标n），n与x长度相等，向量n的第m个元素n(m)表示样值x(m)的位置。位置向量n一般都是单位增向量，产生语句为：n=ns:nf; 其中ns表示序列x(n)的起始点，nf表示序

8、列x(n)的终止点。这样将有限长序列x(n)记为x(n)；n=ns:nf。 第1章时域离散信号和时域离散系统 例如， x(n)0.0000 ，0.5878 ，0.9511，0.9511，0.5878，0.0000，0.5878， 0.9511，0.9511，0.5878，0.0000，相应的 n=5, 4, 3, , 5，所以序列x(n)的MATLAB表示如下:n=5:5; x=0.0000，0.5878，0.9511，0.9511，0.5878， 0.0000，0.5878，0.9511，0.9511，0.5878，0.0000第1章时域离散信号和时域离散系统 这里x(n)的11个样值是正弦

9、序列的采样值，即x(n)=sin(n/5)n=5, 4, , 0, , 4, 5所以，也可以用计算的方法产生序列向量：n=5:5; x= sin(pi*n/5); 第1章时域离散信号和时域离散系统 这样用MATLAB计算产生x(n)并绘图的程序如下：%fig121.m：sin(pi*n/5)信号产生及图1.2.1绘图程序n=5:5; %位置向量n从5到5x=sin(pi*n/5); %计算序列向量x(n)的11个样值subplot(3, 2, 1); stem(n, x, .); line(5, 6, 0, 0)axis(5, 6, 1.2, 1.2); xlabel(n); ylabel(x

10、(n)运行程序输出波形如图1.2.1所示。 第1章时域离散信号和时域离散系统 1.2.1常用的典型序列常用的典型序列1 单位采样序列单位采样序列(n) （1.2.2） 单位采样序列也称为单位脉冲序列，特点是仅在n=0时取值为1，其它均为零。它类似于模拟信号和系统中的单位冲激函数(t)，但不同的是(t)在t=0时，取值无穷大，t0 时取值为零，对时间t的积分为1。单位采样序列和单位冲激信号如图1.2.2所示。 0 00 1)(nnn第1章时域离散信号和时域离散系统 图1.2.2单位采样序列和单位冲激信号第1章时域离散信号和时域离散系统 2 单位阶跃序列单位阶跃序列u(n) (1.2.3) 单位阶

11、跃序列如图1.2.3所示。它类似于模拟信号中的单位阶跃函数u(t)。(n)与u(n)之间的关系如下列式所示： （1.2.4） （1.2.5） 0 00 1)(nnnu) 1()()(nunun0)()(kknnu第1章时域离散信号和时域离散系统 图1.2.3单位阶跃序列第1章时域离散信号和时域离散系统 令nk=m，代入式(1.2.5)得（1.2.6） nmmnu)()(第1章时域离散信号和时域离散系统 3 矩形序列矩形序列RN(n) （1.2.7） 式中，N称为矩形序列的长度。当N=4时，R4(n)的波形如图1.2.4所示。矩形序列可用单位阶跃序列表示，如下式： （1.2.8） nNnnRN其

12、它 010 1)()()()(NnununRN第1章时域离散信号和时域离散系统 图1.2.4矩形序列第1章时域离散信号和时域离散系统 4 实指数序列实指数序列x(n)=anu(n)a为实数如果|a|1，则称为发散序列。其波形如图1.2.5所示。 第1章时域离散信号和时域离散系统 图1.2.5实指数序列 第1章时域离散信号和时域离散系统 5 正弦序列正弦序列式中, 称为正弦序列的数字域频率(也称数字频率)，单位是弧度(rad)，它表示序列变化的速率，或者说表示相邻两个序列值之间变化的弧度数。如果正弦序列是由模拟信号xa(t)采样得到的，那么 )sin()(nnx)sin()sin(| )()()

13、sin()(aannTtxnxttxnTt第1章时域离散信号和时域离散系统 （1.2.9） 因此得到数字频率与模拟角频率之间的关系为 (1.2.9)式具有普遍意义，它表示凡是由模拟信号采样得到的序列，模拟角频率与序列的数字域频率成线性关系。由于采样频率Fs与采样周期T互为倒数，因而有T第1章时域离散信号和时域离散系统 上式表示数字域频率是模拟角频率对采样频率的归一化频率。本书中用表示数字域频率，和f表示模拟角频率和模拟频率。 （1.2.10） sF第1章时域离散信号和时域离散系统 6 复指数序列复指数序列复指数序列用下式表示:式中, 0为数字域频率。设=0，用极坐标和实部虚部表示如下式：nnx

14、)j(0e)(nnx0je)()sin(j)cos()(00nnnx第1章时域离散信号和时域离散系统 由于n取整数，下面等式成立：上面公式中M取整数，所以对数字域频率而言，正弦序列和复指数序列都是以2为周期的周期信号。在以后的研究中，在频率域只分析研究一个周期就够了。nnM00j)2( jee)cos()2cos(00nnM)sin()2sin(00nnM第1章时域离散信号和时域离散系统 7 周期序列周期序列如果对所有n存在一个最小的正整数N，使下面等式成立：（1.2.11） 则称序列x(n)为周期性序列，周期为N。 nNnxnx),()(第1章时域离散信号和时域离散系统 例如：式中数字频率是

15、/4，n取整数，可以写成下式：因此, 是周期为8的周期序列，波形如图1.2.6所示。下面讨论一般正弦序列的周期性。)4sin()(nnx)8(4sin)(nnx)4sin()(nnx第1章时域离散信号和时域离散系统 图1.2.6正弦序列 第1章时域离散信号和时域离散系统 设 那么如果则要求N=(2/0)k。式中，k与N均取整数，且k的取值要保证N是最小的正整数，满足这些条件，正弦序列才是以N为周期的周期序列。)sin()(0nAnx)sin()(sin()(000NnANnANnx)()(nxNnx第1章时域离散信号和时域离散系统 具体正弦序列有以下三种情况：（1） 当2/0为整数时，k=1，

16、正弦序列是以2/0为周期的周期序列。例如, ,， 该正弦序列周期为16。 )8sin(n801620第1章时域离散信号和时域离散系统 （2） 2/0不是整数，是一个有理数时，设2/0=P/Q，式中P、Q是互为素数的整数，取k=Q，那么N=P，则该正弦序列是以P为周期的周期序列。例如, sin(4n/5), 2/0=5/2, k=2, 该正弦序列是以5为周期的周期序列。（3） 2/0是无理数，任何整数k都不能使N为正整数，因此，此时的正弦序列不是周期序列。例如，0=1/4, sin(0n)即不是周期序列。对于复数指数序列的周期性也有和上面同样的分析结果。 ne0j第1章时域离散信号和时域离散系统

17、 以上介绍了几种常用的典型序列，对于任意序列，可以用单位采样序列的移位加权和表示，即 （1.2.12） 这种任意序列的表示方法，在信号分析中是一个很有用的公式。例如, x(n)的波形如图1.2.7所示，可以用（1.2.12）式表示成： mmnmxnx)()()()6()5(2)4()2(5 . 1) 1()(2) 1(5 . 0)2(2)(nnnnnnnnnx第1章时域离散信号和时域离散系统 图1.2.7用单位采样序列移位加权和表示序列第1章时域离散信号和时域离散系统 1.2.2序列的运算序列的运算序列的简单运算有加法、乘法、移位、翻转及尺度变换。1 加法和乘法加法和乘法序列之间的加法和乘法，

18、是指它的同序号的序列值逐项对应相加和相乘，如图1.2.8所示。 第1章时域离散信号和时域离散系统 图1.2.8序列的加法和乘法第1章时域离散信号和时域离散系统 2 移位、翻转及尺度变换移位、翻转及尺度变换序列x(n)如图1.2.9(a)所示，其移位序列x(nn0)(当n0=2时)如图1.2.9(b)所示。当n00时, 称为x(n)的延时序列；当n00 , 序列右移；n0，序列左移。如n=1，得到h(1-m)，如图1.3.2（d）所示。接着将h(m)和h(nm)相乘后，再相加, 得到y(n)的一个值。对所有的n重复这种计算, 最后得到卷积结果，如图1.3.2(f)所示, y(n)表达式为y(n)

19、=1, 2, 3, 4, 3, 2, 1 其实这种图解法可以用列表法代替，上面的图解过程如表1.3.1所示。 第1章时域离散信号和时域离散系统 图1.3.2例1.3.4线性卷积第1章时域离散信号和时域离散系统 表表1.3.1图解法（列表法）图解法（列表法） 第1章时域离散信号和时域离散系统 2) 解析法如果已知两个卷积信号的解析表达式，则可以直接按照卷积式进行计算，下面举例说明。 【例例1.3.5】设x(n)=an(n)，h(n)=R4(n)，求y(n)=x(n)*h(n)。解解 mmnmnuamRnxnhny)()()()()(4第1章时域离散信号和时域离散系统 要计算上式，关键是根据求和号

20、内的两个信号乘积的非零值区间确定求和的上、下限。因为nm时，u(n-m)才能取非零值； 0m3时，R4(m)取非零值, 所以，求和区间中m要同时满足下面两式： mn0m3这样求和限与n有关系，必须将n进行分段然后计算。 第1章时域离散信号和时域离散系统 n0时，y(n)=00n3时，乘积的非零值范围为0mn，因此 n4时，乘积的非零区间为0m3，因此 11011)(aaaanynnnmmn143011)(aaaanynmmn第1章时域离散信号和时域离散系统 写成统一表达式为naaanaaannynnn411301100)(1411第1章时域离散信号和时域离散系统 3） 用MATLAB计算两个有

21、限长序列的卷积MATLAB 信号处理工具箱提供了conv 函数，该函数用于计算两个有限长序列的卷积（或计算两个多项式相乘）。C=conv(A, B)计算两个有限长序列向量A和B的卷积。如果向量A和B的长度分别为N和M，则卷积结果向量C的长度为NM1。如果向量A和B为两个多项式的系数，则C就是这两个多项式乘积的系数。应当注意，conv函数默认A和B表示的两个序列都是从0开始，所以不需要位置向量。第1章时域离散信号和时域离散系统 当然, 默认卷积结果序列C也是从0开始，即卷积结果也不提供特殊的位置信息。例1.3.4中的两个序列满足上述条件，直接调用conv函数求解例1.3.4的卷积计算程序ep13

22、4.m如下：%ep134.m：例1.3.4的卷积计算程序xn=1 1 1 1 ; hn=1 1 1 1; yn=conv(xn, hn); 运行结果：yn=1, 2, 3, 4, 3, 2, 1 第1章时域离散信号和时域离散系统 显然，当两个序列不是从0开始时，必须对conv函数稍加扩展。设两个位置向量已知的序列：x(n)；nx=nxs:nxf，h(n)； nh=nhs:nhf，要求计算卷积：y(n)=h(n)*x(n)以及y(n)的位置向量ny。下面编写计算这种卷积的通用卷积函数convu。第1章时域离散信号和时域离散系统 根据卷积原理知道， y(n)的起始点和终止点分别为：nys=nhs+

23、nxs，nyf=nhf+nxf。调用conv函数写出通用卷积函数convu如下：function y, ny=convu(h, nh, x, nx) %convu 通用卷积函数，y为卷积结果序列向量，%ny是y的位置向量, h和x是有限长序列，%nh和nx分别是h和x的位置向量nys=nh(1)+nx(1); nyf=nh(end)+nx(end); %end表示最后一个元素的下标y=conv(h, x); ny=nys:nyf; 第1章时域离散信号和时域离散系统 如果h(n)=x(n)=R5(N+2)，则调用convu函数计算y(n)=h(n)*x(n)的程序如下： h=ones(1, 5)

24、; nh=2:2; x=h; nx=nh; y, ny=convu(h, nh, x, nx)运行结果：y = 1 2 3 4 5 4 3 2 1ny=4 3 2 1 0 1 2 3 4 第1章时域离散信号和时域离散系统 线性卷积服从交换律、结合律和分配律。它们分别用公式表示如下：（1.3.9） （1.3.8） （1.3.10） )()()()(nxnhnhnx)()()()()()(2121nhnhnxnhnhnx)()()()()()()(2121nhnxnhnxnhnhnx第1章时域离散信号和时域离散系统 以上三个性质请读者自己证明。（1.3.8）式表示卷积服从交换律。（1.3.9）和（

25、1.3.10）式分别表示卷积的结合律和分配律。设h1(n)和h2(n)分别是两个系统的单位脉冲响应，x(n)表示输入序列。按照（1.3.9）式的右端，信号通过h1(n)系统后再通过h2(n)系统，等效于按照（1.3.9）式左端，信号通过一个系统，该系统的单位脉冲响应为h1(n)*h2(n)，如图1.3.3(a)、(b)所示。 第1章时域离散信号和时域离散系统 该式还表明两系统级联，其等效系统的单位脉冲响应等于两系统分别的单位脉冲响应的卷积。按照（1.3.10）式，信号同时通过两个系统后相加，等效于信号通过一个系统，该系统的单位脉冲响应等于两个系统分别的单位脉冲响应之和，如图1.3.3(c)、(

26、d)所示。换句话说，系统并联的等效系统的单位脉冲响应等于两个系统分别的单位脉冲响应之和。第1章时域离散信号和时域离散系统 图1.3.3卷积的结合律和分配律第1章时域离散信号和时域离散系统 需要再次说明的是，关于系统级联、并联的等效系统的单位脉冲响应与原来两系统分别的单位脉冲响应的关系，是基于线性卷积的性质，而线性卷积是基于线性时不变系统满足线性叠加原理。因此, 对于非线性或者非时不变系统，这些结论是不成立的。 第1章时域离散信号和时域离散系统 再考察（1.3.11）式，它也是一个线性卷积式，它表示序列x(n)与单位脉冲序列的线性卷积等于序列本身x(n)，（1.3.11） mnnxmnmxnx)

27、()()()()(第1章时域离散信号和时域离散系统 如果序列与一个移位的单位脉冲序列(nn0)进行线性卷积，就相当于将序列本身移位n0(n0是整常数)，如下式表示：上式中求和项只有当m=nn0时才有非零值，因此得到：（1.3.12） mmnnmxnnnxny)()()()()(00)()()(00nnxnnnx第1章时域离散信号和时域离散系统 【例例1.3.6】在图1.3.4中，h1(n)系统与h2(n)系统级联，设求系统的输出y(n)。 1|),()()4()()()()(21anuanhnnnhnunxn第1章时域离散信号和时域离散系统 图1.3.4例1.3.5框图第1章时域离散信号和时域

28、离散系统 解解先求第一级的输出m(n)，再求y(n)。 )3()2() 1()()3()2() 1()()()()()()()()()4()()4()()()()4()()()()()(3214241nuanuanuanuannnnnuanuanRnhnmnynRnununnunnunnnunhnxnmnnnnnn30)(iininua第1章时域离散信号和时域离散系统 1.3.4系统的因果性和稳定性系统的因果性和稳定性如果系统n时刻的输出只取决于n时刻以及n时刻以前的输入序列，而和n时刻以后的输入序列无关，则称该系统具有因果性质，或称该系统为因果系统。如果n时刻的输出还取决于n时刻以后的输入序

29、列，在时间上违背了因果性，系统无法实现，则系统被称为非因果系统。因此系统的因果性是指系统的可实现性。 第1章时域离散信号和时域离散系统 线性时不变系统具有因果性的充分必要条件是系统的单位脉冲响应满足下式： （1.3.13） 满足（1.3.13）式的序列称为因果序列，因此因果系统的单位脉冲响应必然是因果序列。因果系统条件（1.3.13）式从概念上也容易理解，因为单位脉冲响应是输入为(n)的零状态响应，在n=0时刻以前即n0时，没有加入信号，输出只能等于零，因此得到因果性条件（1.3.13）式。0, 0)(nnh第1章时域离散信号和时域离散系统 所谓稳定系统，是指对有界输入，系统输出也是有界的。系

30、统稳定的充分必要条件是系统的单位脉冲响应绝对可和，用公式表示为（1.3.14） nnh| )(|第1章时域离散信号和时域离散系统 证明证明先证明充分性。 因为输入序列x(n)有界，即kkknxkhnyknxkhny| )(|)(| )(|)()()(8为任意常数BnBnx , | )(|第1章时域离散信号和时域离散系统 因此如果系统的单位脉冲响应满足（1.3.14）式，那么输出y(n)一定也是有界的，即kkhBny| )(| )(| | )(|ny第1章时域离散信号和时域离散系统 下面用反证法证明其必要性。如果h(n)不满足(1.3.14)式，即, 那么总可以找到一个或若干个有界的输入来引起无

31、界的输出，例如： nnh| )(|0)(00)()()()(nhnhnhnhnxkknxkhny)()()(第1章时域离散信号和时域离散系统 令n=0, 有上式说明n=0时刻的输出为无界，系统不稳定，证明了（1.3.14）式条件的必要性。 kkkkhkhkhkhkxkhy| )(| )(|)()()()()0(第1章时域离散信号和时域离散系统 【例例1.3.7】设线性时不变系统的单位系统脉冲响应h(n)=anu(n)，式中a是实常数，试分析该系统的因果稳定性。解解由于n0时，h(n)=0，因此系统是因果系统。|1|1lim|lim| )(|100aaaanhNNnNnnnNn第1章时域离散信号

32、和时域离散系统 只有当|a|1时, 才有因此系统稳定的条件是|a|1；否则，|a|1时，系统不稳定。系统稳定时，h(n)的模值随n加大而减小，此时序列h(n)称为收敛序列。如果系统不稳定，h(n)的模值随n加大而增大，则称为发散序列。 nanh|11| )(|第1章时域离散信号和时域离散系统 【例例 1.3.8】设系统的单位脉冲响应h(n)=u(n)，求对于任意输入序列x(n)的输出y(n)， 并检验系统的因果性和稳定性。解解kknukxnhnxnynunh)()()()()()()(因为当nk0的方向递推，是一个因果解。但对于差分方程，其本身也可以向n0的方向递推，得到的是非因果解。因此差分

33、方程本身不能确定该系统是因果系统还是非因果系统，还需要用初始条件进行限制。下面就是向方向n0递推的例题。 第1章时域离散信号和时域离散系统 【例例1.4.2】设差分方程为求输出序列y(n)。 , 0, 0)(),()( )() 1()(nnynnxnxnayny)()() 1(1nnyany121111) 1(|)1() 1()2(,1)0()0() 1(,00)1 () 1 ()0(,1nany，nnayaynayaynyayn时时时时第1章时域离散信号和时域离散系统 将n1用n代替，得到：这确实是一个非因果的输出信号。用差分方程求系统的单位脉冲响应，由于单位脉冲响应是当系统输入(n)时的零

34、状态响应，因此只要令差分方程中的输入序列为(n)，N个初始条件都为零，其解就是系统的单位脉冲响应。实际上例题1.4.1（1）中求出的y(n)就是该系统的单位脉冲响应，例题1.4.2求出的y(n)则是一个非因果系统的单位脉冲响应。) 1()(nuanyn第1章时域离散信号和时域离散系统 最后要说明的是，一个线性常系数差分方程描述的系统不一定是线性非时变系统，这和系统的初始状态有关。如果系统是因果的，一般在输入x(n)=0(nn0)时，则输出y(n)=0(nn0)，系统是线性非时变系统。下面介绍用MATLAB求解差分方程。 第1章时域离散信号和时域离散系统 MATLAB 信号处理工具箱提供的fil

35、ter函数实现线性常系数差分方程的递推求解，调用格式如下：yn=filter(B, A.xn)计算系统对输入信号向量xn的零状态响应输出信号向量yn，yn与xn长度相等，其中，B和A是（1.4.2）式所给差分方程的系数向量，即 B=b0, b1, , bM, A=a0, a1, ,aN其中a0=1，如果a01，则filter用a0对系数向量B和A归一化。第1章时域离散信号和时域离散系统 yn=filter(B, A.xn，xi)计算系统对输入信号向量xn的全响应输出信号yn。所谓全响应，就是由初始状态引起的零输入响应和由输入信号xn引起的零状态响应之和(在 2.4.3 节介绍)。其中, xi是

36、等效初始条件的输入序列，所以xi是由初始条件确定的。MATLAB信号处理工具箱提供的filtic就是由初始条件计算xi的函数, 其调用格式如下：xi=filtic(B, A, ys, xs)其中，ys和xs是初始条件向量：ys= y(1)，y(2)，y(3)，y(N)，xs= x(1)， x（2），x(3)，x(M) 。如果xn是因果序列，则xs=0，调用时可缺省xs。第1章时域离散信号和时域离散系统 例1.4.1的MATLAB求解程序ep141.m如下：%ep141.m：调用filter解差分方程y(n)ay(n1)=x(n)a=0.8; ys=1; %设差分方程系数a=0.8, %初始状态

37、: y(1)=1xn=1, zeros(1, 30); %x(n)=单位脉冲序列, 长度N=31B=1; A=1, -a; %差分方程系数xi=filtic(B, A, ys); %由初始条件计算等效初始条件的输入序列xi第1章时域离散信号和时域离散系统 yn=filter(B, A, xn, xi); %调用filter解差分方程, 求系统输出信号y(n)n=0:length(yn)-1; subplot(3, 2, 1); stem(n, yn, .)title(a); xlabel(n); ylabel(y(n) 第1章时域离散信号和时域离散系统 程序中取差分方程系数a=0.8时，得到系

38、统输出y(n)如图1.4.1（a）所示，与例1.4.1的解析递推结果完全相同。如果令初始条件y(1)=0 (仅修改程序中ys=0)，则得到系统输出y(n)=h(n)，如图1.4.1（b）所示。第1章时域离散信号和时域离散系统 图1.4.1例1.4.1求解程序输出波形第1章时域离散信号和时域离散系统 1.5模拟信号数字处理方法模拟信号数字处理方法在绪论中已介绍了数字信号处理技术相对于模拟信号处理技术的许多优点，因此人们往往希望将模拟信号经过采样和量化编码形成数字信号，再采用数字信号处理技术进行处理；处理完毕，如果需要，再转换成模拟信号。这种处理方法称为模拟信号数字处理方法。其原理框图如图1.5.

39、1所示。图中的预滤与平滑所起的作用在后面介绍。本节主要介绍采样定理和采样恢复。第1章时域离散信号和时域离散系统 图1.5.1模拟信号数字处理框图第1章时域离散信号和时域离散系统 1.5.1采样定理及采样定理及A/D变换器变换器对模拟信号进行采样可以看做一个模拟信号通过一个电子开关S。设电子开关每隔周期T合上一次，每次合上的时间为T，在电子开关输出端得到其采样信号 。该电子开关的作用等效成一宽度为，周期为T的矩形脉冲串pT(t)，采样信号 就是xa(t)与pT(t)相乘的结果。采样过程如图1.5.2(a)所示。如果让电子开关合上时间0，则形成理想采样，此时上面的脉冲串变成单位冲激串，用p(t)表

40、示。p+(t)中每个单位冲激处在采样点上，强度为1，理想采样则是xa(t)与p(t)相乘的结果，采样过程如图1.5.2(b)所示。用公式表示为 )(atx)(atx第1章时域离散信号和时域离散系统 （1.5.1） 上式中(t)是单位冲激信号，在上式中只有当t=nT时，才可能有非零值，因此写成下式：（1.5.2） nnTttP)()()()()()()(aaanTttxtPtxtxn)()()(aanTtnTxtxn第1章时域离散信号和时域离散系统 图1.5.2对模拟信号进行采样 第1章时域离散信号和时域离散系统 下面研究理想采样前后信号频谱的变化，从而找出为了使采样信号能不失真地恢复原模拟信号

41、，采样速率Fs(Fs=T1)与模拟信号最高频率fc之间的关系。我们知道在傅里叶变换中，两信号在时域相乘的傅里叶变换等于两个信号分别的傅里叶变换的卷积，按照（1.5.2）式，推导如下： 第1章时域离散信号和时域离散系统 设 对（1.5.1）式进行傅里叶变换，得到 （1.5.3） )(FT)j ()(FT)j ()(FT)j (aaaatpPtxXtxXkkkaP)(2)j (s第1章时域离散信号和时域离散系统 式中，s=2/T，称为采样角频率，单位是rad/s。因此2/2/j1de )(1sTttkkTttTakkTP)(2)j (s（1.5.4） 第1章时域离散信号和时域离散系统 （1.5.5

42、） kkkkXTkXTkXTPXX)jj (1d)()j (1d)()j (221)j ()j (21)j (sasasaaa第1章时域离散信号和时域离散系统 上式表明理想采样信号的频谱是原模拟信号的频谱沿频率轴，每间隔采样角频率s重复出现一次，或者说理想采样信号的频谱是原模拟信号的频谱以s为周期，进行周期性延拓而成的。第1章时域离散信号和时域离散系统 在图1.5.3中，设xa(t)是带限信号，最高频率为c，其频谱Xa(j)如图1.5.3(a)所示。p(t)的频谱P(j)如图1.5.3(b)所示，那么按照（1.5.5）式，的频谱如图 1.5.3(c)所示，图中原模拟信号的频谱称为基带频谱。如果

43、满足s2c，或者用频率表示该式，即满足Fs2fc，基带谱与其它周期延拓形成的谱不重叠，如图1.5.3(c)所示情况，可以用理想低通滤波器G(j)从采样信号中不失真地提取原模拟信号，如图1.5.4所示。)(atx)j (aX第1章时域离散信号和时域离散系统 但如果选择采样频率太低，或者说信号最高截止频率过高，使Fs2fc, Xa(j)按照采样频率Fs周期延拓时，形成频谱混叠现象，用图1.5.3(d)表示。这种情况下，再用图 1.5.4 所示的理想低通滤波器对Xa(t)进行滤波，得到的是失真了的模拟信号。下面用公式表示： 第1章时域离散信号和时域离散系统 （1.5.6） ss21| 021| )j

44、 (TGscaascaaa1aaaa21 )()(21 )()()j (FT)()j ()j ()(FT)j (txtytxtyYtyGXtyY第1章时域离散信号和时域离散系统 这里需要说明的是，一般频谱函数是复函数，相加应是复数相加，图1.5.3和图1.5.4仅是示意图。一般称Fs/2为折叠频率，只有当信号最高频率不超过Fs/2时，才不会产生频率混叠现象，否则超过Fs/2的频谱会折叠回来而形成混叠现象，因此频率混叠在Fs/2附近最严重。 第1章时域离散信号和时域离散系统 图1.5.3采样信号的频谱 第1章时域离散信号和时域离散系统 图1.5.4采样恢复 第1章时域离散信号和时域离散系统 总结

45、上述内容，采样定理叙述如下：（1） 对连续信号进行等间隔采样形成采样信号，采样信号的频谱是原连续信号的频谱以采样频率s为周期进行周期性的延拓形成的，用公式（1.5.5）表示。（2） 设连续信号xa(t)属带限信号，最高截止频率为c，如果采样角频率s2c，那么让采样信号通过一个增益为T、 截止频率为s/2的理想低通滤波器，可以唯一地恢复出原连续信号xa(t)。否则, s/T区域有较多的高频分量，表现在时域上，就是恢复出的模拟信号是台阶形的。因此需要在D/AC之后加平滑低通滤波器，滤除多余的高频分量，对时间波形起平滑作用，这也就是在图1.5.1模拟信号数字处理框中，最后加平滑滤波器的原因。虽然这种

46、零阶保持器恢复的模拟信号有些失真，但简单、易实现，是经常使用的方法。实际中，将解码器与零阶保持器集成在一起，就是工程上的D/AC器件。 第1章时域离散信号和时域离散系统 图 1.5.10零阶保持器的频率特性 第1章时域离散信号和时域离散系统 习题与上机题习题与上机题1. 用单位脉冲序列(t)及其加权和表示题1图所示的序列。 第1章时域离散信号和时域离散系统 题1图第1章时域离散信号和时域离散系统 2给定信号：（1） 画出x(n)序列的波形，标上各序列值；（2） 试用延迟的单位脉冲序列及其加权和表示x(n)序列；（3） 令x1(n)=2x(n2)，试画出x1(n)波形；（4） 令x2(n)=2x

47、(n+2)，试画出x2(n)波形；（5） 令x3=x(2n)，试画出x3(n)波形。 其它04061452)(nnnnx第1章时域离散信号和时域离散系统 3 判断下面的序列是否是周期的; 若是周期的，确定其周期。（1） A是常数（2） 4 对题1图给出的x(n)要求：（1） 画出x(n)的波形；（2） 计算，并画出xe(n)波形；,873cos)(nAnx)81( j)(nenx)()(21)(enxnxnx第1章时域离散信号和时域离散系统 （3） 计算，并画出x0(n)波形；（4） 令x1(n)=xe(n)+x0(n), 将x1(n)与x(n)进行比较，你能得到什么结论？ 5 设系统分别用下

48、面的差分方程描述，x(n)与y(n)分别表示系统输入和输出，判断系统是否是线性非时变的。 （1）y(n)=x(n)+2x(n1)+3x(n2) （2）y(n)=2x(n)+3 )()(21)(0nxnxnx第1章时域离散信号和时域离散系统 （3）y(n)=x(nn0)n0为整常数（4）y(n)=(n)（5）y(n)=x2(n)（6）y(n)=x(n2)（7）（8）y(n)=x(n)sin(n)nmmxny0)()(第1章时域离散信号和时域离散系统 6 给定下述系统的差分方程，试判定系统是否是因果稳定系统，并说明理由。（1） （2） y(n)=x(n)+x(n+1)（3） （4） y(n)=x(nn0)（5） y(n)=ex(n)10)(1)(Nkkn