【札记】Poincaré-Dulac变换简介
0. 引子
Poincaré–Dulac 变换是一种用于非线性动力系统的近恒等坐标变换（near-identity transformation）。其核心目标是：
它是研究非线性振荡、模态耦合、分岔理论和 Normal Form理论的重要工具。
1. 从一般非线性动力系统出发
考虑平衡点附近的动力系统：
\[\dot{\mathbf x}=A\mathbf x+\mathbf f_2(\mathbf x)+\mathbf f_3(\mathbf x)+\cdots \tag{1}\]
其中：
- \(A\) 是线性动力学矩阵。
- \(\mathbf f_k\) 是 \(k\) 阶齐次多项式非线性项。
- 假设平衡点为 \(\mathbf x=0\)。
首先对角化线性部分（假设 \(A\) 可对角化）： \[A=\operatorname{diag}(\lambda_1,\ldots,\lambda_N). \tag{2}\] 于是每个模式满足 \[\dot x_j=\lambda_jx_j+ \sum_{|\mathbf m|\ge 2}c_{j,\mathbf m}\mathbf x^{\mathbf m}, \tag{3}\]其中 \[\mathbf x^{\mathbf m}=x_1^{m_1}\cdots x_N^{m_N}. \tag{4}\] 我们的目标是判断哪些非线性耦合项可以通过坐标变换消除。
2. Poincaré–Dulac 变换
引入一个近恒等变换： \[\boxed{\mathbf x=\mathbf y+\mathbf h_2(\mathbf y)+\mathbf h_3(\mathbf y)+\cdots} \tag{5}\] 其中 \(\mathbf h_k\) 为 \(k\) 阶齐次多项式。考虑一个特定的非线性项： \[\dot x_j=\lambda_jx_j+c\,\mathbf x^{\mathbf m}. \tag{6}\] 尝试通过如下变换消去该项 \[x_j=y_j+\alpha\mathbf y^{\mathbf m} \tag{7}\] 消去该项。保留相应阶数，可得 \[\dot y_j = \lambda_jy_j+ \left[ c+(\lambda_j-\mathbf m\cdot\boldsymbol\lambda)\alpha \right]\mathbf y^{\mathbf m}. \tag{8}\] 因此选择 \[\alpha=\frac{c}{\mathbf m\cdot\boldsymbol\lambda-\lambda_j} \tag{9}\] 便能够消去该非线性项。但这里存在一个特殊情况： \[\lambda_j=\sum_{k=1}^N m_k\lambda_k \tag{10}\] 此时分母为零，无法使用上述变换消除该项。这就是 Poincaré–Dulac 共振条件。
3. Normal Form
经过逐阶变换，动力系统被化简为： \[\dot y_j=\lambda_jy_j+ \sum_{\mathbf m\in\mathcal R_j} g_{j,\mathbf m}\mathbf y^{\mathbf m} \tag{11}\] 其中共振集合为 \[\mathcal R_j= \left\{ \mathbf m: |\mathbf m|\ge2,\quad \mathbf m\cdot\boldsymbol\lambda=\lambda_j \right\}. \tag{12}\] 这就是 Poincaré–Dulac normal form。
需要注意：这是形式幂级数层面的结论。无穷阶变换不一定收敛，在有限阶截断下，还需要控制非共振分母过小带来的问题。
4. 物理图像：共振项为什么无法消除？
考虑三个振荡模式： \[\begin{aligned} \dot a_1&=-i\omega_1a_1,\\ \dot a_2&=-i\omega_2a_2,\\ \dot a_3&=-i\omega_3a_3+g a_1a_2. \end{aligned} \tag{13}\] 我们希望消除非线性耦合 \(ga_1a_2\)，对应的变换系数为 \[\alpha=\frac{g}{-i(\omega_1+\omega_2-\omega_3)}. \tag{14}\]
- 当\(\omega_3\ne\omega_1+\omega_2\) 时，可以在形式级数意义下消除这一项。
- 当\(\omega_3=\omega_1+\omega_2\)时，变换奇异。
物理上，这是三波混频的共振条件：两个模式能够有效地向第三个模式传递能量，不能通过近恒等坐标变换将这一物理过程消除。相反，非共振耦合在适当的微扰条件下可以被吸收到坐标的重新定义之中。