matlab 微分方程组参数拟合,常微分方程组参数拟合问题

EDA365欢迎您登录！

您需要 登录 才可以下载或查看，没有帐号？注册

x

代码如下;

' e9 M; r1 r: Ufunction Kinetics4: X% H  p' J- G8 r

clear all

4 ?& z% ?' ^! E2 d- zclc' h% W: G2 G. R* \

k0 = [2.159*10^8 4.734*10^6 174686 149892 0.99 0.7 0.5 0.78];%参数初值，其实是不知道的  ( W; z+ S& H7 `3 z! k

lb = [-inf -inf -inf -inf -inf -inf -inf -inf];                   % 参数下限

% ?- }+ Z$ {- \ub = [inf inf inf inf inf inf inf inf];    % 参数上限0 G9 t' h& @2 j- g5 P

u0 = [0,0]; %y初值

3 y$ W! `# c5 r; j5 R$ C4 a. da=[0;0.06680583;0.192617534;0.301693824;0.419783943;0.527041818;0.598345407;0.65055108;0.6876854;0;0.0295189820000000;0.0760025420000000;0.143652564000000;0.303435528000000;0.419715167000000;0.525496977000000;0.565069320000000;0.743536298000000;0;0.0411725930000000;0.0608472370000000;0.0986899250000000;0.241132264000000;0.324760943000000;0.466871464000000;0.531688237000000;0.683614442000000];. A3 D9 S7 F  Q8 M5 @% d

b=[0;0.2378;0.3288;0.3556;0.39;0.3877;0.3766;0.3580;0.3384;0;0.1373;0.2123;0.2755;0.3561;0.3549;0.3475;0.3486;0.2518;0;0.115879961000000;0.178182567000000;0.219982937000000;0.315058656000000;0.339244190000000;0.348664305000000;0.351963292000000;0.313743618000000];0 y4 D' f9 L! P8 a5 v) d

yexp=[a,b];3 h7 Y: V$ U* p& p) X4 m

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

[& M" y+ ?: {

: \9 L% _3 [! w" k! D[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...% h8 \7 v/ W2 {- a& n

lsqnonlin(@ObjFunc4LNL,k0,lb,ub,[],u0,yexp);

<