数字信号处理6

1、第第4章章 离散时间系统离散时间系统l离散时间系统举例l离散时间系统的分类*l冲激和阶跃响应lLTI离散时间系统的时域特性*l简单互联方案l有限维LTI离散时间系统lLTI离散时间系统的分类*lLTI 离散时间系统的频域表示*l相位延迟和群延迟*离散时间系统离散时间系统l离散时间系统的功能是对给定的输入序列 进行处理得到输出序列 。输出序列应含有更多所需的特性，或可以从中提取关于输入序列的某种信息。l在大多数应用中，用到的离散时间系统是单输入单输出系统。 x n y n离散时间系统离散时间系统l数学上，离散时间系统用运算子 表述，该运算将输入序列变换成系统输出的另外一个序列。l离散时间系统也可

2、有会有多个输入或者多个输出。() 4.1 离散时间系统举例离散时间系统举例双输入、单输出离散时间系统：双输入、单输出离散时间系统： 单输入、单输出离散时间系统：单输入、单输出离散时间系统： 调制器 加法器乘法器 单位延迟 单位超前4.1 离散时间系统举例离散时间系统举例一个更复杂的单输入、单输出离散时间系统：一个更复杂的单输入、单输出离散时间系统： 4.1 离散时间系统举例离散时间系统举例 累加器：累加器： 在时刻在时刻 的输出的输出 是在时刻是在时刻 输入序列输入序列 以及在时刻以及在时刻 的前一个的前一个 输出样本之和输出样本之和，后者是从，后者是从 到时刻到时刻 所有输入样本值的所有输入

3、样本值的和。和。 此系统累加了所有的输入数据。此系统累加了所有的输入数据。 n y nn x n1ny1n1n 1yy1nnnxxx nnx n4.1 离散时间系统举例离散时间系统举例 累加器的输入累加器的输入输出关系的另一种变形是输出关系的另一种变形是 上面形式的累加器用于因果输入序列，此时被上面形式的累加器用于因果输入序列，此时被称为初始条件。称为初始条件。 4.1 离散时间系统举例离散时间系统举例lM点滑动平均滤波器系统点滑动平均滤波器系统用于平滑随机变化的数据。在绝大多数应用中，数据 有界序列。因此M点均值 也是有界序列。 x n y n4.1 离散时间系统举例离散时间系统举例若观测过

4、程没有偏差，则可以简单地通过增加M来提高对噪声数据估计的准确性。M点滑动平均滤波器可以直接实现包括M-1次相加、以值1/M为因子的一次相乘和能存储M-1个过去输入数据样本的存储器。下面将推导滑动平均滤波器的一种更为有效的实现。 11011111y1MMMnx nx nx nMMx nx nx nMx nMM4.1 离散时间系统举例离散时间系统举例l使用上述递归方程计算序列第n时刻的M点滑动均值 ，现在只需要2次相加和1次与1/M相乘。 y n4.1 离散时间系统举例离散时间系统举例 其中其中 带有噪声带有噪声 的信号。的信号。 x ns nd n s n d n4.1 离散时间系统举例离散时间

5、系统举例l指数加权移动平均滤波器 用上式计算的移动平均需要包括1次相加和1次相乘。而且，只需存储前一次移动平均的结果。 不需要存储以前的输入数据样本。不需要存储以前的输入数据样本。 4.1 离散时间系统举例离散时间系统举例 当当 时，指数加权平均滤波器通过对数据时，指数加权平均滤波器通过对数据样本以指数形式加权，对当前数据样本的权重样本以指数形式加权，对当前数据样本的权重较大，而以前数据样本的权重较小。如下所示较大，而以前数据样本的权重较小。如下所示： 014.1 离散时间系统举例离散时间系统举例l线性内插器通常用于估计离散时间序列中相邻的一对样本值之间的样本值的大小。 因子为4的内插4.1

