| 5 | 1/1 | 返回列表 |
| 查看: 2292 | 回復(fù): 18 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
xunking銀蟲 (初入文壇)
|
[求助]
FFT算法的一個問題請教
|
|
|
以下是matlab 的help fft 文檔中的實例,有一個問題不太明白,望不吝指教: Fs = 1000; % Sampling frequency T = 1/Fs; % Sample time L = 1000; % Length of signal t = (0:L-1)*T; % Time vector % Sum of a 50 Hz sinusoid and a 120 Hz sinusoid x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t); y = x + 2*randn(size(t)); % Sinusoids plus noise plot(Fs*t(1:50),y(1:50)) title('Signal Corrupted with Zero-Mean Random Noise') xlabel('time (milliseconds)') NFFT = 2^nextpow2(L); % Next power of 2 from length of y Y = fft(y,NFFT)/L; f = Fs/2*linspace(0,1,NFFT/2); % Plot single-sided amplitude spectrum. plot(f,2*abs(Y(1:NFFT/2))) title('Single-Sided Amplitude Spectrum of y(t)') xlabel('Frequency (Hz)') ylabel('|Y(f)|') 請問這里"Y = fft(y,NFFT)/L;” 為什么是除以L,而不是除以NFFT。 L點x序列補零后使之成為長NFFT點的序列,由IFFT的定義應(yīng)該是除以NFFT才對?但是實驗情況確實除以L,不知道是不是MATLAB里補零后FFT的定義不一樣了。 我實驗中驗證了,確實應(yīng)該除以L而不應(yīng)該是NFFT。 (NFFT=2^15長度的FFT變換,fft除以長度L,幅值才大概是對的大概,除以NFFT反而不對,實驗如下: n=1:1024; LL=1024; Vs=sin(2*pi*100*n/1000); y=fft(Vs,2^15)/LL;%y=fft(Vs,2^15)/2^15 figure;plot((1:2^15)/2^15*1000,2*abs(y)) ) |
新蟲 (著名寫手)
新蟲 (著名寫手)
鐵桿木蟲 (職業(yè)作家)
銀蟲 (初入文壇)
|
謝謝你的回復(fù)。 我大概理解了,這確實是一個好問題,可能很多人在學(xué)習(xí)中忽略了。我覺得應(yīng)該如下理解: 從帕斯瓦爾定理,即sum(Vs(n)^2)=sum(|y(k)|^2)/NFFT,當(dāng)補充零點后每個轉(zhuǎn)換的“幅值”確實應(yīng)該下降的。因為補零后頻域更密集的采樣,如果原來對應(yīng)的幅值仍然不變的話,能量就不守恒了。所以補零后由y(k)計算的信號幅值不是原始信號幅值,其中需要一個比例的變換。 所以正確的理解過程應(yīng)該如下: n=1:1024; LL=1024; Vs=sin(2*pi*100*n/1000); y=fft(Vs,2^15)/2^15 y_real=y*2^15/LL %增加一個比例變換 figure;plot((1:2^15)/2^15*1000,2*abs(y_real)) |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料化工調(diào)劑 總分330 +12 | zhubinhao 2026-03-27 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +4 | 26要上岸 2026-03-28 | 4/200 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 【雙一流院校新能源、環(huán)境材料,材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 8/400 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 10/500 |
|
|
[考研] 打過很多競賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|