MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   资料存档 (https://www.labfans.com/bbs/forumdisplay.php?f=72)
-   -   MATLAB中的FFT:0Hz频率错误 (https://www.labfans.com/bbs/showthread.php?t=23483)

poster 2019-12-10 20:41

MATLAB中的FFT:0Hz频率错误
 
我想在MATLAB中使用fft来分析保存为excell文件的一些实验数据。我的代码:

A=xlsread('Book.xls'); G=A'; x=G(2, : ); N=length(x); F=[-N/2:N/2-1]/N; X = abs(fft(x-mean(x),N)) X = fftshift(X); plot(F,X) 但是它绘制了一个带有0Hz错误分量的图形,我的真实频率约为395Hz,并且未在图形中显示。请告诉我有什么问题。

任何帮助,将不胜感激。



[B]回答:[/B]

假设我们从文件中读取信号:

G = xlsread('Book.xls'); t = G(:,1); x = G(:,2); N = length(x); 首先,我们从时间轴估计采样频率,并建立频率矢量:

Fs = 1 ./ abs( t(2)-t(1) ); F = (-N/2:N/2-1)*Fs/N; 然后计算FFT并绘制:

X = abs( fft(x-mean(x),N) ); X = fftshift(X); stem(F,X) 最终找到峰值和相应的频率:

>> [~,ind] = max(X); >> F(ind) ans = -400 您可能需要放大原点附近以更清楚地看到事物:

xlim([-1000 1000])

[url=https://stackoverflow.com/questions/2930101]更多&回答...[/url]


所有时间均为北京时间。现在的时间是 01:04

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.