回答情况：
提问时间：
问题标签：
我正在尝试实现Runge 4来求解一个微分方程，并需要一些洞察力。我重复一遍:我不是在找解决方案，我在找人帮我指出我的代码出了什么问题，因为我已经盯着这个看了多久了.
代码的工作原理:在函数声明和for循环之间设置初始值(同样，p和q在问
浏览 0提问于2015-01-20得票数 1
回答已采纳
当使用数值方法(如)时，浮子在计算机上的有限精度会影响解(布鲁沃尔定律)。有人能在实践中解释一下这是怎么回事吗？如果A= 3/10，那么如何确定B和C？
浏览 0提问于2018-09-11得票数 1
回答已采纳
我正在对两个方程组实现Runge四阶方法。H是段数，所以T/h是步长。 x1 = [x10] k11 = f1((i-1)*T/h, x1[-1], x2[-1]) k21 = f1((i-1)*T/h + T/h/2, x1[-1] + T/h/2*k11, x2[-1] + T/h/2*k12)
浏览 0提问于2016-12-11得票数 0
回答已采纳
我试图把这个FORTRAN程序(钟摆的运动)转换成CUDA FORTRAN，但是我只能用一个带两个线程的块。有没有办法使用超过两个线程. REAL :: Q,B,W USE CB INTEGER, PARAMETER :: N=10,L=100,M=1 REAL :: PI,H,T,Y1,Y2,G1,G1F,G2,G2F
REAL :: DK1
浏览 5修改于2013-01-29得票数 3
回答已采纳
我用python编写了关于runge方法的代码，但是每次当程序实现任何微积分时，程序都需要微分方程。from math import *n=input("Enter the number of equations n:")x=np.array([]) x0=input("Enter the initial value of x{}:".format(i))
x=np.append(x,[x0])
浏览 2提问于2015-07-30得票数 0
回答已采纳
假设我有一个4 ODE: dX/dt = F(X)的系统，其中X是一个向量(4维)，F: R^4 -> R^4，F称为vectorDE_total_function，我试图用RK-4来计算解。 previous_vector = np.array ([theta_1, omega_1, theta_2, omega_2]); temp_vector = previous_vector;
RK_vector = np.array([0.0,0.0,0.0,
浏览 1修改于2016-12-11得票数 1
回答已采纳
我编写了一个简单的Runge四阶积分程序来求解一个常微分方程组，并使用OpenMP并行化。我不知道这是否是我们所能做的最好的，以最大限度地提高代码的性能，只需付出很少的努力。以下是我的尝试： (链接到gitlab存储库)void xprsys(const int n,const vector<double>& x, vector<d
浏览 0修改于2017-07-06得票数 5
回答已采纳
我已经在MATLAB中编写了一个自适应步长RK4来求解常微分方程组。代码运行时没有错误，但是当我尝试绘制x与y的关系时，它没有生成所需的曲线。我只是得到了一条平坦的直线，而不是环形。从r输出常量值的事实可以明显看出这一点。检查每一行的输出后，它们不是输出常量或错误，也不是inf或NaN，而是输出实部和虚部(复数)。我不知道为什么会发生这种情况，我相信这是我麻烦的根源。parsec = 3.08*10^18;theta_1 = 0.0;g
浏览 2提问于2016-09-16得票数 1
我已经用Java实现了Runge 4算法，正如本文所描述的。protected void integrateRK4(double time) {
final double HO2 = H/2; Vector2[] currentPositions = new Vector2[
浏览 4修改于2017-03-17得票数 1
回答已采纳
我试图用4方法(RK4)来解决一个ODE系统。我正在对下面的算法进行代码测试，发现解决方案不等于解析解(而且误差很大)。下面我包含了I.V.P.dy/dt= f(t，y)的代码测试。我试图在这段代码中找到错误，但找不到它们。任何帮助都是非常感谢的。全球t增长率clear-allset dt .01reset-tickstickask turtles [ set state
浏览 2提问于2015-10-13得票数 0
回答已采纳
我正在用Octave建立一个函数，它可以求解N耦合的常微分方程类型：dy/dt = G(x,y,…,z,t)这三种方法中的任何一种(Euler，Heun和Runge 4)。function sol = coupled_ode(E, dfuns, steps, a, b, ini, method) h=range/steps; columns = size(dfuns)(2)+1;
sol=
浏览 0修改于2017-10-10得票数 2
回答已采纳
我尝试使用Runge-Kutta方法将其与lsode函数进行比较。但它的性能相当差，我使用的其他方法(向前和向后Euler，Heun)都比lsode做得更好，以至于它们与lsode几乎没有区别。下面是我使用的Runge-Kutta方法v(1) = 2;q(1) = 1/sqrt(3);
hold on; plot(
浏览 2修改于2019-01-09得票数 0
我试图在MATLAB中实现。我没有得到，我不确定代码或运行它所用的命令中是否有问题。function RKSystems(a, b, m, N, alpha, f) t = a; w(j) = alpha(j); for i = 1:m
浏览 3提问于2013-02-18得票数 2
回答已采纳
我在这个问题上遇到了一些麻烦。我得到了这组方程dy / dt = x + a * y以及默认值a=0.1、b=0.1、c=14以及龙格-库塔算法： T = np.linspace(0,Tmax,N+1) xv[0] = xvinit for i in range(N):
浏览 0修改于2020-08-04得票数 1
你好，我要编写python函数，用Runge-Kutta 2 2cond级数求解Lorenz微分方程。
浏览 12修改于2020-06-15得票数 1
回答已采纳
我为Runge 4编写了代码，用于求解ODEs系统。让u1 = x和u2 = x' = u1'，然后系统看起来像：u2' = -k*u1u' = f(u, t) u, v = u我的全部代码是：
im
浏览 1修改于2020-09-09得票数 0
回答已采纳
我已经有一段时间没有做任何C/c++了，但是我想使用gsl库编写一个ODE解算器来求解下面的ODE集$$ up'(r)=-(2*(r-1)/(r*(r-2)))*up(r)-((r*r/((r-2)*(r-2)))-(2/r*(r-2)))*u(r) $$
因此，在gsl表示法中，我的y=u、y1==up和上面的RHS定义了f和f1。根据这些定义，可以计算雅可比矩阵和dfdr (通常将它们的'time‘变量称为't’而不是'r')。这样做的原因是因为我在使用Mathematica时遇到了速度问题。我在ODE解算
浏览 5提问于2012-10-25得票数 1
回答已采纳
好的，我试着解决一个初值问题，试着用龙格-库塔2找到f函数给出的速度，我的程序会编译，但当我运行它时，我在任何时候都会得到相同的v值，但我就是找不到哪里做错了。有人能帮帮我吗？implicit none
real(8) a,b,hinteger i,n,j
a=0h=(b-a)/nm=70.0v0=20.0t=a
open(1, file="velocida
浏览 0修改于2017-06-27得票数 0
double xin, double yin[], double h) double yout[] = new double[N];
/* fourth order Runge-Kutta
浏览 3修改于2017-05-18得票数 0
回答已采纳
我有一个耦合方程组:静水平衡方程、质量连续性方程和理想气体的状态方程。在数学语法中，
我想用Runge数值技术求解这些耦合方程，这里我展示了我为解决这个问题而设计的python代码：from pylab i
浏览 10修改于2016-12-28得票数 2
回答已采纳
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号