请教用matlab代码怎么写

Question

 

Answer 1

%将时间离散化，时间间隔为dt
n=500;          %总时刻数
x=zeros(4,n);   %产生4行n列个零来放置坐标矩阵
y=zeros(4,n);
dt=0.02;        %时间间隔
v=10;           %取速度v=10
%四个人初始位置的横,纵坐标:这里取距离a=100
x(1,1)=100; y(1,1)=0;
x(2,1)=0;   y(2,1)=0;
x(3,1)=0;   y(3,1)=100;
x(4,1)=100; y(4,1)=100;

for j=1:n-1
    for i=1:3   %前三个人的位置更新:
        d=sqrt((x(i+1,j)-x(i,j))^2+(y(i+1,j)-y(i,j))^2);    %第i个人和第i+1个人之间j时刻的距离
        dircos=(x(i+1,j)-x(i,j))/d;   %第i个人j时刻运动方向的cos和sin值
        dirsin=(y(i+1,j)-y(i,j))/d;
        x(i,j+1)=x(i,j)+v*dt*dircos;    %第i个人j+1时刻的位置的横,纵坐标
        y(i,j+1)=y(i,j)+v*dt*dirsin;
    end
                %第四个人的位置更新:
    d=sqrt((x(1,j)-x(4,j))^2+(y(1,j)-y(4,j))^2);    %第4个人和第1个人之间j时刻的距离
    dircos=(x(1,j)-x(4,j))/d;   %第4个人j时刻运动方向的cos和sin值
    dirsin=(y(1,j)-y(4,j))/d;
    x(4,j+1)=x(4,j)+v*dt*dircos;    %第4个人j+1时刻的位置的横,纵坐标
    y(4,j+1)=y(4,j)+v*dt*dirsin;
    
                %给第j时刻四个人的位置作图:
    plot(x(1,j),y(1,j),'r*',x(2,j),y(2,j),'y*',x(3,j),y(3,j),'b*',x(4,j),y(4,j),'g*')
    hold on
    pause(0.001)  %暂停0.001秒
end