猫和老鼠的问题.pdf

猫和老鼠的问题 一、基本假设： (1)老鼠由于猫的存在使增长率降低, 假设降低的程度与猫数量成正比; (2)猫由于老鼠为它提供食物的作用使其死亡率降低或使之增长，假定增长的程度与老鼠数 量成正比。 二、符号说明： x——老鼠在 t 时刻的数量； a——老鼠独立生存时的增长率； e——猫掠取老鼠的能力； f——老鼠对猫的供养能力. y——猫在 t 时刻的数量； b——猫独立生存时的死亡率； K——捕获能力系数. 三、 模型(一) 不考虑人工捕获 该模型反映了在没有人工捕获的自然环境中老鼠与猫之间的制约关系, 没有考虑老鼠和 猫自身的阻滞作用, 是最简单的模型. 四、 模型(一) 求解 利用微分方程的相关理论，知原方程组的解是周期解，设周期为 T，则为了解释问题中的数 据，需计算 x 、y 的平均值： 五、模型(二) 考虑人工捕捞 K——捕获能力系数 类似可计算 x 、y 的平均值： 结论：使用灭鼠灵后猫的数量平均值降低，而老鼠平均值增加；进一步捕捞能力系数下降也 导致猫数量上升。 推广：解释杀虫剂的反效果——杀虫剂在杀死害虫的同时也杀死其天敌益虫，这将导致害虫 量的增加。 六、 模型检验 针对一组具体的数据用 Matlab 软件进行计算. 设老鼠和猫的初始数量分别为 对于数据 t 的终值经试验后确定为 15，即模型为： 首先，建立 m-文件 shier.m 如下： function dx=shier(t,x) dx=zeros(2,1); dx(1)=x(1)*(1-0.1*x(2)); dx(2)=x(2)*(-0.5+0.02*x(1)); 其次，建立主程序 shark.m 如下： [t,x]=ode45(shier,[0 15],[25 2]); plot(t,x(:,1), -,t,x(:,2),*) plot(x(:,1),x(:,2)) 数值解如下图： 为实线， 为“*”线 相图 为 由上两图知：x(t)与y(t)都是周期函数 模型（二） 考虑人工捕获 设表示捕获能力的系数为 K，相当于老鼠的自然增长率 由 a 降为a-K，猫的自然死亡率由 b增为 b+K 仍取 设战前捕获能力系数 K=0.3, 战争中降为 K=0.1, 则战前与战争中的模型分别为: 模型求解: 首先,建立 m-文件 shier1.m function dx=shier1(t,x) dx=zeros(2,1); dx(1)=x(1)*(0.7-0.1*x(2)); dx(2)=x(2)*(-0.8+0.02*x(1)); 和 shier2.m function dx=shier2(t,x) dx=zeros(2,1); dx(1)=x(1)*(0.9-0.1*x(2)); dx(2)=x(2)*(-0.6+0.02*x(1)); 定义上述两个方程 建立主程序 shark1.m [t,x]=ode45(shier1,[0 15],[25 2]); plot(t,x(:,2)./(x(:,1)+x(:,2)), -) hold on [t,x]=ode45(shier2,[0 15],[25 2]); plot(t,x(:,2)./(x(:,1)+x(:,2)), g*); 求解两个方程，并画出两种情况下猫数在毛和老鼠总数中所占比例 林奇瀚 F0907103 5090719075