24小时热门版块排行榜

返回列表

_romantic_镇

木虫 (著名写手)

应助: 2 (幼儿园)
金币: 3339
散金: 387
红花: 6
帖子: 1221
在线: 677小时
虫号: 3186948
注册: 2014-05-07
专业: 化学工程及工业化学

[求助] 初学者求助，matlab模仿程序已有1人参与

function KineticsEst5
% 动力学ODE方程模型的参数估计

clear all
clc

k0 = [0  0  0  0];       % 参数初值
lb = [ -inf  -inf  -inf  -inf ];                % 参数下限
ub = [+inf  +inf  +inf  +inf ]; % 参数上限
x0 = [8.5 28.8 27.6];
tspan=[0 0.222 0.333 0.444];
yexp =  [8.5000 28.8000 27.6000
4.8000 23.2000 35.3000
4.2000 21.6000 36.5000
4.0000 21.2000 37.3000];                % yexp: 实验数据[x1 x2 x3]

% 使用函数fmincon()进行参数估计
[k,fval,flag] = fmincon(@ObjFunc4Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp);
fprintf('\n使用函数fmincon()估计得到的参数值为:\n')
fprintf('\tk1 = %.4f\n',k(1))
fprintf('\tk2 = %.4f\n',k(2))
fprintf('\tk3 = %.4f\n',k(3))
fprintf('\tk4 = %.4f\n',k(4))
fprintf('  The sum of the squares is: %.1e\n\n',fval)
k_fmincon = k;

% 使用函数lsqnonlin()进行参数估计
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc4LNL,k0,lb,ub,[],x0,yexp);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n')
Output

% 以函数fmincon()估计得到的结果为初值，使用函数lsqnonlin()进行参数估计
k0 = k_fmincon;
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc4LNL,k0,lb,ub,[],x0,yexp);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n以fmincon()的结果为初值，使用函数lsqnonlin()估计得到的参数值为:\n')
Output

% ------------------------------------------------------------------
function f = ObjFunc4Fmincon(k,x0,yexp)
tspan = [0.00 : 0.01 : 0.44];
[t x] = ode45(@KineticEqs,tspan,x0,[],k);
y(:,1:3) = x(:,1:3);
f = sum((y(:,1)-yexp(:,1)).^2) + sum((y(:,2)-yexp(:,2)).^2) ...
+ sum((y(:,3)-yexp(:,3)).^2) ;

% 以函数fmincon()估计得到的结果为初值，使用函数lsqnonlin()进行参数估计
k0 = k_fmincon;
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc4LNL,k0,lb,ub,[],x0,yexp);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n以fmincon()的结果为初值，使用函数lsqnonlin()估计得到的参数值为:\n')
Output

% ------------------------------------------------------------------
function f = ObjFunc4LNL(k,x0,yexp)
tspan = [0.00 : 0.01 : 0.44];
[t x] = ode45(@KineticEqs,tspan,x0,[],k);
y(:,1:3) = x(:,1:3);
f1 = y(:,1) - yexp(:,1);
f2 = y(:,2) - yexp(:,2);
f3 = y(:,3) - yexp(:,3);
f = [f1; f2; f3];

% ------------------------------------------------------------------
function dxdt = KineticEqs(t,x,k)
dxdt =  ...
[ (k(1)*x(2)- k(1)*x(1))
(k(1)*x(2)+(k(3)-k(2)-k(4))*x(3))
(k(3)*x(2)-k(4)*x(3))
];

回复此楼

» 猜你喜欢

深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组已经有18人回复
论物质与能量的统一模型及物理现象解释已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有95人回复
基于基元I统一理论的数学相关应用推导已经有0人回复
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律已经有1人回复
基元I理论下三大核心空间现象精准推导与细节解析已经有0人回复
基于基元 I 统一理论的反重力理论推导已经有0人回复
基于基元I统一理论的量子力学本源推导已经有1人回复
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费已经有18人回复
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）已经有0人回复
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）已经有0人回复

1楼 2021-12-20 20:54:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

应助: 387 (硕士)
金币: 17803.5
红花: 53
帖子: 2879
在线: 606.3小时
虫号: 1517335
注册: 2011-11-30
性别: GG
专业: 常微分方程与动力系统

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
_romantic_镇: 金币+50, ★有帮助 2022-03-15 11:07:58

程序中的问题：
1、tspan = [0.00 : 0.01 : 0.44]与tspan=[0 0.222 0.333 0.444]没有对应
2、f = sum((y(:,1)-yexp(:,1)).^2) + sum((y(:,2)-yexp(:,2)).^2) + sum((y(:,3)-yexp(:,3)).^2) 后面的代码是多余的，计算的结果不使用
执行结果：

使用函数fmincon()估计得到的参数值为:
k1 = -0.6403
k2 = 1.8125
k3 = 3.6968
k4 = 1.9950
The sum of the squares is: 4.0

赞一下(1人)

回复此楼

QQ：2120156492

2楼2021-12-27 22:26:21