基于动力学理论的复杂系统变点检测方法技术方案

本发明专利技术涉及一种基于动力学理论的复杂系统变点检测方法，方法包括：对时间序列进行相空间重构；根据随机选取的重构吸引子对系统状态进行预测；使用滑动时间窗口对系统连续预测并构建预测准确度序列；对预测准确度序列应用贝叶斯在线变点检测算法，计算运行时长分布；根据运行时长分布推断变点位置。与现有技术相比，本发明专利技术具有较好的灵敏度、鲁棒性和对不同复杂系统的泛用性。复杂系统的泛用性。复杂系统的泛用性。

【技术实现步骤摘要】
基于动力学理论的复杂系统变点检测方法


[0001]本专利技术涉及复杂系统科学领域，尤其是涉及基于动力学理论的复杂系统变点检测方法。

技术介绍

[0002]现实中的许多复杂系统往往会呈现出随时间快速变化的动态结构,而对于这种网络的时变性研究往往富有理论和实践意义。这种结构变化可能是自发的，但在更多情况下是系统本身接收到某种外界输入而诱发的.例如，在人脑功能网络中，神经元之间的连接方式和连接强度可能会因为某种药物输入而呈现完全不同的状态；在生态系统中，某一种群突然的爆发式增长或灭亡将导致其他物种之间的联系也发生变化；在经济系统中，各类新颁布的经济政策会影响整个市场的动态等等。
[0003]在时变复杂系统中，人们首先关注的问题即为对发生变化的时间点进行定位。这就是所谓的变点检测(Change
‑
point detection,CPD)问题。从统计学意义上，变点即系统观测变量的统计分布开始发生变化的时刻，而在变点前和变点后，系统的观测都体现出平稳性。因此，传统的变点检测算法往往依赖于对系统时间序列的概率分布或特征性统计量进行估本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于动力学理论的复杂系统变点检测方法，其特征在于，包括以下步骤：(1)：从系统观测中提取时间序列{x(t)}
i≤N
，其中x＝[x1，...，x
n
]，选定感兴趣的待测序列{x
i
(t)}，初始化t＝t0；(2)：选取指标集L＝{i1，i2，...，i
k
}，对待测序列{x
i
(t)}在t时刻进行预测，得到预测映射ψ
L
(x
i1
(t)，x
i2
(t)，...，x
ik
(t))和预测结果采用相空间重构以确保预测映射的存在性，即当重构维数k满足k＞2d时，d为时间序列所对应吸引子的维数，时间延迟向量[x
i
(t)，x
i
(t
‑
τ)，...，x
i
(t
‑
(k
‑
1)τ)在拓扑上与原系统观测[x
i1
(t)，x
i2
(t)，...，x
ik
(t)]等价,即存在同胚；Ψ：其中：τ为延迟时长，取τ＝1；(3)：随机选取s个不同的指标集L，对每个指标集L重复步骤(2)，形成预测结果的分布，并计算最终预测结果预测结果的标准差和预测偏差即当随机选取不同的重构维数k的指标集L对系统状态进行预测时，所得的预测结果视为一随机变量，并用核密度估计方法对该随机变量的分布进行拟合；当拟合完成后，计算该分布的期望并将其作为最终预测结果，即分布的期望并将其作为最终预测结果，即此外得到预测结果的标准差和相对于真实值的预测偏差(4)：令t＝t+1，返回步骤(2)，最终得到预测准确度序列{Std(t)}和{Err(t)}；(5)：使用BOCD变点检测算法，对序列{Std(t)}或{Err(t)}进行检测，得到运行长度矩阵R
(M+1)
×
(M+1)
；(6)：根据运行长度矩阵得到变点准确位置。2.根据权利要求1所述的一种基于动力学理论的复杂系统变点检测方法，其特征在于，步骤(2)中还包括用低维方法对Ψ和进行估计，所述低维方法采用利用高斯过程回归方法；首先需要根据训练集数据对Ψ进行训练，训练完成后即得到：其中：ψ
L
是同胚Ψ中对应于x
i
(t+1)的分量。3.根据权利要求2所述的一种基于动力学理论的复杂系统变点检测方法，其特征在于，所选取的训练集长度即时间窗口长度根据待...

【专利技术属性】
技术研发人员：侯嘉文，林伟，冷思阳，孙毅，杨伟，
申请(专利权)人：复旦大学，
类型：发明
国别省市：