快速傅立叶变换(fft)算法实验.docVIP

快速傅立叶变换(FFT)算法实验 实验目的： 1. 加深对DFT 算法原理和基本性质的理解； 2. 熟悉FFT 的算法原理和FFT 子程序的算法流程和应用； 3. 学习用FFT 对连续信号和时域信号进行频谱分析的方法。 程序流程图: 傅立叶变换是一种将信号从时域到频域的变换形式，是声学、语音、电信和信号处理等领域中的一种重要分析工具。离散傅立叶变换（DFT）是连续傅立叶变换在离散系统中的表现形式，由于DFT 的计算量很大，因此在很长时间内其应用受到很大的限制。快速傅立叶变换（FFT）是离散傅立叶变换的一种高效运算方法。由于我们在计算DFT 时一次复数乘法需用四次实数乘法和二次实数加法；一次复数加法则需二次实数加法。运算一个X（k）需要4N 次复数乘法及2N+2（N-1）=2（2N-1）次实数加法。所以整个DFT运算总共需要4N^2 次实数乘法和N*2(2N-1)=2N(2N-1)次实数加法。如此一来，计算时乘法次数和加法次数都是和N^2 成正比的，当N 很大时，运算量是可观的，因而需要改进对DFT 的算法减少运算速度。 我们可以将DFT 运算中有些项合并。我们先设序列长度为N=2^L，L 为整数。将N=2^L 的序列x(n)(n=0,1,……，N-1)，按N的奇偶分成两组，也就是说我们将一个N 点的DFT 分解成两个N/2 点的DFT，他们又重新组合成一个如下式所表达的N 点DFT： 一般来说，输入被假定为连续的。当输入为纯粹的实数的时候，我们就可以利用左右对称的特性更好的计算DFT。这样的RFFT 优化算法是包装算法：首先2N 点实数的连续输入称为“进包”。其次N 点的FFT 被连续运行。最后作为结果产生的N 点的合成输出是“打开”成为最初的与DFT 相符合的2N 点输入。使用这一思想，我们可以划分FFT 的大小，它有一半花费在包装输入O（N）的操作和打开输出上。 AD原理图 源程序 #include math.h #include DSP2833x_Device.h // DSP2833x Headerfile Include File #include DSP2833x_Examples.h // DSP2833x Examples Include Fi #include comm.h #include i_cmplx.h #include ext_inf.h /*********************************************************************************/ #if (CPU_FRQ_150MHZ) // Default - 150 MHz SYSCLKOUT #define ADC_MODCLK 0x3 // HSPCLK = SYSCLKOUT/2*ADC_MODCLK2 = 150/(2*3) = 25.0 MHz #endif #if (CPU_FRQ_100MHZ) #define ADC_MODCLK 0x2 // HSPCLK = SYSCLKOUT/2*ADC_MODCLK2 = 100/(2*2) = 25.0 MHz #endif #define ADC_CKPS 0x1 // ADC module clock = HSPCLK/2*ADC_CKPS = 25.0MHz/(1*2) = 12.5MHz #define ADC_SHCLK 0xf // S/H width in ADC module periods = 16 ADC clocks #define AVG 1000 // Average sample limit #define ZOFFSET 0x00 // Average Zero offset #define BUF_SIZE 160 // Sample buffer size /*********************************************************************************/ ///////////////// unsigned int SampleLong; #define SAMPLELONG 3 //////////////////// unsigned int Ad_data[1536]={0}; unsigned int Ad_data1[1536]={0}; unsigned int convcount = 0; volatile unsigned int adcon