数学模型之慢跑者与狗.doc

慢跑者与狗 2、一个慢跑者在平面上沿椭圆以恒定的速率v=1跑步, 设椭圆方程为:x=10+20cost, y=20+5sin t．突然有一只狗攻击他．这只狗从原点出发, 以恒定速率w 跑向慢跑者, 狗的运动方向始终指向慢跑者．分别求出w=20,w=5时狗的运动轨迹． 1．模型建立 设t 时刻慢跑者的坐标为(X(t，Y(t，狗的坐标为(x(t，y(t．则X=10+20cos t, Y=20+15sin t. 狗从(0,0出发, 与导弹追踪问题类似，狗的运动轨迹的参数方程为: 2．模型求解 (1 w=20时, 建立Ｍ文件eq3．m 如下: function dy=eq3(t,y dy=zeros(2,1; dy(1=20*(10+20*cos(t-y(1/sqrt ((10+20*cos(t-y(1^2+(20+15*sin(t-y(2^2; dy(2=20*(20+15*sin(t-y(2/sqrt ((10+20*cos(t-y(1^2+(20+15*sin(t-y(2^2; 取t0=0，tf=10，建立主程序chase3．m 如下： t0=0;tf=10; [t,y]=ode45('eq3',[t0 tf],[0 0]; T=0:0．1:2*pi; X=10+20*cos(T; Y=20+15*sin(T; plot(X,Y,'-' hold on plot(y(:,1,y(:,2,'*' 在chase3．m 中，不断修改tf 的值, 分别取tf=5, 2．5, 3．5, …, 至3．15时, 狗刚好追上慢跑者． t0=0;tf=3.15; [t,y]=ode45('eq3' T=0:0.1:2*pi; X=10+20*cos(T; d 20cos d d 15sin d (00, (00 x t x t y t y t x y ?=+-????=+-???==? Y=20+15*sin(T; plot(X,Y,'-' hold on plot(y(:,1,y(:,2,'*' (2 w=5时 建立M 文件eq4．m 如下: function dy=eq4(t,y dy=zeros(2,1; dy(1=5*(10+20*cos(t-y(1/sqrt ((10+20*cos(t-y(1^2+(20+15*sin(t-y(2^2; dy(2=5*(20+15*sin(t-y(2/sqrt ((10+20*cos(t- y(1^2+(20+15*sin(t-y(2^2; 取t0=0，tf=10，建立主程序chase4．m 如下： t0=0;tf=10; [t,y]=ode45('eq4',[t0 tf],[0 0]; T=0:0．1:2*pi; X=10+20*cos(T; Y=20+15*sin(T; plot(X,Y,'-' hold on plot(y(:,1,y(:,2,'*' 在chase3．m 中，不断修改tf 的值, 分别取tf=20, 40, 80,…, 可以看出, 狗永远追不上慢跑者 t0=0;tf=200; [t,y]=ode45('eq4' ,[t0 tf],[0 0]; T=0:0.1:2*pi; X=10+20*cos(T; Y=20+15*sin(T; plot(X,Y,'b+',y(:,1,y(:,2,'*'