| 查看: 1177 | 回复: 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 |
回复此楼» 猜你喜欢
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有72人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
-
3月 中国石油和化学工业联合会科学技术奖 马上开始，申获奖要点
已经有4人回复
-
中国石油和化学工业联合会科学技术奖的评审标准是什么？
已经有1人回复
-
经费如何使用，花不完怎么办？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
805592836_01
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 823
- 帖子: 10
- 在线: 23.4小时
- 虫号: 2393210
- 注册: 2013-04-01
- 性别: MM
- 专业: 能源化工
5楼2014-03-12 11:31:29
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
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2014-03-12 10:35:53
