24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

northerwood

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 895.4
散金: 6
帖子: 297
在线: 87.3小时
虫号: 1055710
注册: 2010-07-11
性别: GG
专业: 固体力学

[求助] 如何用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论文润色/翻译怎么收费? 已经有290人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

春风化雨乐未央，行健不息须自强

1楼 2011-09-26 20:58:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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];