参数设置
阻尼系数0.20
磁铁强度1.0
摆长1.0
磁铁高度0.4
分辨率60
磁铁1（红）
磁铁2（绿）
磁铁3（蓝）
当前位置 X
—
当前位置 Y
—
速度 |v|
—
最终磁铁
—
点击画布从任意初始位置启动磁摆。吸引盆通过RK4积分计算每个网格点的收敛磁铁后绘制。
将每个网格点按磁摆最终落入的磁铁颜色着色，绘制出美丽的分形吸引盆。体验初始条件的微小差异如何导致完全不同的结果。
点击画布从任意初始位置启动磁摆。吸引盆通过RK4积分计算每个网格点的收敛磁铁后绘制。
磁摆的运动由牛顿第二定律和阻尼项控制，其核心是计算摆锤在三个磁铁引力场下的加速度。每个磁铁对摆锤的引力遵循类似磁偶极子的简化规律，与距离的$n$次方成反比（$n$是可调参数）。
$$ m\ddot{\mathbf{r}}= -c\dot{\mathbf{r}}+ \sum_{i=1}^{3}\frac{\mathbf{r}_i - \mathbf{r}}{|\mathbf{r}_i - \mathbf{r}|^{n+2}}$$其中，$m$是摆锤质量（常归为1），$\mathbf{r}$是摆锤位置向量，$\dot{\mathbf{r}}$和$\ddot{\mathbf{r}}$分别是速度和加速度。$c$是阻尼系数（对应模拟器的“阻尼”滑块），$\mathbf{r}_i$是第$i$个磁铁的位置。$n$是力场指数（对应“力场指数”滑块），决定了引力随距离衰减的快慢。
为了数值求解这个微分方程，模拟器使用了四阶龙格-库塔法（RK4）。这是一种高精度的数值积分方法，通过计算每个时间步内四个斜率的加权平均来更新状态。
$$ \begin{aligned}\mathbf{k}_1 &= f(t_n, \mathbf{y}_n) \\ \mathbf{k}_2 &= f(t_n + \frac{h}{2}, \mathbf{y}_n + h\frac{\mathbf{k}_1}{2}) \\ \mathbf{k}_3 &= f(t_n + \frac{h}{2}, \mathbf{y}_n + h\frac{\mathbf{k}_2}{2}) \\ \mathbf{k}_4 &= f(t_n + h, \mathbf{y}_n + h\mathbf{k}_3) \\ \mathbf{y}_{n+1}&= \mathbf{y}_n + \frac{h}{6}(\mathbf{k}_1 + 2\mathbf{k}_2 + 2\mathbf{k}_3 + \mathbf{k}_4) \end{aligned}$$这里，$\mathbf{y}$是状态向量（包含位置和速度），$f$是计算导数的函数（即上面的运动方程），$h$是时间步长（对应“时间步长”滑块）。RK4是CAE中求解动力学问题的基石之一。
CAE与工程仿真：在汽车碰撞安全仿真中，微小的制造公差（如焊点位置的毫米级差异）可能通过非线性动力学放大，导致仿真预测的变形模式截然不同。理解混沌特性有助于工程师评估设计鲁棒性。
天体轨道与航天：太阳系内小天体（如小行星）的轨道长期演化是典型的N体问题，具有混沌特性。这使得精确预测数百万年后的轨道极其困难，对航天器轨道设计和深空探测任务有重要影响。
加密与信息安全：混沌系统对初始条件的极端敏感性被用于生成伪随机序列和设计加密算法。即使加密密钥只有微小的不同，产生的密文也会完全无法解密，保障了信息的安全。
流体力学与湍流研究：湍流的本质就是混沌。虽然纳维-斯托克斯方程是确定的，但流体中微小的扰动会指数增长，导致流动状态无法长期精确预测。这对飞机机翼设计、气象预报至关重要。
开始使用此模拟器时，有几个需要注意的要点。首先，是“初始位置稍有变化就会导致结果大相径庭，因此计算必须极度精确”这一误解。虽然初始值敏感性确实是混沌的核心，但模拟本身的数值误差也会被视为“另一个初始值”。因此，计算方法（本工具中为RK4）和时间步长（Δt）的选择至关重要。例如，将Δt从0.01秒粗略地改为0.1秒，会导致外观轨迹截然不同，从而观察到的将不是原本的“确定性混沌”，而是“数值误差导致的混沌”。在实际工程中，也存在因网格尺寸或时间步长设置过于粗糙而产生非物理性振荡的陷阱。
其次，是参数设置的平衡。将阻尼系数γ设得接近零可以保持能量守恒，从而绘制出复杂而优美的轨迹，但在现实物理系统中，空气阻力等因素必然导致阻尼存在。反之，若γ设置过大，摆锤会迅速停止运动，使混沌的趣味性大打折扣。例如，在质量m=1、磁力k=1000的系统中，γ=0.1〜0.5通常是适度阻尼与复杂行为之间的平衡点。在现实世界的CAE分析中，将阻尼器的阻尼系数设置为不切实际的值，导致模拟结果与实际试验数据不符，就是一个典型案例。
最后，切勿认为“此模拟完全真实”。这里使用的磁力模型是简化模型，并未完全再现现实磁偶极子间的力。参数ε同时也是为稳定计算而设的“调整项”。这与CAE中先用易于理解的线性弹簧模型近似复杂现象的思路相通。第一步的目标并非完美再现，而是提取现象的本质。
此磁力摆模拟器所使用的技术和概念，直接关联到多个尖端工程领域。首先要提的是汽车和飞机的“操纵稳定性分析”。车辆高速过弯时，轮胎侧向力与侧偏角呈非线性关系。微小的方向盘转角或路面状况差异，有时会使车辆行为变得难以预测，这也是一种混沌行为。模拟器中使用的状态空间表示法和RK4法，正是求解此类车辆运动微分方程的相同手法。
其次是结构物的“非线性振动分析”。例如，承受强风的超高层建筑或大跨度桥梁，当变形增大时，材料特性或几何形状发生变化，导致恢复力不再与位移成正比（呈现非线性）。此时，特定风速下的小扰动可能引发巨大振幅（混沌振动）。与存在三个磁铁的情况类似，结构物也存在多个复杂的模态（振动形态），它们之间会进行能量交换。
再者，“机器人手臂的轨迹规划与控制”也密切相关。多关节机器人手臂的动力学具有强非线性，约束条件也很复杂。到达目标点的轨迹看似无数（如同吸引盆），且表现出初始值敏感性——初始姿态的微小误差会导致最终姿态的巨大偏差。在此培养的“数值追踪非线性系统行为”的直觉，可以说是机器人模拟的基础本身。
如果对此模拟器产生了兴趣，接下来可以按以下步骤深化学习。首先从数学背景入手。推荐通过混沌理论入门书学习“逻辑斯蒂映射”。$$ x_{n+1} = a x_n (1 - x_n) $$ 这个极其简单的递推公式，仅通过改变参数a，就能经历稳定→周期倍分岔→混沌这一与磁力摆相同的过程。由此可以把握“非线性正是混沌之源”的直觉。
接着，学习模拟技术本身。针对本工具的核心“数值积分法”，建议尝试对欧拉法、改进欧拉法（休恩法）以及RK4法进行实际编码，并用简单问题（例如自由落体或弹簧质量系统）进行比较。你将能切身体会RK4为何精度更高，以及其与计算成本（函数求值次数）之间的权衡。这将成为理解CAE软件“求解器”选择的基础知识。
最后，了解“非混沌的系统”也很重要。学习线性系统以及混沌的前一阶段——分岔理论，反而能凸显出混沌的特殊性。在工程上，大多数设计都是在“不发生混沌的稳定区域”内进行的。在此模拟器中调整参数并探索“混沌从何处开始（分岔点）”，这一过程类似于在实际工作中确认安全裕度。下一个挑战主题，便是尝试数值计算此磁力摆的“李雅普诺夫指数”，以量化混沌的强度。