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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » Matlab同时拟合活化能，反应级数，速率常数
51/1返回列表
查看: 5230  |  回复: 26
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

heisekaishui

银虫 (小有名气)

[交流] Matlab同时拟合活化能，反应级数，速率常数

各位，您好!
本人研究的方向是催化化学，最近在评估催化剂效果时想引入动力学的内容，对其中的一些地方有一些疑问。
受月之蓝版主的这个帖子http://muchong.com/bbs/viewthread.php?tid=6425538&target=self&page=1以及其它前辈的启发，对于反应级数以及反应速率常数的拟合有一定的了解。
但是如果我的动力学模型是这种形式：-dC/dt=k0*exp(-E/R*T)*C^m*P^n（只有一个反应物和氢气），其中P与T分别为压力与时间，就是说是一个加氢反应，试验中可以得到不同压力、不同温度下浓度随时间的变化值，应该如何拟合k0,E,m,n。
我的一个想法是先拟合不同压力下的k1（k1=k0*exp(-E/R*T)*P^n)（例如以上述帖子的代码方法），再根据k1=k2*P^n（k2=k0*exp(-E/R*T)）拟合出k2和n，此后拟合出不同温度下的k2，根据k2=k0*exp(-E/R*T)）拟合出E和k0
但是我觉得这种方法较为复杂，一个matlab里有好几个拟合程序，不知道是否有更简洁的CODE。

PS:不知道版里有没有相似的帖子有类似的CODE可以让我学习一下（我没搜到啊）。谢谢
回复此楼

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2014-08-25 17:30:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whqs8426212

铜虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
16楼: Originally posted by heisekaishui at 2014-08-27 11:00:36
非常感谢，对，数据格式是这样的。
主要是数据还没有做出来，如果可以的话可以让我看一下程序吗，我想看一下大概的思路。...

数据造了三组，对于看四个参数背景不了解，所以初值也选的不合理，以至于运行结果比较糟糕
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

heisekaishui
19楼2014-08-27 11:35:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 27 个回答

heisekaishui

银虫 (小有名气)
纯粹的初学者，如果叙述有什么问题希望大家指出。
赞一下(1人) 回复此楼
2楼2014-08-25 17:31:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangjiguo

木虫 (职业作家)
★
heisekaishui(金币+5): 谢谢参与
得金币了，说声谢谢！
赞一下 回复此楼
5楼2014-08-25 17:54:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heisekaishui

银虫 (小有名气)
同学们金币诚可贵，回答更重要啊。
赞一下 回复此楼
9楼2014-08-25 18:17:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
coal-xust13楼
2014-08-26 15:56   回复  
查看全部 27 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