测试信号分析与处理第2版宋爱国刘文波王爱民电子课件第4章节离散傅里叶变换及其快速算法.pptVIP

一、圆卷积的计算方法 注意:式（5-15）表示序列 不动，序列 对 进行了求模运算。为循环矩阵，其元素排列是有规律的，第一行表示的是原点（ ）不动时，序列的反折（倒序）；接下来的各行分别是上一行序列的循环右移。 一、圆卷积的计算方法 4．时域圆卷积定理法 利用DFT的一个重要性质――时域圆卷积定理来计算圆卷积。若 ， 则 二、 圆卷积与线卷积的关系 线卷积具有明确的物理意义，直接计算比较复杂。对于两个有限长序列求线卷积能否用圆卷积来代替，即采用FFT计算线卷积而使两者结果又完全相同呢？答案是肯定的，但需要满足一个条件：就是将进行线卷积的两序列的长度（设两序列的点数分别为 ）均通过补零的办法，加长至 ，然后再进行N点的圆卷积，则圆卷积的结果与线卷积的结果相同。 二、 圆卷积与线卷积的关系 设 分别由点通过补零，加长至N点，其线卷积为，可表示为 计算结果的长度要多出一些零值，但非零值长度仍为 点。 其圆卷积为 ，可表示为 二、 圆卷积与线卷积的关系 而 可得 （5-16） 式中，下标p表示序列的周期化； 是指对线卷积 进行周期为N的延拓后得到的周期序 列； 两序列的圆卷积的结果，是的主值序列。 二、 圆卷积与线卷积的关系 上述过程说明:加长至N点长的 、 两序列的圆卷积 与 线卷积作周期延拓所得到的序列 的主值序列相同。在这个条件下（两序列均加长至N点），就可以通过计算序列的圆卷积来求解线卷积。从式（5-16）的推导过程还可以看出，如果两序列不加长至N，其线卷积的周期延拓序列将发生重叠或混叠现象（因为 线卷积长度为 ），相应计算出的圆卷积也将产生失真，圆卷积的主值序列和线卷积就不相同。 三、 用FFT计算有限长序列的线卷积 根据上述圆卷积与线卷积的关系，可以得出用FFT求解两序列线卷积的原理框图，如图5-9 所示。其计算的具体步骤如下 1）若两序列 、 的长度为N，将序列加长至2N－1，并应修正为2的幂次（基2算法）； 2）计算 、 ； 3）计算 ； 4）计算 。 图5-9 用FFT求线性卷积 三、 用FFT计算有限长序列的线卷积 在MATLAB中直接实现线卷积计算的函数有 conv, conv2, convn。其中conv2和 convn分别用于2维、n维的卷积运算。conv则用于向量卷积与多项式乘的计算，调用的格式为c＝conv（a，b）。式中，a、b表示两个序列，c＝a*b。在MATLAB中，序列可用向量来表示，若向量a的长度为na，向量b的长度为nb，则向量c的长度为na＋nb－1。 四、分段快速卷积——重叠相加法 分段快速卷积的方法:将长序列分成若干小段，每小段分别与短序列作卷积运算，然后将所有的分段卷积结果相叠加，就是线卷积的最后结果，这种方法又称为重叠相加法。 四、分段快速卷积——重叠相加法 设 的长度为M， 为一长序列，将 进行分段，每段的长度为 ，将每一段分别与 进行线卷积，然后将结果重叠相加，如图5-10所示。 图5-10 叠加相加法的分段以及的重叠情况 四、分段快速卷积——重叠相加法 设将 分为 第 段 表示为 （5-17） 则 （5-18） 四、分段快速卷积——重叠相加法 由于 长度为 ， 的长度为M，故 的长度为 ，即的范围为 （5-19） 将式（5-19）与式（5-17） 的范围比较， 显然 比长点，而 的范围是 （5-20） 四、分段快速卷积——重叠相加法 将式（5-19）与式（5-20）比较，可知的