工程数学作业第十四次方健(12.29)(2)

plot(CAc,resid,'*')

xlabel('浓度拟合值（mol/L）') ylabel('残差R (mol/L)') refline(0,0)

% 参数辨识结果

fprintf('Estimated Parameters:\\n')

fprintf('\\tk = %.4f ± %.4f\\n',beta(1),ci(1,2)-beta(1)) fprintf('\\tn = %.2f ± %.2f\\n',beta(2),ci(2,2)-beta(2)) % ------------------------------------------------------------------ function f = OptObjFunc(beta,tspan,CA0,CAm) [tCAc] = ode45(@KineticsEqs,tspan,CA0,[],beta); f = CAc - CAm;

% ------------------------------------------------------------------ functiondCAdt = KineticsEqs(t,CA,beta)

dCAdt = -beta(1)*CA^beta(2); % k = beta(1), n = beta(2) 运行结果如下： beta =

0.0047 1.4555 resnorm = 0.0941 resid =

0 -0.2408 0.1813 0.0303 -0.0169 -0.0435 0.0131 exitflag = 3 output =

firstorderopt: 0.0041 iterations: 3 funcCount: 12 cgiterations: 0

algorithm: 'trust-region-reflective' message: [1x457 char] lambda =

lower: [2x1 double] upper: [2x1 double] jacobian =

(2,1) -394.7680 (3,1) -566.6235 (4,1) -629.8164

(5,1) -588.7875 (6,1) -478.0379 (7,1) -305.6976 (2,2) -4.0422 (3,2) -5.4804 (4,2) -5.7590 (5,2) -4.5555 (6,2) -3.1228 (7,2) -1.3857 ci =

0.0031 0.0063 1.2700 1.6410 Estimated Parameters: k = 0.0047 ± 0.0016 n = 1.46 ± 0.19

10987浓度CA, mol/L ExpModel654321 050100150时间t, s200250300

0.20.150.10.05残差R (mol/L)0-0.05-0.1-0.15-0.2-0.251234567浓度拟合值（mol/L）8910下面将lsqnonlin替换成lsqcurvefit函数观察结果变化：

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库工程数学作业第十四次方健(12.29)(2)在线全文阅读。