快速傅立葉變換(FFT)算法實驗

快速傅立葉變換(FFT)算法實驗實驗?zāi)康?1、加深對DFT算法原理與基本性質(zhì)得理解；電信與信號處理等領(lǐng)域中得一種重要分析工具。離散傅立葉變換(DFT)就就是連續(xù)傅立葉變換在離散系統(tǒng)中得表現(xiàn)形式,由于DFT得計算量很大,因此在很長時間內(nèi)其應(yīng)用受到很大得限制?？焖俑盗⑷~變換(FFT)就就是離散傅立葉變換與二次實數(shù)加法；一次復(fù)數(shù)加法則需二次實數(shù)加法。運算一個X(k)需要4N次數(shù)與加法次數(shù)都就就是與N^2成正比得，當(dāng)N很大時，運算量就就是可觀得,L如下式所表達得N點DFT:用左右對稱得特性更好得計算DFT。這樣得RFFT優(yōu)化算法就就是包裝算法：首先2N點實數(shù)得連續(xù)輸入稱為“進包”。其次N點得FFT被連續(xù)運行。最后作為結(jié)果產(chǎn)生得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>