自适应哈密顿神经网络

摘要

利用机器学习来预测混沌系统的研究的快速增长，最近重新燃起了人们对哈密顿神经网络（HNN）的兴趣，其物理约束由汉密尔顿运动方程定义，代表了一类主要的物理增强神经网络。 我们引入了一类能够对非线性物理系统进行适应性训练预测的 HNN：通过基于目标哈密顿系统的少量分岔参数值的时间序列的神经网络，HNN 可以预测其他参数值下的动态状态，其中网络尚未接触到有关这些参数值的系统的任何信息。 HNN 的架构与之前的架构不同，我们加入了一个输入参数通道，使 HNN 参数可识别。 我们证明，使用范例哈密顿系统，使用最少四个参数值的时间序列训练 HNN，赋予神经机在整个参数区间内预测目标系统状态的能力。 利用集合最大 Lyapunov 指数和对齐指数作为指标，我们表明我们的参数认知 HNN 可以成功预测混沌过渡的路线。 物理增强机器学习是研究的前沿领域，我们的适应性强的 HNN 提供了一种理解具有广泛应用的机器学习的方法。

我简介

机器学习中的一个艰巨挑战是缺乏对人工神经网络内部工作原理的理解。 随着机器学习越来越多地融入到支持现代社会运作的许多重要结构和系统中，必须对底层神经网络的内部齿轮有一个普遍的理解。 例如，前馈神经网络或多层感知器构成了现代深度学习机器的基础，在图像、视频和音频处理领域有着广泛的应用LeCun 等人 (2015)。 这样的神经机器通常由输入层、大量隐藏层和输出层组成。 从输入层开始，同一层的节点彼此不交互，但它们通过一组权重和偏差与下一层的节点连接，其值通过训练确定，其中随机梯度的范式方法血统 (SGD) Goodfellow 等人 (2016) 经常使用。 不同层的网络如何协同工作来解决特定问题仍然未知。 在另一项研究中，储层计算是一类循环神经网络Jaeger (2001)； Mass 等人 (2002);耶格和哈斯 (2004)； Manjunath 和 Jaeger (2013) 自 2017 年以来，作为非线性和混沌动力系统的无模型、完全数据驱动预测的强大范例，获得了相当大的发展势头Haynes 等人 (2015)；较大等人 (2017); Pathak 等人 (2017);卢等人 (2017); Duriez 等人 (2017);卢等人 (2018); Pathak 等人 (2018a, b);卡罗尔（2018）；中井和齐木 (2018)；罗兰和帕利茨（2018）；翁等人 (2019);格里菲斯等人 (2019);蒋和赖（2019）； Vlachas 等人 (2019);范等人 (2020);张等人(2020)。 储层计算机由输入层、单个隐藏层和输出层组成。 与多层感知器的网络结构不同，水库计算机的隐藏层网络具有复杂的拓扑结构，其中节点按照某种概率分布相互耦合。 另一个区别是，在前馈神经网络中，仅连接隐藏层和输出层神经元的权重和偏差由训练确定，而在水库中计算这些参数以及隐藏层中复杂网络的权重是预先定义的。 训练有素的储层机可以对混沌系统的状态演化进行准确的预测，其持续时间通常比使用非线性时间序列分析中的传统方法可以实现的时间长几倍。 考虑到混沌的特点：对初始条件的敏感依赖，这排除了长期预测，这是值得注意的。 然而，人们对水库计算机的内部网络动态如何表现或“管理”以准确地（在一段时间内）复制真实系统的混沌演化知之甚少。

目前，开发一个通用的可解释框架来涵盖各种类型的机器学习是不可行的。 在这方面，一个有前途的追求方向是所谓的物理增强机器学习，其中神经网络被设计来解决特定的物理问题，目的是通过利用潜在的物理原理或约束来提高学习效率。 这个想法是在大约三十年前De Wilde (1993)提出的，当时哈密尔顿力学的原理被纳入神经网络的设计中，产生了哈密尔顿神经网络（HNN） ）最近重新受到关注 Greydanus 等人 (2019)；托特等人 (2019); Bertalan 等人 (2019); Choudhary 等人 (2020a); Garg 和 Kagi (2019)。 与传统的神经网络相比，HNN 中的能量是守恒的。 已经证明，HNN 可以被训练为拥有预测目标哈密顿系统在可积和混沌状态下动态演化的能力，前提是该网络使用取自同一组参数值的数据进行训练，其中预测为 Greydanus 等人 (2019)；托特等人 (2019); Bertalan 等人 (2019); Choudhary 等人 (2020a); Garg 和 Kagi (2019)。 最近，HNN 的原理已被推广到由拉格朗日运动方程Cranmer 等人 (2020) 和一般类型的常微分描述的系统Chen 等人 (2019)方程 Sanchez-Gonzalez 等人 (2019) 或坐标变换 Dulberg 和 Cohen (2020)； Choudhary 等人 (2020b) 在机器人领域的应用 Lutter 等人 (2019)； Havoutis 和 Ramamoorthy (2010)。

