matlab的FFT函數(shù)介紹_第1頁
matlab的FFT函數(shù)介紹_第2頁
matlab的FFT函數(shù)介紹_第3頁
matlab的FFT函數(shù)介紹_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、matlab的 fft 函數(shù)2010-04-26 22:16相關(guān)語法：y = fft(x)y = fft(x,n)y = fft(x,dim)y = fft(x,n,dim)定義如下：相關(guān)的一個例子：fs = 1000;%采樣頻率t = 1/fs;%采樣時間l = 1000;%總的采樣點(diǎn)數(shù)t = (0:l-1)*t;%時間序列（時間軸）%產(chǎn)生一個幅值為0.7 頻率為50hz正弦 +另外一個信號的幅值為1 頻率為120hz的正弦信號x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);y = x + 2*randn(size(t);%混入噪聲信號plot(fs*t(1

2、:50),y(1:50)%畫出前 50 個點(diǎn)title(signal corrupted with zero-mean random noise)xlabel(time (milliseconds)nfft = 2nextpow2(l); %y = fft(y,nfft)/l;%求得最接近總采樣點(diǎn)的2n, 這里應(yīng)該是 210=1024進(jìn)行 fft變換（除以總采樣點(diǎn)數(shù)，是為了后面精確看出原始信號幅值）f = fs/2*linspace(0,1,nfft/2+1);%頻率軸 ( 只畫到fs/2即可，由于y 為實數(shù)，后面一半是對稱的)% 畫出頻率幅度圖形，可以看出 50hz 幅值大概 0.7,120

3、hz 幅值大概為 1. plot(f,2*abs(y(1:nfft/2+1)title(single-sided amplitude spectrum of y(t) xlabel(frequency (hz)ylabel(|y(f)|)ps：前段時間，只是為了自己看明白， 沒有管太多，剛才上空間發(fā)現(xiàn)幾位的留言，感覺慚愧。其實要掌握 matlab 中 fft 的用法，主要有兩點(diǎn)注意的地方：1、從公式上看， matlab 的 fft 序號是從 1 到 n，但是絕大多數(shù)教材上是從0 到n-1。2、y=fft （x）之后，這個 y 是一個復(fù)數(shù)，它的模值應(yīng)該除以（ length(x)2），才能得到各個頻率信號實際幅值。（注：專業(yè)文檔是經(jīng)驗性極

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論