2025-11-25T05:31:17.416859

Dynamics of reservoir computing for crises prediction

Sisodia, Jalan

Reservoir computing has emerged as a powerful framework for time series modelling and forecasting including the prediction of discontinuous transitions. However, the mechanism behind its success is not yet fully understood. This letter elucidates the functioning of reservoir computing by examining its successful prediction of boundary and attractor merging crises. We investigate in detail how reservoirs's internal dynamics mimic the actual system, that enables it to accurately reproduce the scaling exponent near boundary crisis. We establish this across distinct systems, exemplified by the logistic and Gauss maps. The study contributes to the broader understanding of the internal dynamics that enable learning algorithms to anticipate critical transitions.

academic

Dynamics of reservoir computing for crises prediction

基本信息

论文ID: 2510.13246
标题: Dynamics of reservoir computing for crises prediction
作者: Dishant Sisodia, Sarika Jalan
分类: nlin.CD (Nonlinear Sciences - Chaotic Dynamics)
发表时间: October 16, 2025
所属机构: Complex Systems Lab, Department of Physics, Indian Institute of Technology Indore
论文链接: https://arxiv.org/abs/2510.13246

摘要

储层计算(Reservoir Computing, RC)已成为时间序列建模和预测的强大框架，包括对不连续转换的预测。然而，其成功背后的机制尚未完全理解。本文通过研究储层计算对边界危机和吸引子合并危机的成功预测，阐明了储层计算的工作原理。作者详细研究了储层内部动力学如何模拟实际系统，使其能够准确再现边界危机附近的标度指数。这一结论在不同系统中得到验证，以logistic映射和Gauss映射为例。该研究为理解学习算法预测临界转换的内部动力学机制做出了贡献。

研究背景与动机

问题定义

核心问题: 尽管储层计算在预测混沌系统的危机现象方面表现出色，但其成功预测的内在机制仍不清楚
重要性:
- 危机现象在现实系统中可能导致灾难性后果(如生态系统崩溃、金融市场崩盘)
- 理解预测机制有助于开发更可靠的预警系统
- 为机器学习在复杂系统中的应用提供理论基础

现有方法局限性

机制理解不足: 现有研究主要关注RC的预测性能，缺乏从非线性动力学角度的机制分析
理论基础薄弱: 仅有少数研究从纯数学角度分析RC，缺乏动力学系统的物理直觉
验证范围有限: 缺乏在不同类型系统中验证机制普适性的研究

研究动机

本文旨在通过分析RC内部动力学结构，揭示其成功预测危机现象的机制，特别关注：

RC如何复现原系统的不动点结构
RC如何学习并再现临界标度行为
这种机制在不同系统中的普适性

核心贡献

机制阐明: 首次从动力学系统角度深入分析了RC预测危机现象的内在机制
不动点分析: 开发了基于方向纤维(Directional Fibers)的方法来分析高维RC系统的不动点结构
标度行为复现: 证明了训练后的RC能够准确复现原系统在危机点附近的标度指数(γ = 1/2)
普适性验证: 在logistic映射和Gauss映射两个不同系统中验证了机制的普适性
理论贡献: 建立了RC学习二次非线性和不变密度的理论框架

方法详解

任务定义

输入: 混沌时间序列 u(t) 和分岔参数 ε 输出: 系统在参数空间中的动力学行为预测目标: 预测边界危机(BC)和吸引子合并危机(AMC)的发生及其标度行为

模型架构

RC基本结构

RC系统由以下组件构成：

输入权重矩阵 W^in: 将n维输入映射到m维储层空间(m >> n)
储层矩阵 A: 400×400的Erdős–Rényi随机网络邻接矩阵
参数权重矩阵 W^b: 将分岔参数映射到储层空间
输出权重矩阵 W^out: 通过岭回归优化，将储层状态映射回输出空间

动力学方程

储层状态更新规则：

r[i+1] = (1-α)r[i] + α tanh(Ar[i] + W^in u[i] + k_b W^b(ε-ε_b))

训练完成后的自主系统：

r[i+1] = (1-α)r[i] + α tanh(Ar[i] + W^in W^out r[i] + k_b W^b(ε-ε_b))

技术创新点

1. 方向纤维算法

针对高维RNN系统不动点难以求解的问题，采用方向纤维方法：

定义: 对于映射 f(r): R^N → R^N，方向纤维是满足 f(r) = αc 的状态集合
优势: 形成包含所有不动点的一维流形，在 α = 0 处定位不动点
实现: 选择随机方向c，沿对应纤维遍历状态空间寻找 f(r*) = 0 的点

2. 雅可比矩阵分析

通过引入变分 δri+1 = J δri，其中：

J = (1-α)I + α{I - tanh²(Λr* + Ω)}Λ

其中 Λ = A + W^in W^out, Ω = k_b W^b

3. 标度指数理论分析

对于一维映射 x_{n+1} = f(x_n, μ)，通过泰勒展开分析逃逸条件：

∂f/∂x δx + ∂f/∂μ ε + (1/2)∂²f/∂x²(δx)² > 0

对于二次非线性系统，得到 δx ∼ √ε，因此 ⟨τ⟩ = (μ-μc)^{-1/2}

实验设置

数据集

Logistic映射: x_{n+1} = μx_n(1-x_n)
- 训练参数: μ ∈ {3.92, 3.93, 3.94, 3.95}
- 边界危机点: μc = 4
- 吸引子合并危机: μ ≈ 3.68
Gauss映射: x_{n+1} = exp(-νx_n²) + β
- 固定参数: ν = 10
- 训练参数: β ∈ {-0.76, -0.745, -0.73, -0.715, -0.7, -0.685}
- 边界危机点: β ≈ -0.682
- 吸引子合并危机: β ≈ -0.82