在本文中，我们解决了哈密顿神经网络的适应性，这是机器学习中的一个基本问题。 更准确地说，我们考虑目标哈密顿系统可能会经历某些参数缓慢漂移或突然变化的情况。 缓慢的环境变化可能导致绝热参数漂移，而外部干扰可能导致参数突然变化。 我们询问是否可以设计 HNN，使用来自目标系统的少量参数值的数据进行训练，以对不在训练集中的参数值具有预测能力。 受到最近基于水库计算预测耗散动力系统关键转变和崩溃的工作的启发 Falahian 等人 (2015)；塞斯特尼克和阿贝尔（2019）； Kim 等人 (2020); Klos 等人 (2020); Kong等人(2021)，我们阐述了一类HNN，其输入层包含一组通道，专门用于输入神经网络感兴趣的不同参数的值。 参数通道的数量等于目标哈密顿系统中自由变化参数的数量。 最简单的情况是目标系统具有单个分叉或控制参数，因此神经网络只需要一个输入参数通道。 我们证明，通过将这样的参数通道合并到前馈类型的 HNN 中，并使用来自少量分岔参数值（例如四个）的时间序列数据进行训练，我们可以有效地使 HNN 适应参数变化。 也就是说，经过如此训练的 HNN 继承了控制目标哈密顿系统动态演化的规则。 当训练参数集中不存在的感兴趣的参数值通过参数通道输入 HNN 时，机器能够生成与目标系统在该特定参数值下的统计匹配的动态行为。 HNN 因此变得具有适应性，因为它从未接触过目标系统在此参数值下的任何信息或数据，但神经机器可以重现动态行为。 使用 Hénon-Heiles 模型作为原型目标哈密顿系统，我们证明了我们的适应性 HNN 可以成功预测任何与训练参数集中的参数值相当接近的参数值的动态行为，无论是可积的还是混沌的。 值得注意的是，通过将一组系统变化的分岔参数值输入参数通道，HNN 可以成功预测目标哈密顿系统中向混沌的转变，我们使用两种度量来表征：系综最大李雅普诺夫指数和对齐指数。 值得强调的是，现有的 HNN 文献中 Greydanus 等人 (2019)；托特等人 (2019); Bertalan 等人 (2019); Choudhary 等人 (2020a); Garg 和 Kagi (2019)，训练和预测是在目标哈密顿系统的同一组参数值下完成的，但我们的工作不仅如此，还通过增强和扩展的可预测性使 HNN 显着变得更加强大。

我们注意到，在物理学中，机器学习已被用来解决粒子物理学中的难题Guest 等人 (2018); Radovic 等人 (2018)、量子多体系统 Carleo 和 Troyer (2017)、光学系统逆向设计 Peurifoy 等人 (2018)，以及量子信息Dunjko 和 Briegel (2018)； Carleo 等人 (2019). 然而，底层神经网络的工作机制仍然很大程度上未知Iten等人(2020)。 这里研究的物理增强 HNN 与这些应用不同，因为我们专注于利用物理原理使神经网络在参数变化方面具有前所未有的预测能力。

在秒。 II，我们描述了铰接式参数识别 HNN 的架构和训练方法。 在秒。 III，我们提出了在宽参数区域中预测 Hénon-Heiles 系统动力学行为的结果，包括基于计算系综最大 Lyapunov 指数和最小对准指数的混沌过渡。 一般来说，预测精度取决于训练所需参数值与状态的“接近”程度。 在秒。 IV，我们解决了许多相关问题，例如训练参数值的选择、多参数通道以及由一维莫尔斯势定义的哈密顿系统的 HNN。 第二节提供了总结性的讨论和推测。 V。

II 参数认知哈密顿神经网络

物理增强机器学习的中心思想是“迫使”神经网络的动态演化遵循一定的物理规则或约束，例如汉密尔顿运动方程 Greydanus 等人 (2019)；托特等人 (2019); Bertalan 等人 (2019); Choudhary 等人 (2020a); Garg 和 Kagi (2019)、拉格朗日方程 Cranmer 等人 (2020)，或最小作用原理 Karkar 等人 (2020)；杨等人(2020)。 特别是，HNN 的结构使得底层神经动力系统实际上是一个哈密顿系统，在进化过程中能量守恒。 与之前的作品不同Greydanus 等人 (2019)；托特等人 (2019); Bertalan 等人 (2019); Choudhary 等人 (2020a); Garg and Kagi (2019)，目标哈密顿系统的分岔参数通过额外的输入通道作为神经网络的输入“变量”，使 HNN 学习将输入时间序列与特定值关联起来的分岔参数。 使用少量不同分叉参数值的时间序列来训练 HNN，它可以获得通过分叉参数“感知”目标系统的动态（或动态“气候”）变化的能力。

图1： 参数识别 HNN 的结构。 输入通道由蓝色圆圈表示，连接到第一个隐藏层（紫色圆圈）。 “$\alpha$”表示的蓝色圆圈是参数输入通道，输入目标哈密顿系统的分岔参数值，以及时间序列$𝐪(t)$和$𝐩(t)$通过相应的输入通道，进入第一个隐藏层。 有两个隐藏层。 输出变量是目标系统的哈密顿量相对于规范坐标和动量的偏导数，与哈密顿量一起确定下一个时间步的动态状态。

我们的铰接参数认知 HNN 的结构如图 1 所示，其中输入包含三个部分：目标系统的位置和动量变量以及分叉参数。 具体来说，我们使用两个隐藏层，其中每层包含 $200$ 个人工神经元（节点）。 第三层是输出，其中包含单个节点，其动态状态对应于目标系统的哈密顿量。 让 $𝐲$ 表示每层的动态变量集。 第$i$层动态变量到第$(i+1)$层动态变量的变换遵循以下规则：

$${𝐲}^{i+1}={𝝈}^i({𝒲}^i\cdot {𝐲}^i+{𝐛}^i),$$ (1)

其中 ${𝝈}^i$ 是给定的非线性激活函数，${𝒲}^i$ 是权重矩阵，${𝐛}^i$ 是与 $i$ 中的神经元相关的偏差向量第层，通过训练来确定。 我们将输出设置为输入变量的空间导数，以迫使神经网络的动力学遵循汉密尔顿运动方程。 导数是通过反向传播算法计算的。 一旦输出已知，损失函数定义为

$$ℒ=\Vert \frac{\partial H}{\partial 𝐪}+\frac{d{𝐩}_{\text{real}}}{dt}\Vert +\Vert \frac{\partial H}{\partial 𝐩}-\frac{d{𝐪}_{\text{real}}}{dt}\Vert ,$$ (2)