6、离散时间系统举例离散时间系统举例l因子为因子为2的内插器的内插器 (a)大小为512512的原始灰度图像；(b)大小为256256的下抽样后的图像；(c)用双线性内插法得到的放大图像4.1 离散时间系统举例离散时间系统举例l中值滤波器 一组（2K+1）个数的中值定义为这样的一个数：若该组中的K个数值大于该数值，而剩下的K的数小于改数值，则该数据即为中值。可以根据数组中数值的大小对数排序，然后选取位于中间的那个数。 例子：考虑一组数字2，3，10，5，1 排序后为 3，1，2，5，10 因此，中值2，3，10，5，1=2 4.1 离散时间系统举例离散时间系统举例中值滤波器是通过在输入序列 上滑动

7、一个长度为奇数的窗口来实现的，每次滑动一个样本值。在任意时刻n，输出 都为当前窗口中所有输入样本值的中值。中值滤波器常用于去除加性随机冲激噪声，这类噪声将会导致受干扰信号中存在大量突发错误。中值滤波器也通常被用来平滑由冲激噪声干扰的信号。 x n y n4.2 离散时间系统的分类离散时间系统的分类 线性系统移不变系统因果系统稳定系统移不变系统无源无损系统4.2 离散时间系统的分类离散时间系统的分类l线性离散时间系统 定义：若 和 分别是输入序列为 和 的响应，则当输入为 时，其响应为 叠加定理必须对任意常数叠加定理必须对任意常数 和和 以及所有可能以及所有可能输入输入 和和 成立。成立。 1y

8、n 2yn 1xn 2xn 12x nxnxn 12yyynnn 1xn 2xn4.2 离散时间系统的分类离散时间系统的分类 累加器 如果输入为 则输出为 因此，此系统是线性的。 11-ynnx 22-ynnx 12x nxnxn 12-1212-yyynnnnxxxxnn4.2 离散时间系统的分类离散时间系统的分类 当输入为 和 时，得到的输出 和 为 当输入为 时，输出为 1xn 2xn 1yn 2yn 1110yy1nnx 2220yy1nnx 12xnxn 120yy1nnxx4.2 离散时间系统的分类离散时间系统的分类 另一方面，另一方面， 若 则 12112200121200yy=

9、y1+y1=y1 +y1 +nnnnnnxxxx12y -1y1y1 12yyynnn4.2 离散时间系统的分类离散时间系统的分类对于有因果输入的累加器，线性条件 只有当初始条件 、 和 以及所有常量 和 都满足上面的初始条件。若给出的该系统不是处于零初始状态，则这个条件不满足。 12y -1y1y1y11y12y14.2 离散时间系统的分类离散时间系统的分类l非线性离散时间系统 中值滤波器是一个非线性离散时间系统考虑窗长为3的中值滤波器若输入为 ，则其输出是长度为3的序列 。若输入为 ，则其输出是长度为3的序列 。 ，13,4,5 02xnn ，1y3,4,4 02nn ，22,-1,-1

10、02xnn ，2y2,-1,-1 02nn4.2 离散时间系统的分类离散时间系统的分类然而，若输入为 ，则其输出为 。注意： 。因此，中值滤波器是非线性离散时间系统。非零初始状态的累加器的第二种形式是另外一个例子。 12+x nxnxn y3,4,3n 12y+y3,3,3ynnn4.2 离散时间系统的分类离散时间系统的分类l移不变系统 对于移不变离散时间系统，若 是对 的响应，则对输入 的响应可以简单地表示为 ，式中 是任意正或负整数。输入和输出间的这种关系必须对任意输入序列及其输出成立。 若序列或系统的序号n与离散时刻相关联时，上面的限制通常被称为时不变性质。 1yn 1xn 10 x n

11、xnn 10yynnn0n4.2 离散时间系统的分类离散时间系统的分类 时不变性质保证对于一个给定的输入信号，系统相应的输出与输入信号所加入的时刻无关。 举例：上抽样器 输入与输出的关系： /0, 2 ,0ux nLnLLxn，其他4.2 离散时间系统的分类离散时间系统的分类 当输入为 时，输出 为 由上抽样器的定义得 因此，上抽样器是一个时变系统。 10 xnx nn 1,uxn4.2 离散时间系统的分类离散时间系统的分类l线性时不变系统 线性时不变系统（LTI）：既满足线性又满足时不变性的系统。线性时不变系统在数学上容易分析和表示，易于设计。在过去几十年来，用这类系统开发了很多有用的算法。

