| 查看: 943 | 回复: 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)所以需要乘以质量 |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有189人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有4人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有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 (硕士)
- 金币: 17808.5
- 红花: 53
- 帖子: 2879
- 在线: 606.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
4楼2022-03-06 20:23:03
