function fit_nonl
clear all;clc
format long

data=[    8.6667    2.0000
   10.0000    5.0000
   10.6667   10.0000
   12.6667   15.0000
   13.3333   30.0000
   15.0000   60.0000
   15.3333  120.0000
   16.6667  200.0000
   17.3333  300.0000
   18.3333  550.0000
   20.0000  700.0000
];
xspan=data(:,2);    %x的数据，在此输入
texp=data(:,1);    %t的数据，在此输入

k0=[0.1 0.1 1 100];
lb=-[1 1 1 1]*1e9;
ub=[1 1 1 1]*1e9;


%-------------------------------------------------------------------------

% 使用函数lsqnonlin()进行参数估计

options=optimset('maxfunevals',1000);
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@objfunc,k0,lb,ub,options,xspan,texp);

ci = nlparci(k,residual,jacobian);
%residual;
fprintf('\n\n拟合结果:\n')
fprintf('\n\t参数 e1 = %.16f',k(1))
fprintf('\n\t参数 e2 = %.16f',k(2))
fprintf('\n\t参数 e3 = %.16f',k(3))
fprintf('\n\t参数 t = %.16f',k(4))
y=kineticseqs(xspan,k);
r2=1-sum((texp-y).^2)./sum((texp-mean(y)).^2);
fprintf('\n\tr^2 = %.16f',r2);
figure
plot(xspan,kineticseqs(xspan,k),'b',xspan,texp,'or'),legend('计算值','实验值','location','best')


%-------------------------------------------------------------------------

function f = objfunc(k,xspan,texp)
f=kineticseqs(xspan,k)-texp;

%------------------------------------------------------------------------
function xt = kineticseqs(x,k)
a=2;
e1=k(1);e2=k(2);e3=k(3);t=k(4);
xt=a/e1 +(a/e2)*(1-exp(-x/t))+a*x/e3;

16楼: Originally posted by 月只蓝 at 2015-01-21 09:01:30
用1stopt软件做比较方便，该软件参见：http://muchong.com/bbs/viewthread.php?tid=4586018&fpage=1&target=blank

求解该问题的代码：
Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Variable X1,X2,X3,Y ...

4楼: Originally posted by feixiaolin at 2013-11-05 11:56:53
你的 “+a/E2” 再明确一下为好。

2楼: Originally posted by jerkwin at 2013-11-05 11:23:25
不劳matlab, origin, 1stopt足矣

2013-11-04_212307.png

5楼: Originally posted by dshuanger at 2013-11-05 12:32:47
对不起呀，应该是： y=a/E1 +(a/E2)*+a*x/E3...

7楼: Originally posted by 月只蓝 at 2013-11-05 14:09:49
（1） MATLAB代码及结果

function fit_nonl
clear all;clc
format long

data=;
xspan=data(:,2);    %x的数据，在此输入
Texp=data(:,1);    %T的数据，在此输入

k0=;
lb=-*1e9;
ub=*1e9;


%--- ...