我有一组如下形式的微分方程:
x1dot = x3;
x2dot = x2;
x3dot = x1;
x4dot = x2 + integral(x1,t,tend)我有x1,x2在tstart和x3,x4在tend的边界条件。如果没有积分项,它就是一个使用BVP4C的简单实现。
我想知道是否有可能从BVP求解器中得到以前的状态解,它可以用于积分。
一种可能是将ode45和fsolve结合使用来解决边值问题,我可以使用以前的解决方案,但这种方法不如BVP设置那么快。
我还觉得,当我使用以前的积分解决方案x1时,在收敛方面可能会有一些困难。
有没有更好/更快或更简单的方法来解决这个问题?
https://stackoverflow.com/questions/41228805
复制相似问题