Lomb Scargle 阶段答案

【问题标题】：Lomb Scargle phaseLomb Scargle 阶段
【发布时间】：2018-09-26 07:08:10
【问题描述】：

有什么方法可以从 Lomb Scargle 周期图中提取相位？我正在使用来自gatspy 的 LombScargle 实现。

from gatspy.periodic import LombScargleFast

model = LombScargleFast().fit(t, y)
periods, power = model.periodogram_auto()
frequencies = 1 / periods
fig, ax = plt.subplots()
ax.plot(frequencies, power)
plt.show()

权力给了我一个绝对的价值。任何方式我都可以像离散傅立叶变换一样提取每个频率的相位。

【问题讨论】：

标签： python-3.x numpy statistics fft astropy

【解决方案1】：

假设您正在寻找一个特定的频率 fo。那么对应的周期可以由 P = 1/fo 给出。我们可以定义一个函数，如下所示：

def phase_plot(t,period):
    #t is the array of timesteps
    phases = (time/period)%1.

这将为您提供该特定频率的所有阶段。

【讨论】：

【解决方案2】：

Lomb-Scargle 方法生成周期图，即每个频率的功率。与直接最小二乘拟合正弦模型相比，这是为了能够提高性能。我不知道gatspy，但astropy 确实允许您计算特定频率的最佳相位，请参阅http://docs.astropy.org/en/stable/stats/lombscargle.html#the-lomb-scargle-model。我想对许多频率执行此操作比计算周期图要慢很多倍。

-编辑-

文档大纲移至： https://docs.astropy.org/en/stable/timeseries/lombscargle.html

【讨论】：

链接已损坏（截至 2019-08-08）