24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者修竹好衰将赠送您 20 个金币

修竹好衰

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 7.1
帖子: 28
在线: 22.6小时
虫号: 21829020
注册: 2020-04-14
专业: 能源化工

[求助] 初学用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)所以需要乘以质量

回复此楼