- 6
- 0
- 约2.48千字
- 约 7页
- 2019-07-05 发布于北京
- 举报
数学建模实验报告大小圆追击问题
数学建模告实实实实一,实实实目
第一:一个小在大内无滑的,用实实实实实实实实实实实实实实实实实实实Matlab实实程探小上或小内一点的运迹。实实实实实实实实实实实实实实实实
实实分析:
由目条件可知,小在大内无滑,实实实实实实实实实实实实实实实实实实
所以于小与大的初始接触点来,实实实实实实实实实实实实实实实实实
点滑的弧即小的路程。由此便实实实实实实实实实实实实实实实实实
可以系小与大。首先,分析目点。实实实实实实实实实实实实实实实实
若小上一点行察,于任意两点,实实实实实实实实实实实实实实实实实
从它不再与大接触到再一次与大接触,实实实实实实实实实实实实实
其运迹必然是相同的,迹弧的实实实实实实实实实实实实实实实实实实
形状,初末位置不做要求,。所以,我可以以初始实实实实实实
位置的接触点来考察运迹。若小实实实实实实实实实实实实实实实实实
内,非心,一点行考察,其迹必与实实实实实实实实实实实实实实实实
其在小上的点相似。故可以先算出实实实实实实实实实实实实实实实实
小上点的迹,然后放移即可。实实实实实实实实实实实实实实实实实实
若小心行考察,易知其迹必一实实实实实实实实实实实实实实实实实实实
实实实实实实实实实实实实实。分析完成后,下面始迹的算方法。
由于除心外,迹是似于的一条曲实实实实实实实实实实实实实实实实实
所以用参数方程行求解。具体做法是以大实实实实实实实实实实实实
实实实实实实实实实实实实实实实实实心原点,建立一个直角坐系。由
于小的初始位置与小上或小内一点实实实实实实实实实实实实实实实实
的运迹无,所以,我不妨小与实实实实实实实实实实实实实实实实实
大初始位置的接触点在实实实实实实实实实实实X实实实实实实实正半上。之后再
以小的心建立一个直角坐系。于实实实实实实实实实实实实实实实实实
个坐系里的任一点,可以利用向量实实实实实实实实实实实实实实
的算方法,由于大心和小心的实实实实实实实实实实实实实实实实实实实
向量是可求的,所以可将点投影到以大心实实实实实实实实实实
建立的坐系中。在算点的参数方程,可使用公式实实实实实实实实实实实实实实实实实实实实实
l=Ra,a实实实实实实实实实实实实实心角。先求出点在以小心原点
的坐系中的参数方程,然后投影至实实实实实实实实实实实实实实
另一坐系中即可。若察的是上一点,实实实实实实实实实实实实实实实
实实实实实实实实实实实实实只需将原程序中的小半径成点到小心
的距离即可。具体算程如下:实实实实实实实
Matlab代:实实
h=0.1;
q=0:h:20*pi;
R=5;r=3; d=2;
v=-R/r;
xR=R*cos(q);
yR=R*sin(q);
xrh=(R-r)*cos(q);
yrh=(R-r)*sin(q);
xr=xh+d*cos(v*q+q);
yr=yh+d*sin(v*q+q);
plot(xh,yh,xx,yy,xz,yz);
实实实实实实实实实实实果截:色大,半径10,色小心迹,实实实实实实实实实
小半径实实实实3,色小上一点迹。实实实实实实实实实实实
实实实实实实实色大,半径5,色小心迹,小半径实实实实实实实实实实实实实实实3,实
色小内一点迹,到小心距离实实实实实实实实实实实实实实实实实2。
第二:狗追人。人按一固定迹运,实实实实实实实实实实实实实实实实实
狗从某一点追逐人,求狗的运迹。实实实实实实实实实实实实实
实实实实实实实实实实实实实实分析:本采用以直代曲的方法求解。
先假人的迹是一个,并确定人,狗实实实实实实实实实实实实实实实实
的运速度以及人,狗的初始位置。实实实实实实实实实实实实实实
然后假定在某段极小的隔内,人的速度大实实实实实实实实实实实
小和方向均不。由于人的位置可求,狗实实实实实实实实实实实实
的当前位置也已得到,于是可以算出狗的速度方向。之后,
以直代替曲,将实实实实实实实dt内狗的运迹画出。并保实实实实实实实实
存狗在的位置,下一次做准。如实实实实实实实实实实实实实实实实实
此循,直至不足循条件。实实实实实实实实实实实实
Matlab代:实实
d=4;
R=5;
x_dog=10;
y_dog=0;
x_man=5;
y_man=0;
v_dog=1;
w_man=0.5;
dt=0.05;
t=0;
hold on;
axis([-10,10,-10,10]);
grid;
title(×??????????aí?);
while(t=50)
t=t+dt;
x_man=R*cos(w_man*t);
y_man=R*sin(w_man*t);
plot(x_man,y_man,r--,x_dog,y_dog,g--);
hold on;
x_dog=x_dog+v_dog*dt*
原创力文档

文档评论(0)