[求助] 请高手帮忙： 在相同时间内，第一个ODE45的解如何传递给第二个ODE45作为输

dongshanyan

关于ODE45求解微分方程组的问题：（问题见附件）
　　在相同的时间内，第一个ODE45的解是第二个ODE45的输入。我不知如何处理第一个ODE45的解与时间的关系，并将其作为输入传递给第二个ODE45。　　　我的QQ：10181694，手机：13636495360。
　　希望得到大家的帮助，在此先谢了！

caiqd

不知你的这个问题是否已经圆满解决，我觉得，不应该分别建立方程，而应该对整个系统写出一组常微分方程组，基本的变量为[S1,X1,S2,X2],方程为
dS1/dt=Qin*(Sin-S1)/V-32*S1/(20+S1)*X1
dX1/dt=Qin*(Xin-X1)/V+16*S1/(20+S1)*X1-0.4*X1
dS2/dt=Qin*(S1-S2)/V-32*S2/(20+S2)*X2
dX2/dt=Qin*(X1-X2)/V+16*S2/(20+S2)*X2-0.4*X2
初始条件就是[1, 2000，0.1, 1500]
对这个系统，用任何一种RK方法，都可以简单地求解。

dongshanyan

　To caiqd:谢谢你给的方法！
　对于你给的方法，我有一点疑问。ODE45求解时，是否先计算S1、X1的值，然后再将S1和X1的值传递到下面的两个式子。因为在我的系统中，必须体现出第一个反应器的出水是第二个反应器的进水，即第一组解为第二组解的输入。请问MATLAB是按照这个顺序来求解的吗？如果是的话，我想那肯定是正确的。

caiqd

其实你不用特别强调先计算S1和X1这种问题，因为两个反应器其实可以认为是同步的，当写为一个方程组时，顺序不是重要的。