快速傅立叶变换(FFT)算法实验

快速傅立叶变换(FFT)算法实验实验目得:1、加深对DFT算法原理与基本性质得理解；电信与信号处理等领域中得一种重要分析工具。离散傅立叶变换(DFT)就就是连续傅立叶变换在离散系统中得表现形式,由于DFT得计算量很大,因此在很长时间内其应用受到很大得限制。快速傅立叶变换(FFT)就就是离散傅立叶变换与二次实数加法；一次复数加法则需二次实数加法。运算一个X(k)需要4N次数与加法次数都就就是与N^2成正比得，当N很大时，运算量就就是可观得,L如下式所表达得N点DFT:用左右对称得特性更好得计算DFT。这样得RFFT优化算法就就是包装算法：首先2N点实数得连续输入称为“进包”。其次N点得FFT被连续运行。最后作为结果产生得N点得合成输出就就是“打开”成为最初得与DFT相符合得2NAD原理图#include<math、h>Examples、h"／/DSP2833xExamplesIncludeFi/*****＊********＊*＊**＊***＊****＊*******＊＊***＊＊**＊**＊*＊****＊*＊**＊******＊*＊*******＊***/##if(CPU_FRQ_150MHZ)//Default－150MHzSYSCLKOUT#defineADC_MODCLK0x3//HSPCLK＝SYSCLKOUT/2*ADC_MODCLK2＝150／(2＊3)=25、0MHz#if(CPU_FRQ_100MHZ)#defineADC_MODCLK0x2//HSPCLK=SYSCLKOUT/2*ADC_MODCLK2=#defineADC_CKPS0x1/／ADCmoduleclock=HSPCLK/2*A#defineADC_SHCLK0xf/／S/HwidthinADCmoduleperiods＝16ADCclocks#defineAVG1000//Averagesamplelimit#defineZOFFSET0x00／/AverageZerooffset#defineBUF_SIZE160//Samplebuffersize******＊＊******＊*************/unsignedintSampleLong;//////////／//／/／／//////／/////////／／/PLEXDDataBuffer[512]={0};／////／／//unsignedinti；／*＊**＊*************＊****＊*＊*****＊***＊＊**＊**＊**＊*＊**＊＊＊*******＊＊***＊＊*＊****＊＊*＊＊＊＊＊/***＊＊＊＊＊******＊*＊＊*＊*＊＊*＊*****＊***＊**********＊*/voidmain(void){ﻩﻩ#ifSAMPLELONG==1ﻩSampleLong=256;#endif＃ifSAMPLELONG=＝2ﻩSampleLong=512;ﻩ#endif#ifSAMPLELONG==3ﻩSampleLong＝1024;//ThisexamplefunctionisfoundintheDSP2833x_SysCtrl、ce、InitGpio();//Skippedforthisexampletable:DINT;stsarecleared、／/ThisfunctionisfoundintheDSP2833x_PieCtrl、cfile、errupts、//TheshellISRroutinesarefoundinDSP2833x_DefaultIsr、c、cInitPieVectTable();//ISRfunctionsfoundwithinthisfile、isterPieVectTable、ADCINT=&adc_isr；EDIS;//ThisisneededtodisablewritetoEALLOWprote//InitPeripherals();//NotrequiredforthisexampleInitAdc();/／Forthisexample,inittheADC//EnableADCINTinPIEEINT；／／EnableGlobalinterruptINTMERTM；/／EnableGlobalrealtimeinterruptDBGMAdcRegs、ADCTRL1、bit、ACQ_PS=ADC_SHCLK；AdcRegs、ADCTRL3、bit、ADCCLKPS=ADC_CKPS；AdcRegs、ADCTRL1、bit、SEQ_CASC＝1;//0Non-CascadedMode;1CascadedModeAdcRegs、ADCTRL2、bit、INT_ENA_SEQ1=0x1；AdcRegs、ADCTRL2、bit、RST_SEQ1=0x1;AdcRegs、ADCCHSELSEQ1、bit、CONV00＝0x6;AdcRegs、ADCMAXCONV、bit、MAX_CONV1=15;AdcRegs、ADCTRL2、bit、SOC_SEQ1=0x1;ﻩfor(i＝0;i＜(SampleLong/2)；i＋+){ﻩﻩﻩﻩ}ﻩﻩ{ﻩﻩﻩﻩn=(long)p*(long)p＋(long)q*(long)q;mod[m]＝sqrt(n)；ﻩm++;ﻩﻩbreak;ﻩfft512(DDataBuffer，512);ﻩp=DDataBuffer[i]、real;ﻩq=DDataBuffer[i]、imag;ﻩﻩﻩﻩn=(long)p＊(long)p+(long)q*(long)q;ﻩﻩm++;ﻩﻩcase1024：/*1024point*/ﻩﻩfft1024(DDataBuffer,1024);ﻩ{ﻩﻩﻩq＝DDataBuffer[i]、imag;n=(long)p*(long)p+(long)q*(long)q；ﻩﻩﻩmod[m]=sqrt(n);m＋+;}}}}{er{}{Ad_data[convcount]=AdcRegs、ADCRESULT0>