12、4.2 离散时间系统的分类离散时间系统的分类在因果系统中，第 个输出样本 仅由对应于 的输入样本决定，而不依赖于对应于 的输入样本。在因果离散时间系统中，若输入为 和 ，响应分别为 和 。 则 也就是说 简单地说，对于因果系统，输出的变化并不先于输入的变化。简单地说，对于因果系统，输出的变化并不先于输入的变化。0n0y n0nn0nn 1xn 2xn 1yn 2yn 12=xnxnnN， 12y=ynnnN，4.2 离散时间系统的分类离散时间系统的分类因果系统：因果系统：非因果系统：非因果系统： 123401212y=+-1 +-2 +-3y=+-1 +-2+y-1 +y-2y=y-1 +nx

13、 nx nx nx nnb x nb x nb x nanannnx n 1y=+-1121y=+-1232+-213uuuuuuuunxnxnxnnxnxnxnxnxn4.2 离散时间系统的分类离散时间系统的分类稳定性有不同的定义方式。当且仅当系统对于有界输入产生有界输出时，我们说该系统是稳定的。若对于 的响应是序列 ，且对于所有的n值，有 则对于所有的n值有 x n y n4.2 离散时间系统的分类离散时间系统的分类举例M点滑动平均滤波器是BIBO稳定的：对于一个有界输入我们有 101y=-Mknx n kM4.2 离散时间系统的分类离散时间系统的分类l无源无损系统 无源离散时间系统是指，

14、对每个具有有限能量的输入序列 ，输出序列 的能量不超过输入的能量，即 对于每一个输入序列，上面的不等式中的等号成立，则该系统是无损系统。 x n y n4.2 离散时间系统的分类离散时间系统的分类l举例考虑由 定义的离散时间系统，其中N是正整数。其输出能量为 因此，当 时，该系统是无源系统；而当 时，则是无损系统。 y-nx n N11=4.3 冲激和阶跃响应冲激和阶跃响应单位样本响应是指输入单位样本序列 时数字滤波器的输出，或简称为冲激响应，记为 。单位阶跃响应是指输入单位阶跃序列 时离散时间系统的输出，或简称为阶跃响应，记为 。 n h n n s n4.3 冲激和阶跃响应冲激和阶跃响应一

15、个离散时间系统的输入输出关系为 若可得其冲激响应 为可知其冲激响应是长度为4的有限长序列 x nn h n4.3 冲激和阶跃响应冲激和阶跃响应 举例累加器的冲激响应 令 ，则其冲激响应为 x nn4.3 冲激和阶跃响应冲激和阶跃响应 例因子为2的线性内插器的冲激响应 令 ，则其冲激响应为 可以看出，该冲激响应是长度为3的有限长序列，也可记为 x nn 1y=+-112uuunxnxnxn4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 4.4.1 输入输出关系线性时不变性质使得LTI离散时间系统可以由其冲激响应完全描述。即若已知冲激响应，就可以得到系统对任意输入的输出。设 表示LTI

16、离散时间系统的冲激响应 输入为 h n4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 我们计算其输出： 由于该系统是时不变的，因此我们可以计算每一个输入的输出，然后再将其叠加得到 。 由于该系统是时不变的 y n22112255nh nnh nnh nnh n输入输出4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 同样的，由于该系统是线性的 因此，由于线性特性，我们得到 0.5221.511220.7555nh nnh nnh nnh n输入输出0.51.5-0.754.4 LTI离散时间系统的时域特性离散时间系统的时域特性 由上面可知，任意输入序列可以表示为形如 的延

