用matlab写一个简单的模拟四个人在平面内移动的程序

下面是一个简单的模拟四个人在平面内移动的Matlab程序示例：

n = 500;         %总时刻数
x = zeros(4, n); %产生4行n列个零来放置坐标矩阵
y = zeros(4, n);
dt = 0.02;       %时间间隔
v = 10;          %速度

%四个人初始位置的横,纵坐标:
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

在这个程序中，首先定义了总时刻数n、坐标矩阵x和y、时间间隔dt和速度v等变量，并初始化了四个人的初始位置。然后，通过for循环，逐一更新每个人的位置，其中前三个人的位置更新时需要计算和更新运动方向的cos和sin值，而第四个人的位置更新则需要计算和更新与第一个人的距离和运动方向的cos和sin值。最后，通过plot函数将四个人在平面内的位置作图，并暂停一定时间，以模拟四个人的移动过程。

这个程序虽然简单，但是可以通过修改变量的值和增加移动规则等来实现更多的功能，比如将四个人移动的轨迹保存下来并进行分析等。同时，这个程序也展示了Matlab的绘图功能和for循环等基本语法的使用，对于初学者来说是一个不错的练手项目。