作者obason23 (累了)
看板MATLAB
标题[运算] 利用ODE计算反应常数 (Michaelis-Menten)
时间Sat Jan 14 06:54:39 2012
我卡关卡好久了 盼望有好心人士协助...
我现在有几个 (t, S) 资料点
想要利用 Matlab 作非线性回归 并求出反应参数 km 和 Ks
Paper 上只有给微分方程 (其实他就是加了负号的 Michaelis-Menten)
dS/dt=-km*S/(Ks+S)
我试着用Matlab 解成 S(t) 形式, 可是
S = exp((C3 - km*t)/Ks - wrightOmega(log(1/Ks) + (C3 - km*t)/Ks))
^^^^^^^^^^^ 这是个 Lambert W-function
然後我就不知道该怎麽办了....
我也有试过在Matlab Central 里前人写过的 .m file (用ODE解)
唯一成功的是一个把所有 parameter 都定死的
其他的file 都不 work...
网路上其他的方法都是先换成线性再求解
我的指导老师希望我用 Matlab 作出非线性回归 然後求参数
有请大家帮忙了 谢谢!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 155.99.214.1
1F:推 gamer:先用数值微分建立dS/dt DATA,再直接最非线性回归就行了 01/14 11:55