matlab 拟合求函数参数问题

Question

已知：
x=2004:1:2008;
y=[68 78 90 105 122];
函数y=k/(1+exp(a+b*x));
求k,a,b
并作图
matlab程序

Answer 1

function f=fun1(x,tdata);
f=x(1)+(1+exp(x(2)+x(3)*tdata)); %其中x(1)=k，x(2)=a，x(3)=b;

（2）调用函数lsqcurvefit，编写程序如下：
td=2004:1:2008;
cd=[68 78 90 105 122];
x0=[0.2 0.05 0.05];
x=lsqcurvefit(@fun1,x0,td,cd)
plot(td,cd,'o-')