数学实验设计.docVIP

数学实验设计 9、 寄主-寄生现象是自然界中常有的,如黄蜂的幼虫寄生在象鼻虫(寄生)体内,一方面,寄主通常靠自然资源为生,不妨假定它的数量变化用离散形式的阻滞增长模型(45)描述,而寄生生物的存在会减少其增长.显然,寄生物数量越多,寄主的增长率减少得越多,最简单的假设是,寄主的减少率与寄生物数量成正比.另一方面,寄生物完全靠寄主为主,可以自然的假定它的相邻两代数量之比与寄主数量成正比,建立一个寄主-寄生模型研究二者数量变化的规律,讨论时间充分长以后的趋势. 记住第K代寄主的种群数量为Xk,最大容量为N,固有增长率为r ,寄生物的种群数量为yk(k=0.1.2….) a,b(0)是比例系数,a反映寄生物由寄主那里摄取营养,从而阻滞寄主增长能力,b反映寄主供养寄生物,从而使寄生物增长(或减少)的能力. 设定一组参数:寄主Xk的最大容量N=100,固有增长率r=1.5,寄生物最好的条件下每代的数量可以翻一番,既Xk =N时yk+1=2yk,于是?bN=2/N=0.02,a略大于b,设a=0.025,下面令初始值X0=50 , y0=10. 由题意可得： N=100;r=1.5;a=0.025;b=0.02; x=zeros(1,3000);y=zeros(1,3000); x(1)=50;y(1)=10; for i=2:3000 x(i)=x(i-1)+r*(1-x(i-1)/N)*x(i-1)-a*x(i-1)*y(i-1); y(i)=b*x(i-1)*y(i-1); end subplot(1,2,1);plot(x,y,.r);grid gtext(x,y的关系图);pause subplot(1,2,2);plot(x,.b); hold on plot(y,.r); gtext(x(k));gtext(y(k));grid jizhu 所以当k的值相对变大时，x在50的周围震荡，y在30的周围震荡。 令