基于matlab的fft算法實(shí)現(xiàn)(論文)(參考版)
2025-06-30 18:15本頁面
【正文】 grid 。直接運(yùn)算 語音信號FFT頻譜特性39。)。ylabel(39。n39。 %添加網(wǎng)格figure(4)subplot(2,1,1)if y~=0 %判斷指數(shù)是否為0 plot(20*log10(abs(y1)))。)。ylabel(39。n39。)。title(39。 %添加網(wǎng)格subplot(2,1,2)。)。ylabel(39。n39。)。 %繪制原始語音信號的幅頻響應(yīng)圖title(39。 figure(3)subplot(2,1,1)。)。ylabel(39。歸一化頻率39。頻率特性圖39。 %繪制原始語音信號的頻率圖plot(w1/pi,20*log10(abs(x2)))。 % 蝶形運(yùn)算， 注意必須先進(jìn)行減法運(yùn)算，然后進(jìn)行加法運(yùn)算，否則要使用中間變量來傳遞y(k) y(k)=y(k)+t。 % 蝶形運(yùn)算的兩個(gè)因子對應(yīng)單元下標(biāo)的關(guān)系 t=y(kp)*WNp。 WNp=WN^p。 %蝶形運(yùn)算 for L=1:M B=2^L/2。 y=x1。 I=I+1。 K=K/2。 end K=NV2。 x1(J+1)=x1(I+1)。J=0。NM1=N1。if length(x1)N x1=[x1,zeros(1,Nlength(x1))]。 %添加網(wǎng)格M=nextpow2(x1)。)。ylabel(39。時(shí)間/t39。)。 %畫出聲音采樣后的時(shí)域波形title(39。 %計(jì)數(shù)讀取信號的點(diǎn)數(shù)t=(1:N)/fs。sound(x1,fs,bits)。[x,fs,bits]=wavread(‘’,2048)。 參考文獻(xiàn)[1] 張志勇，[s].北京航空航天大學(xué)出版社，2006，08[2] 鄭君里 ,[s].高等教育出版社, 2008, 12 [3] ：清華大學(xué)出版社出版，2001[4] 高西全 丁玉美等 數(shù)字信號處理. 北京：電子工業(yè)出版社，2009[5] 余成波，陶紅艷。MATLAB是在國內(nèi)外廣泛使用的一種數(shù)字信號處理系統(tǒng)，相信除了以上優(yōu)點(diǎn)，還有許多我還未發(fā)現(xiàn)，希望能在以后的學(xué)習(xí)中有更深入的認(rèn)識。只有遇到實(shí)際問題并根據(jù)自己對課堂上獲得的專業(yè)知識的理解來解決才能真正的提高自己的能力。另一方面我也發(fā)現(xiàn)了動手實(shí)踐的重要性。通過查閱書本和以前設(shè)計(jì)的程序我發(fā)現(xiàn)了產(chǎn)生錯(cuò)誤的原因并解決了問題完成了設(shè)計(jì)?，F(xiàn)在仔細(xì)想想，這次課程設(shè)計(jì)使得我對 MATLAB 語言的理解與應(yīng)用能力得到了較大的提升，也讓我認(rèn)識到只要深入學(xué)習(xí)，提升的空間永遠(yuǎn)是存在的。為期兩周的課設(shè)很快接近尾聲了，基于MATLAB的FFT算法實(shí)現(xiàn)設(shè)計(jì)已按計(jì)劃如期全部完成，通過這次DSP課程設(shè)計(jì)，我對課堂上所學(xué)到的理論知識的理解加深了許多， 自己動腦、動手設(shè)計(jì)的能力也得到了較大提高。在這過程中我遇到了所多的難題，通過與老師的交流和學(xué)習(xí)，讓我學(xué)會了很多在課堂上沒有理解的難點(diǎn)。在老師的幫助下我順利的完成了這個(gè)課程設(shè)計(jì)，通過這次數(shù)字信號處理課程設(shè)計(jì)，讓我了解了關(guān)于MATLAB軟件在數(shù)字信號處理方面的應(yīng)用，又一次學(xué)習(xí)了MATLAB軟件的使用和程序的設(shè)計(jì)，MATLAB的仿真使我更加深入的了解了數(shù)字處理的過程 ，對我對數(shù)字信號處理的理解加深了一步。除了指數(shù)的符號相反、并多了一個(gè)1/n的因子，離散傅里葉變換的正變換與逆變換具有相同的形式?？焖俑道锶~變換可以計(jì)算出與直接計(jì)算相同的結(jié)果，但只需要的計(jì)算復(fù)雜度。有各種快速算法。5 設(shè)計(jì)總結(jié)與心得體會快速傅氏變換 英文名是fast fourier transform，快速傅氏變換（FFT）是離散傅氏變換(DFT)的快速算法，它是根據(jù)離散傅氏變換的奇、偶、虛、實(shí)等特性，對離散傅立葉變換的算法進(jìn)行改進(jìn)獲得的。該程序嚴(yán)格按程序框圖編寫，思路清晰、容易理解，程序的運(yùn)行過程在命令窗中一目了然。注：控制第L級順序運(yùn)算控制不同種的旋轉(zhuǎn)因子控制同種旋轉(zhuǎn)因子所對應(yīng)的蝶形運(yùn)算蝶形運(yùn)算的流程圖如下圖5所示：開始開始開始J=0：B1N=2^MB=2^L/2L=1:M輸入x(n)K=J+1:2^L:NP=J*2^(ML)T=x(K)+x(K+B)* W^Np x(K+B)=x(K)x(K+B)* W^Np x(K)=T倒序結(jié)束輸出 1 2 3圖5 蝶形運(yùn)算的流程圖 基于MATLAB的FFT所編寫程序的框圖如圖5所示：語音信號采集完成語音信號時(shí)域圖完成語音信號頻率特性圖編寫fft程序，畫出語音信號頻譜圖實(shí)現(xiàn)輸入信號的倒序?qū)崿F(xiàn)一級中不同種蝶形算運(yùn)實(shí)現(xiàn)一級中相同種蝶形運(yùn)算與maTLab自帶的fft比較
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1