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 ]

回复此楼

» 猜你喜欢

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

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的回帖

zhdgzhdg

至尊木虫 (著名写手)

应助: 32 (小学生)
金币: 22701.4
红花: 6
帖子: 1750
在线: 973.9小时
虫号: 95566
注册: 2005-10-05
专业: 化学反应工程

★
nono2009(金币+1):鼓励应助。欢迎常来。 2010-09-01 06:41:01

用LSQNONLIN最优化方法要先给初值，可以用文献值，如果没有文献值，可以给一个估计相近的数量级就可以了，如果用几个不同的初值，结果都相近，可以认为结果和算法可行。当然最后还是要计算残差的。

赞一下(1人)

回复此楼

5楼2010-09-01 06:16:10

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

查看全部 8 个回答

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.

赞一下(13人)

回复此楼

2楼2009-10-28 09:48:28

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