17、迟和超前单位样本序列的线性加权和。 LTI离散时间系统对于序列 的响应为 。 因此。离散时间系统对于 的响应 为 x knk x k h nk x n y n4.4 LTI离散时间系统的时域特性离散时间系统的时域特性l卷积和 即 举例冲激响应为 的LTI离散时间系统，当输入为 时，其输出 为： 当冲激响应为 时，对于输入 ，可以确定其输出 y=kknx k h nkx nk h n y=nx nh n h n x n y n y=knx k h nkx nh nh nm x n 1yn4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 现在 因此， 性质性质 交换律 结合律 分配率 1

18、yynnm4.4 LTI离散时间系统的时域特性离散时间系统的时域特性l卷积和运算可以作如下解释：卷积和运算可以作如下解释： 先将序列 时间反转得到接着，将 平移 （ 表示向右移 个抽样周期； 表示向左移 个抽样周期）形成序列形成乘积序列把 的全部样本求和得到卷积和 的第 个样本 h khkhknn0nn0nh nk v kx k h nk v k y nn4.4 LTI离散时间系统的时域特性离散时间系统的时域特性在实际中，只要冲激响应序列或输入序列是有限长度的，卷积和就可以用来计算任何时刻的输出样本。若输入序列和冲激响应序列都是有限长度的，则输出序列也是有限长度的。若输入序列和冲激响应序列都是

19、无限长的，则不可以用卷积和的方法计算系统的输出结果。对于此类系统，我们仅分析其另一种时域描述，此时，该系统仅涉及卷积的有限和。4.4 LTI离散时间系统的时域特性离散时间系统的时域特性4.4.2 计算卷积和的列表法 两个有限长序列的卷积和计算可以通过列表法得到。 计算序列 与序列 的 卷积，得到序列 采用传统的乘法，对两个序列的样本相乘，但在列之间没有进位。 , 03g nn , 02h nn4.4 LTI离散时间系统的时域特性离散时间系统的时域特性把每个 样本那一列上方对应的三项相加，可得到由卷积和产生的序列 的样本。 y n y n4.4 LTI离散时间系统的时域特性离散时间系统的时域特性

20、 样本如下： MATLAB实现卷积 conv(a,b) y 0g 00y 1g 10g 01y 2g 20g 11g 02y 3g 30g 21g 12y 4g 31g 22y 5g 32hhhhhhhhhhhh4.4 LTI离散时间系统的时域特性离散时间系统的时域特性l4.4.3 用冲激响应表示稳定条件 界输入有界输出（BIBO）稳定条件 对于所有输入序列 ，系统的输出序列 仍保持有界，则称该离散时间系统是稳定的。当且仅当LTI数字滤波器的冲激响应序列 绝对可和时，即 该系统是（BIBO）稳定的。 x n y n h n4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 证明：设 是

21、一个实序列 由于输入序列 是有界的 因此， 意味着 ，即 是有界的。 h n x n xx nBS yyxnBB S= y n4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 为证明逆命题也成立，假设序列 是有界的，即 。 考虑输入为 其中 是符号函数，定义为 y n yy nB ,0,0sgn hnhnx nKhn sgn c 1,01,0csgn cc1K 4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 注意：因为 ，显然 是有界的。 对于该输入，当 时， 为 因此，由 推出 。 1x n x n0n y n yy nBS4.4 LTI离散时间系统的时域特性离散时间系

22、统的时域特性 举例假定一个因果LTI离散时间系统的冲激响应为 对该系统 因此，当 时，有 ，此时上面的系统是BIBO稳定的。当 时，系统不是BIBO稳定的。 1S1=4.4 LTI离散时间系统的时域特性离散时间系统的时域特性l4.4.4用冲激响应表示因果性条件 设两个输入序列 和 在 时刻，该LTI离散时间系统的相应输出样本为 1xn 2xn0nn4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 若该LTI离散时间系统也是因果的，则 由于 则 所以这两个和式相等的唯一条件就是 和 均等于零 1020yynn 120,x nxnnn 1xn 2xn4.4 LTI离散时间系统的时域特性离

