快速傅里葉變換fft的計(jì)算機(jī)實(shí)現(xiàn)_信號與系統(tǒng)課程設(shè)計(jì)論文-閱讀頁
2024-09-18 16:42本頁面
【正文】 }break。 stardivision()。 stardivision()。}void initfunction() //初始化函數(shù){ stardivision()。 printf( 學(xué)號：U201111932\n)。 stardivision()。 //本程序解說語句，程序總共分為4個(gè)步驟，然后為重復(fù) printf(運(yùn)行此程序后：\ 輸入時(shí)域信號取樣點(diǎn)數(shù)N \n\n \ 依次輸入時(shí)域信號N點(diǎn)離散取樣值，存于數(shù)組xn[N]\n\n \ 進(jìn)行DFT變換并顯示變換結(jié)果\n\n \ 進(jìn)行FFT變換并顯示變換結(jié)果\n\n \ 重復(fù)之前步驟\n)。 stardivision()。 //注意事項(xiàng) stardivision()。}/*********************************************時(shí)域信號取樣函數(shù) 功能：交互輸入時(shí)域諧波信號，輸出時(shí)域取樣值 輸入：取樣點(diǎn)數(shù)N，取樣值存儲數(shù)組指針xn 輸出：取樣值**********************************************/void timedomainsignalsample(int N,double *xn){ /*變量說明： w：基波角頻率 T:基波周期 t：時(shí)間變量 maxharmanicorder：波形所含的最大諧波次數(shù)，必須為奇數(shù) 例如：對于y=a[1]sin(wt)+a[3]sin(3wt)+a[5]sin(5wt)+……+a[11]sin(11wt)，maxharmanicorder=11 a[maxharmanic]：存儲諧波分量系數(shù)的數(shù)組，僅取其奇數(shù)項(xiàng) y:t變量的因變量，暫存取樣值 */ double y=0,a[maxnum],w=0,T=0,t=0。 printf(\n 下面請初始化時(shí)域信號\n 表達(dá)式舉例：y=a[1]sin(wt)+a[3]sin(3wt)+a[5]sin(5wt)……+a[11]sin(11wt)，請依據(jù)提示輸入\n)。 scanf(%lf,amp。 T=2*PI/w。 scanf(%d,amp。 while((maxharmanicorder=0)||(maxharmanicorder%2)==0) { printf( 錯(cuò)誤：最大諧波次數(shù)應(yīng)為正奇數(shù)，請重新輸入\n maxharmanicorder=)。maxharmanicorder)。 for(i=1。i+=2) { printf( %d次諧波分量系數(shù)a[%d]=,i,i)。(a[i]))。jN。 for(i=1。i+=2) { y+=a[i]*sin((i)*w*t)。 }} /****************************************************顯示函數(shù) 輸入：時(shí)域信號取樣點(diǎn)數(shù)N，頻域信號取樣數(shù)組指針 輸出：以兩種形式在命令行顯示出頻域信號取樣點(diǎn)的值， 分別是“實(shí)部+虛部”“絕對值*角度值”****************************************************/void display(int N,struct XKstruct *XK) { int i。iN。 //以“實(shí)部+虛部”形式顯示頻域信號取樣點(diǎn) printf( XK[%d] = (%f)exp(%fPIi):\n,i,XK[i].absolutevalue,XK[i].radianmeasureangle)。 for(i=0。i++) { XK[i].real=um[i].real。 XK[i].absolutevalue=sqrt(XK[i].real*XK[i].real+XK[i].image*XK[i].image)。 XK[i].radianmeasureangle=XK[i].phaseangle/PI。 =*+*。 return result。 //輔助變量 double *xnDFTtemp=xnDFT。 for(k=0。k++) { (*XKDFTtemp).real=0。 for(n=0。n++) //變換公式參見報(bào)告DFT定義 { (*XKDFTtemp).real+=(*xnDFTtemp)*cos(2*PI*k*n/NDFT)。 xnDFTtemp++。 XKDFTtemp++。 int flag2,flag1,Ntemp,level,step,length,lenght1,ip。 for(i=0。i++) { XKtemp[i].real=xnFFT[i]。 } //變址運(yùn)算，即把自然順序變成倒位序 flag2=NFFT/2。 for(i=0,j=0。i++) { if(ij) //如果ij,即進(jìn)行變址 { plextemp=XKtemp[j]。 XKtemp[i]=plextemp。 //求j的下一個(gè)倒位序 while(k=j) //如果k=j,表示j的最高位為1 { j=jk。 //k/2，比較次高位，依次類推，逐個(gè)比較，直到某個(gè)位為0 } j=j+k。 for(level=1。level++)。step=level。 //length蝶形結(jié)距離，即第m級蝶形的蝶形結(jié)相距l(xiāng)ength點(diǎn) lenght1=length/2。 //coefficient為蝶形結(jié)運(yùn)算系數(shù)，初始值為1 =0。 //coquotient為系數(shù)商，即當(dāng)前系數(shù)與前一個(gè)系數(shù)的商 =sin(PI/lenght1)。j=lenght11。i=NFFT1。 //i，ip分別表示參加蝶形運(yùn)算的兩個(gè)節(jié)點(diǎn) plextemp=plexmul(XKtemp[ip],coefficient)。 XKtemp[ip].image=XKtemp[i].。 XKtemp[i].image=XKtemp[i].image+。 //改變系數(shù)，進(jìn)行下一個(gè)蝶形運(yùn)算 } } for(i=0。i++) //將XKtemp暫存的頻域取樣信號轉(zhuǎn)存在XKFFT數(shù)組中 { XKFFT[i]=XKtemp
高考資料相關(guān)推薦
鄂ICP備17016276號-1