基于while循环的Runge-Kutta算法我目前正在尝试使龙格库塔4集成工作,但它输出以下内容: Runge Kutta value: inf 虽然它应该给出一个在以下范围内的值: ^{pr2}$ 我使用了以下代码,但似乎无法输出正确的近似 ...2024-10-17 已阅读: n次
在Python脚本中定义全局函数我是Python新手。我正在写一个脚本,将用Runge-Kutta方法数值积分一组常微分方程。由于Runge-Kutta方法是一种有用的数学算法,我将它放在自己的.py文件rk4.py中。 def r ...2024-10-17 已阅读: n次
SciPy中的三维插值我有一个形状[544, 256, 256]数组的三维数据。在每个网格点,都有一个函数的值,比如F(x,y,z)。我需要在网格点之间插值,作为Runge-Kutta四阶格式的一部分。在 例如,我希望能够 ...2024-10-17 已阅读: n次
定义更复杂的静态数组数值方法中的系数往往是固定的。我想知道Cython/C中设置这样的数组或变量的最佳方法是什么。在 在我的例子中,除了系数和级数外,Runge-Kutta积分方法基本相同。现在我正在做类似(简化) # ...2024-10-17 已阅读: n次
计算rk方法/无打印的函数import math import matplotlib import numpy as np from numpy import linspace tmax=10.0 n=2000 G=4 D ...2024-10-17 已阅读: n次
如何在Python中将函数传递给函数?我是Python的初学者/中级者。我已经将4th-order Runge-Kutta method (RK4)编码到Python中。它基本上是解一个摆,但这不是重点。在 我想通过以下方式改进RK4方法 ...2024-10-17 已阅读: n次
在RK4算法中使用lambda函数在Python-showhere中有两种实现经典Runge-Kutta格式的方法。第一个使用lambda函数,第二个没有lambda函数。在 哪一个会更快?为什么? ...2024-10-17 已阅读: n次
如何将数学函数放入python函数中? 我用Runge和Kutta算法来求解dy/dx=f (x,y)型微分方程。你知道吗 我不想在我的代码中有def f(x,y),我想把它输入我的Runge和Kutta函数的参数中 RK4 (x0,y0 ...2024-10-17 已阅读: n次
“无法在不同大小的向量值之间转换”我试着用C语言运行Runge Kutta 4。代码如下: __kernel void updade_state( __global float8 *q,__global float8 *qm, __g ...2024-10-17 已阅读: n次
附加到NumPy(Python)数组首先,我在三自由度纽比阵列上做一个龙格库塔。我的阵列是这样的: 状态=[[X]、[Vx]、[Y]、[Vy]、[Z]、[Vz]] 我运行Runge-Kutta,得到四个K值,然后用[newaxis]进行 ...2024-10-17 已阅读: n次
如何在python中数值求解具有奇点的非线性常微分方程(BVP)?因此,我必须用数值方法求解以下常微分方程y'+f(y)*(y')^2=0:ODE 最初在[y_i,y_f]之间,其中y=y(t)和y(0)=y_i。主要问题是函数f(y)有一个(规则的)奇点(无法发布 ...2024-10-17 已阅读: n次
在python中使用matplotlib打印?我一直试图用四阶Runge-Kutta method来模拟一阶微分方程,但在绘制它时遇到了问题。你知道吗 #simulation of ode using 4th order rk method d ...2024-10-17 已阅读: n次