在天体力学中，数值方法被广泛用于求解微分方程。 本代码中，根据牛顿万有引力定律，利用Runge-Kutta四阶方法对轨道运动方程进行数值积分，模拟物体绕地球运动的轨迹。
输入：位置和速度向量 (x,y,z,vx,vy,vz) 或者开普勒元素 (a, e, i, Omega, w, M) h = 步长步数 = 步数输出：在ECI参考系中传播的卫星PV矢量
调用：[X_RK] = RK_4(X,h,steps)
2021-09-29 19:06:33
4KB
1
这些文件提供了用于求解倾斜冲击关系和Taylor-Maccoll方程的数值程序。 采用四阶Runge-Kutta数值格式隐式求解Taylor-Maccoll方程。
使用了反向方法（JD Anderson，现代可压缩流，第10.4节）
为以下内容计算流的属性-超音速马赫数-零俯仰和偏航-粘稠的完美气体
注-所产生的冲击波本质上是3D的，但是由于冲击是局部平面的，因此可以通过使用2D斜向冲击理论来对其进行局部处理
如果需要，可以将图形注释掉。
2021-09-10 18:56:28
11KB
1
光的光子量子统计可以通过高阶相干来显示。 基于双重Banbury Brown Twiss（HBT）方案，研究了各种量子态的四阶相干，包括Pock态，相干态，热态和压缩真空态。 分析结果是通过综合考虑整体效率和背景而得出的。
2021-08-27 10:33:35
150KB
1
微分方程组数值求解四阶显式龙格库塔法自编程实现-龙格-库塔法求解常微分方程-Matalab自编程序实现.pdf
常微分方程（组）的数值解法有多种，其中Runge-Kutta方法是最为流行算法之一，
几乎可以解决所有常微分方程组的数值解问题。
论坛中不少朋友执着于利用MATLAB自置函数solve求解微分方程。其实solve函数存在局限性
（比如系统输入是时变的，solve函数就不易处理），也不利于大家掌握算法。
本人在此提供两个微分方程组的数值算例及程序代码，供大家参考（附件）。算法采用显式4阶Runge-Kutta算法。
关于Runge-Kutta算法，大家可以参考一般 数值方法 的教材即可。
2021-08-17 11:51:54
81KB
1