编辑注意：这个文件被选为MATLAB中心挑选本周
对于动画，程序使用MATLAB的普通绘图命令与raiveO命令相结合。它基于相应微分方程的ODE45解决方案。
提供了全面的文档，包括如何衍生运动方程的最重要步骤的草图。简单的Mathematica笔记本包含所有操作。
通过>> double_pendulum_init启动程序
您可以在double_pendulum_init中调整诸如棒子长度、质量、初始角位置和速度等参数。m文件。
在未来的版本中实现一些特性可能会很有趣，比如角时间依赖图、相空间画像或Poincaré部分。
玩得开心观察这种简单的丰富的动态行为，但主要是混乱的系统:-)
引用
Alexander Erlich（2021年）。动画双边形(//www.tatmou.com/matlabcentral/fileexchange/27212-animated-double-pendulum)， MATLAB中央文件交换。检索。
意见及评分(16.）
您还可以从以下列表中选择一个网站：
没有运行。
刘。是的，您是对的，内置Matlab ode求解器不保守汉密尔顿系统中的能量。您需要使用不同的求解器，例如基于VERRET方法
有谁知道为什么这个模拟在较大的时间内失败超过100秒（例如1000秒）？似乎总能量不保守。鉴于拉格朗日表格是否正确，它可能是一个matlab问题吗？
非常好的例证。然而，对于为什么摆在初始化上方的均衡时，钟摆开始移动的原因很好奇
（phi1 = pi; dtphi1 = 0; phi2 = pi; dtphi2 = 0;）
这是由于matlab的近似错误，数字pi是用有限个小数表示的吗?
令人惊叹的是，这个程序可能是进一步作品的基本参考。
结果与其他项目和在线模拟器非常不一致，但它们似乎使用迭代计算而不是IVP DEQ的计算。M-File中的方程也与PDF文档不一致。请问你能详细说明吗？
精彩的！
10 !