在我的工作中,我有一个观察到的时间序列和模拟的时间序列。我想比较光曲线并检查相似性以找出最适合的模拟曲线,哪些参数最好地模拟光曲线。
目前我使用 numpy 的互相关函数来完成。但我不确定这是否是最佳选择,因为与具有较低 CC 系数的其他模拟相比,具有最高互相关系数的光曲线并不总是看起来像最佳拟合/模拟。是否有另一种方法来衡量相似性?我读了一些关于卡方统计的东西,但我不确定它是如何工作的,以及它如何应用于我的问题。
我使用的观察数据不是均匀分箱的,所以我使用了 Scipy 的插值函数。我是否也应该平滑观察数据,还是会丢失数据的真实特征。我考虑过使用 savitzky golay 平滑。
目前我正在使用蛮力方法尝试所有可能的参数并模拟相应的光变曲线。问题是这需要大量时间来处理 20 个参数。这些参数或多或少地相互依赖。所以我不能使用最小二乘拟合方法,因为有多个可能的最小值。有没有我忽略的简单方法。还是受限蛮力适合我的最佳选择?
在下图中,您将看到一个带有模拟和观测数据的图。
感谢所有建议。
![红色:模拟曲线,蓝色:观察到的插值数据