jiangkingdom
2009-09-25, 14:54
M文件
function yy=model(beta0,x)
a=beta0(1);
b=beta0(2);
c=beta0(3);
yy=a*b*exp(c*x)/(b-a+a*exp(c*x))
程序
x =
0 3 6 9 12 15 18 21 24 27 30 33 36 39 42 45
beta0 =
4.4300 82.0000
y= 4.4300 8.4000 13.3100 23.1600 41.6300 54.7800 66.1200 75.7900 77.9500 77.8100 80.5100 81.4600 81.9000 82.0500 80.1900 79.6700
betafit = nlinfit(x,y,'model',beta0)
如何求C也就是beta(3）
希望详细点，我对这方面很。。。
感谢了。
function yy=model(beta0,x)
a=beta0(1);
b=beta0(2);
c=beta0(3);
yy=a*b*exp(c*x)/(b-a+a*exp(c*x))
程序
x =
0 3 6 9 12 15 18 21 24 27 30 33 36 39 42 45
beta0 =
4.4300 82.0000
y= 4.4300 8.4000 13.3100 23.1600 41.6300 54.7800 66.1200 75.7900 77.9500 77.8100 80.5100 81.4600 81.9000 82.0500 80.1900 79.6700
betafit = nlinfit(x,y,'model',beta0)
如何求C也就是beta(3）
希望详细点，我对这方面很。。。
感谢了。