文件名称：欧拉公式求圆周率的matlab代码-hw8:到期至16.01.16
文件大小：85KB
文件格式：ZIP
更新时间：2024-06-12 06:53:01
系统开源
欧拉公式求长期率的matlab代码作业8 到期至16.01.16 考虑两个引力相互吸引的物体。 我们选择一个物体作为坐标系的原点。 由于运动将发生在一个平面上，因此我们可以通过二维矢量q =（q_1，q_2）描述第二个物体的位置。 运动的哈密顿量是 从解析解到问题，我们知道，如果H <0，则运动沿椭圆发生。 令e为椭圆的偏心率。 让我们考虑一下初始条件 这意味着H ＝ -1 / 2。 解的周期为2。下图显示了tEnd = 20的数值解，是在dt = 0.0005处经过约125 000欧拉前移之后获得的。 显然，解决方案并不令人满意。 辛欧拉方法是解决该问题的一种更好的算法。 它在时间上像Euler一样是一阶的，并且属于分区方法的组。 我们可以写成 编写一个程序，使用辛欧拉方法对哈密顿两体问题进行积分。 使用上述初始条件，并生成与上述正向Euler方法相同的图。 在第二个图中检查H的守恒性。 选择dt = 0.0005和dt = 0.05。 最后，提交您的代码和四个图表（.png图像）。 分别显示两个dt的轨道和H。
【文件预览】：
hw8-master
----plot2.png(5KB)
----plot1.png(5KB)
----.travis.yml(145B)
----README.md(2KB)
----Makefile(55B)
----kepler.cxx(1KB)
----plot3.png(7KB)
----stuffy_stuff()
--------se.png(20KB)
--------f2.png(16KB)
--------f1.png(15KB)
--------f3.png(5KB)
--------euler.png(8KB)
----plot4.png(5KB)