2回答
我有多个常微分方程问题，我正在解决。我需要解(u)和解的导数(du)。对于较小的ODEs，我可以执行以下操作
du[1]=@. u[2]
du[2]=@*u[1,:]*cos(2*pi*p[8]*t)*2*pi*p[8]/(p[2]*p[4]))/((1-u[2,:]/p[4])*u[1,:]+4*p[3
浏览 28提问于2019-07-13得票数 1
引言x''=-D(y) * x' * sqrt(x'^2 + y'^2)
y''=-D(y) * y' * sqrt(x'^2 + y'^2) + g(y)上述两
浏览 4提问于2017-01-24得票数 1
1回答
我想集成一个控制问题，即形式为dx/dt = A.f(t)的常微分方程，其中x(t)是R^3中的函数，f(t)是R^4中的函数，A是矩阵3x4。在我的特殊情况下，f(t) = F'(t)，即函数F的时间导数。此外，F是1周期的。因此，在区间0，1上积分ODE应该再次得到起始位置。有没有一种特殊的常微分方程解算器，可以将这种周期性的精度提高到很高的程度？
浏览 19提问于2021-01-22得票数 0
回答已采纳
1回答
在Python中，我可以只传递一个函数作为这个RK4函数作为输入获取的参数之一。如果我在Julia中这样做，会有性能损失吗？我的函数如下所示： Vs_k1 = rhs(Vs,Ps) Vs_k2 = rhs((Vs+(1/2)*Vs_k1),Ps) Vs_k3 = rhs((Vs+(1/2)*Vs_k2),Ps)
Vs_k3 = Ps.<e
浏览 2提问于2016-06-27得票数 2
考虑以下与常微分方程组的解相对应的代码片段db/dt = (b-a)^2/(b^2+c^2+1)
dc/dt = t^2*c^2/(a^2+c^2+如何包含z(0) =2的附加微分方程dz/dt = 5*da/dt？我不确定如何将它包含在函数句柄中，因为它是右侧的一个状态变量的导数，而
浏览 0提问于2018-04-04得票数 0
Introduction单d是一阶导数双d是二阶导数
通过数值求解控制LaGrange方程，用Matlab对一些动态系统进行了数值模拟。基本上是一组二阶常微分方程。我从Mathworks (下面的教程链接)上找到了一个关于如何求解基本的二阶常微分方程组的很好的教程。然而，我遇到了一些方程组，在这些方程组中，变量不能得到明确的解，就像弹簧摆的例子一样。对于一种情况，我有以下一组<e
浏览 4提问于2017-01-25得票数 0
1回答
我正在尝试使用julia提供的DifferentialEquations.jl，它工作正常，直到我尝试在二阶常微分方程上使用它。例如，考虑二阶常微分方程x'(0) = 0，x(0) =1using DifferentialEquations;
f
浏览 5提问于2021-05-04得票数 3
这个问题涉及到Python和Mathematica代码。A、B、C和D是我想要求解的四个耦合的复常微分方程。值得注意的是，abs(A)^2 + abs(B)^2 + abs(C)^2 + abs(D)^2 =1总是成立的(由于问题的性质)。实际上，在NDSolve的解决方案中总是如此
。因此，我的问题是，在NDSolve或solve_ivp的定义<em
浏览 3提问于2020-07-12得票数 1
1回答
我有一个函数中的4-4微分方程组(subsystem4)，我用odeint函数解决了它。我设法画出了系统的结果。我的问题是我想要绘制和一些其他方程(例如x，y，vcxdot...)它们都包含在相同的函数(subsystem4)中，但我得到了NameError:没有定义名称'vcxdot‘。此外，我想使用其中的一些方程(不仅是方程系统的结果)作为下面的微分方程系统的输入，并绘制同一时间段(t)内的所有方程。如何访问和绘制函数(
浏览 0提问于2019-01-08得票数 1
2回答
如何使用scipy.integrate.ode找到方程的一阶导数等于0的点？ # This function sets up the following relations
# dx/dt= v , dv/dt = -(Cp
浏览 4提问于2011-03-11得票数 2
回答已采纳
我在拉格朗日力学中遇到了下面的微分方程系统。你能推荐一种数值方法吗，有相关的链接和参考资料，关于我如何解决它。另外，在Matlab或Mathematica上有没有更短的实现？mx (y点)^2+ mgcosy - Mg - (M=m)(x双点) =0其中(x点)或(y点)= dx/dt或dy/dt，双点表示二次导数wrt
浏览 0提问于2012-09-07得票数 0
1回答
我需要帮助解决以下微分方程使用函数ode45从matlab。我的问题是方程上的一个时间相关函数。d(C(t)*W(t))/dt = F1-F2
F1和F2是我想传递给函数的常量。为了求解常微分方程，我使用了欧拉的方法，但为了比较结果，我想使用ODE45或适当的常微分方程求解器。Euler方法中的错误值可能导致错误信息。T = 0.0125; %Duration
浏览 1提问于2014-11-20得票数 1
我在excel电子表格中得到了这些数据。所以导入之后，我只需要做velocity = cumtrapz(t,y)，然后做position = cumtrapz(velocity)
浏览 0提问于2014-04-25得票数 0
回答已采纳
我的问题与有关，它解释了如何用VexCL实现boost::odeint，以实现并发(完整的代码可以找到)。下图显示了我如何看待ODEINT的迭代：
现在我扪心自问，到底是什么/或者它的哪一部分是VexCL中的并行化的？我的印象是，ODE部分是一个单一的任务，因为在给定的示例中，所有ODE方程都在一个块内。也许集成部分运行在三个并行任务中。这导致了四个任务，其中(我认为) ODE任务是一个瓶颈(因为方程可能变得非常大)。
如果这是正确的，我想知道，如何改进这个并发。我认为横
浏览 2提问于2020-02-11得票数 0
回答已采纳
2回答
我正在对一台桥式起重机进行建模，并获得了以下方程：
当涉及到Scilab时，我是新手，到目前为止，我只模拟了(使用ODE)不超过两个自由度的线性系统，这是一些简单的系统，我可以很容易地将其转换为am但是这个系统我不知道如何模拟它，不是因为sin和cos函数，而是因为我不知道如何把它放在状态空间矩阵中。
浏览 0提问于2015-11-16得票数 1
我花了一周多的时间来解决问题：我必须使用牛顿-拉夫森方法来得到u的近似解。我有这样做的脚本，但我需要“线性化”这个非线性常微分方程。k1-k4不是常量。
浏览 1提问于2013-03-21得票数 0
+1, retstep=True);CFL = 0.3
for k in range(m): k1 = derivatives(t,u,c,dx)*dt; k3 =
浏览 79提问于2020-01-05得票数 1
回答已采纳
如果我必须在modelica中使用偏导数，该如何使用。我不确定偏导数是否可以在modelica中求解，但我想知道，如果它可以使用，它应该如何实现。
浏览 0提问于2012-08-07得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
