| 查看: 3087 | 回复: 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 ] |
回复此楼» 猜你喜欢
生物学308分求调剂（一志愿华东师大）
已经有3人回复
-
求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过
已经有3人回复
-
环境工程 085701，267求调剂
已经有14人回复
-
308求调剂
已经有12人回复
-
求调剂：085600材料与化工，考材科基，总分319
已经有21人回复
-
311求调剂
已经有8人回复
-
294分080500材料科学与工程求调剂
已经有13人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有19人回复
-
288资源与环境专硕求调剂，不限专业，有学上就行
已经有23人回复
-
一志愿085600中科院宁波所276分求调剂
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- matlab拟合拟均相动力学参数
已经有6人回复
- 吸附的动力学方程拟合
已经有18人回复
- 【求助】拟合动力学方程求助
已经有13人回复
- 腐蚀动力学曲线如何通过origin拟合呢？
已经有4人回复
- 【求助】关于催化反应中的活化能
已经有24人回复
- origin，sigmaplot等软件自定义曲线拟合（主要是一些动力学方程的模拟）
已经有2人回复
- 【求助】光催化动力学方程拟合
已经有7人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
5楼2010-09-01 06:16:10
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+2,VIP+0):谢谢热心帮助。欢迎常来交流。若愿意成为本版的专家顾问，请PM我。 10-28 09:53
小木虫(金币+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
4楼2010-09-01 06:11:52
iris19860908
实习版主
- 应助: 0 (幼儿园)
- 金币: 2.9
- 散金: 1152
- 红花: 1
- 帖子: 132
- 在线: 94.9小时
- 虫号: 859638
- 注册: 2009-09-29
- 专业: 催化化学
6楼2010-09-02 20:51:55
