j**u
发帖数: 6059 | 1 公式很简单
dM1/dt = -qM1-c1M1+p1c2M2
dM2/dt = -qM2-c2M2+c1M1+c3M3
dM3/dt = -qM3-c3M3+p2c2M2
initial condition: c1M1+c2M2=c3M3, p1+p2=1
S(q)=M1+M2+M3
现在已有数据系列数据S(q)，请问如何做data fitting 得到参数c1,c2,c3和p1,p2? 本
来打算用lsqcurvfit和ode45，但是总是搞不定，大家有什么好办法吗？ | j**u
发帖数: 6059 | 2 仔细看了一下，原来基本上是因为对初值太敏感。
【在 j**u 的大作中提到】
: 公式很简单
: dM1/dt = -qM1-c1M1+p1c2M2
: dM2/dt = -qM2-c2M2+c1M1+c3M3
: dM3/dt = -qM3-c3M3+p2c2M2
: initial condition: c1M1+c2M2=c3M3, p1+p2=1
: S(q)=M1+M2+M3
: 现在已有数据系列数据S(q)，请问如何做data fitting 得到参数c1,c2,c3和p1,p2? 本
: 来打算用lsqcurvfit和ode45，但是总是搞不定，大家有什么好办法吗？
|
|
