| 查看: 1065 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者修竹好衰将赠送您 20 个金币 | ||
[求助]
初学用MATLAB求解动力学模型，但是在非线性二乘法处出现问题 ，麻烦各位大神帮忙看看 已有2人参与
|
||
|
出现的问题是： 错误使用 lsqncommon (line 64) The Levenberg-Marquardt algorithm does not handle bound constraints and the trust-region-reflective algorithm requires at least as many equations as variables; aborting. 出错 lsqnonlin (line 262) lsqncommon(funfcn,xCurrent,lb,ub,options,defaultopt,optimgetFlag,caller,... 出错 kinetic (line 16) lsqnonlin(@fun,k0,lb,ub,[],yexp); 其中反应数据是引用文件：‘reaction.xls’中每7行为一个反应温度下的各个时间点数据，总共3个温度下的数据，共计21行，已在后面图片里面显示 程序： clc clear all; global a b tspan=[0 1 3 5 7 9 11]; %各个反应时间点，0h，1h，3h…… yexp=xlsread('reaction.xls') %y的赋值从reaction这个excel文件中调取 k0=[2*10^8 60000 6.75*10^10 72000 3*10^11 90000]; %每个反应速率常数k所对应的A、E，总计6个参数初值设定（正反应、逆反应、副反应） lb=[0 0 0 0 0 0]; ub=[+inf +inf +inf +inf +inf +inf]; [k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@fun,k0,lb,ub,[],yexp); %===================================================================================== function f=fun(k,yexp) %这个就是拟合的目标函数,得到对应最合适的k f=[]; for i=1:3 %对应三个温度下的数据， tspan=[0 1 3 5 7 9 11]; a=7*(i-1)+1; x0=yexp(a,  ;[t,x]=ode45(@func,tspan,x0,[],k,i); d=a+6; yc1=x(:,1); yc2=x(:,2); yc3=x(:,3); yc4=x(:,4); f11=yexp(a:d,1)-yc1; f12=yexp(a:d,2)-yc2; f13=yexp(a:d,3)-yc3; f14=yexp(a:d,4)-yc4; ff=[f11 f12 f13 f14]; f=[f:ff]; end %===================================================================================== function dxdt=func(t,x,k,i) switch i case 1 m=0.6; V=0.2431; T=303.15; case 2 m=0.6; V=0.2448; T=313.15; case 3 m=0.6; V=0.2491; T=323.15; end r2=-(k(1)*exp(-k(2)/8.314/T)*x(1)/V*x(2)/V-k(3)*exp(-k(4)/8.314/T)*x(3)/V); r3=-r2; r4=k(5)*exp(-k(6)/8.314/T)*x(1)/V*x(1)/V; r1=r2-2*r4; dxdt=[m*r1;m*r2;m*r3;m*r4]; %因为这里r的单位是mol/(g min L)所以需要乘以质量 |
回复此楼» 猜你喜欢
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有523人回复
-
物理学I论文润色/翻译怎么收费?
已经有85人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
修竹好衰2楼2022-03-05 19:23:00
3楼2022-03-06 13:56:17
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17854.5
- 红花: 53
- 帖子: 2880
- 在线: 606.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
4楼2022-03-06 20:23:03
