据我所见,你的代码基本上是好的,但是缺少一些细节。我认为你的问题主要是解释。因此,模拟数据现在是最好看的,下面是一个示例,其中包含我在注释中建议的模拟数据(我添加了关于重要行的注释,以及##用于更改):import numpy as np

import scipy as sy

import scipy.fftpack as syfp

import pylab as pyl

dt = 0.02071

t = dt*np.arange(100000) ## time at the same spacing of your data

u = np.sin(2*np.pi*t*.01) ## Note freq=0.01 Hz

# Do FFT analysis of array

FFT = sy.fft(u)

# Getting the related frequencies

freqs = syfp.fftfreq(len(u), dt) ## added dt, so x-axis is in meaningful units

# Create subplot windows and show plot

pyl.subplot(211)

pyl.plot(t, u)

pyl.xlabel('Time')

pyl.ylabel('Amplitude')

pyl.subplot(212)

pyl.plot(freqs, sy.log10(abs(FFT)), '.') ## it's important to have the abs here

pyl.xlim(-.05, .05) ## zoom into see what's going on at the peak

pyl.show()

如您所见,有两个峰值,在+和-输入频率(.01赫兹),如预期。

编辑:

我不明白为什么这种方法不适用于OP的数据,我也看了一眼。问题是采样时间不是均匀分布的。这是时间的柱状图(代码如下)。

因此,样本之间的时间在短时间和长时间之间大致平均分配。我在这里快速寻找了一个模式,没有什么明显的。

要进行FFT,需要均匀分布的时间样本,所以我插值得到以下结果:

这是合理的(直流偏移,主峰和小谐波)。代码如下:data = np.loadtxt("data.dat", usecols=(0,1))

t = data[:,0]

u = data[:,1]

tdiff = t[1:]-t[:-1]

plt.hist(tdiff)

new_times = np.linspace(min(t), max(t), len(t))

new_data = np.interp(new_times, t, u)

t, u = new_times, new_data

FFT = sy.fft(u)

freqs = syfp.fftfreq(len(u), dt)

# then plot as above

python计算频率_如何计算给定波的频率和时间相关推荐

  1. Python计算两日期之间排除节假日与非上班时间的工作时间

    Python计算两日期之间排除节假日与非上班时间的工作时间 前言 一.基本思路 二.代码示例 总结 文章目录 前言 一.基本思路 二.代码示例 总结 前言 工作中遇见需要写UDF计算事项办理时间的需求 ...

  2. 用python计算准确率_分词结果准确率、召回率计算-python(示例代码)

    使用python计算分词结果的准确率,召回率和F值 测试文件output.txt格式如下: 团 B B 圆 E E 是 BE BE 春 B B 节 E E 千 B B 年 E E 不 B B 变 E ...

  3. python计算方位角_实例讲解:用python 计算方位角(根据两点的坐标计算)记得收藏哦...

    今天为大家分享一篇python 计算方位角实例(根据两点的坐标计算),具有很好的参考价值,希望对大家有所帮助.一起来看看吧! 知道两点坐标,怎么计算两点方向的方位角? 答:首先计算坐标增量dx,dy( ...

  4. 如何用python计算圆周率_如何用python计算圆周率?

    如何用python计算圆周率? python计算圆周率的方法: 圆周率没有精确的计算公式,所以只能用近似的方式计算它的近似值. 我们运用蒙特卡罗方法,思路很简单,在下面图形中随机抛置大量的点, 计算落 ...

  5. 用python计算圆周率_用python计算圆周率π

    原博文 2020-03-22 15:32 − 用python计算圆周率π 一.要求 1.要起义能计算到圆周率后面越多位越好. 2.用进度条显示计算的进度,能给出越多种进度条越好. 3.要求给出圆周率π ...

  6. python计算方位角_[宜配屋]听图阁

    知道两点坐标,怎么计算两点方向的方位角? 答:首先计算坐标增量dx,dy(两个对应坐标分量相减,终点的减始点的). 若dx,dy中有一个为零时,根据另一个的正负决定方位角(0,90,180,270这四 ...

  7. 核心频率个加速频率_流处理器、核心频率、 位宽……这些显卡参数你知道吗?—— 电脑硬件科普篇(八)...

    判断一张显卡的好坏,依然要使用参数 然而我认为显卡是不需要知道参数的 因为在知道性能差距的情况下,用相同分辨率下的游戏帧数进行对比就能选购自己的显卡 但是既然是科普,也就讲一讲吧 来自http://w ...

  8. cpu频率_新版win10修复CPU频率显示不正确的bug

    在此前win10 2020年5月更新后,有用户反馈win10任务管理器无法正确显示CPU使用率和频率,比如会出现CPU频率为0.00GHz. 现在有外媒报道,微软在win10测试版19042.608中 ...

  9. 计算时间差_小王子是怎么用四步法进行时间计算并成功避开日界线的

    特别申明:本文为原创,转载请注明来源 试题引入: 法国作家圣埃克絮佩里的<小王子>当中有如下一段对话.据此回答10-11题. 小王子:"我喜欢看日落.我们去看一回日落吧--&qu ...

最新文章

  1. 英特尔cpu发布时间表_英特尔第11代桌面CPU将会支持PCIe4.0,Z490主板或可支持PCIe4.0...
  2. 图像处理之霍夫变换(直线检測算法)
  3. 2013ACM多校联合(2)
  4. C# WinForm获取程序所在路径方法
  5. bzoj3456: 城市规划
  6. oracle删除所有触发器的命令,Oracle删除当前用户的所有对象(表、视图、触发器、存储过程、函数)...
  7. SAP License:职业生涯与发展之路——写给李婷小妹妹的信
  8. Unity3D基础11:Rigidbody物理类组件
  9. 对象锁 区别 类锁
  10. php有没有能无限遍历,php无限遍历目录代码
  11. Excel如何导出自定义架构xml
  12. Python 复制文件并重命名
  13. 分析图书管理系统的5W1H
  14. pca图解读_主成分分析pca图解读,主成分分析散点图解读
  15. 通信协议:分辨率和偏移量的使用
  16. 如何用计算机解锁苹果手机软件,屏幕显示连接到iTunes,如何解锁iPhone? -
  17. 《失业七个月,面试六十家公司》的深圳体验
  18. Android向服务器发送图片(一)
  19. 基于单片机的超市收银机
  20. 初涉QQ飞车第一天-创建角色

热门文章

  1. 最大公约数(GCD)算法与最小公倍数(LCM)算法-C++实现
  2. fufu笔记之多线程
  3. 金牛能源管理信息化案例
  4. (半成品)Excel文件按要求导出至模板excel中(含配套数据)
  5. 网易秘密启动“D计划” 引爆新一轮电子邮箱大战
  6. 12.3 call_soon、call_at、call_later、call_soon_threadsafe
  7. 启用DWA R8.02 Lite模式
  8. JAVA 基本语法 Scanner类的使用 实现步骤 三步走(Scanner)
  9. 关于本地时间和网络时间校验
  10. Tableau用多线图展示多个指标