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 ]

回复此楼