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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » Matlab同时拟合活化能，反应级数，速率常数
51/1返回列表
查看: 6103  |  回复: 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

铜虫 (正式写手)
引用回帖:
20楼: Originally posted by heisekaishui at 2014-08-28 19:50:16
非常感谢，其实我理解没错的话，其实Matlab解此类问题的方法就是同时拟合不同条件的下的多组数据，得到最优解？
以及我还请问程序里lsqnonlin函数中的数据dataTP1等，我看其它程序中有的会设置yexp = dataTP1(:,2 ...

没区别，一样的

[ 发自小木虫客户端 ]
回复此楼
24楼2014-09-09 12:31:40
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