[數(shù)學(xué)]快速傅里葉變換程序設(shè)計(jì)-資料下載頁(yè)
2025-01-18 13:24本頁(yè)面
【正文】 淚始干，我用這兩句來(lái)說(shuō)我們的呂老師實(shí)在是在恰當(dāng)不過(guò)了，在這一周的課設(shè)中，呂老師每天都要到實(shí)驗(yàn)室來(lái)看看同學(xué)們的設(shè)計(jì)過(guò)程，如果哪位同學(xué)有什么不能解決的問(wèn)題，呂老師都會(huì)不厭其煩的引導(dǎo)，給出思路。他每天都很忙，很多時(shí)候他都在辦公室忙著搞科研，同學(xué)們有問(wèn)題了就會(huì)到他的辦公室去找他，每次去辦公室找老師時(shí)心里都會(huì)有些害怕，擔(dān)心會(huì)因?yàn)榇驍_他他會(huì)批評(píng)我們，但是每次找呂老師時(shí)，他的臉上總是面帶著微笑，然后給同學(xué)們認(rèn)真的講解。而這次，呂老師不但要帶著我們搞課設(shè)，還要指導(dǎo)那些參加飛思卡爾比賽的同學(xué)，呂老師這一周過(guò)的真的很累，老師，您辛苦了在這里我也同樣要感謝雷老師在最后一天對(duì)我們的答辯，因?yàn)閰卫蠋熍R時(shí)有事，故雷老師給我答辯，這本不是雷老師的事，雷老師本可對(duì)我們寬松些，這樣自己也會(huì)輕松許多，然而雷老師卻非常的認(rèn)真，總是問(wèn)的詳詳細(xì)細(xì)的，好多同學(xué)都被問(wèn)的無(wú)言以對(duì)，而且因?yàn)榇蠹叶技敝疝q，害的雷老師很晚才去吃午飯，炎熱的夏天，悶熱的教室，豆大的汗珠順著雷老師的臉頰流下來(lái)，在次，我怎么忘記對(duì)雷老師的感謝同時(shí)，也感謝和自己做相同課設(shè)題目的黃玉亮、石菊等人，我每當(dāng)越到不懂的地方都是先和大家一切探討，所以，這次課設(shè)的成功同樣離不開我的同學(xué)們，謝謝你們?cè)诖?，也感謝圖書館為我們提供便利的條件查詢資料，謝謝參考文獻(xiàn)[1] 俞一彪，孫兵. ：東南大學(xué)出版社，2005[2] 李行一. ：重慶大學(xué)出版社，2002[3] 、結(jié)構(gòu)及應(yīng)用基礎(chǔ)——TMS320F28x[M].北京：機(jī)械工業(yè)出版社，2007.[4]姜沫岐，許涵，俞鵬，段國(guó)強(qiáng) . DSP原理與應(yīng)用—從入門到提高. 北京: 機(jī)械工業(yè)出 版社， 2007[5] 俞卞章. 數(shù)字信號(hào)處理. 西安: 西北工業(yè)大學(xué)出版社， 2002[6] 附錄A1 程序清單include //數(shù)字函數(shù)頭文件define PI define N 128 //采樣次數(shù)Nvoid InitForFFT()。 //FFT初始化函數(shù)void MakeWave()。 //波形發(fā)生函數(shù)void finv(int N1,float *xr,float *xi)。 //倒序運(yùn)算函數(shù)f（N1，Xr,Xi）,對(duì)輸入序列倒序int INPUT[N],DATA[N]。float fWaveR[N],fWaveI[N],w[N]。float sin_tab[N],cos_tab[N]。 //正余弦函數(shù)int Mum。 //Mun為蝶形圖運(yùn)算的級(jí)數(shù)void FFT(float Xr[N],float Xi[N]) //時(shí)間抽取法FFT程序，要求采樣點(diǎn)數(shù)為2的整數(shù)冪次方， //X[]和X[]分別為輸入序列的實(shí)部和虛部{ int S,B。 //S為旋轉(zhuǎn)因子的冪數(shù)，B為蝶形圖運(yùn)算輸入數(shù)據(jù)的距離， //也即各級(jí)旋轉(zhuǎn)因子的個(gè)數(shù) int m,j,k。 float X,Y。 finv(N,Xr,Xi)。 //倒序運(yùn)算函數(shù)，對(duì)輸入序列倒序 for(m=1。m=Mum。m++){ B=(int)(pow(2,m1)+)。 //B=2^(m1) for(j=0。jB。j++) //每級(jí)需要進(jìn)行B種蝶形運(yùn)算 { S=j*(int)pow(2,Mumm)+。 for(k=j。k=N1。k+=(int)(pow(2,m)+)) { X=Xr[k+B]*cos_tab[S]+Xi[k+B]*sin_tab[S]。 Y=Xi[k+B]*cos_tab[S]Xr[k+B]*sin_tab[S]。 Xr[k+B]=Xr[k]X。 Xi[k+B]=Xi[k]Y。 Xr[k]=Xr[k]+X。 Xi[k]=Xi[k]+Y。 } } }for(m=0。mN/2。m++){ w[m]=sqrt(Xr[m]*Xr[m]+Xi[m]*Xi[m])。 //計(jì)算功率普}}main(){int i。InitForFFT()。 //FFT初始化函數(shù)MakeWave()。 //波形發(fā)生函數(shù)for(i=0。iN。i++){ fWaveR[i]=INPUT[i]。 fWaveI[i]=。 w[i]=。}Mum=(int)(+log(N)/log(2))。 //Mum為蝶形運(yùn)算的級(jí)數(shù)，N=2^MunFFT(fWaveR,fWaveI)。for(i=0。iN。i++)DATA[i]=w[i]。while(1)。}void InitForFFT() //FFT初始化函數(shù)，建立正余弦函數(shù)表{int i。for(i=0。iN。i++) { sin_tab[i]=sin(PI*2*i/N)。 cos_tab[i]=cos(PI*2*i/N)。 }}void MakeWave() //波形發(fā)生函數(shù){int i。for(i=0。iN。i++){ INPUT[i]=sin(PI*2*i/N*3)*1024。 //f=3Hz,正弦函數(shù)} }//倒序運(yùn)算函數(shù)finv(N1,Xr,Xi),對(duì)輸入序列倒序//N1為序列長(zhǎng)度；Xr[],Xi[]分別為輸入序列的實(shí)部和虛部//倒序原理：倒序列的加1實(shí)在最高位加1，滿2向次高位進(jìn)1，最高位變0，依次往下//從當(dāng)前倒序值可求得下一倒序值void finv(int N1,float*xr,float*xi) //倒序運(yùn)算函數(shù)f(N1,Xr,Xi),對(duì)輸入序列倒序{ int m,n,N2,k。 //m為正序列；n為到序列；k為各個(gè)權(quán)值；//N2為最高位的權(quán)值 float T。 //臨時(shí)變量T N2=N1/2。 //最高位加1相當(dāng)于十進(jìn)制加上最高位的權(quán)N1/2 n=N2。 //第一個(gè)倒序值 for(m=1。m=N12。m++) //第0個(gè)和最后一個(gè)不倒序 { if(mn) //為了避免再次調(diào)換，只需對(duì)mn的部分調(diào)換順序 { T=xr[m]。xr[m]=xr[n]。xr[n]=T。 T=xi[m]。xi[m]=xi[n]。xi[n]=T。 } k=N2。 //最高位權(quán)值 while(n=k) { n=nk。 //次高位位1，繼續(xù)上下進(jìn)位，滿2置0 k=(int)(k/2+)。 //向下權(quán)值依次比上級(jí)減半 } n=n+k。 //得到下一倒序值 }} 25
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1