四只狗相向而行动态模拟matlab程序

% 例：四只狗相向而行
clear;clc;clf
a=[1,1];b=[0,1];c=[0,0];d=[1,0];         %四狗的位置（x,y）
t=0;               %初始时刻
dt=0.02;        %时间步长
v=1/60;  
k=0;              %第k个位置
d1=inf;            %相隔距离

%画出初始位置
hold on;
h1=plot(1,1,'r.');   
h2=plot(0,1,'.y');  
h3=plot(0,0,'k.');    
h4=plot(1,0,'.b'); 

grid off
axis([-0.5,1.5,0,1]);
set(h1,'erasemode','none','markersize',5);
set(h2,'erasemode','none','markersize',5);
set(h3,'erasemode','none','markersize',5);
set(h4,'erasemode','none','markersize',5);
while d1>1e-3

    k=k+1;               
    t=t+dt;               
    %b狗的位置
    w1=a-b;               
    d1=norm(w1);        
    w1=w1/d1;