| 查看: 1116 | 回复: 5 | |||
[交流]
使用MATLAB估计SCR动力学参数，运行有错误，请大神指点一下。（附件中红色为目标值） 已有3人参与
|
|
function K1K2K3K4 format long clear all clc tspan = [ ]; x0 = [0];%初始状态 k0 = [1 1 1 1 1 1]; %定义指前因子K1K3活化能K2K4 lb = [0 0 0 0 0 0];%下限 ub = [+inf +inf +inf +inf +inf +inf];%上限 data=... %第一列为时间，第二列NH3为输出值 [ ]; yexp = data(:,2); %NH3输出值 [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\tk1 = %.9f ± %.9f\n',k(1),ci(1,2)-k(1)) fprintf('\tk2 = %.9f ± %.9f\n',k(2),ci(2,2)-k(2)) fprintf('\tk3 = %.9f ± %.9f\n',k(3),ci(3,2)-k(3)) fprintf('\tk3 = %.9f ± %.9f\n',k(4),ci(4,2)-k(4)) fprintf(' The sum of the squares is: %.9e\n\n',resnorm) function f = ObjFunc(k,tspan,x0,yexp) % 目标函数 [t Xsim] = ode45(@reaction,tspan,x0,[],k); Xsim1=Xsim(:,1); ysim(:,1) = Xsim1(2:end); size(ysim(:,1)); size(yexp(:,1)); f = (ysim(:,1)-yexp(:,1)); function dCdt = reaction(t,C,k) % ODE模型方程 R=8.3145;%气体常数 Rseg=288; Pamb=101325; n=2;%（微元个数） f=0.81; Vc=0.005788730358;%(体积） Cs=155.88;%储氨能力（mol/m3) mf=225/3600;%废气流量(kg/s) Tin=;%输入温度（K) nNH3=400; dC1dt=K5*mf*(Tin-C1)-K6(C1^4-298.15^4); T=(C1+Tin)/2; dC2dt=K1*exp(-K2/(R*T))*C3*(1-C2)-K3*exp(-K4/(R*T))*C2; dC3dt=n/(f*Vc)*nNH3-C3*((Rseg*n/(Pamb*f*Vc))*mf*T+Cs*K1*exp(-K2/(R*T))*(1-C2))+Cs*K3*exp(-K4/(R*T))*C2; dCdt = [dC1dt;dC2dt;dC3dt]; Error in K1K2K3K41>reaction (line 902) dC2dt=K1*exp(-K2/(R*T))*C1*(1-C2)-K3*exp(-K4/(R*T))*C2; Error in odearguments (line 88) f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ode45 (line 114) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in K1K2K3K41>ObjFunc (line 884) [t Xsim] = ode45(@reaction,tspan,x0,[],k); Error in lsqnonlin (line 189) initVals.F = feval(funfcn{3},xCurrent,varargin{:}); Error in K1K2K3K41 (line 872) [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... Caused by: Failure in initial user-supplied objective function evaluation. LSQNONLIN cannot continue. |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Book1.xlsx
- 附件 2 : Book1.xlsx
2016-09-26 14:24:42, 380.43 K
2016-09-26 14:27:55, 402.16 K
» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有222人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急求：使用MATLAB编写遗传算法的程序，求解非线性函数的参数
已经有2人回复
- 急求：使用MATLAB编写遗传算法的程序，求解非线性函数的参数
已经有1人回复
- 【MATLAB书籍053】《分类,参数估计和状态估计：一个使用MATLAB的工程方法》(英文版)
已经有57人回复
- matlab用fmincon函数求极大似然估计
已经有0人回复
- 《Visual+C++_MATLAB图像处理与识别实用案例精选》
已经有1人回复
- Matlab指导
已经有0人回复
- 求助Matlab
已经有0人回复
- matlab 线性插值问题
已经有2人回复
- Matlab求助
已经有0人回复
- 求Matlab高手指导
已经有12人回复
- [诚心求问！]利用matlab做动力学参数估计时遇到问题，求高手赐教！
已经有1人回复
- 未完成的matlab
已经有5人回复
- matlab和isight的诡异问题
已经有9人回复
- matlab求助
已经有13人回复
- MATLAB关于lsqcurvefit程序几个错误，求大神指导一下
已经有6人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- 【软件】Mathworks Matlab R2011a ISO 正式注册版
已经有127人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.6小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2016-09-26 20:08:52
3楼2016-09-26 21:47:39
4楼2016-09-27 07:45:55
5楼2016-11-21 14:34:33
grangergaga
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 2875.4
- 散金: 751
- 红花: 3
- 帖子: 327
- 在线: 84.6小时
- 虫号: 2511069
- 注册: 2013-06-17
- 专业: 多相流热物理学
6楼2021-10-22 22:47:08