可以计算。 通过训练过程，我们优化了方程中的权重和偏差。 (1) 通过最小化损失函数。 这是通过使用标准 SGD 方法Goodfellow 等人 (2016) 来完成的。 从网络搭建、训练到进行预测的整个过程都是使用开源包Tensorflow和KerasAbadi等人（2016）完成的； Chollet等人(2015)。

表 1 总结了我们的 HNN 的结构和参数。 它有大约40,000个未知参数需要通过训练来优化确定。 即使不使用并行或 GPU 加速，计算也可以非常高效。 预期的是，在使用来自少量不同分岔参数值的时间序列数据进行训练后，HNN 可以在较宽的参数区间内预测目标哈密顿系统的动态行为，其中参数变化是通过输入参数通道实现的到 HNN。

表格1： HNN规格

Description	Values
Number of hiden layers	$2$
Neurons per layer	$200$
Optimizer	Adam
Epochs	$500$
Activation functioin	tanh

图2： Hénon-Heiles 系统中不同类型的动力学行为。 能量值为$E=1/6$。 (a-d) 在由 $q_1=0$ 定义的截面的庞加莱曲面上，分别为 $\alpha =0$、0.7、0.9 和 1.0 的可积、混合、小和很大程度上混沌动力学。

III 用于预测混沌过渡的自适应哈密顿神经网络

为了测试我们的参数认知 HNN 在预测哈密顿系统中的状态演化和动态转换方面的适应性，我们使用范例 Hénon-Heiles 模型 Hénon 和 Heiles (1964)。 它是一个二自由度系统，用于研究不同类型的哈密顿动力学，包括可积、混合和混沌行为以及它们之间的转换。 它起源于引力三体系统Hénon and Heiles (1964)，并在分子动力学Waite and Miller (1981)等领域得到应用；费特和小弗莱克 (1984)； Vendrell 和 Meyer (2011)。

III.1系统说明

Hénon-Heiles 哈密顿量是

$$H=\frac{1}{2}\left(p_1^2+p_2^2\right)+\frac{1}{2}\left(q_1^2+q_2^2\right)+\alpha \left(q_1^2{q}_2-\frac{1}{3}{q}_2^3\right),$$ (3)

其中$q_1$和$q_2$表示坐标，$p_1$和$p_2$是相应的动量，$\alpha >0$是分叉设置非线性势函数大小的参数，描述例如分子中的解离能 Waite 和 Miller (1981)；费特和小弗莱克 (1984)； Vendrell 和 Meyer (2011)。 系统动力学方程。 (3)不仅取决于$\alpha$，还取决于系统在动态演化过程中守恒的能量$E$。 势函数的最大值为$E_{\text{max}}=1/(6{\alpha }^2)$。 对于$\alpha =0$,$E_{\text{max}}$发散，因此所有轨迹都是有界的。 对于$\alpha >0$，如果粒子能量超过$E_{\text{max}}$，哈密顿系统将变得开放，散射轨迹可以逃逸到无穷大。 为了训练适应性强的 HNN，需要有界轨迹，因此我们设置 $0\le \alpha \le 1$ 和 $E\le 1/6$。 （对于高于阈值的粒子能量，可能会出现混沌散射动力学和分形几何de Moura 和 Letelier (1999)；Seoane 等人 (2006, 2007)。） 随着 $\alpha$ 的值从 0 增加到 1，可能会出现具有不同特征的动态行为，例如可积、混合和混沌。 特别是，对于$\alpha =0$，方程中的非线性项。 (3) 消失，系统变成谐振子——可积系统。 此时，整个相空间仅包含周期轨道和准周期轨道，如图2(a)所示。 随着 $\alpha$ 从零增加，系统变得非线性，相空间中会出现 Kol'mogorov-Arnol'd-Moser (KAM) 岛屿中的混沌海洋，从而产生混合动力学，如图所示图2(b)为$\alpha =0.7$和$E=1/6$。 对于$\alpha =0.9$、$\alpha =1$和$E=1/6$，相空间中的大部分轨迹是混沌的，如图1和图2所示。 2(c) 和 2(d)。

III.2适应性训练与测试

我们训练的目标是向HNN“灌输”一定的适应性能力。 为了实现这一点，我们选择分叉参数 $\alpha$ 的多个不同值。 对于每个 $\alpha$ 值，我们随机选择能量低于逃逸阈值 $E_{\text{max}}$ 的初始条件，并对哈密顿运动方程进行数值积分，以生成相空间中的粒子轨迹。 由于混合动力学，训练数据包含可积轨道和混沌轨道。 具体来说，轨迹的时间间隔为$0\le t\le 1000$，包含数百个振荡周期，我们使用采样时间步长$dt=0.1$收集训练数据。 与训练数据相关的能量在 ${10}^{-6}$ 范围内保持恒定。

一般来说，SGD 方法确定的自适应 HNN 的权重和偏差取决于训练数据集。 为了减少预测误差，可以使用 HNN 集合Kong 等人 (2021)。 具体来说，对于 $\alpha$ 的每个值，我们为训练生成 $20$ 不同的数据集，从而形成 20 个 HNN 的集合。 训练的参数设置如表 1 所示。 2.

表2： Hénon-Heiles 系统训练参数列表

Description	Values
Neural network ensembles	$20$
Energy samples	$7$
Orbit per energy	$1$
Orbit length	$1000$
Time step	$0.1$
Training parameter set	$\alpha \in \{0.2,0.4,0.6,0.8\}$

