24小时热门版块排行榜    

查看: 3065  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

apolloking

金虫 (小有名气)

[交流] 【求助】使用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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhdgzhdg

至尊木虫 (著名写手)


nono2009(金币+1):鼓励应助。欢迎常来。 2010-09-01 06:41:01
用LSQNONLIN最优化方法要先给初值,可以用文献值,如果没有文献值,可以给一个估计相近的数量级就可以了,如果用几个不同的初值,结果都相近,可以认为结果和算法可行。当然最后还是要计算残差的。
5楼2010-09-01 06:16:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

Rick

★ ★ ★
小木虫(金币+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.
2楼2009-10-28 09:48:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhdgzhdg

至尊木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
目标函数正确,先用用四阶的龙格库塔法沿床层积分求出xcal,在用LSQNONLIN(马夸特)方法拟合曲线。
4楼2010-09-01 06:11:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iris19860908

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
我也是初步设计反应动力学的,想请问楼主是如何消除内外扩散的呢?
6楼2010-09-02 20:51:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856调剂,是学校就去 +6 sllhht 2026-03-19 7/350 2026-03-19 19:50 by 制度的
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-19 3/150 2026-03-19 13:36 by houyaoxu
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
信息提示
请填处理意见