资源描述:
《毕于慧作业---酶促反应》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、酶促反应一问题分析:为了研究酶促反应中嘌呤霉素对反应速度与底物浓度之间关系的影响,我们建立数学模型,反映该酶促反应的速度与底物浓度以及经嘌呤霉素处理与否之间的关系。酶促反应的基本性质是:底物浓度较小时,反应速度大致与浓度成正比,底物浓度较大或渐近饱和时,反应速度趋于固定值。设计了两个试验:酶经过嘌呤霉素处理,酶未经嘌呤霉素处理。所得试验数据如下:底物浓度(ppm)0.020.060.110.220.561.10反应速度处理764797107123139159152191201207200未处理6751848698115131124144158160/二模型假设:反应过
2、程中除底物浓度以外的其他条件(如温度,PH值等)保持不变。三符号说明:y---酶促反应的速度;x---底物浓度;β,β---待定系数;四模型建立:将数据绘成散点图:酶经嘌呤霉素处理后,图为:酶未经嘌呤霉素处理,图为:由图知:数据符合非线性模型五模型求解:利用matlab进行运算:(—)计算回归系数的初值,利用经嘌呤霉素处理后的数据,y=转化为线性模型:即代码:>>x=[0.020.020.060.060.110.110.220.220.560.561.101.10];>>y=[764797107123139159152191201207200];>>m=1./x;>>
3、n=1./y;>>t=[ones(12,1),m'];>>formatlong>>[b,bint,r,rint,stats]=regress(n',t)结果:b=0.005107181641590.00024722096110bint=0.003538575527240.006675787755930.000175687508860.00031875441334r=-0.004310334959730.003808366048110.001081747357270.000118263399280.00077543637650-0.000160400319990.000
4、058394893150.00034803408547-0.00031304554931-0.00057352326544-0.00050101009546-0.00033192796986rint=-0.00495843977711-0.003662230142340.002174938141440.00544179395478-0.003091325958440.00525482067298-0.004130539015760.00436706581432-0.003425295345150.00497616809816-0.004398754505510.004
5、07795386552-0.004138435153950.00425522494025-0.003841052530410.00453712070136-0.004463521051620.00383742995300-0.004708624222100.00356157769122-0.004626247284860.00362422709393-0.004466570589660.00380271464994stats=0.8556947438456859.297545124106190.00001642169065所以,θ1=0.00510718164159,
6、置信区间为[0.003538575527240.00667578775593];θ2=0.00024722096110,置信区间为[0.000175687508860.00031875441334]。β1=1/θ1=195.8027,β2=θ2/θ1=0.04841(二)求解非线性模型:针对要拟合的非线性模型y=,建立一个m-文件(f1.m)如下:functionk=f1(beta,x)k=beta(1)*x./(beta(2)+x);运算时输入代码:>>x=[0.020.020.060.060.110.110.220.220.560.561.101.10];>>y=
7、[764797107123139159152191201207200];>>beta0=[195.80270.04841]';>>[beta,r,J]=nlinfit(x',y','f1',beta0);>>betaci=nlparci(beta,r,J);>>beta,betaci结果:beta=1.0e+002*2.126818382332920.00064118312791betaci=1.0e+002*1.972028439682332.281608324983510.000456679140520.00082568711530所以,β1的估