| 查看: 884 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
805592836_01铁虫 (初入文坛)
|
[求助]
未完成的matlab 已有2人参与
|
|
|
有实验测得的一系列时间点 t=（0,180,300,420,600,900,1200,1800,2400） CA=（1.845，1.414,1.237,1.065,0.873,0.629,0.440,0.226,0.124） 微分方程：k(1)*k(2)*k(3)*CA0*(CB.^3)./(1+k(2)*CA+k(3)^(1/6)*CB.^(1/2)).^7 CB=0.4为定值 希望由此得到三个参数k(1),k(2),k(3)，以及实验值与拟合值的比较图。现已经编了部分代码，但不是很确定对否，希望高手可以修改一下 clear all; clc k0=[1,1,1]; lb=[0,0,0]; ub=[+inf,+inf,+inf]; tspan=[0,180,300,420,600,900,1200,1800,2400]; yexp=[1.414,1.237,1.065,0.873,0.629,0.440,0.226,0.124]'; y0=1.845; % 使用fmincon()进行参数估计 [k,fval,flag]=fmincon(@ObjFunc4Fmincon,k0,lb,ub,option,tspan,y0,yexp); fprintf('\tk1=%.4f\n',k(1)), %本征动力学参数 fprintf('\tk2=%.4f\n',k(2)), % A的吸附平衡常数KA fprintf('\tk3=%.4f\n',k(3)), % H2的吸附平衡常数KB fprintf('The sum of the squares is:%.1e\n\n',fval), k_fmincon=k; % 使用函数Isqnonlin()进行参数估计 yy=[y0 yexp']; options = optimset('MaxFunEvals',100000) [k,resmorm,residual,exitflag,output,lambda,jacobian]=... Isqnonlin(@ObjFunc4LNL,k0,lb,ub,options,tspan,y0,yexp); ci=nlparci(k,residual,jacobian); fprintf('n\n使用Iaqunonlin()估计得到的参数值:\n'),Output % 以函数fmincon()估计得到的结果为初值，使用Isqnonlin()进行参数估计 k0=k_fmincon; [k,resmorm,residual,exitflag,output,lambda,jacobian]=... Isqnonlin(@ObjFunc4LNL,k0,lb,ub,options,tspan,y0,yexp); ci=nlparci(k,residual,jacobian); fprintf('n\n以fmincon()结果为初值，使用函数Isqnonlin()估计得到的参数值为:\n') Output %------------------------------------------------------------ function f=ObjFunc4Fmincon(k,tspan,y0,yexp) [t y]=ode45(@KineticsEqs,tspan,y0,yexp,k); %---------------------------------------------------------- function f=Objfunc4LNL(k,tspan,yexp) [t y]=ode45(@KineticEqs,tspan,tspan,[],k); %---------------------------------------------------------- function dydt = KineticsEqs(t,y,k) y0=1.845; CB=y; CB0=0.4; CA0=y0; k1=k(1); KA=k(2); KB=k(3); dxdt=k(1)*k(2)*k(3)*CA0*(CB.^3)./(1+k(2)*CA0+k(3)^(1/6)*CB.^(1/2)).^7 |
回复此楼» 猜你喜欢
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有6人回复
-
物理学I论文润色/翻译怎么收费?
已经有270人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有30人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有63人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab中的计算问题？
已经有8人回复
- MATLAB求助，关于importdata的问题
已经有6人回复
- 用matlab制作的M文件保存在哪里才能直接调用？
已经有9人回复
- matlab一运行就出现的问题
已经有7人回复
- matlab的m文件关联问题
已经有17人回复
- 求助一个matlab关于分类讨论的问题
已经有9人回复
- matlab编程中遇到的问题
已经有4人回复
- 帮忙修改Matlab中出现的错误
已经有4人回复
- matlab读取图像的问题
已经有12人回复
- 分享matlab的一些代码（入门贴）
已经有19人回复
- 求MATLAB清华版的教学视频
已经有14人回复
- 关于Matlab作图的问题
已经有14人回复
- MATLAB初学者，这样的式子怎么就是画不了图，总是出错!求指导!!!!!!!!!!!!!
已经有8人回复
- matlab程序将图中的数据记录下来
已经有4人回复
- Matlab在地学中的应用
已经有24人回复
- matlab中dydt的编写要求
已经有12人回复
- 求助svm在matlab中的应用
已经有7人回复
- 关于matlab 2010b 的一个小问题
已经有7人回复
- MATLAB安装后出现一堆java的问题
已经有12人回复
- 【求助】关于matlab的高维数组
已经有6人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2014-03-12 10:35:53
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2014-03-12 09:56:13
805592836_01
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 823
- 帖子: 10
- 在线: 23.4小时
- 虫号: 2393210
- 注册: 2013-04-01
- 性别: MM
- 专业: 能源化工
3楼2014-03-12 10:02:43
805592836_01
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 823
- 帖子: 10
- 在线: 23.4小时
- 虫号: 2393210
- 注册: 2013-04-01
- 性别: MM
- 专业: 能源化工
5楼2014-03-12 11:31:29
