快速傅里葉變換fft的計(jì)算機(jī)實(shí)現(xiàn)_信號(hào)與系統(tǒng)課程設(shè)計(jì)論文(參考版)
2024-09-02 16:42本頁面
【正文】 }。iNFFT。 } coefficient=plexmul(coefficient,coquotient)。 XKtemp[i].real=XKtemp[i].real+。 //蝶形運(yùn)算，詳見報(bào)告中定義 XKtemp[ip].real=XKtemp[i].。i=i+length) //控制同一蝶形結(jié)運(yùn)算，即計(jì)算系數(shù)相同蝶形結(jié) { ip=i+lenght1。j++) //控制計(jì)算不同種蝶形結(jié)，即計(jì)算系數(shù)不同的蝶形結(jié) { for(i=j。 for(j=0。 =cos(PI/lenght1)。 //同一蝶形結(jié)中參加運(yùn)算的兩點(diǎn)的距離 =1。step++) // 控制蝶形結(jié)級(jí)數(shù) { //step表示第step級(jí)蝶形 length=2(step1)。 //計(jì)算level的值，即計(jì)算蝶形總級(jí)數(shù)level=log（2）N for(step=1。(Ntemp=Ntemp/2)!=1。 //把0改為1 } //使用蝶形運(yùn)算完成FFT運(yùn)算 Ntemp=NFFT。 //把最高位變成0 k=k/2。 } k=flag2。 XKtemp[j]=XKtemp[i]。iflag1。 flag1=NFFT1。 XKtemp[i].image=0。iNFFT。 int i,j,k。 }}/**********************************************************************************//*FFT函數(shù) 輸入：時(shí)域信號(hào)取樣點(diǎn)數(shù)NFFT，時(shí)域信號(hào)取樣數(shù)組指針xnFFT，頻域信號(hào)取樣數(shù)組指針XKFFT 功能：進(jìn)行DFT，離散時(shí)域信號(hào)值存儲(chǔ)在xn[N]，離散頻域信號(hào)值存儲(chǔ)在XK[N]*//***********************************************************************************/void FFTfunction(int NFFT,double *xnFFT,struct plexnum *XKFFT){ /*變量說明：coefficient指蝶形結(jié)運(yùn)算系數(shù) coquotient指運(yùn)算系數(shù)商，當(dāng)前系數(shù)與之前一個(gè)系數(shù)的商 plextemp為輔助變量，臨時(shí)存儲(chǔ)XKtemp用 XKtemp：為不改變XKFFT指針值所設(shè)置的替代數(shù)組 flagflag2：倒位序算法的變址標(biāo)志 Ntemp：為不改變NFFT值所設(shè)置的替代數(shù) i、j、k:輔助變量*/ struct plexnum coefficient,coquotient,plextemp,XKtemp[maxnum]。 } xnDFTtemp=xnDFT。 (*XKDFTtemp).image+=(*xnDFTtemp)*sin(2*PI*k*n/NDFT)。nNDFT。 //初始化為0 (*XKDFTtemp).image=0。kNDFT。 //為了不改變?cè)羔樀闹担x了暫時(shí)指針 struct plexnum *XKDFTtemp=XKDFT。}/**********************************************************************************//*DFT函數(shù) 輸入：時(shí)域信號(hào)取樣點(diǎn)數(shù)NDFT，時(shí)域信號(hào)取樣數(shù)組指針xnDFT，頻域信號(hào)取樣數(shù)組指針XKDFT 功能：進(jìn)行DFT，離散時(shí)域信號(hào)值存儲(chǔ)在xnDFT[NDFT]，離散頻域信號(hào)值存儲(chǔ)在XKDFT[NDFT]*//***********************************************************************************/void DFTfunction(int NDFT,double *xnDFT,struct plexnum *XKDFT){ int k,n。 =**。 }}/****************************************************復(fù)數(shù)乘法函數(shù) 輸入：兩個(gè)復(fù)數(shù)結(jié)構(gòu)體，不分先后 返回值：兩個(gè)復(fù)數(shù)的乘積的復(fù)數(shù)結(jié)構(gòu)體****************************************************/struct plexnum plexmul(struct plexnum mul1,struct plexnum mul2){ struct plexnum result。 XK[i].phaseangle=atan(XK[i].image/XK[i].real)。 XK[i].image=um[i].image。iN。 //以“絕對(duì)值*角度值”形式顯示頻域信號(hào)取樣點(diǎn) }}/**********************************************************************格式轉(zhuǎn)換函數(shù) 功能：將復(fù)數(shù)結(jié)構(gòu)體數(shù)組轉(zhuǎn)化成復(fù)頻域信號(hào)取樣點(diǎn)結(jié)構(gòu)體數(shù)組 輸入：時(shí)域信號(hào)取樣點(diǎn)數(shù)N，復(fù)數(shù)結(jié)構(gòu)體指針um，頻域信號(hào)取樣點(diǎn)指針XK**********************************************************************/void ComplexnumToXKStruct(int N,struct plexnum *um,struct XKstruct* XK){ int i。i++) { printf( XK[%d] = %f + i(%f)\n,i,XK[i].real,XK[i].image)。 for(i=0。 } xn[j]=y。i=maxharmanicorder。j++) { t=T/N*j。 } for(j=0。 scanf(%lf,amp。i=maxharmanicorder。 } printf( 請(qǐng)依次輸入諧波分量系數(shù)\n)。 scanf(%d,amp。maxharmanicorder)。 printf( 請(qǐng)輸入最大諧波次數(shù)maxharmanicorder,例如上面舉例中的最大諧波次數(shù)maxharmanicorder
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1