训练结束后，方程中的所有权重和偏差都得到了。 (1) 被确定。 HNN 集合中每个网络的哈密顿量及其导数可以针对任何输入进行评估，从而得出平均导数值。 为了表征不同分岔参数值的预测精度，我们使用均方根误差 (RMSE)，该误差可以根据 HNN 预测轨道与真实轨道之间的差异计算得出。 对于$\alpha =0$，运动是可积的，因此预测轨道总是接近某个真实轨道，从而导致误差极小。 在这种情况下，我们利用了 HNN 的一个特性，即它直接产生哈密顿函数，从中可以计算势函数。 因此，可以方便地使用预测势函数与真实势函数之间的相对误差来表征 HNN 性能，其定义为

$$⟨\mathrm{\Delta }V⟩\equiv \frac{\overline{|V_{\text{pred}}-V_{\text{real}}|}}{\overline{V_{\text{real}}}},$$ (4)

其中平均值是在 区域中获取的。 预测的电势分布由 $V_{\text{pred}}=H_{\text{pred}}-C$ 给出，其中 $C=\text{min}(H_{\text{pred}})$ 使得 $V_{\text{pred}}$ 的最小值为零。 请注意，方程中的平均值。 (4)是根据物理空间中的二维域中的积分计算的，需要指定域的边界。 设置边界的标准的自然选择是，但有时预测的轨道会发散。 从数字上看，有不同的方法可以克服这个困难。 例如，如果根据标准设置边界：，则几乎所有轨道都是有界的，从而可以计算误差$⟨\mathrm{\Delta }V⟩$。

我们证明 HNN 可用于重建目标系统的哈密顿量。 考虑 的相空间点并使用泰勒级数展开关于原点的哈密顿量：

$$H_{\text{pred}}(\alpha )=\sum _{i_1,i_2,i_3,i_4}{\beta }_{i_1{i}_2{i}_3{i}_4}(\alpha )q_1^{i_1}{q}_2^{i_2}{p}_1^{i_3}{p}_2^{i_4},$$ (5)

其中$\beta$为展开系数，根据$0\le \text{sum}(i_1,i_2,i_3,i_4)\le 3$求和，共包含$35$项。 与真实的哈密顿方程相比。 (3)，只有六项是非零的。

图3： (a) 预测 Hénon-Heiles 系统势函数的相对误差。 显示的是单位间隔内的误差与分叉参数 $\alpha$ 的关系。 HNN 使用 $\alpha$ 四个值的时间序列数据进行训练：$\alpha \in \{0.2,0.4,0.6,0.8\}$（四个绿色五边形）。 在包含$\alpha$这四个值的阴影区域中，相对误差小于$2\%$，这证明了HNN在预测目标哈密顿系统中参数值不在训练集。 适应性甚至延伸到阴影区域之外，但误差较大（但仍在 $8\%$ 之内）。 (b) Coefficients of the Taylor expansion for $H_{\text{pred}}$ versus the bifurcation parameter $\alpha$, where ${\beta }_{2000}$, ${\beta }_{0200}$, ${\beta }_{0020}$ and ${\beta }_{0002}$ correspond to the first four square terms in Eq. (3) whose true value is $1/2$, and ${\beta }_{2100}$ and ${\beta }_{0300}$ correspond to the two cubic terms that are proportional to $\alpha$. 展开式中的其他项没有出现在原始哈密顿量中，其中前两个最大的项是${\beta }_{3000}$和${\beta }_{1200}$，对应于其他三次势项。

我们在分叉参数的四个值处训练 HNN：$\alpha \in \{0.2,0.4,0.6,0.8\}$。 对于每个 $\alpha$ 值，我们选择七个随机初始条件，其能量低于阈值。 图3(a)显示了预测$\alpha \in (0,1)$势函数的相对误差。 $\alpha$训练中的间隔可以分为两部分：阴影区域$\alpha \in [0.2,0.8]$包含了$\alpha$中使用的四个值，以及两者上的空白区域阴影区域的一侧。 在阴影区域内，相对误差小于$2\%$，但远离阴影区域误差增大。 图3(b) 显示了预测的哈密顿量的展开系数。 与真实哈密顿量中的项相比，我们的 HNN 可以准确地预测线性项。 对于非线性项，HNN 再现了分岔参数 $\alpha$ 变化的行为，其中图 3(b) 中阴影区域的误差很小，但相对而言区外较大。

图4： 测试 HNN 对两个训练点之间参数值的适应性。 (a,b) 分别为 $\alpha =0.7$ 势函数的真实和预测等高线图，其中后者是通过在不同位置提取具有恒定动量输入并归一化的哈密顿量而获得的。 (c,d) 分别是来自初始条件 $[q_1,q_2,p_1,p_2]=[0,0,1/\sqrt{6},1/\sqrt{6}]$ 的真实轨道和预测轨道，它们是准周期的。

为了更详细地检查我们的参数识别 HNN 的适应性，我们以两个训练点 $\alpha =0.6$ 和 $\alpha =0.8$ 之间的 $\alpha =0.7$ 为例，对于其中绝大多数能量为$E=1/6$的轨道都是准周期的，如图2（b）所示。 图4(a)和4(b)分别显示了的真实势函数和预测势函数，它们本质上是无法区分的。 图4(c)和4(d)显示了从相同初始条件开始的一些具有代表性的真实轨道和预测轨道，它们在定性上彼此一致，但在细节上有所不同。 特别值得强调的是，对于预测的准周期轨道，能量可以保持在恒定值。 事实上，我们已经测试了储层计算的方法Haynes等人(2015)；较大等人 (2017); Pathak 等人 (2017);卢等人 (2017); Duriez 等人 (2017);卢等人 (2018); Pathak 等人 (2018a, b);卡罗尔（2018）；中井和齐木 (2018)；罗兰和帕利茨（2018）；翁等人 (2019);格里菲斯等人 (2019);蒋和赖（2019）； Vlachas 等人 (2019);范等人 (2020);张等人(2020)用于预测轨道，发现虽然它通常在短时间内（例如几个周期）产生更准确的轨道，但从长远来看，能量不守恒，并且预测误差变大。 总的来说，由于测试分叉参数值 $\alpha =0.7$ 夹在两个训练点之间，因此我们的参数认知 HNN 表现出很强的适应性。

