python – 使用scipy.signal.spectral.lombscargle进行句点发现

内容导读

互联网集市收集整理的这篇技术教程文章主要介绍了python – 使用scipy.signal.spectral.lombscargle进行句点发现，小编现在分享给大家，供广大互联网技能从业者学习和参考。文章包含1293字，纯文字阅读大概需要2分钟。

内容图文

新的Scipy v0.11提供了一个用于光谱分析的软件包.不幸的是,文档很少,并且没有很多可用的示例.

作为一个婴儿的例子,我试图做一个正弦波的周期发现.不幸的是,它预测的时间为1而不是预期的2pi.有任何想法吗？

# imports the numerical array and scientific computing packages
import numpy as np
import scipy as sp
from scipy.signal import spectral

# generates 100 evenly spaced points between 1 and 1000
time = np.linspace(1, 1000, 100)

# computes the sine value of each of those points
mags = np.sin(time)

# scales the sine values so that the mean is 0 and the variance is 1 (the documentation specifies that this must be done)
scaled_mags = (mags-mags.mean())/mags.std()

# generates 1000 frequencies between 0.01 and 1
freqs = np.linspace(0.01, 1, 1000)

# computes the Lomb Scargle Periodogram of the time and scaled magnitudes using each frequency as a guess
periodogram = spectral.lombscargle(time, scaled_mags, freqs)

# returns the inverse of the frequence (i.e. the period) of the largest periodogram value
1/freqs[np.argmax(periodogram)]

这将返回1而不是预期的2pi~ = 1 / 0.6366.有任何想法吗？

解决方法:

请注意,spectrum.lombscargle的最后一个参数是根据docstring的角频率：

Parameters
----------
x : array_like
Sample times.
y : array_like
Measurement values.
freqs : array_like
Angular frequencies for output periodogram.

内容总结

以上是互联网集市为您收集整理的python – 使用scipy.signal.spectral.lombscargle进行句点发现全部内容，希望文章能够帮你解决python – 使用scipy.signal.spectral.lombscargle进行句点发现所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错，欢迎将互联网集市网站推荐给程序员好友。

内容备注

版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 gblab@vip.qq.com 举报，一经查实，本站将立刻删除。

内容手机端