数学建模-猎狗追兔子问题

1、数学建模论文数学建模（2014春）课程期末论文题目猎狗追兔子问题题号A成员姓名学号班级学院联系学生1学生2摘要（一）对于问题一：自然科学中存在许多变量，也有许多常量，而我们要善于通 过建立合适的模型找到这些变量之中的不变量。猎狗追赶兔子的问题是我们在生活中常见的实例，而题目把我们生活中的普通的例子抽象成为高等数学中微分方程的例子，通过对高阶微分方程的分析，建立微分方程模 型，并用数学软件编写程序求解，得出结论，解决生活中常见的实际问题。（二）对于问题二：学习使用matlab进行数学模型的求解，掌握常用计算机软件的使用方法。关键词微分方程导数的几何意义猎狗追兔子数学建模数学软件一、问题重述如图1

2、所示，有一只猎狗在 B点位置，发现了一只兔子在正东北方距离它250m的地方O处, 此时兔子开始以8m/s的速度正向正西北方向，距离为 150m的洞口 A全速品&去.假设猎狗在追赶兔子的时候，始终朝着兔子的方向全速奔跑。请回答下面的问题：猎狗能追上兔子的最小速度是多少？在猎狗能追上兔子的情况下，猎狗跑过的路程 是少？ 假设猎狗在追赶过程中，当猎狗与兔子之间的 距离为30m寸，兔子由于害怕导致奔跑速度每秒减半， 而狗却由于兴奋奔跑速度每秒增加0.1倍，在这种情况下回答前面两个问题。二、问题分析与假设在猎狗追赶兔子的时候猎狗一直朝着兔子的方向追赶，所以可以建立平面直角坐标 系，通过导数联立起

3、猎狗运动位移，速度和兔子的运动状态。1,假设兔子的运动是匀速的。2.假设猎狗的运动轨迹是一条光滑并且一阶导数存在的曲线。3,猎狗的运动时匀速或者匀变速的。4.猎狗运动时总是朝向兔子。三、模型的建立及求解3.1 符号规定1. (x, y)：猎狗或者兔子所在位置的坐标。2. t :从开始到问题结束经过的时间。3. a:猎狗奔跑的路程。4. v:猎狗的奔跑速度。3.2 模型一的建立与求解猎狗能够抓到兔子的必要条件：猎狗的运动轨迹 在。酸有交点以0朋y轴，以。的x轴建立坐标系，则由图有 0(0, 0) , A(0,150),B(250 , 0),兔子的初始位置 0 点，而猎狗初始位置是B点，t (s)

4、后猎狗到达了 C(x, y),而兔子到达了 D (0, 8t),则有CD勺连线是猎狗 运动轨迹的一条切线，由导数的几何意义有：dy y 8tdx xdav dtda Jdx2 dy2三式联立消去t,得到；设：若猎狗可以追上兔子则有当兔子在 OA,猎狗在OB之间运动时此方程有解，设:得到：得到：两式联立相加得到：1 .如果q=1即v=8 m/s得到所以此情况无交点，所以v=8m/s猎狗无法追上兔子；2 .如果q<1即v>8m/s得到此情况有交点，所以有可能能够追上兔子，如果要追上兔子需要y<=150;解得至1J:即所以这种情况下能够追上的最小速度是.3 .如果q>1利用上

5、式得到，所以这种情况不能追上兔子。综上讨论，猎狗可以追上兔子的最小速度为。3.3 模型二的建立与求解如果猎狗可以追上兔子那么猎狗的轨迹和兔子的轨迹必相交与一点，此时兔子的路程 ，所用放的时间，那么猎狗的的路程a=tv;带入数值解得a=。3.4 模型三的建立与求解模型三利用matlab 试验，得到代码如下：a=8;dogxa=;dogya=;rabbitxa=;rabbitya=;d=1;dogx=250;dogy=0;rabbitx=0;rabbity=0;t=0;dt=0.001;for b=0:100dogx=250;dogy=0;rabbitx=0;rabbity=0;t=0;c=b;a

6、=8;while(sqrt(dogx-rabbitx)A2+(dogy-rabbity)A2)>d&rabbity<150)if(sqrt(dogx-rabbitx)A2+(dogy-rabbity)A2)<=30)b=b*1.1Adt;a=a*0.5Adt;endt=t+dt;dogx=dogx+b*dt*(rabbitx-dogx)/sqrt(dogx-rabbitx)A2+(dogy-rabbity)A2);dogy=dogy+b*dt*(rabbity-dogy)/sqrt(dogx-rabbitx)A2+(dogy-rabbity)A2);rabbitx=r

