数学建模实验报告4酵母培养物的离散阻滞增长模型..docVIP

成绩 实验名称 酵母培养物的离散阻滞增长模型 院系 理学院 专业 理科实验 班级 姓名 学号 日期 实验 目的 简述本次实验目的：掌握酵母培养物的阻滞增长差分方程模型的分析与建模方法，并会利用Matlab作数值计算与误差分析. 实验 准备 你为本次实验做了哪些准备：学习离散阻滞增长模型，学习如何用差分方程建立离散阻滞增长模型，学习用MATLAB做数值计算与误差分析 实验 进度 本次共有 1 个练习，完成 1 个。 实验 总结 日 本次实验的收获、体会、经验、问题和教训： 本次试验我收获颇丰，学会了用差分方程建立离散阻滞增长模型，离散阻滞增长模型就是一阶非线性差分方程，离散阻滞增长模型中一个关键的要素是什么因素导致增长率的递减，我们学习离散阻滞增长模型，不仅要学会计算结果，作图分析，更重要的是我们要学会分析其机理。 教师 评语 一．实验题目： 已知从测量酵母培养物增长的实验收集的数据如表： 时刻/h 0 1 2 3 4 5 6 7 8 9 生物量/g 9.6 18.3 29.0 47.2 71.1 119.1 174.6 257.3 350.7 441.0 时刻/h 10 11 12 13 14 15 16 17 18 生物量/g 513.3 559.7 594.8 629.4 640.8 651.1 655.9 659.6 661.8 二．实验要求 1、作图分析酵母培养物的增长数据、增长率、与相对增长率. 2、建立酵母培养物的增长模型. 3、利用线性拟合估计模型参数，并进行模型检验，展示模型拟合与预测效果图. 4、利用非线性拟合估计模型参数，并进行模型检验，展示模型拟合与预测效果图. 5、请分析两个模型的区别，作出模型的评价. 三．实验内容 (1)对于此问，可直接根据数据作图 先求相对增长率随时间的变化，程序如下： k=[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]; x=[9.6,18.3,29.0,47.2,71.1,119.1,174.6,257.3,350.7,441.0,513.3,559.7,594.8,629.4,640.8,651.1,655.9,659.6,661.8]; n=1; for n=1:18 dx(n)=x(n+1)-x(n); end r=dx./x(1:18); plot(0:17,r,kv) xlabel(时间k（小时）),ylabel(增长率 (%)) title(增长率与时间) 模拟效果图如下： 再求增长量随时间的变化，程序如下： k=[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]; x=[9.6,18.3,29.0,47.2,71.1,119.1,174.6,257.3,350.7,441.0,513.3,559.7,594.8,629.4,640.8,651.1,655.9,659.6,661.8]; n=1; for n=1:18 dx(n)=x(n+1)-x(n); end plot(0:17,dx,ko) xlabel(时间k（小时） ),ylabel(增长量 (克)) title(增长量与时间) 模拟效果图如下： （2）建立酵母培养物的模型 k---时刻（小时）； x(k)---酵母培养物在第k小时的生物量（克）； r(k)---用前差公式计算的生物量在第k小时的增长率； r---生物量的固有增长率； N---生物量的最大容量。 在营养有限的环境下，假设用前差公式计算的增长率r(k)随着生物量x(k)的增加而线性递减，即 r_k=(x_(k+1)-x_k)/x_k=r*(1-x_k/N),k=0,1,2… 根据以上模型假设，即可建立离散阻滞增长模型 x_(k+1)=x_k+r*x_k*(1-x_k/N),k=0,1,2… （3）首先，根据r_k和x_k的数据多项式拟合出（2）问中的r,N；然后根据生物量的观测数据直接取x_0=9.6，用(2)问中的循环语句进行迭代计算，算出0~18小时酵母生物量的模拟值，并计算误差平方和，绘制模拟效果图和模拟误差图。 程序如下：t=[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]; x=[9.6,18.3,29.0,47.2,71.1,119.1,174.6,257.3,350.7,441.0,513.3,559.7,594.8,629.4,640.8,651.1,655.9,659.6,661.8]; r=(x(2:19)-x(1:18))./x(