基于while循环的Runge-Kutta算法我目前正在尝试使龙格库塔4集成工作,但它输出以下内容: Runge Kutta value: inf 虽然它应该给出一个在以下范围内的值: ^{pr2}$ 我使用了以下代码,但似乎无法输出正确的近似 ...2024-05-29 已阅读: n次
在Python脚本中定义全局函数我是Python新手。我正在写一个脚本,将用Runge-Kutta方法数值积分一组常微分方程。由于Runge-Kutta方法是一种有用的数学算法,我将它放在自己的.py文件rk4.py中。 def r ...2024-05-29 已阅读: n次
SciPy中的三维插值我有一个形状[544, 256, 256]数组的三维数据。在每个网格点,都有一个函数的值,比如F(x,y,z)。我需要在网格点之间插值,作为Runge-Kutta四阶格式的一部分。在 例如,我希望能够 ...2024-05-29 已阅读: n次
在双_标量中遇到RungeKutta方法溢出我不是编程专家,但我试图通过尝试边界条件的错误来模拟龙格-库塔方法。这对我来说没问题,只是如果一个程序试图扩展时间域(即通过增加xStop向向量添加更多项),程序会突然爆炸并溢出 你能告诉我为什么它这 ...2024-05-29 已阅读: n次
定义更复杂的静态数组数值方法中的系数往往是固定的。我想知道Cython/C中设置这样的数组或变量的最佳方法是什么。在 在我的例子中,除了系数和级数外,Runge-Kutta积分方法基本相同。现在我正在做类似(简化) # ...2024-05-29 已阅读: n次
Python中将一组函数作为参数传递给另一个函数我想用Python实现我自己的基本Runge-kutta4集成器。格式应该是这样的: ---- EXAMPLE set of equations ---- f1 = lambda x: x**2 f2 ...2024-05-29 已阅读: n次
计算rk方法/无打印的函数import math import matplotlib import numpy as np from numpy import linspace tmax=10.0 n=2000 G=4 D ...2024-05-29 已阅读: n次
如何在Python中将函数传递给函数?我是Python的初学者/中级者。我已经将4th-order Runge-Kutta method (RK4)编码到Python中。它基本上是解一个摆,但这不是重点。在 我想通过以下方式改进RK4方法 ...2024-05-29 已阅读: n次
在RK4算法中使用lambda函数在Python-showhere中有两种实现经典Runge-Kutta格式的方法。第一个使用lambda函数,第二个没有lambda函数。在 哪一个会更快?为什么? ...2024-05-29 已阅读: n次
如何将数学函数放入python函数中? 我用Runge和Kutta算法来求解dy/dx=f (x,y)型微分方程。你知道吗 我不想在我的代码中有def f(x,y),我想把它输入我的Runge和Kutta函数的参数中 RK4 (x0,y0 ...2024-05-29 已阅读: n次
“无法在不同大小的向量值之间转换”我试着用C语言运行Runge Kutta 4。代码如下: __kernel void updade_state( __global float8 *q,__global float8 *qm, __g ...2024-05-29 已阅读: n次
附加到NumPy(Python)数组首先,我在三自由度纽比阵列上做一个龙格库塔。我的阵列是这样的: 状态=[[X]、[Vx]、[Y]、[Vy]、[Z]、[Vz]] 我运行Runge-Kutta,得到四个K值,然后用[newaxis]进行 ...2024-05-29 已阅读: n次
tigramite#Tigramite–时间序列数据集的因果发现3.0版,如http://arxiv.org/abs/1702.07007(python包)[Github](https://github.com/jak ...2024-05-29 已阅读: n次
pyradau13radau13是一个隐式的13阶自适应runge-kutta实现 Ernst Haier和Gerhard Wanner的Radau IIA方法。Fortran代码是 可从安永海尔的网站获得 http ...2024-05-29 已阅读: n次