求MATLAB二元指数拟合，在线等答案。。。好的话再加分，非常感谢！

Question

cs=10;c0=20；ux=4;L0=20;t=[0,8,28,36,56]';y=[10,8.5,7.0,6.1,7.2]';
y=cs-(cs-c0)*exp(a*t/-ux)+d*L0/(a-d)*(exp(a*t/-ux)-exp(d*t/-ux)），求a和d，c0的值
呃，改一下条件，c0=10，然后最好有具体程序，顺便每句帮解释下，还是万分感谢！

Answer 1

你先把这段code保存为testfun.m，是一个函数。

function y=testfun(x,xdata)
n=size(xdata,2);
y1=zeros(1,n);
cs=10;c0=20;ux=4;L0=20;
for i=1:n
y1(i)=cs-(cs-c0)*exp(x(1)*xdata(i)/-ux)+x(2)*L0/(x(1)-x(2))*(exp(x(1)*xdata(i)/-ux)-exp(x(2)*xdata(i)/-ux));
end
y=y1;
end

然后运行下面这段code，也是在editor里面写好再运行。

x0=[0.01,0.03];
xdata=[0,8,28,36,56];
ydata=[10,8.5,7.0,6.1,7.2];

xf=lsqcurvefit(@testfun,x0,xdata,ydata);

x=0:.1:56;
y=testfun(xf,x);

plot(xdata,ydata,'ro',x,y,'k');

xf(1)是拟合后的a, xf(2)是拟合后的d。

Answer 2

a=0.24847477145119115093940820315765

d=0.30938097165974515201369854366931

c0 不是已知了吗？