| 24小時(shí)熱門版塊排行榜 |
| 2 | 1/1 | 返回列表 |
| 查看: 900 | 回復(fù): 1 | |||
小卒wwt銀蟲 (小有名氣)
|
[求助]
matlab中fft(x)與fft(x,n)的區(qū)別 已有1人參與
|
|
matlab中的fft函數(shù)主要有兩種形式,fft(x)和fft(x,n),當(dāng)n是2的m次方時(shí),兩者的計(jì)算結(jié)果一致,但有個(gè)疑問是,當(dāng)n不是2的m次方時(shí),對(duì)信號(hào)進(jìn)行fft變化應(yīng)該采用哪種方式?我在網(wǎng)上看到是采用fft(x,n)居多,matlab幫助中采用的也是這種形式,但也發(fā)現(xiàn)過采用第一種形式的,那么到底該采用哪種形式呢?還是要具體問題具體分析?這是程序代碼 %程序1 clear;clc;fs=100; t=[0:1/fs:83]; p=1.3*sin(0.48*2*pi*t)+0.5*sin(1.8*2*pi*t)+0.9*sin(2.2*2*pi*t); figure(1); plot(t,p); grid on; title('信號(hào) p(t)'); xlabel('t'); ylabel('p'); Y=fft(p); mag=abs(Y)*2/length(t); f=(0:length(t)-1)'*fs/length(t); figure(2); plot(f,mag); grid on; title('頻譜圖 (理想值:[0.48Hz,1.3]、[1.8Hz,0.5]、[2.2Hz,0.9]) '); xlabel('f (Hz)') ylabel('幅值') axis([0,10,0,1.4]); %程序2 clear;clc;fs=100; t=[0:1/fs:83]; p=1.3*sin(0.48*2*pi*t)+0.5*sin(1.8*2*pi*t)+0.9*sin(2.2*2*pi*t); figure(1); plot(t,p); grid on; title('信號(hào) p(t)'); xlabel('t'); ylabel('p'); nfft=2^nextpow2(length(t)); Y=fft(p,nfft); mag=abs(Y)*2/length(t); f=(0:nfft-1)'*fs/nfft;%頻率序列 figure(2); plot(f,mag); grid on title('頻譜圖 (理想值:[0.48Hz,1.3]、[1.8Hz,0.5]、[2.2Hz,0.9]) '); xlabel('f (Hz)') ylabel('幅值') axis([0,10,0,1.4]); |
鐵桿木蟲 (正式寫手)
|
兩種形式的效果是一樣的。 程序2中“nfft=2^nextpow2(length(t));”改為“nfft=2^nextpow2(length(t)/2);”(由采樣定理決定),可知兩種方式下的結(jié)果是一樣的。 X為信號(hào),n為變換點(diǎn)數(shù)。 Y = fft(X) 是對(duì)信號(hào)X進(jìn)行快速傅里葉變換; Y = fft(X,n)就是對(duì)信號(hào)X的前n個(gè)點(diǎn)進(jìn)行快速傅里葉變換,如果n大于x的點(diǎn)數(shù),則直接取前n個(gè)點(diǎn),若小于n,則X先進(jìn)行補(bǔ)零擴(kuò)展為n點(diǎn)序列再求N點(diǎn)FFT。 一般情況下,n要取最接近X長(zhǎng)度的2的整數(shù)冪,這樣可以實(shí)現(xiàn)更快的FFT,提高計(jì)算效率。 |

| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9 | Ymlll 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +18 | rare12345 2026-03-18 | 18/900 |
|
|
[考研] 328求調(diào)劑,英語六級(jí)551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|