matlab动力学参数,用matlab求动力学参数

本人MATLAB小白，最近审稿人要求做个动力学模型，我不会用MATLAB拟合求速率常数，希望木虫的大神可以教我，谢谢。

(dC_A)/dt=-k1*CA-k7*CA

(dC_B)/dt=k1*CA-k2*CB- k5*CB

(dC_C)/dt=k2*CB-k3*CC- k6*CC

(dC_D)/dt=k3*CC-k4*CD

(dC_E)/dt= k4*CD

(dC_F)/dt=k5*CB+k6*CC+ k7*CA

此外还采用Nelder-Mead 方法通过最小化目标函数优化速率常数

F_obj=∑_m▒{∑_(n=1)^(n_p)▒[C_(m,n,calc)-C_(m,n,exp) ] ^2 }

C_(m,n,calc)是给定的时间n下组分的实验浓度；C_(m,n,exp)是给定的时间n下组分的计算浓度

第一组：

T= (0.5, 1, 1.5, 2, 3, 4, 5, 6)

A= (18.7, 6.5, 0.7, 2.1, 1.7, 1.5, 1.2, 0.2)

B= (48.9, 31.3, 15.3, 6.9, 2.2, 1.8, 1.7, 0.5)

C= (19.3, 34.7, 41.2, 40.9, 37.2, 31.7, 26.1, 19.4)

D= (9.7, 19.6, 28.7, 34.5, 38.4, 38.2, 36.3, 33.4)

E= (0.2, 2.6, 5.7, 9.3, 17.1, 24.8, 30.5, 36.7)

第二组：

T= (0.5, 1, 1.5, 2, 3, 4, 5, 6)

A&#