《计算物理(本科)》[第5章].pptVIP

第五章 快速傅立叶变换(FFT) §5.1 离散傅立叶变换（DFT） 一、离散傅立叶变换 在有限多个等距、离散数据的基础上的傅立叶变换称为离散傅立叶变换。 设已知函数f(x)在区间0≤x2π上N个等分点 2π l/N，(l=0,1,2,…,N-1)处的值已知fl 。实际上，N分点的值fl常常是采样的数据。 设采样周期T=2π/N，采样值则为fl =f(xl )= f(lT ) ，(l=0,1,2,…,N-1) §5.2 快速傅立叶变换（FFT） 快速傅立叶变换（FFT）算法 傅立叶变换公式中具有 e- i2klπ/N 的周期性，在N个Fk中，实际只有N个不同的项，即1, e- i2π/N ,…, e- i2(N-1)π/N ，或ω0 , ω1,ω2,…, ω (N-1) ，这些都是不相等的。如果把Fk公式中各项先按ω归类，合并同类项，则可以大大减少运算量。特别是当N=2m （m为正整数）时，运算次数减少到Nlog2N。这种方法就称为快速傅立叶变换（FFT）算法 采样点条件：N=2m （m为正整数）。这一条件不难满足，可以大大简化编程和提高计算效率。 §5.3 快速傅立叶变换的应用 一、 MATLAB的实现 快速傅立叶变换（FFT）算法不需要自编程序，MATLAB自身提供了若干个函数。最常用的函数有快速傅立叶变换函数fft()、快速傅立叶逆变换函数ifft(