额～
早上敲错一个步长～
所以算错了～
果然不是笨狗～
很快可以咬到人了～
运行下面程序可知：
clear
t=0;
d0=0.1;
r=0.01;
s=20;
gx=0;gy=0;
for i=1:1000
rx(i)=10+20*cos(t+r*(i-1));
ry(i)=20+5*sin(t+r*(i-1));
d(i)=sqrt((rx(i)-gx(i))^2+(ry(i)-gy(i))^2);
if d(i)<d0;
break
end
x=rx(i)-gx(i);
y=ry(i)-gy(i);
gx(i+1)=r*s*x/sqrt(x^2+y^2);
gy(i+1)=r*s*y/sqrt(x^2+y^2);
gx(i+1)=gx(i)+gx(i+1);
gy(i+1)=gy(i)+gy(i+1);
end
plot(gx,gy,rx,ry)
|