有限离散傅里叶变换的循环褶积.PPT

第五章 有限离散傅里叶变换 一、有限离散傅里叶变换的循环褶积 ——循环褶积定理1 二、 循环褶积与普通褶积的关系 ——循环褶积定理2 三、利用快速傅里叶变换（FFT）计算褶积 * §3 有限离散傅里叶变换的循环褶积 一、有限离散傅里叶变换的循环褶积 ——循环褶积定理1 在第四章， 我们知道， 两个离散信号的连续频谱 相乘，对应于两个离散信号的褶积运算。 这就是滤波 理论的基本关系式。 为了能用快速傅里叶变换算法在频谱域上实现滤波， 我们必须分析两个有限离散频谱相乘所对应的信号究 竟是什么？ 这就是本节要讨论的问题。 设有两个有限离散信号 则有 这里 但注意到上式右边表达式算出来的是以 为周期 函数， 因此从理论上讲上述 可以取任何整数。 也就是说， 我们可以将 扩充称以 为 周期的函数， 这种扩充便于下面的讨论。 令 则 是一个新的有限离散频谱， 它 对应的有限离散信号为 问题1 与 有何关系？ 即 称之为有限离散信号 和 的循环褶积或周期褶积， 记作 即 这里 也以 为周期， 且 循环褶积定理1 设 和 是长度为 的有限离散 信号， 将 扩充成以 为周期的函数， 则有限离 散频谱 所对应的有限离散信号 注 （1）若令 则有 据此可知， 循环褶积与普通褶积的几何含义类似， 只是要注意这里的 已扩充成以 为周期的函数. （2）循环褶积也可以用两个同心圆来表示： 将两个同心圆的圆周 等分， 用分点代表离散信号， 并将 顺时针排列， 逆时针排列。 若计算 则 将 与 相对应， 然后各对应项相乘，再相加即 可得 参见下图。 例1 设 和 分别是长度为 和 的信号： 将 扩充成长度为 的信号： 求 和 的循环褶积 解 我们可以用同心圆，或类似于普通褶积的图形 来计算， 也可以直接用循环褶积定义计算。 总之， 二、循环褶积与普通褶积的关系 ——循环褶积定理2 设有两个有限离散信号 则 与 的普通褶积 为 因为当 时 当 时 所以当 或 时 故 另一方面， 在（2）式中， 当 即 时， 而 故当 时所有 从而 又当 即 时， 而 故当 时所有 从而 总之， 即 是长度为 的有限离散信号。 设 与 是以 为周期 的信号， 并且 与 的关系为： 当 时， 记 与 的循环褶积为 问题2 1）在 范围内， 与 有何 关系？ 2）在什么条件下，有 下面分别讨论这两个问题： 1）当 时， 因为当 时， 所以 从而 又因为当 时， 所以 从而 总之， 2）当 中的每一项 时， 因为当 时， 从而 所以当 时， 才可能非零。 故 又因为 而当 即 时， 而 故当 时， 从而每一项 此时 总结以上讨论结果，可得 循环褶积定理2 设有两个有限离散信号 与 是以 为周期的信号， 并且 当 时， 则当 满足 时， 褶积 与循环褶积 相等， 即 注 （1）若 即 时， 于是有 这时， 循环褶积 可以完全反映褶积 故滤波 （2）若 即 时， 不成立， （3）当 时， 与 部分相等。 例2 设 与 的普通褶积 以 为周期的信号 的循环褶积 循环褶积定理2中的条件 这里 故 这与计算结果相同。 当 时， 作业：P131 8 三、利用快速傅里叶变换（FFT）计算褶积 利用快速傅里叶变换计算褶积的理论根据是循环 褶积定理2， 根据这个定理， 选取适当的 就可把褶积运算转化为循环褶积运算。 再根据循环 褶积定理1， 循环褶积运算又可转换为有限离散频谱 相乘， 而有限离散频谱可利用快速傅里叶变换算法 来计算。 下面具体说明之。 1. 在一般情况下利用FFT计算褶积的方法 设 在一般情况下， 利用FFT计算褶积 的步 骤如下： 1）选择 使之满足 2）构造以 为周期的信号 使之满足当 时， 3）利用FFT分别计算 的有限离散频谱 4）将 与 相乘，得 5）利用IFFT对 作反变换，得 则 就是我们所要计算的褶积 2. 在特殊情况下利用FFT计算褶积的方法 在特殊情况下， 即 时， 点， 我们可以分若干段来计算褶积。 根据循环褶积的特 具体步骤参见 P123. 实验一：滤波的MATALAB实现 一、实验目的 （1）掌握滤波在时域的表现形式和频域的表现形式。 （2）掌握褶积（也叫卷积）的MATALB实现。 （3）掌握快速傅里叶变换（FFT）的MAT