LOMB(T,H,OFAC,HIFAC) 计算 Lomb 归一化周期图(光谱功率作为频率的函数)的一系列 N 个数据点 H, 在时间 T 采样,这不一定是均匀间隔的。 T 和 H 必须是大小相等的向量。 该例程将计算频谱功率对于增加的频率序列(以频率的倒数单位时间阵列 T) 高达 HIFAC 乘以平均奈奎斯特频率,具有OFAC 的过采样因子(通常 >= 4)。
返回值是所考虑的频率数组 (f), 相关的频谱功率 (P) 和估计的功率显着性值(概率)。 注:返回的意义是误报零假设的概率,即数据由独立的高斯随机变量。 低概率值表明在相关的周期信号中具有高度的重要性。
尽管此实现基于 Press 中描述的实现, Teukolsky 等人。 数字配方在 C,第 13.8 节,而不是使用三角函数,这利用了 MATALB 的数组运算符来计算方程中定义的精确频谱功率第 577 页的 13.8.4。
2021-06-18 19:51:32
2KB
matlab
1