[论文解读] MIARMA: An information preserving method for filling gaps in time series. Application to CoRoT light curves

J. Pascual-Granado, R. Garrido+1 more|arXiv (Cornell University)|Oct 3, 2014

Stellar, planetary, and galactic studies参考文献 21被引用 8

一句话总结

MIARMA 是一种针对时间序列的新型插补方法，利用前向-后向自回归移动平均（ARMA）建模技术，保留了恒星光 light curves 的原始频谱特性，尤其适用于 CoRoT 数据。与线性插值相比，该方法可消除功率谱中的虚假混叠，揭示此前未被检测到的微结构信号，尤其在低频脉动星（如 Be 星和 δ Scuti 星）中表现显著。

ABSTRACT

The method here presented intends to minimize the effect of the gaps in the power spectra by gap-filling preserving the original information, that is, in the case of asteroseismology, the stellar oscillation frequency content. We make use of a forward-backward predictor based on autoregressive moving average modelling (ARMA) in the time domain. The method MIARMA is particularly suitable for replacing invalid data such as those present in the light curves of the CoRoT satellite due to the pass through the South Atlantic Anomaly, and eventually for the data gathered by the NASA planet hunter Kepler. We select a sample of stars from the ultra-precise photometry collected by the asteroseismic camera on board the CoRoT satellite: the δ Scuti star HD 174966, showing periodic variations of the same order as the CoRoT observational window, the Be star HD 51193, showing longer time variations, and the solar-like HD 49933, with rapid time variations. We showed that in some cases linear interpolations are less reliable to what was believed. In particular: the power spectrum of HD 174966 is clearly aliased when this interpolation is used for filling the gaps; the light curve of HD 51193 presents a much more aliased spectrum than expected for a low frequency harmonic signal; and finally, although the linear interpolation does not affect noticeably the power spectrum of the CoRoT light curve of the solar-like star HD 49933, the ARMA interpolation showed rapid variations previously unidentified that ARMA interprets as a signal. In any case, the ARMA interpolation method provides a cleaner power spectrum, that is, less contaminated by spurious frequencies. In conclusion, MIARMA appears to be a suitable method for filling gaps in the light curves of pulsating stars observed by CoRoT since the method preserves their frequency content, which is a necessary condition for asteroseismic studies.

研究动机与目标

• 解决天文时间序列中因数据缺失而引入的虚假频率混叠问题，尤其针对 CoRoT 光曲线。
• 克服线性插值在低频脉动星中无法保留原始频谱内容的局限性。
• 开发一种能保持真实信号内容的插补方法，以支持可靠的星震学分析。
• 确保该方法适用于多种恒星类型，包括 δ Scuti 星、Be 星和类太阳振荡星。
• 提供一种稳健的现有插补技术替代方案，避免将噪声误认为信号进行重建。

提出的方法

• 采用基于自回归移动平均（ARMA）模型的前向-后向预测器，估算时间序列中的缺失值。
• 将有效数据段建模为 ARMA 过程，以捕捉潜在的时间动态特性，而无需假设解析形式。
• 通过迭代应用该方法填充缺失段，同时保持原始信号的谱特性。
• 使用赤池信息准则（AIC）评估模型拟合度，以区分噪声与真实信号。
• 通过在合成白噪声数据上测试验证该方法，确认其不会引入人工功率。
• 与 k-插补和线性插值方法进行比较，评估信号保留能力和谱保真度。

实验结果

研究问题

• RQ1 线性插值在低频脉动星的有缺失光曲线中，其对真实频谱内容的保留程度如何？
• RQ2 基于 ARMA 的插补方法能否消除 CoRoT 观测孔径函数引起的虚假混叠？
• RQ3 MIARMA 是否能揭示传统插值方法所掩盖的、此前未被检测到的光曲线微结构？
• RQ4 在合成白噪声数据中，MIARMA 与 k-插补在保持白噪声统计特性方面表现如何？
• RQ5 ARMA 插补能否可靠地区分有缺失时间序列中的本征信号与噪声？

主要发现

• 线性插值在 δ Scuti 星和 Be 星（如 HD 174966 和 HD 51193）的光曲线上引入了显著的功率谱混叠，扭曲了真实的振荡频率。
• 对于类太阳星 HD 49933，线性插值对功率谱无明显影响，但 MIARMA 揭示了此前未被识别的快速变化，可解释为真实物理信号。
• 无论恒星类型如何，MIARMA 均成功消除了所有功率谱混叠，即使在长周期变星中，线性插值失效的情况下亦然。
• 在纯随机白噪声时间序列上应用 MIARMA 时，插补值的标准差为 σ = 0.0753，接近零，表明未引入人工信号。
• 相比之下，k-插补算法在相同白噪声数据上产生 σ = 0.2730 的标准差，表明其将噪声误认为信号进行重建，违背了信号保真原则。
• 赤池信息准则（AIC）确认，由 MIARMA 插补的 HD 49933 光曲线中的微结构并非白噪声，支持其为真实物理信号的解释。