24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ronghengyidu

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2854.5
帖子: 218
在线: 62.2小时
虫号: 1003427
注册: 2010-04-22
专业: 生物化工与食品化工

[求助] matlab拟合反应动力学参数已有1人参与

本人是matlab大白，最近仿照黄华江编著的“实用化工计算机模拟”例题7-5拟合动力学参数，
         动力学方程为dc1/dt=k1*c1/(k2+k3*c1+k4*c2+k5*c1*c2)，拟合k1,k2,k3,k4,k5,
         matlab程序和动力学数据见附件，运行出现如下错误，请各位大拿帮助。


Error: File: KineticsEst.m Line: 47 Column: 60
Unbalanced or unexpected parenthesis or bracket.

Error in run (line 64)
evalin('caller', [script ';']);

Error: File: KineticsEst.m Line: 47 Column: 60
Unbalanced or unexpected parenthesis or bracket.

Error in run (line 64)
evalin('caller', [script ';']);

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : KineticsEst.m

2016-10-03 16:39:52, 1.98 K

附件 2 : KineticsData.m

2016-10-03 16:40:16, 365 bytes

1楼 2016-10-03 15:44:36

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

ronghengyidu

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2854.5
帖子: 218
在线: 62.2小时
虫号: 1003427
注册: 2010-04-22
专业: 生物化工与食品化工

送红花一朵

引用回帖:

2楼: Originally posted by lu_yu_lan at 2016-10-03 21:21:39
楼主的问题，没有看明白，找到了黄华江编著的“实用化工计算机模拟”例题7-5拟合动力学参数

7个微分方程，5个拟合参数k1,k2,k3,k4,k5，初值t=0时，x为0.1883, 0.2507, 0.0467, 0.0899, 0.1804, 0.1394, 0 ...

非常感谢您的回复，我只有一个微分方程，dc1/dt=k1*c1/(k2++k3*c1+k4*c2+k5*c1*c2), 拟合5个参数k1,k2,k3,k4,k5。我的t，c1，c2数据在附件中，matlab程序也在附件中，这两个文件都可以用记事本打开。在这里我把程序和实验数据再贴出来。

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

clear all
clc

k0 = [0.5  0.5  0.5  0.5  0.5];       % 参数初值
lb = [0  0  0  0  0];                % 参数下限
ub = [+inf  +inf  +inf  +inf  +inf]; % 参数上限
x0 = [150  1500];
KineticsData;
yexp = ExpData(:,2:3);                % yexp: 实验数据[x1 x2]

% 使用函数fmincon()进行参数估计
[k,fval,flag] = fmincon(@ObjFunc4Fmincon,k0,[],[],[],[],lb,ub,[],[],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('\tk5 = %.4f\n',k(5))
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 = ExpData(:,1);
[t x] = ode45(@KineticEqs,tspan,x0,[],k);
y(:,1:2) = x(:,1:2);
f = sum((y(:,1)-yexp(:,1)).^2) + sum((y(:,2)-yexp(:,2)).^2));

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

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

实验数据如下：

动力学数据:
t    c1       c2
0    150       1500
10 137.3 1487.3
30 107.7 1457.7
60 72.8    1422.8
90 48.8    1398.8
120 22.9    1372.9
150 11.8    1361.8
180 5.4    1355.4
220 1.9    1351.9
240 0.9    1350.9

期望楼上能够指点，如果您能够给出openLU程序，我也非常感谢，因为我第一次听说这个程序，怕仿照编程又要出错。非常感谢您在节假日能够回答我的问题。