作者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