最新信号处理导论中文版.docx

PAGE PAGE # PAGE PAGE # 第三章离散系统 本章的讨论重点是离散系统，尤其是 离散线性时不变系统。线性时不变系统的输入输出 (I/O)方 程可以用输入信号与系统冲激响应的离散卷积来表示。 根据系统的冲激响应是否是有限延时还是无限延时可以分为 有限冲激响应 (FIR )和无限冲激响 应(IIR)两种。本章的主要目的是为 FIR滤波器设计算法。FIR滤波算法可以分为按块(Block to Block) 和样值处理(Sample to Sample)算法两种。 分批处理算法中，输入信号视为一次抽样的块。将这一块信号与滤波器冲激响应卷积得到一个输 出块。 如果输入序列时限非常长或者是无限延时，这种方法需要做些改进，比如说可以将输入信号分成 多个块，每一块的长度都可以分别处理，可以一次滤波一块，然后再把输出拼凑在一起。 样值处理算法中，一次只处理一个抽样。滤波器可以看作是一台状态机器，也就是说，把输入抽 样与滤波器当前的状态结合起来计算当前的输出抽样， 同时也更新滤波器的内部状态为下一次处理作准 备。 当输入信号特别长的时候，这种方法对于实时运算特别有效。滤波器自身特性变化的自适应滤波 就适合于使用这种算法。目前的 DSP芯片对这种算法也很有效。 §3.1输入输出规则 离散系统所实现的就是将输入的离散抽样序列 x(n),根据一定的输入/输出(I/O)规则转换成输 出序列的运算。I/O规定了怎样由已知的输入计算输出。 样值处理方法，我们可以认为其 I/O规则就是一次处理一个输入抽样。 [xi, x2, x 丄，x n,L] - ' [ y , y , y 丄，y 丄] 按块处理的方法，输入序列划分成块，每次处理一块。 因此其I/O规则也就是将输入向量根据某种函数映射成输出向量。 y = H〔x 1 对于线性系统，这种映射就是用矩阵 H作线性变换。线性定常系统，其变换矩阵 H根据系统的冲 激响应有特定的结构。 例 3.1.1 例 3.1.2 y(n) = 2x(n) 3x(n-1) 4x(n-2) n时刻的输出是此前连续三个输入抽样的加权和。也就是说， n时刻，线性系统必须记住前两个时 刻的抽样 x(n-1)、x(n-2)。 例3.1.3将长度为L=4的输入抽样‘X。, xi, X2 , X3 /视为一块，例3.1.2所示的线性系统将其转换成长 度为6的输出序列。 「y。1 2 1 0 0 01 I y1 1 3 2 0 0 jX0 1 y = y2 4 3 2 0 |x =Hx y3 0 4 3 2 |X2 y I 0 0 4 1 3 lx _ y4 J P 0 0 4 一 输出序列的长度比输入序列长度大 2，因为系统必须保存两个抽样，最后的两个输出可以认为是 输入消失后(input-off)的过渡状态。如果输入的抽样为 L=5,那么，输出的序列为： ]y° 2 0 0 0 01 y1 3 2 0 0 0 1 y2 4 3 2 0 0 1 X 1 y = y3 ■ 0 4 3 2 0 X2 =Hx y4l 0 0 4 3 2 X 3 y5 1 0 0 0 4 3 一X4 一 Li _y6 一 1 0 0 0 0 4 例3.1.4、例3.1.2的输入输出方程也可以用下列样值处理的算法来实现： y( n)=2x( n)+3w 1(n )+4w2 (n) W2( n+1)=w 1(n) W1( n+1)=x (n) 附加的W1(n)、W2(n)可以视为系统的内部状态。当前的输入结合当前的内部状态足以计算当前的 输出。由有下一个输入x(n +1)所产生的输出y(n +1)要求我们知道已经更新的内部状态。而此时的内部状 态(n+1时刻的内部状态)已经更新。也就是说， n+1时刻，我们有： y(n +1)= 2x(n+1)+3w 1(n+1)+4w 2(n+1) W2( n+2)=w 1(n+1) W1( n+2)=x( n+1) 这样的计算是从某个时刻开始并且不断重复，我们可以归结为以下算法： 一旦内部状态的当前值在计算输出 y的时候使用过以后，他 for each new in put x do: y ： = 2x+3w1+4w2W2： =W1w for each new in put x do: y ： = 2x+3w1+4w2 W2： =W1 w1: =x W2｝更新的次序非常重要，也就是首先更新 W2，接下来更新 W1, 以避免把正确的值覆盖。 例3.1.2、例3.1.3、例3.1.4是同一个离散系统的等效描述方式。究竟是采用哪一种形式取决于应 用的场所，也就是要看输入序列是有限长还是无限长、输入抽样是否在接收到以后应该立刻处理还是可 以延缓处理。 上面的例子实际上是用下述 I/O方程描述的、具