伦疤相

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。我想这对许多频率这样做比计算周期图慢很多倍。