对于$\alpha =1$，能量为$E=1/6$时，大部分轨道都是混沌的，相空间中KAM岛的部分变得相对微不足道，如图2(d)。 此时，真实势函数的等高线图呈三角形，如图5(a)所示。 预测的电位等值线图如图5(b)所示，与真实情况吻合得相当好。 图5(c)和5(d)显示了相同初始条件下的真实混沌轨道和预测混沌轨道。 虽然细节不同，但 HNN 正确预测了轨道是混乱的。 事实上，正如第 2 节所示。 III.3，对于表征轨道统计行为的两个量，例如最大李雅普诺夫指数和对齐指数，预测轨道得到的结果与真实轨道的结果相同。 一般来说，测试参数值越接近训练点之一，预测精度越高。

图5： 测试 HNN 对训练间隔之外的参数值的适应性。 (a,b) $\alpha =1$ 势函数的真实和预测等值线图。 (c,d) 来自相同初始条件 $[q_1,q_2,p_1,p_2]=[0,0,1/\sqrt{6},1/\sqrt{6}]$ 的真实轨道和预测轨道，它们在细节上有所不同，但都是混沌的。

III.3 哈密顿系统的适应性预测

在典型的哈密顿系统中，随着非线性参数的增加，遍历性的转变路径如下Lichtenberg and Lieberman (1992)。 在弱非线性状态下，系统是可积的，其中运动是准周期的，并且发生在由不同初始条件生成的圆环上，如图2(a)中的Hénon-Heiles系统所示。 随着非线性参数$\alpha$的增加，出现各种大小的混沌海，导致混合相空间，如图2(b)所示。 在强非线性区域，例如$\alpha =1$，大部分相空间构成混沌海，只有一小部分仍然被KAM岛占据，如图2（ c）和（d）。 在这里，我们通过证明参数识别 HNN 的适应性提供了强有力的证据，证明它可以准确预测转换场景，并且基于仅少数非线性参数值的时间序列进行训练。

与耗散系统不同，在耗散系统中，吸引子（周期性或混沌）吸引盆中的随机初始条件导致轨迹最终都位于同一吸引子中，在哈密顿系统中，不同的初始条件通常会导致不同的动态不变集。 由于哈密顿系统的这一特征，为了研究转变场景，需要从整个相空间中的初始条件进行计算，从而对所得轨道进行统计评估和表征。 我们关注两个统计量：最大李雅普诺夫指数和最小对齐指数，前者表征无限接近轨迹的指数分离率，后者衡量两个任意向量沿轨迹的相对“接近度”（例如，它们是否成为平行、反平行或两者都不是）Skokos (2001)。 对于混沌轨迹，无穷小向量分别沿着不稳定或稳定方向呈指数拉伸或收缩。 结果，一个随机向量将沿着轨迹接近不稳定方向，并且两个随机向量将快速彼此对齐。 特别是，给定两个初始向量 ${𝐮}_1^0$ 和 ${𝐮}_2^0$，在 $i$ 个时间步之后，它们变为 ${𝐮}_1^i$ 和 ${𝐮}_2^i$ 分别。 最小对齐索引定义为

$${\gamma }^i\equiv \text{min}(\Vert {𝐮}_1^i+{𝐮}_2^i\Vert ,\Vert {𝐮}_1^i-{𝐮}_2^i\Vert ).$$ (6)

当混乱出现时，${\gamma }^i$ 的值将随着时间迅速趋近于零。

图6： 测试参数认知 HNN 在预测 Hénon-Heiles 系统中向混沌转变的适应性。 (a-c) 集合最大李亚普诺夫指数 ${\lambda }_M$、集合最小对齐索引 ${\gamma }_m$ 以及阈值 ${10}^{-8}$ 和分数 $f_c$相空间中的混沌分别与非线性参数$\alpha$的关系。 大约 $\alpha \gtrsim 0.7$ 发生向混沌的转变。 橙色和蓝色分别对应于真实结果和 HNN 预测结果。 预测行为与真实行为之间存在合理的一致性，证明了 HNN 的适应性预测能力。

对于确定了权重和偏差的经过适当训练的 HNN，输出包含预测的哈密顿量，其相对于坐标和动量向量的偏导数可以直接基于神经网络的架构计算。 这些偏导数构成了底层动力系统的速度场，然后可以确定其雅克比矩阵，然后可以根据该矩阵计算机器预测的李亚普诺夫指数和对齐指数（参见附录A)。 Lyapunov指数和最小对齐指数的真实值可以直接从目标系统的原始哈密顿量(3)计算出来。

在我们的计算中，我们在单位间隔$\alpha \in [0,1]$内取100个等距的分叉参数值。 对于每个 $\alpha$ 值，我们选择 200 个随机初始条件，并为每个初始条件计算最大 Lyapunov 指数和最小对齐索引的值。 如果最大指数为正并且最小对齐索引小于${10}^{-8}$，则轨迹被视为混沌Zotos (2015)。 我们将 200 个轨迹集合中的最大李雅普诺夫指数和最小对齐索引分别表示为 ${\lambda }_M$ 和 ${\gamma }_m$，它们是 $\alpha$ 的函数。 另一个有趣的量是混沌轨迹的分数，表示为 $f_c$，它也取决于 $\alpha$。 特征量的三元组 ${\lambda }_M$、${\gamma }_m$ 和 $f_c$ 可以根据 HNN 和原始哈密顿量作为 $\alpha$ 然后可以进行比较来评估我们的参数识别 HNN 的预测适应性能力。

