猎狗奔跑的仿真图.m

来自「构图追击问题构图追击问题构图追击构图追击问题问题」· M 代码 · 共 22 行

clear all
close all
clc
 
a=8;b=19 %è á 1·μ  ù è ab=19￡¨ ú× μúò  ê ù èòa ó￡ 
jstxb=[];jstyb=[];zscxb=[];zscyb=[];
d=0.1;dt=0.05;t=0;
 
jstx=200;jsty=0;zscx=0;zscy=0;
while (sqrt((jstx-zscx)^2+(jsty-zscy)^2)>=d)
    t=t+dt;
    jstx=jstx-b*dt*jstx/sqrt(jstx^2+(a*t-jsty)^2);
    jstxb=[jstxb,jstx];
    jsty=jsty+b*dt*(a*t-jsty)/sqrt(jstx^2+(a*t-jsty)^2);
    jstyb=[jstyb,jsty];
    zscy=a*t;
    zscyb=[zscyb,zscy];
end
zscxb=zeros(length(zscyb));
plot(jstxb,jstyb,'o',zscxb,zscyb,'.');
size(zscyb)