23、散时间系统的时域特性 如果 也满足上面的条件。 于是，当且仅当LTI离散时间系统的冲激响应序列 是因果序列时，该系统才是因果的。 举例离散时间系统为因果系统 则其冲激响应为 h n 1234y=123nx nx nx nx n 1234h n 4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 举例离散时间累加器为因果系统 则其冲激响应为 举例因子为2的内插器不是一个因果系统 y=nlnln =nlh nln 1y112uuun = xnxnxn4.4 LTI离散时间系统的时域特性离散时间系统的时域特性 则其非因果冲激响应为 注意：可以对具有有限长冲激响应的非因果离散时间系统中插入适当

24、数量的延迟来实现因果系统。例如，对离散时间因为为2的线性内插器的输出延迟一个样本周期，可以得到对应的因果形式，新的输入输出关系为 = 0.5 1 0.5h n4.5 简单互联方案简单互联方案 l4.5.1级联 若两个系统的冲激响应分别为 和 ，则级联后的系统的冲激响应 表示为两者的线性卷积，即注意：一般来说，由于卷积满足交换律，因而级联中的滤波器的顺序对整个冲激响应没有影响。两个稳定系统的级联仍是稳定系统。类似地，两个无源（无损）系统的级联仍是无源（无损）系统。 1h n 2hn h n4.5 简单互联方案简单互联方案 级联方案可用于生成逆系统。 如果级联系统满足 则LTI系统成为LTI系统

25、的逆系统，反之亦然。 逆系统的一个应用就是将从变形的结果中恢复出来。 1h n4.5 简单互联方案简单互联方案 举例离散时间累加器的冲激响应是单位阶跃序列 。 因此，逆系统必须满足条件 可以推出当 时， 且 因此，该系统的冲激响应为 n0n 20hn 4.5 简单互联方案简单互联方案l4.5.2 并联 两个并联的LTI离散时间系统的冲激响应分别为 和 ，则整个系统的冲激响应为 1h n 2hn 12h nh nhn4.5 简单互联方案简单互联方案 下面的离散时间系统由四个简单的离散时间系统互联而成，他们对应的冲激响应分别为 4.5 简单互联方案简单互联方案 总冲激响应为 其中4.6 LTI离散

26、时间系统的分类离散时间系统的分类 l4.6.1基于冲激响应长度的分类基于冲激响应长度的分类 若 具有有限长度，即 则它是一个有限冲激响应（FIR）离散时间系统，此时，卷积和简化为 上面的卷积和是有限的，可以直接用于计算 h n y n4.6 LTI离散时间系统的分类离散时间系统的分类滑动平均系统以及线性内插器都是FIR离散时间系统。若 无限长，则称其为无限冲激响应（FIR）离散时间系统。本书讨论的一类FIR滤波器是可以由线性常系数差分方程描述的因果系统。举例离散时间累加器 就是一个IIR系统。 h n yy1nnx n4.6 LTI离散时间系统的分类离散时间系统的分类4.6.2 基于输出计算过

27、程的分类非递归离散时间系统：若仅仅知道当前和过去时刻的输入样本就可以顺序计算出输出样本，则该系统为非递归离散时间系统。递归离散时间系统：若计算输出时除了需要知道当前和过去时刻的输入样本外，还需要知道过去时刻的输出样本，则该系统为递归离散时间系统。 4.6 LTI离散时间系统的分类离散时间系统的分类l4.6.3 基于冲激响应系数的分类实离散时间系统：冲激响应是实值序列的离散时间系统。复离散时间系统：冲激响应是复值序列的离散时间系统。4.7 LTI离散时间系统的频域表示离散时间系统的频域表示 在实际中遇到的多数离散时间信号，都可以表示为很多个甚至无限多个具有不同角频率的正弦离散时间信号的线性组合。

