s实验六_FFT变换.docVIP

实验四 FFT实验 一 实验目的 了解FFT的原理； 了解在DSP中FFT的设计及编程方法； 熟悉对FFT的调试方法 二 实验内容 本实验要求使用FFT变换求一个时域信号的频域特性，并从这个频域特性求出该信号的频域值。使用DSP汇编语言实现对FFT的DSP编程。 三 实验原理 对于有限长离散数字信号{x[n]},0，它的频谱离散数学值{X(K)}可由离散傅立叶变换（DFT）求得。DFT定义为： 也可以方便的把它改写成如下形式： 式中（有时简写为）代表。不难看出，是周期性的，且周期为N，即 m,l=0, 的周期性是DFT的关键之一。为了强调起见，常用表达式取代以便明确地给出地周期为N。 由DFT地定义可以看出，x[n]为复数序列地情况下，完全可以直接运算N点DFT需要次复数乘法和N(N-1)次复数加法。因此，对一些相当大的N值（如1024点）来说，直接计算它的DFT所需要的计算量是很大地。一个优化的实数FFT算法是一个组合后的算法。原始的2N个点的实输入序列组合成一个N点的复序列，然后对复序列进行N点的FFT运算，最后再由N点复数输出拆散成2N点的复数序列，这2N点的复数序列与原始的2N点的实数输入序列的DFT一致。 FFT的基本思想在于：将原来的N点序列分成两个较短的序列，这些序列的DFT可很简单地组合起来得到原序列的DFT。例如，若N为偶数，将原有的N点系列分成两个（N/2）点序列，那么计算N点的DFT将只需要[（N/2）2*2]=N2/2次复数乘法，即比直接计算少作一半乘法。因子(N/2)2表示直接计算（N/2）点DFT所需要地乘法次数，而乘数2代表必须完成两个DFT。上述处理方法可以反复使用，即（N/2）点的DFT计算也可以化成两个（N/4）点的DFT（假定 N/2为偶数），从而又少作一半乘法，使用这种方法，在组合输入和拆散输出地操作中，FFT地运算量减半。这样，利用实数FFT算法来计算实输入序列地DFT的速度几乎是一般复FFT算法的两倍。 假定序列{x[n]}的点数N是2的幂，按上述处理方法，定义两个分别为x[n]的偶数项和奇数项的(N/2)点序列x1[n]和x2[n],即： x1[n]=x[2n] n=0,1,…,(N/2)-1 x2[n]=x[2n+1] n=0,1,…,(N/2)-1 {x[n]}的N点DFT可写成： + ＝＋ 因考虑到可写成： ＝[]== 故X(k)可写为： + =X1(k)+ X2(k) 式中X1(k)和X2(k)是x1(n)和x2(n)的（N/2）点DFT。上式表明，N点DFT 可分解为按上式地规则加以组合地两个（N/2）点DFT。 四FFT的DSP编程 FFT地实现步骤主要有以下四步： 将输入序列压缩和位倒序。 N点的复数FFT。 奇数号部分和偶数号部分分离。 产生最后的输出数据。 最初，将原始输入的2N点的实序列a(n),存储在4N大小的数据处理缓冲区的下半部分。 .asg AR1,FFT_TWID_P LD #FFT_DP,DP STM #SYSTEM_STACK,SP rfft_task: STM #FFT_ORIGIN,AR3 ; to data STM #data_input,AR4 RPT #K_FFT_SIZE*2-1 ；K_FFT_SIZE=128 MVDD *AR4+,*AR3+ 0C00h 0C01h … 0CFEh 0CFFh 0D00h a(0) 0D01h a(1) … … 0DFFEh a(254) 0DFFFh a(255) 然后通过子程序调用，有以下几步计算 ①驱除原始数据的直流分量。 由于采样的需要，信号发生器输出的为叠加了直流分量的正弦信号，所以有此一步。 .asg AR3,FFT_INPUT .def data_ave data_ave: STM #FFT_ORIGIN,FFT_INPUT ;AR3 1st originalinput LD #0,A RPT #K_FFT_SIZE*2-1,A ADD *FFT_INPUT+,A ;累加原始数据到A NOP LD A,#-1-k_LOGN ;A右移8