Matlab 用一系列微分方程描述的模型拟合实验数据

我有一系列的实验动力学数据，其形式为浓度=ftime，即：

时间=时间 时间=时间 Btime=f2time Ctime=F3时间 Dtime=f4time Ptime=f5time Atime+…+PTime=已知值 我有一个表示为微分方程组的模型：

dL/dt = -k_tot*MR3*(L)^2, where k_tot = k_LA + k_LB + k_LC + k_LD
dMR3/dt = -k_tot*MR3*L
dP/dt = k_AP*A + k_BP*B + k_CP*C + k_DP*D
dA/dt = 2*k_LA*MR3*(L)^2 - k_AP*A
dB/dt = 2*k_LB*MR3*(L)^2 - k_BP*B
dC/dt = 2*k_LC*MR3*(L)^2 - k_CP*C
dD/dt = 2*k_LA*MR3*(L)^2 - k_DP*D

我想使用ODE模型ode45，以便：

查找/估计k_XY k_LA、k_LB、…、k_AP等 计算数据与实验数据拟合曲线。 其思想是建立基于ODE的函数，该函数允许在给定的方程组中拟合指定微分方程的指定数据

有人能帮我看一下Matlab代码的例子吗

---

非常感谢

这可以表述为最优控制问题，目标函数是ODE解到数据的距离，ODE是约束，参数是控制变量。其主要思想是向前求解ODE，保留该计算的轨迹，并累积目标函数wrt的梯度。跟踪后的参数向后。然后将其插入梯度下降方案，如CG。非常感谢！然而，我需要几个月的时间来编写Matlab代码并调试它。。。