第四章快速傅里叶变换(FFT)汇编.pptxVIP

第四章 快速傅里叶变换（FFT）;Contents;从DFT到FFT的缘由;从哪个方面减少运算量呢？;浅谈FFT实现原理;算法理论推导： 设序列 的长度为N，且满足N=2M ，M为自然数。 按n奇偶把 分解为两个N/2点的子序列 则x(n)的DFT为 ; 因为 所以 ;所以X(k)又可表示为 ;8点DFT一次时域抽取分解运算流图; 完成一次蝶形运算需要一次复数乘法和两次复数加法运算。因此在8点DFT一次时域抽取分解中，共需要计算两个N/2点DFT运算和N/2个蝶形运算。 所以按照上图的计算DFT时，总的复数乘法次数为 复数加法次数为 ; 类似地，我们将 按奇偶分解成两个N/4点子序列 ，其表达式分别如下： ;式中 同理，由 的周期性和 的对称性 最后得到：;用同样的方法可以计算出 其中;8点DFT二次时域抽取分解运算;8点DIT-FFT运算流图; 从上面的蝶形算法，当 时，其运算应该有M级蝶形，每一级都由N/2个蝶形运算构成。因此每一级运算都需要N/2次复数乘法和N次复数加法，所以M级蝶形运算的乘法次数为： 加法次数： ;一个简单的算例 计算序列 的8点DFT。分别用基本的DFT算法和FFT算法实现，体会计算过程中的时间复杂度。 当然针对计算机来说，计算乘法所需要的时间远大于加法。一般的计算机，差不多相差十倍左右。 用计算机产生随机的1024个点构成序列，然后取N=1024.此时计算时间差距就会加大。N=2048时，时间差距会更大。 ;为了更进一步的体会FFT的物理意义，引入一个算例： 假设对某信号经过ADC之后，得到一个序列 ，此时我们不知道其具体的函数表达式。但是我们可以对其做FFT运算。现在我们做一个验证： 假设我们有一个信号，它含有2V的直流分量，频率为50Hz、相位为-30度、幅度为3V的交流信号，以及一个频率为75Hz、相位为90度、幅度为1.5V的交流信号。用数学表达式就是如下： S=2+3*cos(2*pi*50*t-pi*30/180)+1.5*cos(2*pi*75*t+pi*90/180) 现在对其进行变换，取样点为1024，采样频率为1024Hz 注意这里的取样频率只要满足原始频率的2倍即可，且取样点和取样频率根据频率分辨率来选取。 ; 从图中的结果可以得出，当频率为0、50、75时，对应的幅度值依次为2、3、1.5，相位依次为0、90、-30。当然，这看上去几乎是没有误差的，因为我们取样频率和所取点数比较大，当N=Fs=256时，会存在误差，但是这个误差完全不影响我们对信号函数的分析与判断。从而进一步的验证了DFT与FFT算法的正确性。 ;上图是从模拟信号到频域离散信号的完整的过程。这其中对应有几个概念容易混淆。因此对此做出区分：; 数字频率、模拟频率与采样频率：模拟频率通常用 表示，数字频率用 表示。此时的数字频率主要是针对序列而言，因此没有采样频率 就不会有数字频率的概念，所以数字频率与采样频率和模拟频率一定满足某种关系。我们知道有如下过程： 因此 其中 为取样间隔。对应的采样频率; DFT取样点N与采样频率和频率分辨率（步进值）的关系 首先我们根据奈奎斯特定律得到： ， 为连续信号的截至（最高）频率（因为它可能有很多频率成分）。现在问题就是取样点N的选取？现在我们从序列的FT说起！ FT完成的任务就是对时域连续信号抽样之后通过序列的傅里叶变换得到的频谱关系，我们知道这个谱是连续的。此时的谱是一个周期谱，那么周期是多少呢？是不是与时域采样频率有关系呢？答案是肯定的，此时的周期就是 （注意此时是频域，频谱图中的横坐标表示的频率）。 ; 对于周期连续的谱，计算机还是无法操作，因此我们还得对FT变换之后的谱进行频域采样（DFT的实质）。类似于时域采样，频域FT谱线的采样频率为多少合适呢？显然此时只需一个周期即可，即在一个周期中取N个点，满足 ： 式中 为频率分辨率， 为采样频率。 因此关于频率分辨率和取样点N的关系就得出来了。 ;首先必须得画出序列对应的FT谱线图，横坐标为 ，纵坐标为频率。此时的FT谱线应该是对上图的周期延拓。这就是我们一般只取其主值[-pi,pi]得原因。;得到FT谱