24小时热门版块排行榜

返回列表

apolloking

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1083.1
散金: 35
红花: 2
帖子: 156
在线: 60.5小时
虫号: 75702
注册: 2005-06-20
专业: 能源化工

[交流] 【求助】使用Matlab拟合反应动力学方程问题已有4人参与

求助大家一个问题

在固定床积分反应器上测定催化剂催化性能
条件如下
给定不同的温度和空速测定相应的转化率
如
温度空速转化率
350 5000 0.224067
375 5000 0.308978
400 5000 0.418433
425 5000 0.583276
450 5000 0.712133
475 5000 0.834063
500 5000 0.927528
525 5000 0.971279
550 5000 0.989578
350 10000 0.167333
375 10000 0.26823
400 10000 0.359565
425 10000 0.492762
450 10000 0.605102
475 10000 0.706325
500 10000 0.870149
525 10000 0.936043
550 10000 0.972126
350 15000 0.161258
375 15000 0.229717
400 15000 0.305041
425 15000 0.445914
450 15000 0.562024
475 15000 0.702027
500 15000 0.847299
525 15000 0.901723
550 15000 0.956323
350 20000 0.103751
375 20000 0.143094
400 20000 0.227758
425 20000 0.297666
450 20000 0.456384
475 20000 0.592814
500 20000 0.743624
525 20000 0.825191
550 20000 0.897261
350 25000 0.089297
375 25000 0.149173
400 25000 0.213542
425 25000 0.304452
450 25000 0.433972
475 25000 0.563372
500 25000 0.703383
525 25000 0.801953
550 25000 0.870735
也就是一个温度 对于不同的空速有不同的转化率
在试验中消除了内外扩散影响
反应器按照等温操作时理想置换积分反应器处理
ra=k0*exp（-E/R/T）*CA0/FA0*（1-x）
ra=dx/dW

两式联立
dx/dW=k0*exp（-E/R/T）*CA0/FA0*（1-x）
文献中说是用四阶的龙格库塔法沿床层积分上式采用Nelder-Mead单纯形法进行参数优化参数优化目标函数如附图所示

最终求得k0=7.32E9
E为64550J/mol

请问这个用matlab应该怎么去做啊？
谢谢

[ Last edited by apolloking on 2009-10-27 at 12:36 ]

回复此楼

» 猜你喜欢

2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
孩子确诊有中度注意力缺陷已经有6人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab 求指点动力学方程拟合过程中导数的获取已经有10人回复
matlab拟合拟均相动力学参数已经有6人回复
吸附的动力学方程拟合已经有18人回复
【求助】拟合动力学方程求助已经有13人回复
腐蚀动力学曲线如何通过origin拟合呢？已经有4人回复
【求助】关于催化反应中的活化能已经有24人回复
origin，sigmaplot等软件自定义曲线拟合（主要是一些动力学方程的模拟）已经有2人回复
【求助】光催化动力学方程拟合已经有7人回复
【求助】使用Matlab预估动力学方程问题已经有13人回复

1楼2009-10-27 12:22:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Rick

应助: (幼儿园)
在线:
虫号: 44235

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+2,VIP+0):谢谢热心帮助。欢迎常来交流。若愿意成为本版的专家顾问，请PM我。 10-28 09:53

There are multiple things involved:
1. Numerical Integration of dX/dW. Check 'help ode15s'
2. Optimization. Simplex method is coded in MATLAB function 'fminsearch'. however, I would suggest you try 'nlinfit' in your case.

3. Pose your problem properly. I would say this is the hardest step and it might take a while.

Anyway, 1 and 2 are good starting points in my opinion.