离散阻滞增长模型及其应用-精品教育文档.pptVIP

31 3.4.2 酵母培养物的增长 5. 模型求解和模型检验 方法一 首先，根据表 3.2 中的 k r 和 k x 的数据多项式拟合 出 (3.4.9) 式的参数 r 和 N ； 然后根据生物量的观测数据直接取 0 9.6 x ? ，用 循环语句按照 (3.4.2) 式进行迭代计算，算出第 0 至 18 小时的酵母生物量的模拟值，并计算误差平方和，绘 制模拟效果图和模拟误差图（图 3.11 ） . 长率 k r 随着生物量 k x 的增加而线性递减，即 1 1 , 0,1, 2, k k k k k x x x r r k x N ? ? ? ? ? ? ? ? ??? ? ? ? ? (3.4.9) k x ? ? 32 3.4.2 酵母培养物的增长 5. 模型求解和模型检验 方法一计算的 MATLAB 脚本： t=0: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(1:18); a1=polyfit(x(1:18),r,1); r1