数学建模酶促反应matlab求解,数学建模——酶促反应论文.doc

PAGE

word文档 可自由复制编辑

word文档 可自由复制编辑

数学建模

摘要

本文针对嘌呤霉素在某项酶促反应中对反应速度和底物浓度之间的关系的影响的问题，根据实际可知符合底物浓度与反应速度的模型有两种，即Michaelis-Menten模型和指数增长模型。对于Michaelis-Menten模型，例题中已经详细分析，不再详细讨论。本论文旨在建立指数模型对实际数据进行拟合分析。

由酶促反应的基本性质知，酶浓度x和反应速度y之间满足当底物浓度较小时，反应速度大致与浓度成正比；在底物浓度很大，渐进饱和时，反应速度将趋于一个固定值，由此建立一个指数增长模型 并使用Matlab中nlinfit函数对给出数据进行非线性回归，用cftool函数对结果进行验证，确定出

此时 。

为使模型更加准确，改进模型为 ，用同样的方法进行拟合与分析，得出 ， 和 ，此时 。

同过两个对模型进行预测与做残差图等方法，我们发现第二个模型相比第一个有所改进。

我们通过对实际问题的仔细分析，把实际问题转化成为数学上求解线性回归的问题，并建立了广为大家所熟悉的数学模型指数模型。通过数学软件的求解，得出模型中变量的系数。由于模型中的有些参数是估计的，考虑到实际与理论的差距，为了是使理论分析更贴近生活实际，我们从简略模型到优化模型进行了进一步分析，通过计算机利用数学软件MATLAB对问题进行了求解分析，得到了比较客观的分析结果。最后我们还根据模型的特点，对模型进行了推广，使其更具有一般性，能够解决更多实际问题。

关键词：指数模型 非线性回归 MATLAB nlinfit cftool 残差图

一、问题提出

某生化系学生为了研究嘌呤霉素在某项酶促反应中对反应速度和底物浓度之间的关系的影响，设计了两个实验，一个实验中使用的酶是经过嘌呤霉素处理的，而另一个是未经过嘌呤霉素处理的，所的实验数据见下表

底物浓度/ppm

0.02

0.06

0.11

0.22

0.56

1.10

反应速度

处理

76

47

97

107

123

139

159

152

191

201

207

200

未处理

67

51

84

86

98

115

131

124

144

158

160

——

对实验处理结果实际数据做非线性回归分析，其结果如何？试做模型的残差图进行比较。

二、基本假设

1、假设1. 当底物浓度较小时，反应速度大致与浓度成正比(即一级反应)；

2、假设2. 当底物浓度很大时，渐进饱和时，反应速度将趋于一个固定值—最终反应速度(即零级反应)

3、假设3. 反应速度与底物浓度成正比关系，即线性关系；

4、假设4. 反应速度不受外界温度的影响，而酶的活性在反应中一直保持不变。

三、符号说明

符号

意义

单位

备注

酶促反应速度

ppm/h

底物浓度

ppm

系数

/

a

系数

b

系数

四、问题分析

酶促反应动力学简称酶动力学，主要研究酶促反应的速度和底物浓度以及其他因素的关系。在底物浓度低时，酶促反应是一级反应；在底物浓度高时，酶促浓度是零级反应。

当反应浓度低时，反应速度大致与浓度成正比；当浓度很大时，渐进饱和时，反应速度趋近于一个固定值——最终反应速度。

以下指数增长模型满足这个性质

即 与

下面分别对这两个模型进行分析求解。

五、模型的建立与求解

5.1 模型一建立与求解

5.1.1 模型一的分析

由给出模型，先用Matlab中的nlinfit函数可求出系数，此处需先给出系数的初始值进行迭代，根据函数意义，为最终反映速度， 反映该酶促反应达最终速度的快慢，粗略估其值为220与10.

5.1.2 模型一模型的建立

5.1.3 模型一模型的求解

利用MATLAB统计工具箱中的nlinfit命令代入初值进行求解，将得到的结果作为初值再次代入到模型中求解，得到稳定的系数(见下表)。

参数

参数估计值

参数置信区间

192.0945

[173.8772 210.3117]

11.3854

[7.7571 15.0137]

rmse(剩余标准差)=17.4400

5.1.4 模型一结果的分析及验证

用MATLAB中cftool函数进行验证，得出结果如下：

f(x) = a*(1-exp((-b)*x))

Coefficients (with 95% confidence bounds):

a = 192.1 (173.9, 210.3)

b =