种群logistic增长模型.docVIP

种群logistic增长模型 生命科学院 09科五 卢春燕 20092501092 一、实验原理 logistic增长模型 ： 种群在有限环境下的“S”型增长曲线拟合的方程称为logistic方程： 其积分式为： K——环境容纳量；N——种群的数量； r——种群的瞬时增长率；t——时间。 实验步骤 1、制备草履虫培养液； 2、确定培养液中草履虫的初始密度； 3、定期观测和记录； 4、方程参数的估计 （1）K值的估计（均值法） K=111 （2）a、r的估计 求出K值后，将logistic方程的积分式变形为： 两边取对数，即为： 设 , b=-r，x=t， 则logistic方程的积分式变为： y=a+bx， 利用直线回归方法求得a、b（则r = -b） (3) 曲线的拟合 1） 将求得的K、a和 r代入logistic方程，建立logistic增长模型。 2） 计算得到各个增长时间种群大小的理论估计值，依照理论估计值绘制logistic方程的理论曲线。 3）可以进一步将理论估计值与实验观测值进行显著性检验，确定无显著性差异。 实验结果与讨论 表1 草履虫在培养液中增长实验数据统计分析表 天数 重复1只mL) 重复2 只/mL 重复3 只/mL) 平均值只/mL) (K-N)/N ln[(K-N)/N)] a-rt exp logistic 0 3 3 3 36 3.583519 1.346 3.8 22.92429 1 10 7 10 9 11.33333 2.427748 1.2981 3.7 23.80783 2 19 11 28 19.33333 4.741379 1.556328 1.2502 3.5 24.71587 3 27 16 31 24.66667 3.5 1.252763 1.2023 3.3 25.64836 4 5 61 81 49 1.265306 0.235314 1.1544 3.2 26.60518 5 66 179 87 110.6667 0.003012 -5.80513 1.1065 3.0 27.58616 6 35 40 15 30 2.7 0.993252 1.0586 2.9 28.59106 7 12 13 28 17.66667 5.283019 1.664498 1.0107 2.7 29.61956 8 11 10 19 13.33333 7.325 1.991293 0.9628 2.6 30.67129 9 13 8 23 14.66667 6.568182 1.882237 0.9149 2.5 31.7458 10 7 3 18 9.333333 10.89286 2.388107 0.867 2.4 32.84256 K的估计值为111(只/mL) 将logistic方程的积分式变形为： 两边取对数，即为： 设 , b=-r，x=t， 则logistic方程的积分式变为： y=a+bx， 利用直线回归方法求得a、b（则r = -b） 求得a=1.3460 ,b=-0.0479,代入逻辑斯蒂方程 111 求得 N = 1+e1.3460-0.0479r 图1 草履虫观察值散点图及拟合增长曲线图 表2 草履虫实验数据理论估计值与实验观测值显著性检验分析表 天数 观察值（只/mL） 理论值（只/mL) X2 X21，0.01 显著性 0 3 23 16.45865 6.63 极显著差异 1 9 24 8.5986 6.63 极显著差异 2 19 25 0.96453 6.63 无差异 3 25 26 0.009047 6.63 无差异 4 49 27 18.01841 6.63 极显著差异 5 111 28 247.2087 6.63 极显著差异 6 30 29 0.028896 6.63 无差异 7 18 30 4.428451 6.63 显著差异 8 13 31 9.24372 6.63 极显著差异 9 15 32 8.658394 6.63 极显著差异 10 9 33