28、若已知LTI系统对单个正弦信号的响应，就可以利用系统的叠加性，求出对更复杂信号的响应。 4.7 LTI离散时间系统的频域表示离散时间系统的频域表示l4.7.1 频率响应 特征函数（Eigen Function）对于某些特定的输入信号，输出信号是输入信号乘以一个复常量。LTI系统的重要性质 冲激响应为 的LTI离散时间系统如图所示 h n4.6 LTI离散时间系统的分类离散时间系统的分类 输入输出关系由卷积和给出 若输入为复指数序列 则输出信号表示为 4.7 LTI离散时间系统的频域表示离散时间系统的频域表示 还可表示为 其中 对于复指数输出信号 ，LTI离散时间系统的输出是具有相同频率的复指数

29、信号乘以复常量 。 就是该系统的特征函数。 j ne()jH ej ne4.7 LTI离散时间系统的频域表示离散时间系统的频域表示l 是LTI离散时间系统的频率响应l 提供了系统的频域描述。l 是系统冲激响应 的傅里叶变换 是变量 的周期 的复函数，可用如下形式表示其中 称为LTI离散时间系统的幅度（幅频）响应 ()jH e()jH e()jH e h n()jH e2()jH e4.7 LTI离散时间系统的频域表示离散时间系统的频域表示 称为LTI离散时间系统的相位（相频）响应在很多实际应用中，离散时间系统的设计指标以幅度响应或相位响应或两者同时给出。在某些情况下，幅度函数会以分贝形式给出

30、其中 称为增益函数，称为 衰减函数或损失函数。 ( ) ( )g( )( )g 4.7 LTI离散时间系统的频域表示离散时间系统的频域表示l例4.31 滑动平均滤波器的频率响应考虑M点的滑动平均滤波器，其冲激响应表示为 频率响应表示为4.7 LTI离散时间系统的频域表示离散时间系统的频域表示 相应的幅度响应和相位响应为 其中 是 的阶跃函数，定义为 ( ) 4.7 LTI离散时间系统的频域表示离散时间系统的频域表示 Matlab中freqz(h,w)可用来求指定冲激响应量h在一组给定频率点w上的频率响应值。由这些频率响应值，还可以用函数real和imag计算实部和虚部，或用函数abs和angl

31、e计算幅度和相位。 4.7 LTI离散时间系统的频域表示离散时间系统的频域表示幅度响应在 处具有最大值1，而在 处为0；相位响应在 的每个零点表现处相差 的不连续，而在其他地方呈线性且斜率为 ；由于计算了正反切函数，当计算相位响应时可能出现多个 的跳变，可通过在跳变点增加的 倍数来展开相位从而将相位函数写为 的连续函数。Matlab中unwrap函数可用于上述目的。02/,1,2,/ 2k M kM()jH e(1)/ 2M224.7 LTI离散时间系统的频域表示离散时间系统的频域表示l4.7.5 稳态响应和瞬态响应 一个因果稳定LTI系统，当以振幅为常量的正弦序列作为输入时产生的输出，将在某

32、一个时刻之后具有稳态输出，它也是振幅为常量并与输入具有相同角频率的正弦序列。 具有实冲激响应 的LTI离散时间系统的频率响应为 ，LTI系统的输入为 h n()jH e4.7 LTI离散时间系统的频域表示离散时间系统的频域表示 将 表示为两个复指数序列之和的形式，即 其中 输入 的响应 为 输入 的输出是 的复共轭，即 x n01 2jnjg nAe e g n v n g n v n4.7 LTI离散时间系统的频域表示离散时间系统的频域表示 输出 的表示为 由于 是实的，则有 也可表示为 y n h n()()jjH eH e( )() y n4.7 LTI离散时间系统的频域表示离散时间系统的频域表示l4.7.7 滤波的概念数字滤波器：让输入序列中的某些频率分类没有任何失真的通过，同时阻止其他频率分类通过的LTI离散时间系统滤波的关键将任意一个输入序列表示为无穷多个指数序列的线性加权和，通过适当选择LTI数字滤波器在对应于输入信号正弦分量的频率上的幅度函数值，相对于其它频谱分量，某些频谱就会被选择性地严重削弱或滤除。 4.7 LTI离