我正在研究使用最优控制方法(轨迹生成和优化)的无限维优化算法。我想要做这个轨迹优化的系统是非线性的,形式是$\dot{x}(t) = f(x(t),u(t),t)$。换句话说,我有时变的非线性输入。
是否有可能使用boost::odeint解决这样的ODE问题?我在文档中没有发现任何提示,但我可能没有看到它。
发布于 2014-04-09 21:32:49
是的,你可以用odeint来解决这类问题。显式步进者期望传递给odeint的系统函数( ODE)具有签名。
ode( x , dxdt , t );其中x是当前状态的输入参数,dxdt是r.h.s的输出参数。对于ODE,t是时候了。例如,驱动振荡器的实现方式如下
typedef std::array< double , 2 > state_type;
struct oscillator
{
double driving_strength;
double dribving_frequency;
void operator()( state_type const &x , state_type &dxdt , double t ) const
{
dxdt[0] = x[1];
dxdt[1] = -x[0] + driving_strength * sin( driving_frequency * t );
}
};
state_type x;
oscillator osc;
// initialize x and osc
runge_kutta4< state_type > stepper;
integrate_const( stepper , osc , x , t_start , t_end , dt );https://stackoverflow.com/questions/22973864
复制相似问题