[討論] nlinfit 遇到 ode
function nonlinefit
clc;clear;
t=[0 8.33 16.67 58.33 116.67 241.67 416.667];
y=[0 0.05 0.0571 0.1 0.171 0.229 0.321];
beta=nlinfit(t,y,@ode,[1 1])
a=beta(1)
b=beta(2)
plot(t,y,'o')
[t,y] =ode45(@ode,t,0)
function y1=ode(t,y)
a=const(1);
b=const(2);
y1=10.581*((((1-y)./(3.76+y)).^2)*a
-(y./(3.76+y).*((1.76+2*y)./(3.76+y)).^2)*b);
這是我的code
我上網也查不到相關的資料,不然就是寫得不夠清楚
想請問版上的各位大大
如何用nlinfit去fit出 ode的參數呢?
(例如 d(Ca)/d(t)=kCa,我們知道一些(t,Ca)的數據點,然後用nlinfit去找出參數k)
感謝!!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.37.104.70
※ 編輯: ghost17612 來自: 114.37.104.70 (04/15 21:57)
推
04/16 01:59, , 1F
04/16 01:59, 1F
→
04/19 17:59, , 2F
04/19 17:59, 2F