我正在尝试使用ODEINT(它是Boost软件包之一)以矩阵形式求解一组线性常微分方程(ODE)。
我要解决的问题具有以下形式(以MATLAB格式):
dxdt(1:10) = A1*x(1:10) + B1;
dxdt(11:20) = A2*x(11:20) + B2;
其中A1和A2是常数矩阵,而B1和B2是常数向量; “ x”是要求解的状态,“ dxdt”是状态“ x”的时间导数。我的问题是如何使用ODEINT在c ++中对其进行编程。我尝试了一些在互联网上找到的东西,但是找不到任何合适的方法。任何帮助将不胜感激。