并导出频域采样定理.PPTVIP

由此可见，循环卷积既可以在时域直接计算，也可以按照图3.4.1所示的计算框图在频域计算。由于DFT有快速算法，当L很大时，在频域计算循环卷积的速度快得多，因而常用DFT(FFT)计算循环卷积。图3.4.1用DFT计算循环卷积的原理框图在实际应用中，为了分析时域离散线性时不变系统或者对序列进行滤波处理等，需要计算两个序列的线性卷积。与计算循环卷积一样，为了提高运算速度，也希望用DFT(FFT)计算线性卷积。而DFT只能直接用来计算循环卷积，因此，下面先导出线性卷积和循环卷积之间的关系以及循环卷积与线性卷积相等的条件，最后得出用图3.4.1线性卷积的条件。假设h(n)和x(n)都是有限长序列，长度分别是N和M。它们的线性卷积和循环卷积分别表示如下：其中所以(3.4.1)(3.4.2)L对照(3.4.1)式可以看出，上式中即(3.4.3)(3.4.3)式说明，yc(n)等于yl(n)以L为周期的周期延拓序列的主值序列。我们知道，yl(n)长度为N＋M－1，因此只有当循环卷积长度L≥N＋M－1时，yl(n)以L为周期进行周期延拓时才无时域混叠现象。此时取其主值序列显然满足yc(n)=yl(n)。由此证明了循环卷积等于线性卷积的条件是L≥N＋M－1。图3.4.2中画出了h(n)、x(n)、h(n)*x(n)以及L分别取6、8、10时h(n)Lx(n)的波形。由于h(n)长度N=4，x(n)长度M=5，N＋M－1=8，因此只有L≥8时，h(n)Lx(n)波形才与h(n)*x(n)相同。图3.4.2线性卷积与循环卷积波形图综上所述，取L≥N＋M-1，则可按照如图3.4.1所示的计算框图用DFT(FFT)计算线性卷积。其中DFT和IDFT通常用快速算法(FFT)来实现，故常称其为快速卷积。实际上，经常遇到两个序列的长度相差很大的情况，例如MN。若仍选取L≥N＋M－1，以L为循环卷积区间，并用上述快速卷积法计算线性卷积，则要求对短序列补很多零点，而且长序列必须全部输入后才能进行快速计算。因此要求存储容量大，运算时间长，并使处理延时很大，不能实现实时处理。况且在某些应用场合，序列长度不定或者认为是无限长，如电话系统中的语音信号和地震检测信号等。显然，在要求实时处理时，直接套用上述方法是不行的。解决这个问题的方法是将长序列分段计算，这种分段处理方法有重叠相加法和重叠保留法两种。下面只介绍重叠相加法。设序列h(n)长度为N，x(n)为无限长序列。将x(n)等长分段，每段长度取M，则于是，h(n)与x(n)的线性卷积可表示为式中（3.4.4a）（3.4.4b）（3.4.4c）(3.4.4b)式说明，计算h(n)与x(n)的线性卷积时，可先计算分段线性卷积yk(n)=h(n)*xk(n)，然后把分段卷积结果叠加起来即可如图3.4.3所示，每一分段卷积yk(n)的长度为N＋M－1因此相邻分段卷积yk(n)与yk＋1(n)有N－1个点重叠，对于重叠点的序列值必须把重叠部分的yk(n)与yk＋1(n)相加，才能得到正确的卷积序列y(n);当第二个分段卷积y1(n)计算完后，叠加重叠点便可得输出序列y(n)的前2M个值；同样道理，分段卷积yi(n)计算完后，就可得到y(n)第i段的M个序列值。显然,每一段可用图3.4.1所示的快速卷积法计算分段卷积yk(n)，其中L=N＋M－1。图3.4.3用重叠相加法计算线性卷积时域关系示意图因此，这种方法不要求大的存储容量，且运算量和延时也大大减少，最大延时TDmax=2MTs+To，Ts是系统采样间隔，To是计算1个分段卷积所需时间，一般要求ToMTs。这样，就实现了边输入边计算边输出，如果计算机的运算速度快，可以实现实时处理。用DFT计算分段卷积yk(n)的方法如下：(1)i=0；L=N＋M－1；计算并保存H(k)=DFT[h(n)]L;(2)读入xk(n)=x(n)RM(n－kM)，构造变换区间［0，L－1］上的序列