基于matlab的fft算法程序設(shè)計(jì)word格式-wenkub.com
2024-11-08 15:27 本頁(yè)面
【正文】 如有侵權(quán)，請(qǐng)告知，我看到會(huì)立刻處理。 6 設(shè)計(jì)體會(huì)與建議 設(shè)計(jì)體會(huì) 通過(guò)此次課程設(shè)計(jì)， 讓我對(duì)快速傅里葉變換即 FFT有了更深刻的了解，不再局限于書(shū)上所學(xué)的傅里葉變化的公式，它不僅僅是一種計(jì)算方法，更是在通信系統(tǒng)中起到很大的作用，通過(guò)傅里葉變換我們可以畫(huà)出幅頻圖，從中得到信號(hào)的頻率組成。 0 20 40 60 80 100020406080頻率 ( H z )幅值正弦信號(hào)幅頻譜圖 N = 1 2 80 20 40 60 80 100020406080頻率 ( H z )幅值正弦信號(hào)加噪聲后的幅頻譜圖 N = 1 2 8 圖 4 原始信號(hào)幅頻圖 圖 5 加噪聲后信號(hào)幅頻圖 通過(guò) 幅頻曲線分 析，圖 4 中可看出只含有一種頻率的信號(hào)， f=10Hz；而從圖5 中可看到除了原始信號(hào)的頻率外還有許多的位置頻率在干擾著信號(hào)。通過(guò) IFFT 轉(zhuǎn)換 加噪聲后的 信號(hào)波形 39。 ylabel(39。 plot(ti,magx)。 magx=real(xifft1)。正弦信號(hào) 加噪聲后的 幅頻譜圖 N=12839。 ylabel(39。%做頻譜圖 axis([0,100,0,80])。*fs/length(g)。 grid。)。t39。 figure(1)。)。y39。 xlabel(39。 ti=[0:length(xifft)1]/fs。)。幅值 39。 xlabel(39。%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換 figure(2)。 %進(jìn)行 FFT 變換并做頻譜圖 y=fft(x,N)。 title(39。)。 subplot(121)。 開(kāi)始 設(shè)置抽樣頻率 fs、頻率抽樣點(diǎn) N、信號(hào)頻率 f0 生成正弦信號(hào) x=sin(2*pi*f0*t) 對(duì)該信號(hào)進(jìn)行 FFT 變換 畫(huà)出該信號(hào)的幅頻譜圖 對(duì)該信號(hào)加入一個(gè)噪聲 畫(huà)出該信號(hào)的幅頻圖 再進(jìn)行 FFT 變換 用 IFFT 恢復(fù) 加噪聲后的 信號(hào) 結(jié)束 利用 IFFT 將原始正弦信號(hào)還原 t=n/fs。 （ 3） 設(shè)計(jì) 程序 利用 分析出的 取樣頻率以及頻率取樣點(diǎn)來(lái) 畫(huà)出該信號(hào)在加噪聲與不加噪聲的情況下的幅頻圖 ，并進(jìn)行比對(duì) 。如果要提高頻率分辨力，則必須增加采樣點(diǎn)數(shù)，也即采樣時(shí)間。而每個(gè)點(diǎn)的相位呢，就是在該頻率下的信號(hào)的相位。每一個(gè)點(diǎn)就對(duì)應(yīng)著一個(gè)頻率點(diǎn)。 N 個(gè)采 樣點(diǎn)，經(jīng)過(guò) FFT 之后，就可以得到 N個(gè)點(diǎn)的 FFT 結(jié)果。 3 設(shè)計(jì)過(guò)程 設(shè)計(jì)原理 （ 1） FFT 變換原理 N點(diǎn)序列的 DFT 和 IDFT 變
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1