[论文解读] Wavelet analysis: a new significance test for signals dominated by intrinsic red-noise variability

Paweł Lachowicz|ArXiv.org|Jun 23, 2009

Image and Signal Denoising Methods参考文献 28被引用 4

一句话总结

本文提出了一种基于小波的显著性检验方法，适用于以红噪声变异为主导的信号，特别适用于具有 $1/f^\alpha$ 型功率谱密度的X射线光 light curves。通过将真实数据的小波系数与匹配源本征功率谱密度（PSD）的模拟信号小波系数进行比较，该方法能够在 99.73% 和 99.99% 的置信水平下可靠识别准周期振荡（QPOs），成功验证了 XTE J1550–564 和 RE J1034+396 中的 QPOs。

ABSTRACT

We develop a new statistical test for the wavelet power spectrum. We design it with purpose of testing signals which intrinsic variability displays in a Fourier domain a red-noise component described by a single, broken or doubly-broken power-law model. We formulate our methodology as straightforwardly applicable to astronomical X-ray light curves and aimed at judging the significance level for detected quasi-periodic oscillations (QPOs). Our test is based on a comparison of wavelet coefficients derived for the source signal with these obtained from the averaged wavelet decomposition of simulated signal which preserves the same broad-band model of variability as displayed by X-ray source. We perform a test for statistically significant QPO detection in XTE J1550--564 microquasar and active galaxy of RE J1034+396 confirming these results in the wavelet domain with our method. In addition, we argue on the usefulness of our new algorithm for general class of signals displaying 1/f^alpha-type variability.

研究动机与目标

• 解决在以内在红噪声变异为主导的信号中，小波功率谱缺乏可靠显著性检验的问题。
• 开发一种方法，以考虑X射线源典型的宽带 $1/f^\alpha$ 型功率谱密度。
• 在具有非平稳、类似红噪声变异的时间序列中实现准周期振荡（QPOs）的稳健检测。
• 提供一种统计上可靠的、基于模拟的显著性检验方法，同时保留原始信号的PSD特性。

提出的方法

• 该方法使用时间核（Time-Kernel, TK）方法生成合成时间序列，使其匹配观测源的 $1/f^\alpha$ 功率谱密度模型。
• 利用Morlet小波变换，计算真实X射线光 light curve 和模拟信号的小波系数。
• 通过计算多次模拟中波小波功率的分布，建立显著性阈值，形成置信包络。
• 将真实信号的小波功率谱与模拟背景进行比较，识别出超过 99.73% 和 99.99% 置信水平的峰值。
• 归一化因子 $\xi = 2\Delta t \bar{x}^{-2}$ 确保全局小波谱的单位为 (rms/mean)² Hz⁻¹。
• 该方法使用 dyadic 网格尺度 $a_m = a_0 2^{m\Delta m}$ 以实现与尺度相关的时频分辨率。

实验结果

研究问题

• RQ1 能否开发一种小波显著性检验方法，能够正确处理X射线光 light curves 中的内在红噪声变异？
• RQ2 如何利用小波分析在具有 $1/f^\alpha$ 型功率谱的信号中可靠检测准周期振荡（QPOs）？
• RQ3 与现有显著性检验相比，该新方法在识别真实QPOs与噪声波动方面有多大程度的改进？
• RQ4 该方法能否推广至其他具有 $1/f^\alpha$ 型变异的信号，而不仅限于X射线天文学？

主要发现

• 该方法以 99.99% 的置信水平成功确认了在 RE J1034+396 中存在一个频率为 $f \simeq 2.7 \times 10^{-4}$ Hz 的 QPO，与先前基于傅里叶的分析结果一致。
• 在 RE J1034+396 中，该 QPO 在整个 601 点观测段内保持相干性，其 70% 的功率在 3σ 水平上显著。
• 在 XTE J1550–564 中，该方法确认了在 37–94 秒 timescale 范围内的 QPO 检测结果，与先前发现一致。
• 真实信号的小波功率谱在 QPO 持续时间的显著部分超过了 99.73% 的置信水平，验证了其统计显著性。
• 模拟背景在 $10^{-4}$–$10^{-3}$ Hz 频带内与预期的 5.5% 均方根（rms）值高度匹配，证实了模型的保真度。
• 该方法在存在红噪声变异的情况下，仍能稳健检测到相干且持续时间较长的 QPOs，优于标准的小波显著性检验方法。