2回答
我使用Numerical Recipes 3，但我经常在编译它们的库时遇到麻烦。我想知道其他的选择。
浏览 6提问于2011-10-02得票数 14
回答已采纳
1回答
我用四阶龙格-库塔法求解两个耦合的常微分方程。由于应用此方法，我在打印z的值时遇到了问题。以下是源代码以供参考。请通过打印z和theta的更新值来帮助我修复此代码。谢谢。k4) * dt / 6) z, dz, theta, dtheta = np.asarray(u).T
pri
浏览 33提问于2021-04-17得票数 1
我正在尝试用python中的四阶龙格-库塔法求解椭圆型微分方程。在执行之后，我得到了应该获得的实际绘图的一小部分，并伴随着一个错误，说明： "RuntimeWarning: invalid value encountered in double_scalars" importx[i]+dx#print(y)plt.plot(x,
浏览 16提问于2021-01-21得票数 2
回答已采纳
3回答
我有一个二维向量场，也就是说，对于每一点(x, y)，我有一个向量(u, v)，而u和v是x和y的函数。📷
我需要一个算法，它能有效地计算给定向量场的一些轨迹。轨迹必须满足平面上的某种最小密度(对于平面上的每一点，我们必须有一个“附近的”轨迹)，或者其他一些条件来得到一个合理的
浏览 0提问于2012-09-08得票数 2
回答已采纳
我用c++写了一个可以解决微分方程的程序。问题是，它似乎不能很好地与ROOT一起工作。它编译得很好，但是当我执行的时候，我得到的是：
#6 0x0000000000403b8a in main ()这是我的程
浏览 0提问于2016-01-22得票数 0
我是Matlab编程的初学者，也是Runge方法的初学者。这是我的问题..。使用odearguments错误(第93行)初始条件向量为5。由@(N，X)K1;K2;K3;K4;K5返回的向量和初始条件向量必须具有
相同数量的元素。ode45中的错误
浏览 2提问于2015-03-28得票数 1
1回答
我正在为我的一门课做作业，我应该用我选择的程序(我选择了Matlab)写一段代码，用四阶龙格-库塔法求解贝塞尔函数微分方程。我可以把它分成两个耦合的一阶DEs：(Z_n)'+(1/x)*Z_n+(x^2-n^2)/x^2*J_n=0。在这项任务之前，我没有使用Matlab或任何其他编程语言的
浏览 2提问于2013-04-14得票数 0
回答已采纳
我用四阶龙格-库塔法求解带有NumPy阵列的duffing振子的微分方程，但是我收到了一个错误。RuntimeWarning: overflow encountered in double_scalars 有谁知道可能导致溢出错误的错误来源，以及我可以尝试解决它的方法吗？n-1] + h/6 * (k
浏览 41提问于2021-04-27得票数 0
1回答
我知道我可以使用Sympy，但这里的目标是实现这个方法，代码已经准备好用Fortran语言编写，所以基本上我有一个数据库，其中有正确的解决方案值，我必须在代码中得到类似的解决方案。然而，我们有一些问题，我用线性方程组(一阶和二阶)做了同样的几次，但是这是牛顿万有引力定律的二阶非线性方程。代码没有错误，我的问题是我的代码做错了什么，我无法得到正确的结果。下面我将展示一些预期的值和我正在得到的<
浏览 1提问于2020-11-17得票数 1
回答已采纳
所以我用加速度计把x和y加速度计的值积分，然后把球往那个方向移动一步。所以我可能不得不改变我的策略，用某种时间解决方案来代替。现在，速度越快，我移动的步伐就越大。相反，也许我应该有一个定时器，每毫秒移动一个像素，如果速度是高的，或者只有10毫秒，如果
浏览 2提问于2014-07-23得票数 1
回答已采纳
我已经得到了一个步长h，用于实现高阶常微分方程的Runge-Kutta方法。我的问题出现在划分从开始时间t_0到最终时间t_f的时间间隔时。t = linspace(t_0, t_f, N)t = t_0:h:t_f
但这并不总是包括端点t_f。
浏览 1提问于2015-08-07得票数 0
你好，我要编写python函数，用Runge-Kutta 2 2cond级数求解Lorenz微分方程。
初始条件为(x，y，z)=(0,1,0)from
浏览 12提问于2020-06-14得票数 1
回答已采纳
我正在尝试使用四阶龙格库塔方法来求解Lorenz system，其中 dx/dt=a*(y-x)dx/dt=x*y-c*z 由于这个系统不依赖于时间的显式，它可能会忽略迭代中的这一部分result[time-1])
result[time]=result[time-1]+(1/6)*(k[0]+2*k[1]+2*k[2]+k[3]) 结果应该是Lorenz吸引子，但
浏览 23提问于2019-04-28得票数 1
回答已采纳
2回答
我在试着预测篮球的抛射运动。请注意，我没有考虑任何空气阻力的影响，只考虑了G的影响。为此，我使用了欧拉方法。def euler(euler_center, euler_velocity):(339, 167)(325, 151
浏览 0提问于2018-06-18得票数 1
我已经在网上找了一段时间了，为一个物理引擎寻找集成的方法，我正试图编写一个有趣的代码(我喜欢那里的书呆子:P)。我找到了欧拉法、RK4和Verlet (以及时间校正版本)。令人惊讶的是，在我所有的研究中，我没有找到我想要做的事情的技术名称！也许描述我的具体问题会让我的问题更清晰。假设我想要模拟一个球在圆形(或球形，一旦我实现了3d)引力场中移动。这个球将遇到力向量，它可以用来计算球在该特定刻度上所处的</em
浏览 0提问于2012-02-24得票数 1
回答已采纳
3回答
我有一个使用龙格库塔算法求解微分方程系统的问题。到目前为止，我已经将二阶PDE重写为一组两个耦合的方程，其中 g(L1,L2) = A*(B*L1-C*L2-D) k1 = f(L1,L2)
k2 = f(L1 + 0.5 * dt
浏览 3提问于2013-11-09得票数 0
1回答
我正在写一个获得谐振子系统真实速度的仿真程序，如下所示其中，当P(0)=1 0；0 1时，P=p1 p2;p2 p3可以用龙格-库塔积分法求解function z=getPos(x)
%% Function to get velocity
浏览 2提问于2014-11-17得票数 0
2回答
我读了费曼的物理讲座第九章，并尝试自己的模拟。我用Riemann积分来计算速度和位置。虽然所有的入门都是一样的，但我的轨道看起来就像双曲线。
浏览 5提问于2020-06-14得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
