《毕于慧作业---酶促反应.doc》由会员分享，可在线阅读，更多相关《毕于慧作业---酶促反应.doc（5页珍藏版）》请在金锄头文库上搜索。
1、酶促反应一 问题分析：为了研究酶促反应中嘌呤霉素对反应速度与底物浓度之间关系的影响，我们建立数学模型，反映该酶促反应的速度与底物浓度以及经嘌呤霉素处理与否之间的关系。酶促反应的基本性质是：底物浓度较小时，反应速度大致与浓度成正比，底物浓度较大或渐近饱和时，反应速度趋于固定值。设计了两个试验：酶经过嘌呤霉素处理，酶未经嘌呤霉素处理。所得试验数据如下：底物浓度(ppm)0.020.060.110.220.561.10反应速度处理764797107123139159152191201207200未处理6751848698115131124144158160/二 模型假设：反应过程中除底物浓度以外的其
2、他条件（如温度，PH值等）保持不变。三 符号说明：y - 酶促反应的速度；x - 底物浓度；， - 待定系数；四 模型建立：将数据绘成散点图：酶经嘌呤霉素处理后，图为:酶未经嘌呤霉素处理，图为：由图知：数据符合非线性模型五 模型求解：利用matlab进行运算：（）计算回归系数的初值，利用经嘌呤霉素处理后的数据， y = 转化为线性模型： 即 代码： x=0.02 0.02 0.06 0.06 0.11 0.11 0.22 0.22 0.56 0.56 1.10 1.10; y=76 47 97 107 123 139 159 152 191 201 207 200; m=1./x; n=1./
3、y; t=ones(12,1),m; format long b,bint,r,rint,stats=regress(n,t)结果：b = 0.00510718164159 0.00024722096110bint = 0.00353857552724 0.00667578775593 0.00017568750886 0.00031875441334r = -0.00431033495973 0.00380836604811 0.00108174735727 0.00011826339928 0.00077543637650 -0.00016040031999 0.0000583948931
4、5 0.00034803408547 -0.00031304554931 -0.00057352326544 -0.00050101009546 -0.00033192796986rint = -0.00495843977711 -0.00366223014234 0.00217493814144 0.00544179395478 -0.00309132595844 0.00525482067298 -0.00413053901576 0.00436706581432 -0.00342529534515 0.00497616809816 -0.00439875450551 0.00407795
5、386552 -0.00413843515395 0.00425522494025 -0.00384105253041 0.00453712070136 -0.00446352105162 0.00383742995300 -0.00470862422210 0.00356157769122 -0.00462624728486 0.00362422709393 -0.00446657058966 0.00380271464994stats = 0.85569474384568 59.29754512410619 0.00001642169065所以,1= 0.00510718164159,置信
6、区间为0.00353857552724 0.00667578775593；2= 0.00024722096110,置信区间为0.00017568750886 0.00031875441334。1=1/1=195.8027，2=2/1= 0.04841（二）求解非线性模型：针对要拟合的非线性模型y = ，建立一个m-文件(f1.m)如下：function k=f1(beta,x)k=beta(1)*x./(beta(2)+x);运算时输入代码： x=0.02 0.02 0.06 0.06 0.11 0.11 0.22 0.22 0.56 0.56 1.10 1.10; y=76 47 97 10
7、7 123 139 159 152 191 201 207 200; beta0=195.8027 0.04841; beta,r ,J=nlinfit(x,y,f1,beta0); betaci =nlparci(beta,r,J); beta,betaci结果：beta = 1.0e+002 *2.12681838233292 0.00064118312791betaci = 1.0e+002 * 1.97202843968233 2.28160832498351 0.00045667914052 0.00082568711530所以，1的估计值为212.6819，置信区间为197.202
8、8 228.1608;2的估计值为0.0641,置信区间为0.0457 0.0826。代码： a=0.02;0.01;1.10; b=beta(1)*x./(beta(2)+x); plot(x,y,+,a,b)结果：代码： nlintool(x,y,f1,beta0)结果：所以，最终反应速度为1=212.6819，半速度点为2=0.0641。在同一模型中考虑嘌呤霉素处理的影响，需要用混合反应模型，即将y = 转化为y = , 其中，x1为底物浓度； x2为一示性变量； x2=1表示经过处理，x2=0表示未经处理；1是未经处理的最终反应速度；1是经处理后最终反应速度的增长值；2是未经处理的反应的半速度点；2是经处理后反应的半速度点的增长值。基于对数据的分析, 1=170,2=0.05,1=60,2=0.01代码：x1=0.02 0.02 0.06 0.06 0.11 0.11 0.22 0.22 0.56 0.56 1.10 1.10; x2=1,0; y=76 47 97 107 123 139 159 152 191 201 207 200;beta0=170,60,0.05,0.01; beta,r,J=nlinfit(x1,x2,y,f2,beta0);结果：
