![](/img/trans.png)
[英]How can I integrate the Ordinary differential equation using odeint from scipy.integrate?
[英]Solution for differential equation is completly different in boost::odeint and scipy.integrate
我正在尝试将快速原型从python移植到C ++。 我尝试使用一个简单的微分方程来测试表示法,但是对于起始值[2,0],结果却大不相同。 Python正在下降,而C ++解决方案正在迅速上升。
它适用于以下示例: 如何将查找表中的时变参数合并到boost :: odeint,c ++
但这对我的例子不起作用
TransferF::TransferF(const double& deltaT) : dt(deltaT), t(0.0), y(2)
{
// initial values
y[0] = 2.0; // x1
y[1] = 0.0; // x2
}
void TransferF::ode(const state_type &x, state_type &y, double t)
{
y[0] = x[0];
y[1] = x[1];
y[2] = (-2*y[1] - y[0] + 1) / (pow(y[0],2));
}
和在py中相同:
def modelt(x,t):
y = x[0]
dydt = x[1]
dy2dt2 = (-2*dydt - y + 1)/ (y **2)
return [dydt,dy2dt2]
x3 = odeint(modelt,[2,0],timev)
我希望时间序列的结果相同,但是pythons解决方案正在下降,C ++正在上升。
C ++代码存在细微的不一致。 输出向量y
应该只包含导数y', y"
,而不能包含函数y
本身:
void TransferF::ode(const state_type &x, state_type &y, double t)
{
y[0] = x[1];
y[1] = (-2*x[1] - x[0] + 1) / (pow(x[0],2));
}
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.