[论文解读] Irregular Hamiltonian Systems

Olivera Mišković, Jorge Zanelli|ArXiv.org|Jan 31, 2003

Control and Stability of Dynamical Systems参考文献 3被引用 7

一句话总结

本文提出了一套系统性框架，用于分析具有不规则约束的哈密顿系统——在这些系统中，由于约束线性相关，标准狄拉克程序失效。该框架将约束分类为三类：线性、多重线性和高阶约束。它提供了一种修改后的哈密顿方法，以保持与原始拉格朗日系统的动力学等价性，表明只有当结果约束为第二类时，高阶约束才能被线性约束替代；否则，原始拉格朗日系统在动力学上不一致。

ABSTRACT

Hamiltonian systems with linearly dependent constraints (irregular systems), are classified according to their behavior in the vicinity of the constraint surface. For these systems, the standard Dirac procedure is not directly applicable. However, Dirac's treatment can be slightly modified to obtain, in some cases, a Hamiltonian description completely equivalent to the Lagrangian one. A recipe to deal with the different cases is provided, along with a few pedagogical examples.

研究动机与目标

• 解决在约束线性相关（即不规则系统）时狄拉克标准哈密顿程序失效的问题。
• 根据约束在约束曲面上的行为对不规则约束进行分类，区分线性、多重线性和高阶类型。
• 提供一种系统化的修改哈密顿分析的方法，以保持与原始拉格朗日描述的动力学等价性。
• 阐明在何种情况下将高阶约束替换为线性约束可得到物理等价系统，以及在何种情况下不能。

提出的方法

• 将不规则约束分为三类：线性（秩亏恒定）、多重线性（在交点处具有简单零点）和高阶（无线性近似，例如 $\phi = C(z^s - \bar{z}^s)^k$）。
• 对于多重线性约束，将其替换为在失效点处等价的正则线性约束集合 $\varphi^i \approx 0$，以确保与原始约束曲面的一致性。
• 对于高阶约束，建议将 $\phi \approx 0$ 替换为 $\chi = \phi^{1/k} \approx 0$，但仅当 $\chi$ 为第二类约束时成立。
• 使用雅可比矩阵的秩 $\Re(J)$ 作为判断正则性的标准，当在约束曲面 $\Sigma$ 上有 $\Re(J) < R$ 时判定为失效。
• 将修改后的方法应用于示例 $L = \dot{x}\dot{z} + yz^2$，表明将 $z^2 \approx 0$ 简单替换为 $z \approx 0$ 会改变物理自由度的数量。
• 通过检查线性化约束 $\chi$ 是第一类还是第二类来判断动力学等价性：仅当 $\chi$ 为第二类时，才能保持原始动力学。

实验结果

研究问题

• RQ1 在何种条件下，高阶约束 $\phi = C(z^s - \bar{z}^s)^k \approx 0$ 可被替换为线性约束 $\chi = \phi^{1/k} \approx 0$ 而不改变物理动力学？
• RQ2 为何在哈密顿系统中将 $z^2 \approx 0$ 替换为 $z \approx 0$ 会导致与原始拉格朗日描述不同的物理自由度数量？
• RQ3 如何根据约束在约束曲面附近的局部行为对不规则约束进行系统分类？
• RQ4 在何种情况下一组不规则约束在几何上等价于一组正则约束，而在何种情况下这种等价性在动力学上失效？
• RQ5 决定从高阶约束导出的线性化约束是第一类还是第二类的因素是什么？为何这一区别对动力学至关重要？

主要发现

• 多重线性约束（如 $\phi = qp \approx 0$）可在失效点替换为等价的正则线性约束集合 $\{q \approx 0, p \approx 0\}$，从而恢复狄拉克程序的适用性。
• 对于高阶约束（如 $\phi = z^2 \approx 0$），将其替换为 $\chi = z \approx 0$ 会导致系统物理自由度为零，与原始拉格朗日描述中存在一个物理自由度相矛盾。
• 仅当 $\chi$ 为第二类约束时，替换 $\chi = \phi^{1/k} \approx 0$ 才能获得动力学等价系统；若 $\chi$ 为第一类，则动力学与原始拉格朗日系统不同。
• 雅可比矩阵的秩 $\Re(J)$ 决定正则性：当在 $\Sigma$ 上 $\Re(J) < R$ 时判定为失效，这种失效可能发生在子流形上或全局范围内。
• 当正则性在全局范围内失效时（例如在 $L = \dot{x}\dot{z} + yz^2$ 示例中 $\Re(J) = 1 < R = 2$），无法通过标准线性化方法对系统进行正则化而不改变其动力学。
• 高阶约束 $\phi \approx 0$ 仅当 $\chi$ 为第二类时才与 $\chi \approx 0$ 动力学等价；否则，原始拉格朗日系统被视为动力学不一致。