图 6(a-c)分别显示了粒子能量 $E=1/6$时，${\lambda }_M$、${\gamma }_m$ 和 $f_c$ 相对于 $\alpha$ 的机器预测值和真实值。 可以看出，$\alpha \gtrsim 0.7$出现了混乱，${\lambda }_M$变为正值，${\gamma }_m$减小为${10}^{-8}$,$f_c$ 从零开始增加。 在图6(a)中，的${\lambda }_M$的真实值基本上为零，但HNN预测的${\lambda }_M$为略微积极。 显着的特征是，对于 $\alpha >0.7$，两种类型的 ${\lambda }_M$ 值都开始明显增加。 事实上，${\lambda }_M$ 的真实行为与预测行为之间存在合理的一致性。 ${\gamma }_m$ 和 $f_c$ 与 $\alpha$ 的行为中存在类似的特征，如图 2 和 3 所示。 分别为6(b)和6(c)。 这些结果有力地证明了我们的参数认知 HNN 能够对哈密顿系统中的不同动态行为进行适应性预测。

IV 与哈密顿神经网络适应性相关的问题

我们通过提出以下三个问题来解决 HNN 的适应性。 首先，HNN的适应性是否可以通过增加分叉参数的训练值数量来增强？ 第二，多参数通道能否实现适应性？ 第三，适应性对于不同的目标哈密顿系统是否成立？

IV.1 训练参数值数量的影响

图7： 增加训练参数值的数量$N$对自适应预测误差的影响。 显示的是预测文本中解释的三个模拟设置的势函数时的整体误差$⟨\mathrm{\Delta }V⟩$。 将 $N$ 增加到超过 4 并不会导致错误显着减少，这表明 HNN 已经获得了在四个不同分叉参数值下训练所需的适应性。

