2013082 9猎狗追击问题.doc

猎狗追击问题 邹奥东3110101152问题描述：一只兔子从原点以速度匀速向轴正方向运动，在外部有一只猎狗在追击兔子，猎狗的追击速度是匀速的，在坐标为的地方是兔子的洞穴所在地，兔子只要跑到了洞穴就算安全。第一问是要求画出追击轨迹，第二问是要求画出可以追到兔子的猎狗的出发地的集合（理论上面验证是一个椭圆，而事实画出来的也是）。第一问问题解决：1. 这是一个追击问题：假如当猎狗的初始出发地点是，理论上面解出来的猎狗的运动轨迹是：公式 1。没错，这个方程极其复杂，而这还只是一个区域上取定的一个初始出发点的解答，而我们要解出这个区域上理论上所有解的样子是几乎不可能的。所以我们要采用计算方法中的欧拉方法来近似模拟猎狗的运动。2.欧拉方法的描述：我们把时间划分为一段一段非常小的时间段，那么在的时间内猎狗的运动可以近似的看作是匀速直线运动了一小段，假设在现在这个时刻兔子的位置坐标为，狗的位置坐标是，那么在时间过后狗的位置坐标变成了。（其中，公式 2）那么我们把总时间都分成这么一小段一小段时间后并且如此处理狗的运动坐标，那么轨迹就可以很好的模拟出来。程序代码：x1=0;y1=0; %兔子的初始位置；x2=15;y2=7; b=19; %狗的初始位置；XX1=x1;y1; %存储兔子的初始位置；XX2=x2;y2; %存储狗的初始位置；delt_t=0.005;a=0.7; %兔子和狗的速度比为a； w=22; %狗的速度为22；i=0;n=350;for i=0:n x1=0; t=i*delt_t; if(x1-x2)2+(y1-y2)20.01)|(y1-b)20.01) break; endy1=a*w*t; %兔子的瞬时位置 d=sqrt(x1-x2)2+(y1-y2)2); x2=x2+w*delt_t*(x1-x2)/d; y2=y2+w*delt_t*(y1-y2)/d; %狗的瞬时位置XX1=XX1,x1;y1; XX2=XX2,x2;y2; %不断更新存储位置 plot(XX1(1,:),XX1(2,:),XX2(1,:),XX2(2,:),linewidth,2); axis(-20,40,0,35); %设置观摩所坐标轴位置 pause(0.02) %形成动态图 i=i+1;enddisp(总时间)i*delt_t 程序截图（是动态过程）：图表 1 图表 2第二问问题解答：无需多说的是第二问的解答肯定是在第一问的基础上面进行，第一问已经很好的模拟出来了定常猎狗初始位置的情况下的追击曲线。那么我们第二只要在一个区域内充分的选取到所有点来做初始位置并且判断是否能够到达即可。然后如果可以达到，我们就在图画中标出这点。理论上的边界是：。我们用数值方法和理论值做比较的结果是：程序代码：x1=0;y1=0;w=22;p=100; %x2=15;y2=7; b=5;a=0.7;n=20; %delt_t=b/(p*a*w);for e=-2*b:0.2:2*b for j=-2*b:0.2:4*b x20=e;x2=x20; y20=j;y2=y20;w=22;%5k=1; for i=0:p x1=0; t=i*delt_t; y1=a*w*t;% d=sqrt(x1-x2)2+(y1-y2)2); x2=x2+w*delt_t*(x1-x2)/d; y2=y2+w*delt_t*(y1-y2)/d;% if(x1-x2)2+(y1-y2)20.01) break; end if(abs(y1-b)0.01) break; end endif (x1-x2)2+(y1-y2)20.01) plot(x20,y20,.);pause(0.001); hold on; axis(-40,40,-40,40); end end end程序截图：图表 3图表 4图表 5图表 6图表 7图表 8图表 9程序总结：这个程序和初值的关系非常大，之前取得初值点的范围是：10000个以上的点，然后每个点都要算很久，然后有没有设置动态图像，所以要等10分钟以上还不出结果。而且精度一开始不适合设置过于精确，因为我一开始设置的是0.0001以上的精度，基本上估计30分钟都不够他算的，所以适合一开比较粗糙一点的取点运算，然后再精确一点点地去多取一点。然后再通过和文献所给的标准椭圆做对比的时候发现总会在上部多上一点点，而我这个时候的精度是设置的0.01，而多出的那一点点从直观上面来看已经超出了0.01.我觉得有两个解释：第一个是plot命令在