快速傅里葉變換fft的計(jì)算機(jī)實(shí)現(xiàn)_信號(hào)與系統(tǒng)課程設(shè)計(jì)論文(編輯修改稿)
2024-10-04 16:42 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 ine maxnum 128 //宏定義最大數(shù)：用以初始化數(shù)組，且限制了時(shí)域信號(hào)取樣點(diǎn)數(shù)須小于128define PI //宏定義圓周率struct XKstruct //頻域信號(hào)采樣點(diǎn)結(jié)構(gòu)體{ double real。 //實(shí)部 double image。 //虛部 double absolutevalue。 //絕對(duì)值 double phaseangle。 //相位角 double radianmeasureangle。 //用圓周率表示的角度}。struct plexnum //復(fù)數(shù)結(jié)構(gòu)體{ double real。 //實(shí)部 double image。 //虛部}。struct plexnum plexmul(struct plexnum ,struct plexnum )。 //函數(shù)聲明語(yǔ)句，在定義處進(jìn)行功能介紹void ComplexnumToXKStruct(int,struct plexnum*,struct XKstruct*)。 void DFTfunction(int ,double*,struct plexnum*)。void FFTfunction(int ,double *,struct plexnum *)。void xndistribute(int N,double *xn,double *an,double *bn)。void initfunction()。void display(int ,struct XKstruct*)。void stardivision()。void timedomainsignalsample(int ,double *)。void main(){ int i,N。 //i為輔助變量，N為時(shí)域信號(hào)取樣點(diǎn)數(shù) double xn[maxnum]。 //xn[N]存放時(shí)域信號(hào)采樣值 struct plexnum XK[maxnum]。 //XK[N]存放頻域信號(hào)采樣值 struct XKstruct XK1[maxnum],XK2[maxnum]。 //XK1[N]存放DFT變換后的頻域信號(hào)采樣值，XK2[N]存放FFT變換后的頻域信號(hào)采樣值 //double w。 //基波角頻率 //double a[maxnum]。 //諧波分量系數(shù) initfunction()。 //初始化函數(shù)，在定義出進(jìn)行功能介紹 while(1) { printf(開始運(yùn)行:\n)。 printf(輸入時(shí)域信號(hào)取樣點(diǎn)數(shù)N \n N=)。 //步驟1：輸入時(shí)域信號(hào)取樣點(diǎn)數(shù)N scanf(%d,amp。N)。 printf(\n輸入時(shí)域信號(hào)N點(diǎn)離散取樣值，存于數(shù)組xn[N]\n)。 //步驟2：依次輸入時(shí)域信號(hào)N點(diǎn)離散取樣值，存于數(shù)組xn[N] printf( 選項(xiàng)1：直接輸入離散取樣值\n 選項(xiàng)2：輸入正弦諧波分量信息，讓計(jì)算機(jī)進(jìn)行取樣\n 選擇（1/2）)。 scanf(%d,amp。i)。 while((i!=1)amp。amp。(i!=2)) { printf(請(qǐng)重新選擇：（1/2）)。 scanf(%d,amp。i)。 } switch(i) { case 1: for(i=0。iN。i++) { printf( xn[%d]=,i)。 scanf(%lf,amp。(xn[i]))。 }break。 case 2: timedomainsignalsample(N,xn)。break。 } printf(\n進(jìn)行DFT變換并顯示變換結(jié)果\n )。 //步驟3：進(jìn)行DFT變換并顯示變換結(jié)果 DFTfunction(N,xn,XK)。 //DFT函數(shù)，將時(shí)域信號(hào)采樣數(shù)組xn[N]進(jìn)行DFT，在定義處進(jìn)行功能介紹 ComplexnumToXKStruct(N,XK,XK1)。 //將復(fù)數(shù)結(jié)構(gòu)體XK轉(zhuǎn)化成復(fù)頻域信號(hào)取樣點(diǎn)XK1 display(N,XK1)。 //顯示函數(shù)，在定義處進(jìn)行功能介紹 printf(\n進(jìn)行FFT變換并顯示變換結(jié)果\n )。 //步驟4：進(jìn)行FFT變換并顯示變換結(jié)果 switch(N) //根據(jù)N的值判定能否進(jìn)行FFT，僅當(dāng)N是以2為底的正指數(shù)數(shù)時(shí)才可進(jìn)行，否則直接跳出 { case 1: return。 case 0: case 2: case 4: case 8: case 16: case 32: case 64: case 128: { FFTfunction(N,xn,XK)。 //FFT函數(shù)，將時(shí)域信號(hào)采樣數(shù)組xn[N]進(jìn)行FFT，在定義處進(jìn)行功能介紹 ComplexnumToXKStruct(N,XK,XK2)。 //將復(fù)數(shù)結(jié)構(gòu)體XK轉(zhuǎn)化成復(fù)頻域信號(hào)取樣點(diǎn)XK1 display(N,XK2)。 //顯示函數(shù)，在定義處進(jìn)行功能介紹 printf()。 //格式控制語(yǔ)句 stardivision()。 printf(\n)。 stardivision()。 }break。 default: { printf( 錯(cuò)誤：無(wú)法進(jìn)行FFT變換!!\n 原因：時(shí)域信號(hào)取樣點(diǎn)數(shù)N不是以2為底的正指數(shù)，或者N大于128\n\n)。 stardivision()。 //格式控制語(yǔ)句 printf(\n)。 stardivision()。 } } }}void stardivision() //格式控制函數(shù)：輸出一整排星號(hào)作為分隔符{ printf(*******************************************************************************\n)。}void initfunction() //初始化函數(shù){ stardivision()。 //輸出一排星號(hào)，以下為學(xué)生信息 printf( 班級(jí)：1106班\n)。 printf( 學(xué)號(hào)：U201111932\n)。 printf( 姓名：曾超\n)。 stardivision()。 stardivision()。 //本程序解說(shuō)語(yǔ)句，程序總共分為4個(gè)步驟，然后為重復(fù) printf(運(yùn)行
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1