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