matlab 估计动力学参数

小管

我现在做的是一个可逆的水解反应，方程式为A+B生成C+D，故其动力学方程应为：
-dCA/dt=k1CACB-k2CCCD;
-dCB/dt=dCC/dt=dCD/dt=-dCA/dt；
为了提高反应的转化率，在反应过程中通过精馏将反应生成的产物C及时移走，但是在移走C的过程中，B也有一部分被移除，所以要对上述方程式进行修改：
dCB/dt=k1CACB-k2CCCD-s/v*CC    （1）
dCC/dt=k1CACB-k2CCCD-s/v*CD      （2）
其余两个的不变，又因为在反应过程中反应液的体积v是不断变化的，所以s/v也是一个变量将该变量设为CE，并将CE代入方程（1）和（2）中，就可以得到一个新的动力学方程式，




问题是matlab 程序如何编写，有木有可逆均相反应的动力学拟合matlab 程序参考下，求助各位马友鼎力帮助

                               
登录/注册后可看大图

该贴已经同步到 小管的微博