| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2294 | 回復(fù): 18 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
xunking銀蟲(chóng) (初入文壇)
|
[求助]
FFT算法的一個(gè)問(wèn)題請(qǐng)教
|
|
|
以下是matlab 的help fft 文檔中的實(shí)例,有一個(gè)問(wèn)題不太明白,望不吝指教: 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)|') 請(qǐng)問(wèn)這里"Y = fft(y,NFFT)/L;” 為什么是除以L,而不是除以NFFT。 L點(diǎn)x序列補(bǔ)零后使之成為長(zhǎng)NFFT點(diǎn)的序列,由IFFT的定義應(yīng)該是除以NFFT才對(duì)?但是實(shí)驗(yàn)情況確實(shí)除以L,不知道是不是MATLAB里補(bǔ)零后FFT的定義不一樣了。 我實(shí)驗(yàn)中驗(yàn)證了,確實(shí)應(yīng)該除以L而不應(yīng)該是NFFT。 (NFFT=2^15長(zhǎng)度的FFT變換,fft除以長(zhǎng)度L,幅值才大概是對(duì)的大概,除以NFFT反而不對(duì),實(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)) ) |
鐵桿木蟲(chóng) (職業(yè)作家)
新蟲(chóng) (著名寫(xiě)手)
銀蟲(chóng) (初入文壇)
|
謝謝你的回復(fù)。 我大概理解了,這確實(shí)是一個(gè)好問(wèn)題,可能很多人在學(xué)習(xí)中忽略了。我覺(jué)得應(yīng)該如下理解: 從帕斯瓦爾定理,即sum(Vs(n)^2)=sum(|y(k)|^2)/NFFT,當(dāng)補(bǔ)充零點(diǎn)后每個(gè)轉(zhuǎn)換的“幅值”確實(shí)應(yīng)該下降的。因?yàn)檠a(bǔ)零后頻域更密集的采樣,如果原來(lái)對(duì)應(yīng)的幅值仍然不變的話,能量就不守恒了。所以補(bǔ)零后由y(k)計(jì)算的信號(hào)幅值不是原始信號(hào)幅值,其中需要一個(gè)比例的變換。 所以正確的理解過(guò)程應(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 %增加一個(gè)比例變換 figure;plot((1:2^15)/2^15*1000,2*abs(y_real)) |
銀蟲(chóng) (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 安徽大學(xué)專(zhuān)碩生物與醫(yī)藥專(zhuān)業(yè)(086000)324分,英語(yǔ)已過(guò)四六級(jí),六級(jí)521,求調(diào)劑 +10 | 美味可樂(lè)雞翅 2026-03-26 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)308分調(diào)劑 +11 | 你好明天你好 2026-03-30 | 12/600 |
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +11 | 秋二十二 2026-03-30 | 11/550 |
|
|
[考研] 285求調(diào)劑 +3 | FZAC123 2026-03-30 | 3/150 |
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-30 | 11/550 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +9 | 李李不服輸 2026-03-25 | 9/450 |
|
|
[考研] 289求調(diào)劑 +3 | Acesczlo 2026-03-29 | 4/200 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
[考研] 調(diào)劑求院校招收 +7 | 鶴鯨鴿 2026-03-28 | 7/350 |
|
|
[考研] 本科211總分289,08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚(yú) 2026-03-27 | 12/600 |
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚(yú) 2026-03-30 | 5/250 |
|
|
[考研] 283求調(diào)劑(080500) +14 | A child 2026-03-27 | 14/700 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專(zhuān)業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|