龙格-库塔RK4不比弗雷特好吗?我只是在游戏中测试几个轨道动力学的积分方案。 我在这里采用了RK4的恒定自适应步骤 http://www.physics.buffalo.edu/phy410-505/2011/topic2/app1 ...2024-09-27 已阅读: n次
Euler/RK4-python轨道仿真我开始为8个行星的轨道模拟编写代码(太阳在(0,0)中) 我得到了X和Y的加速度,现在我需要用它们来得到X和Y的速度,先用欧拉,然后是RK4,我被卡住了,我不知道从哪里开始。这是我的密码: ...2024-09-27 已阅读: n次
我的rk4实现的加权平均值中的不同近似值都是sam好吧,我想写一个微分方程数值逼近的四阶龙格库塔方法的实现,这是我数学课程的一部分,也是为了学习一些编程,但问题是它在逼近的每一步中都使用这些系数。我们从一个x值和一个y值开始,它想为后面的x值找到y值 ...2024-09-27 已阅读: n次
Python中ODE集成的Runge–Kutta方法,带有附加约束考虑到一阶导数的附加约束,我有一个关于使用RK4求解二阶微分方程的问题。我正在做一些修改后的here示例 θ′′(t) + b θ′(t) + c sin(θ(t)) = 0 附加约束是: wh ...2024-09-27 已阅读: n次
Streamplot python模块“ImportError:没有名为Streamplot的模块”我试图运行这段代码,但问题是"ImportError: No module named streamplot" 请帮助修复此错误 import matplotlib.pyplot as plt imp ...2024-09-27 已阅读: n次
如何解决IndexError: list assignment index out of range错误?我用rk4方法编写了一个求解ode的代码。 代码如下: from __future__ import division import math import numpy as np from math ...2024-09-27 已阅读: n次
阻尼振荡器动画我能够模拟阻尼振动下的质量弹簧系统。然而,我想添加一个子地块的位置对时间和另一个子地块的速度对位置(相位路径),这样我将有三个动画。如何添加它们?我使用的源代码如下所示 更新:我尝试添加第一个子批次位 ...2024-09-27 已阅读: n次
python中的Runge-Kutta方法以下是我编写的函数: def rk4(f, t0, y0, h, N): t = t0 + arange(N+1)*h y = zeros((N+1, size(y0))) y ...2024-09-27 已阅读: n次
在Python脚本中定义全局函数我是Python新手。我正在写一个脚本,将用Runge-Kutta方法数值积分一组常微分方程。由于Runge-Kutta方法是一种有用的数学算法,我将它放在自己的.py文件rk4.py中。 def r ...2024-09-27 已阅读: n次
使用Python中的ctypes访问数组我正在用C编写一个ode解算器,导出到一个WindowsDLL和一个DLL的Python包装器。我非常习惯Python,但我对C和ctypes也是一个完全的初学者。在 受公认答案here启发的修改后的 ...2024-09-27 已阅读: n次
用RK4求解大型耦合微分方程组我想用RK4方法解一个微分方程组,如下所示 对于这种情况,我要解3个微分方程- d/dt(x0)=-2*x1 d/dt(x1)=-2*x0-8*x1-2√2*x2和 d/dt(x2)=-2√2*x1- ...2024-09-27 已阅读: n次
如何使用张量作为函数的输入?(张量流)我试着用张量流来构造我的模型来解一个微分方程,例如 dX/dt=f(\mu,X,t) 这里,\ mu是一个依赖于X的函数,它很复杂,所以我想用神经网络来预测\mu(X)。你知道吗 首先,我的输入X通 ...2024-09-27 已阅读: n次