7、abbitx+0;rabbity=rabbity+a*dt;endif(rabbity<=150)b=c;break;endendfprintf(' 猎狗的最小速度是：:%2f',b);a=8;b=16;d=1;dogxb=;dogyb=;rabbitxb=;rabbityb=;dogx=250;dogy=0;rabbitx=0;rabbity=0;t=0;dt=0.001;s=0;while(sqrt(dogx-rabbitx)A2+(dogy-rabbity)A2)>d)t=t+dt;if(sqrt(dogx-rabbitx)A2+(dogy-rabbity)A

8、2)<=30)b=b*1.1Adt;a=a*0.5Adt;enddogx0=dogx;dogy0=dogy;dogx=dogx+b*dt*(rabbitx-dogx)/sqrt(dogx-rabbitx)A2+(dogy-rabbity)A2)dogy=dogy+b*dt*(rabbity-dogy)/sqrt(dogx-rabbitx)A2+(dogy-rabbity)A2)dogxb=dogxb,dogx;dogyb=dogyb,dogy;rabbitx=rabbitx+0;rabbity=rabbity+a*dt;rabbitxb=rabbitxb,rabbitx;rabbityb

9、=rabbityb,rabbity;s=s+sqrt(dogx0-dogx)A2+(dogy0-dogy)A2);endfprintf('最短路程是：%1f',s);得到猎狗的最小速度是：16m/s猎狗此时的路程是：312.5m四、模型的检验使用 matlab 进行计算机模拟实验检验模型的可行性： 问题一的检验：h=250;a=8;v=16;dogxb=;dogyb=;rabbitxb=;rabbityb=;d=0.01;dt=0.1;t=0;dogx=h;dogy=0;rabbitx=0;rabbity=0;while(sqrt(dogx-rabbitx)A2+(dogy-r

10、abbity)A2)>d&&t<=19.3) t=dt+t;dogx=dogx-v*dt*dogx/sqrt(dogxA2+(a*t-dogy)A2);dogy=dogy+v*dt*(a*t-dogy)/sqrt(dogxA2+(a*t-dogy)A2); dogxb=dogxb,dogx;dogyb=dogyb,dogy;rabbity=a*t;rabbityb=rabbityb,rabbity;end rabbitxb=zeros(length(rabbityb);plot(dogxb,dogyb,rabbitxb,rabbityb,'*')问题

11、二的模拟：n=250;a=8;v=16;d=0.1;dt=0.1;t=0;dx=n;dy=0;rx=0;ry=0;while(sqrt(dx-rx)A2+(dy-ry)A2)>d&&t<19.3) plot(dx,dy,rx,ry,'y*')pause(0.00001)hold ont=dt+t;dx=dx-v*dt*dx/sqrt(dxA2+(a*t-dy)A2);dy=dy+v*dt*(a*t-dy)/sqrt(dxA2+(a*t-dy)A2); ry=a*t;plot(dx,dy,rx,ry,'y*')endAle Edit V

12、i&w Insert Tools Desktop Window Help 4 口昌母国艮的夏 EIH H 口160i1111140 ；-120 -I100-吗80,-%60-40，.-20""-"" '' E.t|, -,.0£0100150200250五、模型的评价5.1 模型的优缺点模型的优点。(1)模型的使用围比较广泛，可以类推到其他许多模型中。(2)模型具有很高的使用价值。(3)模型对题目中的问题解决合适，模型使用得当。这里写模型的缺点。(4)题目中增加了一些理想化的假设，致使模型的波动比较大(5)不同兔子和猎狗的情况会有差异。5.2 模型的改进可使用仿生学原理，建立我们更加准确的模型。六、参考文献1书来，MATLA弱程与最优化问题，：电子工业，2013。2邹学军，周凯，宋军全，数学建模竞赛辅导教程，大学， 2009。3志林，欧宜贵，数学建模及其典型案例分析，化学工业，2006.4 Matlab 入门教程,wenku.baidu./view/daf8592fff00bed5b9f31d5d.htl 2014.06附录1: Matlab