Lomb Scargle阶段

Question

有什么方法可以从 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()

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

Answer 1

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

Answer 2

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

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

这将为您提供特定感兴趣频率的所有相位。

Lomb Scargle阶段

Lomb Scargle phase

statistics

numpy

fft

python-3.x

astropy