内容发布更新时间 : 2024/11/6 3:53:34星期一 下面是文章的全部内容请认真阅读。
西北农林科技大学实验报告
学院名称:理学院 专业年级:2011级信计1班 姓 名:xxx 学 号:xxx 课 程:数学模型与数学建模 报告日期:2013年12月24日
1 实验题目:
生物种群动态模拟 2 实验问题陈述:
在这个问题中我们采用较为复杂的增长模型g(P)?(r/c)(P?c)(1?P/K)P, 其中参数C表示种群水平的最小值,低于这个水平种群将出现负增长,假设a?10 ,且最小值的种群水平是蓝鲸3000条,长须鲸15000条。
(a)两种鲸鱼种群能否共存?建立一个动力系统模型描述他们的动态。 (b)画出这个模型的向量场。确定每个平衡态是否稳定。
(c)假设存在5000条蓝鲸和70000条长须鲸,关于这两种鲸鱼的将来这个模型能预测些什么?
(d)假设我们低估了蓝鲸的最小有效种群值,这个值实际上接近10000.于是两个种群会发生什么变化? 3 实验目的:
(1)了解生物种群的模型的建立过程, 建立蓝鲸和长须鲸的交互作用种群动态模型,并回答两种群的共存问题、平衡态稳定与否问题、预测问题及调整最小有效种群值两种群变化问题。
(2)熟悉解决类似鲸群问题模型的设计过程并说明数学模型上的动态模拟过程。 (3)本次实习中主要用到的软件主要是MATLAB软件,所以本次实习目的还有了解MATLAB软件的使用来求解生物种群动态模拟。 4 实验内容:
(1)模型分析:
假设:(1)只考虑相互竞争的两个生物种群:蓝鲸与长须鲸的动态; (2)当种群将没有竞争时,这两个种群都为logistic增长; (3)竞争种群的存在将降低对方种群的增长率;
(4)由竞争导致的增长率的减少量与当时竞争的种群数量成正比; (5)假设两个种群间的竞争力是一样的。 (2)符号表示
表1:一些符号说明
分别记蓝鲸和长须鲸的内禀增长率为:rB,rF ?7第 1 页 共 4 页
环境的承载力为:KB,kF 记蓝鲸对长须鲸的竞争力为CB 长须鲸的竞争力为CF (3)模型分析与建立 据题意有:
cB=3000,cF?15000
由假设②,对于单个种群的增长,有:
dB?(rB/cB)(1?B/KB)(B?cB)B…………………………….① dtdF ?(rF/cF)(1?F/KF)(F?cF)F………………………....②
dt
当竞争种群存在时,模型为:
dB?[(rB/cB)(1?B/KB)?CFF](B?cB)B…………..③ dt
dF?[(rF/cF)(1?F/KF)?CPP](F?cF)F…………………④ dt以上为这两种群的lotka-volterra竞争模型。
且已知:rB=0.05,rF=0.08;KB?15?10,KF?40?10, 由假设⑤有:CB?CF??; 由此可得模型为:
44dBB?[0.05(1??10?4)(B?3000)/3000??F]B……………….⑤ dt15dFF?[0.08(1??10?4)(F?15000)/15000??B]F……………..⑥ dt40其为一个二维的自治系统,相空间应该是二维的平面N=(B,F),针对实际问题的要求,应该在这个平面的第一象限S={(B,F):B>=0,F>=0}中讨论系统的动态。
将模型记为:B'?f1(B,F)B;F'?f2(B,F)F,其中: f1(B,F)?(rB/cB)(1?B/KB)(B?cB)?CFF
?(0.05/3000)?(1?B/150000)(B?3000)??F f2(B,F)?(rF/cF)(1?F/KF)(F?cF)?CBB
?(0.08/15000)?(1?F/400000)(F?15000)??B
5 结果分析与讨论:
第 2 页 共 4 页
(1)取x?B?10?4,y?F?10?4:
则程序运行后的平衡点有: [0 , 0] [3000 , 0] [0 , 15000] [0 , 400000] [150000 , 0] [5490 , 399893] [147527 , 397104] [149890 , 17942] [3092 , 15060]
6x 105-...- (x2 (x2/400000 - 1) (x2 - 15000))/187500 = 054长须鲸321000.511.5蓝鲸22.5x 1035
图一 鲸鱼竞争种群动态的向量场
从图一可以看出,两个生物种群将呈现的动态结果是蓝鲸将会灭绝,长须鲸将会稳定于它的饱和状态kF。由上面的图像可以看出不是每个平衡点都是稳定的,只有到达点(5490,39 9893)才能达到最终平衡,此时蓝鲸数量约为5490,长须鲸数量约为399 893。
当蓝鲸为5 000,长须鲸为70 000时,从上图我们可以看出在一段时间内,蓝鲸和长须鲸的数量都会增长。当蓝鲸的最小有效种群为10 000时,程序运行后有:
平衡点为(下面的数据为略去了小数的数据): [0 , 0] [10000 , 0] [0 , 15000] [0 , 400000] [150000 , 0] [19163 , 399626] [140895 , 397235] [149614 , 17936] [10326 , 15201]
第 3 页 共 4 页