应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】使用Matlab拟合反应动力学方程问题
51/1返回列表
查看: 3480  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2009-10-27 12:22:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhdgzhdg

至尊木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
目标函数正确，先用用四阶的龙格库塔法沿床层积分求出xcal，在用LSQNONLIN（马夸特）方法拟合曲线。
赞一下(1人) 回复此楼
4楼2010-09-01 06:11:52
已阅   回复此楼   关注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.
赞一下(13人) 回复此楼
2楼2009-10-28 09:48:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhdgzhdg

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

iris19860908

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
我也是初步设计反应动力学的，想请问楼主是如何消除内外扩散的呢？
赞一下(1人) 回复此楼
6楼2010-09-02 20:51:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 87/4350 2026-08-18 06:35 by 医学老男孩
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 时间戳变了，能看出什么问题？ +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
关闭
请填处理意见
关闭