| 5 | 1/1 | 返回列表 |
| 查看: 2293 | 回復(fù): 18 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
xunking銀蟲 (初入文壇)
|
[求助]
FFT算法的一個問題請教
|
||
|
以下是matlab 的help fft 文檔中的實(shí)例,有一個問題不太明白,望不吝指教: 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點(diǎn)x序列補(bǔ)零后使之成為長NFFT點(diǎn)的序列,由IFFT的定義應(yīng)該是除以NFFT才對?但是實(shí)驗(yàn)情況確實(shí)除以L,不知道是不是MATLAB里補(bǔ)零后FFT的定義不一樣了。 我實(shí)驗(yàn)中驗(yàn)證了,確實(shí)應(yīng)該除以L而不應(yīng)該是NFFT。 (NFFT=2^15長度的FFT變換,fft除以長度L,幅值才大概是對的大概,除以NFFT反而不對,實(shí)驗(yàn)如下: 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)) ) |
銀蟲 (初入文壇)
|
謝謝你的回復(fù)。 我大概理解了,這確實(shí)是一個好問題,可能很多人在學(xué)習(xí)中忽略了。我覺得應(yīng)該如下理解: 從帕斯瓦爾定理,即sum(Vs(n)^2)=sum(|y(k)|^2)/NFFT,當(dāng)補(bǔ)充零點(diǎn)后每個轉(zhuǎn)換的“幅值”確實(shí)應(yīng)該下降的。因?yàn)檠a(bǔ)零后頻域更密集的采樣,如果原來對應(yīng)的幅值仍然不變的話,能量就不守恒了。所以補(bǔ)零后由y(k)計(jì)算的信號幅值不是原始信號幅值,其中需要一個比例的變換。 所以正確的理解過程應(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)) |
新蟲 (著名寫手)
鐵桿木蟲 (職業(yè)作家)
銀蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北化085600材料專碩275|有文章專利|求調(diào)劑 +11 | Micky11223 2026-03-25 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-30 | 8/400 |
|
|
[考研] 291求調(diào)劑 +8 | Y-cap 2026-03-29 | 11/550 |
|
|
[考研] 309求調(diào)劑 +15 | 誰不是少年 2026-03-29 | 15/750 |
|
|
[考研] 310求調(diào)劑 +10 | 爭取九點(diǎn)睡 2026-03-30 | 10/500 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +7 | nhj_ 2026-03-25 | 7/350 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +9 | 天天好運(yùn)來上岸?/a> 2026-03-24 | 10/500 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個好覺_16 2026-03-24 | 4/200 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|