| 查看: 3380 | 回复: 20 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
northerwood金虫 (小有名气)
|
[求助]
如何用MATLAB 实现化学反应方程式（写程序代码）?
|
||||
|
这个属于代码交流分类，已帮你转移 如何用MATLAB 实现化学反应方程式（写程序代码）? 例如：S + (ω+1)M ---------> L + (ω)D [ Last edited by dbb627 on 2011-10-16 at 13:10 ] |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有235人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
|
\code function f=enzym %生物化学米氏方程的仿真 %S+E<=>ES->E+P,反应动力学常数k1(S+E->ES)、k1r(ES->E+S)和k2(ES->E+P) %---- %其中S代表底物，E代表酶，ES表示酶和底物相合在一起的中间产物， %一个底物分子结合一个酶分子，然后生成一个分子新的物质P，并将酶释放出来， %酶本身并没有发生任何变化。 %已知反应动力学常数k1=2、k1r=1和k2=1.5，系统的初值S(0)=30 mol，E(0)=5 mol， ES(0)=0，P(0)=0 [ts,data] = ode45(@enzyme,0:0.01:10,[30 5 0 0]); plot(ts,data(:,1),'r',ts,data(:,2),'b',ts,data(:,3),'g',ts,data(:,4),'m') legend('[S]','[E]','[ES]','[P]') f=[ts,data]; function dmf=enzyme(t,X) S=X(1); E=X(2); ES=X(3); P=X(4); dS=-10*S+2*ES*S+ES; dE=2.5*ES-10*S+2*ES*S; dES=10*S-2*ES*S-2.5*ES; dP=1.5*ES; dmf=[dS;dE;dES;dP];\code 这段代码直接存为m文件 run即可 |
11楼2011-10-13 13:55:08
2楼2011-09-27 16:37:35
northerwood
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 895.4
- 散金: 6
- 帖子: 297
- 在线: 87.3小时
- 虫号: 1055710
- 注册: 2010-07-11
- 性别: GG
- 专业: 固体力学
3楼2011-09-28 08:53:31
4楼2011-09-28 12:35:36