实验二FFT与DFT计算时间的比较及圆周卷积代替线性卷积的.pptVIP

实验二 FFT 与 DFT 计算时间的比 较及圆周卷积代替线性卷积的有 效性实验 一 实验目的 二 实验内容及要求 三 预做实验 一 实验目的 ? 1 ：掌握 FFT 基 2 时间（或基 2 频率）抽选法， 理解其提高减少乘法运算次数提高运算速度 的原理。 ? 2 ：掌握 FFT 圆周卷积实现线性卷积的原理 二 实验内容及要求 ? 1. 对 N=2048 或 4096 点的离散时间信号 x(n) ，试用 Matlab 语言编程分别以 DFT 和 FFT 计算 N 个频率样值 X(k), 比较两者所用 时间的大小。 ? 2. 对 N/2 点长的 x(n) 和 N/2 点长的 h(n) ，试 用 Matlab 语言编程实现以圆周卷积代替线 性卷积，并比较圆周卷积法和直接计算线 性卷积两者的运算速度。 三预做实验 ? 1.FFT 与 DFT 计算时间的比较 （ 1 ） FFT 提高运算速度的原理 （ 2 ）实验数据与结论 ? 2. 圆周卷积代替线性卷积的有效性实验 （ 1 ）圆周卷积代替线性卷积的原理 （ 2 ）实验数据和结论 FFT 提高运算速度的原理 ? FFT 算法将长序列的 DFT 分解为短序列的 DFT 。 N 点的 DFT 先分解为 2 个 N/2 点的 DFT ， 每个 N/2 点的 DFT 又分解为 N/4 点的 DFT ，等 等。最小变换的点数即所谓的“基数”。因 此，基数为 2 的 FFT 算法的最小变换（或称蝶 型）是 2 点的 DFT 。一般地，对 N 点 FFT ，对 应于 N 个输入样值，有 N 个频域样值与之对应。 ? 以基 2 时间按抽选法为例，分解递推公式为： 1 2 1 2 ( ) ( ) ( ) 0, 1 2 ( ) ( ) ( ) 2 k N k N X k X k W X k N k N X k X k W X k ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? X1(K) 和 X2 （ K ）分别对应时间序列 x(n) 的的偶、奇序 列 N/2 点长的 DFT ，以此方法分解下去，得到两点长 N/2 为组的离散时间分组结果。在一次由两点 DFT 复合 产生 4 、 8 、 16…… 到 N 点的最后结果。运算量由 N*N 次 降为 次乘法。 2 log 2 N N 实验数据与结论 ? 令 N 为不同长度的序列时， DFT 和 FFT 的耗时比较： ? N=512 点时： dft_cost_time =0.3750s fft_cost_time =0.0150s ? N=1024 点时： dft_cost_time =2.9220s fft_cost_time =0.1250s ? N=2048 点时： dft_cost_time =23.5000s ； fft_cost_time =0.4690s ? N=4096 点时： dft_cost_time =305.9530s fft_cost_time =103.9220s ? 有上面数据表明，同样长度的信号， DFT 耗时要比 FFT 耗时要少，表明 FFT 算法的有效性。 圆周卷积代替线性卷积的原理 ? 线性卷积的长度及运算量 ； ? 用 FFT 算的步骤； ? 比较两者乘法运算量 ； 线性卷积的长度及运算量 ? 设一离散线性移不变系统的冲激响应为 , 其输入信号 为 . 其输出为 . 并且的长度为 L 点 , 的长度为 M 点 , 则： 1 0 ( ) ( ) ( ) ( ) ( ) L m y n x n h n x m h n m ? ? ? ? ? ? ? 1 0 ( ) ( ) ( ) ( ) ( ) L m y n x n h n x m h n m ? ? ? ? ? ? ? ? ? ? ? 1 h n h M n ? ? ? ? 则 y(n) 的长度为 L+M-1 。 线性卷积的乘法运算量为 LM 次，又由于 FIR 滤波 器 h(n) 的特性是对称序列即 因而，乘法运算次数减半，即为 LM/2 。 用 FFT 算的步骤 1. ( ), ( ) 1 ; x n h n N M L ? ? ? 将 补零点，至少为 点 ? ? 2. ( ) ( ) H k FFT h n ? 求 ； ? ? 3. ( ) ( ) X k FFT x n ? 求 ； 1. ( ), ( ) 1 ; x n h n N M L ? ? ? 将 补零点，至少为 点 4. ( ) ( ) ( ) Y k X k H k ? 求 ； ? ? 5. ( ) ( ) y n IFFT Y k ? 求 。 流程图 ( ) x n ( ) h n ( ) X k IFFT FFT FFT ( ) X k (