作者ghost17612 (就是要ROCK)
看板MATLAB
标题[讨论] nlinfit 遇到 ode
时间Mon Apr 15 21:55:33 2013
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)
1F：推 Raymond0710:换字形 04/16 01:59
2F：→ roxxot:nlinfit 的第三个input写错，当然算出不来.... 04/19 17:59