发布网友 发布时间:2022-11-19 15:47
共1个回答
热心网友 时间:2024-10-27 11:22
复杂反应体系的动力学方程参数实际上是反常微分方程组系数拟合问题。可以按下列步骤进行:
1、根据动力学方程,定义关于微分方程组的自定义函数。
2、构造其数据误差目标函数为最小,即
min F(ki)=Σ[S(ti)-Si]²
3、然后最小二乘法求出系数ki
4、ode45函数求出x1,x2,x3,。。。拟合值
5、计算实验值与拟合值的相关系数R²是否接近于1
通过编程计算,得到如下结果:
k1=0.17701;k2=0.013736;k3=0.013395;k4=0.078272;k5=0.0098657;k6=0.15358;k7=0.45656; k8=0.00072802;k9=0.076499;k10=0.27148;k11=0.026889;k12=0.010795;k13=0.27413; k14=0.055203;k_1=0.16205;