到目前为止，我们已经使用了分叉参数的四个不同值来训练我们的参数识别 HNN。 我们现在研究是否可以通过增加训练参数值的数量来增强适应性。 这里的“增强”是指在包含不在训练集中的值的参数区间中预测哈密顿量的总体误差的减少。 为了测试这一点，我们进行了以下数值实验。 我们选择 $N\ge 3$ 训练参数值，对于每个参数值，我们使用从低于 $M$ 能量值收集的时间序列集合来训练 HNN $M$ 次。阈值（能量低于阈值的 10 个随机初始条件的 10 个时间序列）。 为了使比较有意义，我们选择 $M$ 和 $N$ 的值，使得 $NM$ 近似恒定。 特别是，对于模拟$\mathrm{\#}$1，我们设置$N=3$:$\alpha =0.25$、0.5和0.75以及$M=9$。 对于模拟 $\mathrm{\#}$2，我们选择 $N=4$:$\alpha =0.2$、0.4、0.6 和 0.8 以及 $M=7$。 对于模拟 $\mathrm{\#}$3，我们有 $N=5$:$\alpha =0.1$、0.3、0.5、0.7 和 0.9 以及 $M=5$。 对于每次模拟，我们计算预测势函数时的集成误差$⟨\mathrm{\Delta }V⟩$，如方程（1）中定义的那样。 (4) 为 $0\le \alpha \le 1$。 结果如图7所示。 可以看出，$N=3$ 的误差通常比 $N>3$ 的误差大，但是两种情况($N=4$ 和 5）的误差大约是同样，表明将$N$增加到4以上不会导致适应性预测的误差显着减少。 也就是说，通过对分叉参数的四个值的多个时间序列进行训练，HNN 已经获得了预测其他附近参数值下的系统行为所需的适应性。

IV.2 具有两个参数通道的 HNN

图8： 具有两个输入参数通道的 HNN 的预测性能。 目标哈密顿系统由不对称 Hénon-Heiles 系统给出，如方程 1 所定义。 (7)。 显示的是 $({\alpha }_1,{\alpha }_2)$ 平面中颜色编码的整体预测误差 $⟨\mathrm{\Delta }V⟩$。 训练在绿色五边形所示的16个点进行。 中心区域$({\alpha }_1,{\alpha }_2)\in [0.2,0.8]$的预测误差较小()。

我们构建了具有多个参数通道的参数识别 HNN。 为此，我们修改了 Hénon-Heiles 哈密顿方程。 (3) 到

$$H=\frac{1}{2}\left(p_1^2+p_2^2\right)+\frac{1}{2}\left(q_1^2+q_2^2\right)+{\alpha }_1{q}_1^2{q}_2-\frac{{\alpha }_2}{3}{q}_2^3,$$ (7)

其中${\alpha }_1$和${\alpha }_2$是两个独立的分叉参数，需要两个独立的参数输入通道到HNN。 有界运动的能量阈值可以用数字来评估。 我们针对 ${\alpha }_1$ 和 ${\alpha }_2$ 值的多种组合进行训练：${\alpha }_1,{\alpha }_2\in \{0.2,0.4,0.6,0.8\}$。 训练数据的生成方式如下：对于每个参数对，我们选择低于阈值的五个能量值，并且对于每个能量值，收集单个时间序列。 训练完成后，我们在参数变化的每个方向上以间隔$0.1$预测${\alpha }_1,{\alpha }_2\in [0,1]$的势函数。 图8显示了$({\alpha }_1,{\alpha }_2)$平面中颜色编码的整体预测误差$⟨\mathrm{\Delta }V⟩$。 对于一些数值相差较大的${\alpha }_1$和${\alpha }_2$的组合，阈值能量小于$1/6$。 对于这种情况，方程中的积分域。 (4)根据阈值进行相应修改。 可以看出，在参数区域$({\alpha }_1,{\alpha }_2)\in (0.2,0.8)$内，预测误差约为$5\%$，而该区域外的误差则呈增大趋势。 在两个非对角线的拐角处，由于电位分布的强烈不对称性，误差最大。 图8展示了具有两个参数通道的HNN可以被训练以适应预测。

IV.3 双原子分子系统的 HNN

图9： 针对一维莫尔斯系统训练的参数认知 HNN。 (a) 显示了 $a=1.0$、1.5 和 2.0（实线）的预测电势分布，以及相应的真实分布（虚线）。 $a=1.5$ 的预测势函数不是 $a=1.0$ 和 $a=2.0$ 的插值，证明了 HNN 的适应性预测能力。 (b) 从 $[x,p_x]=[1,1]$ 开始的相空间中的真实轨道和预测轨道。

我们考虑一个不同的目标哈密顿系统，该系统由一维莫尔斯势定义，描述双原子分子 Morse (1929) 之间的相互作用。 该系统之前曾用于混沌散射的研究Lai (1999)；赖等人(2000). 哈密​​顿量由下式给出

$$H=\frac{p_x^2}{2}+V(x)\equiv \frac{p_x^2}{2}+{[1-\exp (-a(x-x_0))]}^2-1,$$ (8)

其中势函数 $V(x)$ 在 $x=x_0$、$V(x_0)=-1$ 和 $V(x\to \mathrm{\infty })\to 0$ 处具有最小值。 以最小势能值作为能量参考点$E$，所有轨道都以为界。 我们设置$x_0=1$并选择$a$作为分叉参数。 训练数据是根据 $a$ 的四个不同值生成的：$a=0.5$、1.0、2.0 和 4.0，其中对于每个训练参数值，使用五个能量值的集合，产生 20 个独立时间序列的集合。 每个时间序列的时间跨度为$0\le t\le 100$，采样时间步长为$\mathrm{\Delta }t=0.1$。

图9显示了$a=1.5$的预测潜力概况，以及两个训练点$a=1$和$a=2$的预测潜力概况。 对于最小势点附近的$x$，结果是准确的，但当位置远离最小点时，误差较大。 一个可能的原因是，对于较大的 $x$ 值，势能变化缓慢，导致动量变化较小。 因此，时间序列的相应部分表现出较小的变化，导致 HNN 的预测误差较大。 经过训练的 HNN 显然获得了一定的适应性，因为 $a=1.5$ 的预测结果不是 $a=1.0$ 和 $a=2.0$ 预测结果的插值。

V讨论

一般来说，开发适应性强的机器学习对于当前感兴趣的关键问题具有广泛的应用。 例如，最重要的问题是根据当前可用的信息来预测当系统的一些关键参数可能发生漂移时系统将来的行为方式。 举个例子，假设一个生态系统当前处于正常状态。 由于环境恶化，其一些参数，例如承载能力和/或物种腐烂率将在未来发生漂移。 当发生一定量的参数漂移时，当系统方程未知并且唯一可用的信息是可以测量之前但包括现在的时间序列数据时，是否可以预测系统是否会崩溃t0>? 适应性机器学习提供了一种可能的解决方案。 例如，最近的研究表明，Kong 等人 (2021) 将参数识别机制纳入储层计算机中，可以针对任何给定数量的参数预测未来可能的关键转变和系统崩溃漂移。 然而，目前正在深入研究的最先进的储层计算方案Haynes 等人（2015）；较大等人 (2017); Pathak 等人 (2017);卢等人 (2017); Duriez 等人 (2017);卢等人 (2018); Pathak 等人 (2018a, b);卡罗尔（2018）；中井和齐木 (2018)；罗兰和帕利茨（2018）；翁等人 (2019);格里菲斯等人 (2019);蒋和赖（2019）； Vlachas 等人 (2019);范等人 (2020);张等人(2020)没有考虑能量守恒等物理约束，因此它们适用于耗散动力系统。

将物理定律与传统机器学习相结合有可能显着提高神经网络的性能和预测能力。 最近的研究表明，在传统前馈神经网络中执行哈密顿运动方程可以提高哈密顿系统在可积和混沌体系中的预测精度Greydanus等人(2019)；托特等人 (2019); Bertalan 等人 (2019); Choudhary 等人 (2020a); Garg 和 Kagi (2019)。 在这些研究中，训练和预测是针对目标哈密顿系统的同一组参数值进行的，因此底层哈密顿神经网络不具有适应性，因为它们无法预测系统在不同条件下的动态行为。参数设置。

我们开发的适应性强的HNN有什么实际意义吗？ 从预测受参数漂移影响的哈密顿系统的未来状态的角度来看，答案也许是否定的。 主要原因是 HNN 需要所有坐标和动量时间序列。 例如，人们可能有兴趣预测一个复杂的多体天体物理系统是否会失去稳定性并在未来变得非常混乱，其中唯一可用的信息是系统运行之前或当前的位置和动量测量。仍处于基本可积的状态。 由于该系统的物理定律是已知的，因此训练所需的数据并不比了解哈密顿量本身更轻。 尽管如此，我们的工作还是对 HNN 的工作产生了深入的了解，如下所示。

我们的参数认知、适应性强的 HNN 具有到标准多层网络的参数输入通道，其损失函数由哈密顿运动方程规定，并且能够成功预测哈密顿系统中向混沌的转变。 特别是，通过使用来自四个不同分叉（非线性）参数值的坐标和动量时间序列进行训练，机器获得了适应性，其成功预测了包含训练参数的整个参数区间中目标系统的动态行为就证明了这一点价值观。 也就是说，训练的好处在于，HNN 不仅学习了目标哈密顿系统的动态“气候”，而且还学习了“气候”如何随分叉参数变化。 因此，机器学习可以被视为神经网络自我调整其动态演化规则以合并目标系统的动态演化规则的过程。

当系统地改变分岔参数的值输入 HNN 时，它可以预测从大部分可积状态到混沌的转变，这由系综最大 Lyapunov 指数和最小对齐指数以及混沌分数作为函数来确定分岔参数。 对于单个参数通道，只要训练参数集包含至少三个或四个不同值，就可以实现自适应预测能力。 对于具有双工参数通道的 HNN，所需训练参数集的大小应至少为 4 x 4。 对于由莫尔斯势函数定义的不同哈密顿系统也完成了适应性预测。 我们期望本文设计的参数识别 HNN 的原理和训练方法适用于一般哈密顿系统。

问题之一是能量表面对分岔参数的依赖性。 随着参数的不断变化，能量面也会随之演化。 如果我们打算预测固定能量值的分岔参数的某些特定值的系统动力学，则训练数据集应包含从较大能量值收集的时间序列，以覆盖所需能量值处的相关相空间区域。

还应该指出的是，使用 HNN 来预测 Hénon-Heiles 系统中从可积动力学到混沌的转变最早是 Choudhary 等人 (2020a) 报道的，它依赖于使用能量 $E$作为非线性参数固定值的控制参数（例如，$\alpha =1$)。 在这里，我们研究了使用 $\alpha$ 作为固定能量值（例如 $E=1/6$)的分叉参数的转变。 这两条路线是等价的，因为 Hénon-Heiles 系统方程。 (3) 在配置空间中具有三重对称性。 这种等价也出现在势函数包含非线性平方项的系统中，例如经典的${\varphi }^4$或FPU模型Fermi等人（1955）； Caiani等人(1998)。 然而，对于双参数哈密顿方程。本文研究的(7)，三重对称性被破缺，破坏了改变非线性参数和能量之间的等价性。 事实上，对于诸如莫尔斯和双摆系统之类的哈密顿系统，等价性并不成立 Morse (1929)；赖等人(2000); Levien 和 Tan (1993)。 我们的适应性强的 HNN 不依赖于任何此类等价性，并且可以有效地预测任何类型的哈密顿系统中向混沌的转变。

