基于matlab的fft算法程序設(shè)計(jì)word格式-在線瀏覽
2025-01-15 15:27本頁(yè)面
【正文】 號(hào)的幅度有什么關(guān)系呢？假設(shè)原始信號(hào)的峰值為A，那么 FFT 的結(jié)果的每個(gè)點(diǎn)（除了第一個(gè)點(diǎn)直流分量之外）的模值就是 A 的 N/2倍。而每個(gè)點(diǎn)的相位呢，就是在該頻率下的信號(hào)的相位。例如某點(diǎn) n所表示的頻率為： nF ( 1) /sn F N? ? ? （ 6） 由公式 （ 26） 可以看出， Fn所能分辨到頻率為為 Fs/N，如果采樣頻率 Fs為 1024Hz，采樣點(diǎn)數(shù)為 1024 點(diǎn)，則可以分辨到 1Hz。如果要提高頻率分辨力，則必須增加采樣點(diǎn)數(shù)，也即采樣時(shí)間。 假設(shè) FFT 之后某點(diǎn) n用復(fù)數(shù) a+b*i表示，那么這個(gè)復(fù)數(shù)的模就是 An=根號(hào)a*a+b*b，相位就是 ： tan 2( , )P a b a? （ 7） 根據(jù)以上的結(jié)果，就可以計(jì)算出 n點(diǎn)（ n≠1 ，且 n=N/2）對(duì)應(yīng)的信號(hào)的表達(dá)式為： / ( / 2) c os( 2 ) 2 / c os( 2 )n n n n n nA N pi F t P A N pi F t P? ? ? ? ? ? ? ? ? ? ? （ 8） 對(duì) 于 n=1 點(diǎn)的信號(hào)，是直流分量，幅度即為 A1/N。 （ 3） 設(shè)計(jì) 程序 利用 分析出的 取樣頻率以及頻率取樣點(diǎn)來(lái) 畫(huà)出該信號(hào)在加噪聲與不加噪聲的情況下的幅頻圖 ，并進(jìn)行比對(duì) 。%設(shè)定采樣頻率 N=128。 開(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。%設(shè)定正弦 信號(hào)頻率 %生成正弦信號(hào) x=sin(2*pi*f0*t)。 subplot(121)。%作正弦信號(hào)的時(shí)域波形 xlabel(39。)。y39。 title(39。)。 %進(jìn)行 FFT 變換并做頻譜圖 y=fft(x,N)。%求幅值 f=(0:length(y)1)39。%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換 figure(2)。 plot(f,mag)。 xlabel(39。)。幅值 39。 title(39。)。 %用 IFFT 恢復(fù)原始信號(hào) xifft=ifft(y)。 ti=[0:length(xifft)1]/fs。 subplot(121)。 xlabel(39。)。y39。 title(39。)。 %對(duì)信號(hào)加信噪比為 16 的高斯白噪聲 SNR=16。 figure(1)。 plot(t,y)。t39。 ylabel(39
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1