力学/混沌
双摆混沌模拟器
使用四阶龙格-库塔法实时数值积分双摆运动方程。绘制初始角度相差仅 0.001° 的两条轨迹，直观体验混沌的『蝴蝶效应』。
初期条件
预设
纵轴（对数尺度）为两条轨道末端距离 d(t)。斜率约等于 Lyapunov 指数 λ。
理论与主要公式
$\ddot\theta_1 = \dfrac{-r\dot\theta_2^2\sin\Delta - 2g\sin\theta_1 - \cos\Delta\bigl(\dot\theta_1^2\sin\Delta - g\sin\theta_2\bigr)}{2-\cos^2\Delta}$
$\ddot\theta_2 = \dfrac{2\dot\theta_1^2\sin\Delta - 2g\sin\theta_2 + \cos\Delta(r\dot\theta_2^2\sin\Delta + 2g\sin\theta_1)}{r(2-\cos^2\Delta)}$
$\Delta = \theta_1 - \theta_2$
積分：Runge-Kutta 4次法、dt=0.005 s
🙋 “确定性的却不可预测”是什么意思？
🙋
双摆明明有严格的力学方程，为什么被称为“混沌”？有方程的话，未来应该可以计算吧？
🎓
理论上确实如此，但问题在于“即使相同的初始条件，观测误差也必然存在”。用这个模拟器试试看吧。同时观察从θ₁=120°和120.001°（差0.001°）出发的两条轨迹，最初是重合的，但几秒后就会变成完全不同的运动。这就是“初值敏感性”——混沌的本质。
🙋
在“轨迹发散”选项卡中，距离越来越大。这就是指数增长吗？
🎓
是的。纵轴是对数刻度，所以如果呈直线增长，则d(t) ∝ exp(λt)，即指数增长。这个斜率λ就是李雅普诺夫指数，双摆的典型值约为λ ≈ 1〜3 s⁻¹。例如，λ=2时，1秒后差异变为e²≈7倍。即使位置精度达到毫米量级，几秒后也会变成米量级的误差。
🙋
但是，如果使用“小振幅（线性）”预设，两条轨迹一直保持重合。不会变成混沌吗？
🎓
重大发现。在小振幅（θ≪1 rad）下，sin θ ≈ θ的线性化成立，变成具有两个固有振动模式的“耦合振动”。这是可积系统，可以进行长期预测。混沌只出现在能量较大（摆幅较大）时。在物理学中，“线性→可积→可预测”“非线性→可能混沌→预测困难”是常见模式。
🙋
天气预报说“超过一周就很难预测”，也是同样的原因吗？
🎓
正是如此。洛伦兹在1960年代发现的“洛伦兹吸引子”是大气简化模型，它表明混沌本质上使气象预测变得困难。由于大气的初始观测误差大约在两周内增长到整个系统的尺度，因此超过两周的精确预报在原理上被认为是困难的。这就是蝴蝶效应名称的由来——“巴西的蝴蝶扇动翅膀引发得克萨斯的龙卷风”这个比喻。
🎓
实际上应用很多。利用混沌流进行流体混合，可以用更少的能量实现均匀混合（如微流控芯片）。混沌加密利用初值敏感性作为加密密钥。此外，分析心电图的混沌特性可以区分健康心跳和心律失常。混沌不是随机的，而是“具有结构的复杂性”，因此有控制和利用的空间。
常见问题
混沌是确定性的——从相同的初始条件出发，一定会遵循相同的轨迹。随机是概率性的（真正的不确定性），即使初始条件相同，每次也会得到不同的结果。双摆由完整的方程描述，因此完全不涉及随机数。但重要的是，由于观测精度有限，实际上变得无法预测。
单摆的运动方程 θ̈ + (g/L)sinθ = 0 是保守系统（能量守恒）的单自由度系统，根据刘维尔定理，相空间体积守恒。根据阿诺德-刘维尔定理，它是可积系统（可通过积分求解），并且根据KAM定理，运动停留在环面上，不会产生混沌。在二自由度以上的非线性系统中，可积性被破坏，可能产生混沌。
并行积分从初始状态 x₀ 和微小偏移 x₀+δ₀ 出发的两条轨迹，测量时刻 t 的距离 δ(t)。λ = lim_{t→∞} (1/t) * ln(δ(t)/δ₀) 就是李雅普诺夫指数。实际上，当距离变得过大时，需要重新归一化（将δ重置为δ₀并继续测量δ）。本模拟器显示的是无重新归一化的简单发散，因此经过足够长的时间后精度会下降。
常微分方程的数值解法，每一步进行四次梯度评估，并使用类似辛普森法则的加权方式计算下一步。由于四阶精度（局部误差∝dt⁵），比简单的欧拉法（一阶）精度更高，即使在双摆这样的非线性系统中也能相对稳定地进行积分。在混沌中，虽然无法追踪真实轨迹，但可以准确计算统计特性（如李雅普诺夫指数）。
在无摩擦（保守系统）设置下，总机械能 E = T + V 应保持恒定。数值积分的误差会导致能量略微漂移，但使用RK4和dt=0.005时，即使在长时间内误差也能被抑制得很小。本模拟器在统计卡片中显示“总机械能”，因此可以检查能量是否随时间保持恒定。
什么是双摆混沌模拟器？
双摆混沌模拟器是CAE和应用物理中的重要基础课题。本交互式模拟器允许您直接调节参数并观察实时结果，从而理解关键规律和变量之间的关系。
通过将数值计算与可视化反馈相结合，本模拟器有效地弥合了抽象理论与物理直觉之间的鸿沟，既是学生的高效学习工具，也是工程师进行快速验算的实用手段。
物理模型与关键公式
本模拟器基于双摆混沌模拟器的核心控制方程构建。理解这些方程有助于正确解读计算结果，并判断参数变化对系统行为的影响。
方程中的每个参数都对应控制面板中的一个滑块。移动滑块时，方程的解会实时更新，帮助您直观建立数学表达式与物理行为之间的对应关系。
实际应用场景
工程设计：双摆混沌模拟器相关概念可用于工程初步估算、参数灵敏度分析和教学演示。在开展更完整的CAE分析之前，可借助本工具快速把握主要物理量级与趋势。
教育与科研：在工程教学中，本工具可将理论与数值计算有效结合。在科研阶段，也可作为假设验证的第一步工具使用。
CAE工作流集成：在运行有限元（FEM）或计算流体力学（CFD）仿真之前，工程师通常先用简化模型评估物理量级、识别主导参数，并确定合理的边界条件，本工具正是为此目的而设计。
常见误解与注意事项
模型假设：本模拟器所用数学模型基于线性、均质、各向同性等简化假设。在将计算结果直接用于设计决策之前，务必确认实际系统是否满足这些假设。
单位与量纲：许多计算错误源于单位换算错误或数量级判断失误。请时刻注意各参数输入框旁标注的单位。
结果验证：始终将模拟器输出结果与物理直觉或手算结果进行核对。若结果出乎意料，请检查输入参数或采用独立方法进行验证。