致谢

我们要感谢万尼瓦尔·布什教职人员奖学金计划的支持，该计划由负责研究和工程的助理国防部长基础研究办公室赞助，并由海军研究办公室通过拨款号 2011 资助。 N00014-16-1-2828。

附录A计算哈密顿神经网络李亚普诺夫指数和对齐指数的算法

给定一个动力系统$d𝐱/dt=𝐟(𝐱)$，雅可比矩阵由$𝒥=\partial 𝐟/\partial 𝐱$给出。 对于哈密顿系统，动态变量为 $𝐱\equiv {[𝐪,𝐩]}^T$ 和

$$𝐟(𝐪,𝐩)=[\frac{\partial H}{\partial 𝐩},-\frac{\partial H}{\partial 𝐪}].$$ (9)

对于 HNN，原则上，哈密顿量 $H$ 由神经网络的一系列操作给出，其权重和偏差如式（1）所示。 (1) 由训练决定。 计算雅可比矩阵 $𝒥$ 的另一种有效方法是有限差分法。 特别是，对于给定的初始条件，我们生成一个由 $N$ 点组成的轨道，时间间隔为 $dt$，并在每个时间步计算 $𝒥$。 令雅可比矩阵序列表示为$𝒥(t_0),𝒥(t_1),\mathrm{\cdots },𝒥(t_N)$，并令$𝒴(t_0)$为单位矩阵$ℐ$。 如果目标哈密顿系统的相空间为$D$维（对于Hénon-Heiles系统为$D=4$)，则存在$D$李亚普诺夫指数。 令$𝝀$为$D$李亚普诺夫指数的向量：$𝝀\equiv {({\lambda }_1,{\lambda }_2,\mathrm{\dots },{\lambda }_D)}^T$并将指数的初始值设置为零：$𝝀(t_0)={(0,0,\mathrm{\dots },0)}^T$。经过 $N$ 步骤后，我们有

$$𝒴(t_N)=𝒥(t_N)\cdot 𝒴(t_{N-1}).$$ (10)

对矩阵 $𝒴(t_N)$ 进行 QR 分解，得到的矩阵分别表示为 $𝒬$ 和 $ℛ$，我们有

${\lambda }_j(t_N)$ $=$ ${\lambda }_j(t_{N-1})+\log |R_{jj}|,$ (11)

$𝒴(t_N)$ $=$ $𝒬,$ (12)

其中$R_{jj}$是矩阵$ℛ$的第$j$对角线元素。 李亚普诺夫指数由下式给出

$${\lambda }_j=\underset{N\to \mathrm{\infty }}{lim}\frac{{\lambda }_j(t_N)}{Ndt},j=1,\mathrm{\dots },D.$$ (13)

最大李亚普诺夫指数为${\lambda }_M={\text{max}}_j({\lambda }_j)$。

为了计算对齐索引，我们引入一个矩阵$ℳ$并将其设置为初始时刻的单位矩阵：$ℳ(t_0)=ℐ$。 令${𝐮}_1(t_0)={[1,0,0,0]}^T$和${𝐮}_2(t_0)={[0,1,0,0]}^T$为初始时刻的两个线性无关向量。 在 $N$ 步骤之后，我们有

(14)

将向量 ${𝐮}_{1,2}(t_N)$ 按其各自的大小归一化为单位长度，我们获得最小对齐索引为

$${\gamma }_m=\underset{N\to \mathrm{\infty }}{lim}\text{min}(\Vert {𝐮}_1(t_N)+{𝐮}_2(t_N)\Vert ,\Vert {𝐮}_1(t_N)-{𝐮}_2(t_N)\Vert ).$$ (15)