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

下載本文檔

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

文檔簡介

1、matlab的 fft 函數(shù)2010-04-26 22:16相關語法：y = fft(x)y = fft(x,n)y = fft(x,dim)y = fft(x,n,dim)定義如下：相關的一個例子：fs = 1000;%采樣頻率t = 1/fs;%采樣時間l = 1000;%總的采樣點數(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 個點title(signal corrupted with zero-mean random noise)xlabel(time (milliseconds)nfft = 2nextpow2(l); %y = fft(y,nfft)/l;%求得最接近總采樣點的2n, 這里應該是 210=1024進行 fft變換（除以總采樣點數(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 的用法，主要有兩點注意的地方：1、從公式上看， matlab 的 fft 序號是從 1 到 n，但是絕大多數(shù)教材上是從0 到n-1。2、y=fft （x）之后，這個 y 是一個復數(shù)，它的模值應該除以（ length(x)2），才能得到各個頻率信號實際幅值。（注：專業(yè)文檔是經(jīng)驗性極

溫馨提示

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

最新文檔

評論

0/150

提交評論