1. 课程导入与耦合系统认知

在工程与科学计算中，我们面对的绝大多数真实系统都不是单一变量能够描述的。一个弹簧-质量-阻尼系统，每个质量块的位移都会通过弹簧和阻尼器影响相邻质量块；一个 RLC 电路网络，各节点的电压与支路电流通过基尔霍夫定律相互约束；一个捕食者-猎物生态系统，物种数量之间通过食物链形成非线性反馈。这类由两个或两个以上相互作用的变量（或子系统）共同决定演化过程的动力系统，就是本课程要研究的多变量耦合动力系统。

所谓“耦合”，指的是子系统之间不是简单叠加，而是存在能量、力、电流、物质或信息的传递通道。数学上，一个 n 维耦合系统通常写成状态空间形式：

dx/dt = f(x, t, u),   x ∈ R^n,  u 为外部输入
y = g(x, t)

其中状态向量 x 的各个分量之间存在交叉项，例如 x₁ 的导数中出现 x₂，x₂ 的导数中出现 x₁，这种交叉依赖就是耦合的数学标志。若 f 对 x 是线性的，可写成 dx/dt = A·x + B·u，矩阵 A 的非对角元素即刻画耦合强度；若含 x₁x₂、sin(x₁) 等项，则为非线性耦合，需要数值方法求解。

三个典型实例：从物理到生态

机械振动：考虑两个质量块 m₁、m₂ 由弹簧 k₁、k₂、k₃ 串联连接，中间质量块同时受两侧弹簧作用。其运动方程为 m₁ẍ₁ = −k₁x₁ + k₂(x₂−x₁)，m₂ẍ₂ = −k₂(x₂−x₁) − k₃x₂。可以看到 x₁ 与 x₂ 通过 k₂ 相互耦合，系统的固有频率不再是单个弹簧-质量系统的 √(k/m)，而是由耦合矩阵的特征值决定，会出现“拍振”和能量在两个质量块间来回传递的现象。

电路网络：两个通过互感 M 耦合的 RLC 回路，或一个多节点电路，节点电压方程中导纳矩阵的非对角项就是耦合项。耦合系数 k = M/√(L₁L₂) 直接决定两个回路的谐振频率分裂程度，这是无线电调谐和滤波器设计的基础。

生态系统：经典的 Lotka-Volterra 捕食者-猎物模型 dx/dt = αx − βxy，dy/dt = δxy − γy，其中 xy 项就是捕食者与猎物相遇造成的双向耦合，系统会围绕平衡点形成周期振荡，而不是单调收敛。

计算任务全景与课程路线图

面对一个耦合系统，计算工作通常沿一条主线展开：先建模、再求解、后分析、终判读。建模阶段要确定状态变量、写出耦合方程、识别参数；求解阶段要根据系统刚性、非线性程度选择显式或隐式积分器；分析阶段要计算平衡点、特征值、模态、分岔与灵敏度；判读阶段则要结合能量守恒、相位关系、收敛阶等指标判断结果是否可信。

本课程后续将按这条主线逐步深入。先建立线性耦合系统的矩阵描述与特征值分析，理解“耦合如何改变系统固有频率与稳定性”；再进入非线性耦合的数值积分，掌握 RK4、隐式 Euler、刚性方程处理等方法的适用边界；随后讨论耦合强度、参数扰动对系统行为的影响，以及分岔与混沌的识别；最后通过机械、电路、生态三类案例完成从建模到结果判读的完整闭环。

学习本课程时，建议始终带着三个问题阅读每一节：这个系统的耦合通道在哪里？用什么数学结构描述它？计算结果中哪些量能反映耦合的作用？ 把这三个问题贯穿始终，就能把零散的方法串成一条可操作的计算链路，而不是记住一堆孤立的公式。