谐振子的Lagrange运动方程构建求解
过冷水诚挚邀请你加入Matlab仿真秀官方交流群进行Matlab学习、问题咨询、 Matlab相关资料下载，群号：927550334
在力学的学习过程中经典分析力学是最基本的入门知识，过冷水之前和大家一起学习了两个小车通过弹簧链接起来的做来回摆动运动的运动轨迹学习。推文中直接给了一个微分方程组，然后解出微分方程组就得到了小车的演化轨迹。本期过冷水带从零开始构建一个微分方程组，而不是单纯解微分方程。
假设有一个一维的谐振子在三维空间中做运动。m=5 kg；kx=2 N/m;ky=7 N/m;kz=4 N/m。根据三维空间受力分析可得：
谐振子的势能函数为：
过冷水在这里需要给大家补充一点点初中弹簧做功的知识。W=Fr是做功的函数，已知F=kr;则W=kr2，而这里势能却是W=0.5kr2，这是因为错误使用W=Fr，该公式是恒力做功，而弹簧力是变力做功，需要使用积分来算功。
谐振子的动能为：
谐振子的总能量为：
谐振子的Lagrange函数为：
已知Lagrange方程方程形式：
Lagrange方程的证明较为复杂，在此过冷水只告诉打算使用了该方程而不证明方程的由来，感兴趣的可以找一本力学方面的书进行探究。Lagrange是一个微分方程，针对于我们的情况，可以很容易的就解出原函数。Lagrange方程的通解为：
设初始条件t=0时:
解得：
这样我们就求得了一个谐振子在三维空间的运动方程。可以得到如下的运动轨迹图
根据轨迹图像可以看出谐振子在三维空间中运动轨迹很复杂。这只是动力学上比较简单一个运动问题。过冷水要带大家做的是有很多个相互作用的空间运动问题。不积跬步无以至千里，本期先和大家分享这个简单的案例，下期带大家了解两个谐振子的一维运动图像。想要深入了解学习请持续关注过冷水的推文。
过冷水发表于 仿真秀 平台原创文章，未经授权禁止私自转载，如需转载请需要和作者沟通表明授权声明，未授权文章皆视为侵权行为，必将追责。如果您希望加入Matlab仿真秀官方交流群进行Matlab学习、问题咨询、 Matlab相关资料下载均可加群：927550334。
精品回顾
过冷水和你分享 matlab读取存储各种文件的方法 文末有独家金曲分享
