1回答
我正在学习如何在python上编写代码，并且我正在试图弄清楚如何在特定时间从ODE系统中找到解决方案的总和。例如，这是来自SciPy Cookbook的例子被称为“为僵尸启示录建模”https://scipy-cookbook.readthedocs.io/items/Zombie_Apocalypse_ODEINT.html以下是该网站的部分代码： # zombie apoc
浏览 9提问于2020-03-27得票数 1
回答已采纳
1回答
我正在尝试解决ODE的dF/dt = A*F, F_initial = eye(9)的一个特定系统。作为Matlab的新手，我正在尝试使用实现的ode45函数，我在网上找到了一些有用的建议。然而，它们都假设A是常数，而在我的例子中，矩阵A是t的函数，换句话说，A随着时间步长而变化。我已经分别求解了A并将其存储在一个9x9xN数组中(因为我的</e
浏览 8提问于2018-08-31得票数 1
回答已采纳
我有一个系统看起来dh/dt=g(h,v)
我想在流形F(v,n,h)=0上解这个方程，这是v中的一个非线性函数。我试图在每个时间步骤中使用类似于v=fzero(@(x) F(x,n,h),0)的方法来求解流形上v的值。但这是难以置信的慢和ode15s (我的系统是一个松弛振荡器)不能满足集成的容限。如何在F(v,n,h)=0定义<e
浏览 6提问于2016-10-03得票数 5
回答已采纳
如果可能的话，我需要您的建议来正确地并行实现以下场景。
假设我有一个用于解决某些问题的ODE系统。我们可以用ode(...)来解决上述问题。deSolve包的。我的问题从这里开始。我想用20个不同的初始条件(即20个这样的系统并行运行)来模拟上面的系统，这样在每个时间步(dt=0.001)之后，我需要从当时第i个系统的解中减去第(i+1)和(i-
浏览 3提问于2018-06-19得票数 0
1回答
如何计算曲线到达特定x坐标所需的时间(在Matlab中)。假设我们有: dx/dt = x^2 + y^2和dy/dt = 5.x.y，曲线从点(a，b)开始。在ode45的帮助下，我能够得到曲线的图形。我需要计算曲线到达x= c，(c>a)所需的时间。有人告诉我，这可以通过插值来完成，但我不知道如何编写代码。
浏览 1提问于2013-12-07得票数 0
我正在尝试求解一个常微分方程组，我的输入激励是时间的函数。
我一直在集成函数中使用interp1，但这似乎不是一种非常有效的方法。我知道它不是，因为一旦我将输入激励更改为sin函数，这不需要在函数内部调用interp1，我就会得到快得多的结果。但是每一步做插值大约需要10-20倍的时间才能收敛。那么，有没有一种更好的方法来求解任意时间相关激励的常微分方程，而不需要做插值或其他一些技
浏览 1提问于2016-04-03得票数 0
我有一个带有微分方程的函数文件，我正在对函数执行一个标准形式的ode23s，即t，m=ode23s('DE_ function '，tspan，mA pA mB pB mC pC mD，optionsDE我为每个mA获得大约150个值/结果/输出，依此类推。我的ode23s很好用。我有相同的mA的实验数据集，所以我必须用它来计算最小平方误差。在我的实验数据中，根据20个
浏览 1提问于2016-05-17得票数 2
回答已采纳
这是我在这里的第一个问题，所以请对我宽宏大量。我想知道是否有一种方法可以集成ODE系统，直到找到指定变量的局部最大值。以下是更多细节：
让我们称我们的ODE系统为dX/dt = F(X) where X(t) = [x1(t), x2(t), ... , xn(t)]。进展:使用scipy.integrate.odeint我能够找到极限环C和它的周期T的解，但是从任意点开
浏览 0提问于2015-06-02得票数 4
与时间相关的哈密顿量的薛定谔方程是：我试图实现一个求解薛定谔方程的时间相关哈密顿量在ode45。然而，由于哈密顿$H(t)$依赖于时间。我不知道如何在ode45中进行插值。psi0 = [0 1];
H = [1 0;0 1]*cos(t); %this is wrong, I do not know how to implement this and pass it to ode45h
浏览 1提问于2016-06-23得票数 0
回答已采纳
1回答
在没有敏感性的情况下解决问题效果很好，但现在我想在finaltime、初始值和控件(我的等式中使用的参数)中使用敏感性。我只在assimulo页面上找到了一个例子(这只是初始敏感度)，我没有正确地理解它。非常感谢
浏览 1提问于2013-04-04得票数 1
我有一个常微分方程系统，我想要解决，但有一个棘手的部分，当系统达到稳定状态时，我想改变一个(或多个)参数的值。params(2); dydt(1) = F*x(1) - G*x(1)*x(2);end
我希望我的代码能够工作，例如，当系统达到稳定状态时，F的值更改为10，G的值更改为2。我正在考虑使用以下命令来检测dydt(1)和dy
浏览 0提问于2016-04-17得票数 2
1回答
我试着用ODEINT来解一个一阶的ODE。下面是代码。我希望得到3个时间点的y值。我正在努力解决的问题是，能否传递mt和nt的n个值来计算dydt。我认为ODEINT通过了mt和nt的所有3个值，而只是第0、1或2，这取决于迭代。该系统类似于一个双罐系统，第一个坦克以未知的速度接收特定的激素，但是我们的传感器将
浏览 1提问于2017-01-25得票数 3
1回答
我很难解出光学bloch方程，它是一个具有复值的一阶ODE系统。我发现scipy可以解决这样的系统，但是他们的网页提供的信息太少，我很难理解。我有8个耦合的一阶ODE，我应该生成一个函数，类似于： compute the time dervative of elements in y那就做complex_ode</
浏览 4提问于2012-02-25得票数 3
回答已采纳
1回答
当因变量超过特定值时，我想终止ODE求解器。考虑一个简单的ODE模型dModel<- function(t, y, parms) {list(c(y))}
out<-ode(c(1),seq(0,100,1我知道有一个ODE事件，但在阅读了文档后，我不知道如何在我的情况下实现。你能提供这个问题的解决方案吗？
浏览 1提问于2018-07-10得票数 0
我想强制ode45采取恒定的步长，所以在求解方程时，它总是在T轴上递增0.01。我该怎么做呢？%clf;tspan=[to tf];
[T,Y]=ode45
浏览 0提问于2010-03-31得票数 6
回答已采纳
我想在数值格式的每一步提取状态变量的值。
deSolve R包中的ode()函数，使用实现的ODE求解器之一对一个常微分方程组进行数值求解。为了做到这一点，它使用了一个动态调整的集成步骤，该步骤基于每个集成末尾的本地截断错误的值。用户基本上指定了所需的输出时间，它定义了时间步长网格，它可以有一个等于或小于输出时间的</e
浏览 2提问于2019-03-31得票数 2
回答已采纳
我改变了我的模型中的参数，但有时它会导致模拟停止，因为病态雅可比，我不知道这个错误意味着什么。
浏览 3提问于2020-07-12得票数 4
回答已采纳
我有一个用ode15s解决的僵硬的颂歌系统。时间导数的类型为： yt = A*y - u.*y; 其中yt、y和u是列向量，A是稀疏矩阵。u的计算必须在每次求值时进行，而且计算量很大(u是非线性的，这是刚度的原因)。我不想预先计算u，因为我想保留求解器在以后选择任何所需时间的自由。然而，我感兴趣的不是解y(t)，而是已经计算过的计算yt的u(t)
浏览 70提问于2020-03-10得票数 0
我正在使用刚性求解器(ode15s)对常微分方程系统进行时间积分。它是有效的，但我想加快速度。方程系统以状态空间的形式给出： % Mx'' + Cx' + Kx = F(t)
% BUNCH我认为(除了ϕ=0之外)的最佳初始条件是在ode求解器(ϕ(k)_initial=ϕ(k-1))的前一次迭
浏览 0提问于2019-04-08得票数 0
3回答
我在Octave中使用OdePkg来解决一个僵硬的系统，例如通过ode5r： yprime = [ param(1) * Y(1) - param(2) Y(2) * Y(3); # ODE(@myODEs, time_span, Y0, ..., param);
求解器将因变量Y存储在
浏览 4提问于2009-07-22得票数 0
回答已采纳
云服务器
ICP备案
云直播
对象存储
实时音视频
腾讯云开发者
