[傅里叶分析与数据库]
[傅里叶分析与数据库]
提问于 2017-05-17 07:29:00
我有一份水面高程档案，我想知道它的周期和振幅。所以我想用python做傅立叶分析，但我没有到达。
我的档案是这样的
- Col 0:时间，
- Col 2:海拔
我执行了以下代码：
df=pd.read_csv(file,sep=';',header=None)
df.columns=['Time','Nothing','HS','Other1','Other2']
N=len(df)
T=df.loc[0,'Time']-df.loc[len(df)-1,'Time'] #Intervalle de temps
freq=np.linspace(0.0, 1, N) #Vecteur fréquence
HSf=scipy.fft(df['HS']) #Transformée de Fourier
plt.plot(freq,abs(HSf))问题是，通过这样做，我获得了与峰值对应的HSf的第一个元素。很明显，它将对应于向量freq (0)的第一个元素。
我知道周期必须在5-8秒左右，但我想用这种分析找到它。
你能帮我吗？
谢谢
回答 1
Stack Overflow用户
回答已采纳
发布于 2017-05-17 11:39:19
我不确定我是否正确地理解了你的问题，但如果你的第0次频率分量是最强的，但你对此不感兴趣，我看到两个选择：
(1)不看全谱，而忽略第一要素：
plt.plot(freq[1:],abs(HSf)[1:])(2)在进行FFT之前减去基线：
HSf=scipy.fft(df['HS']-np.mean(df['HS']))编辑:我也很确定你的频率轴在你的应用中是否有任何意义。您可能想看看numpy.fft.fftfreq。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/44018293
复制相关文章
相似问题
腾讯云开发者
