用matlab实现在各向同性行走模型的基础上， 模拟各向异性的细胞迁移， 主要通过改变持 续时间P和迁移速度S实现。具体的参数设置情况参考如下：（1）Px=10，Sx=0.5，Py=5，Sy=0.2（2）Px=20，Sx=1.0，Py=5，Sy=0.2（3）Px=30，Sx=1.5，Py=5，Sy=0.2（4）Px=50，Sx=2.0，Py=5，Sy=0.2。矩阵维度必须一致。

以下是使用 Matlab 实现各向异性细胞迁移的代码，其中使用了一个名为 `walk` 的函数来模拟各向同性行走模型：

% 定义模拟参数
P = [10 20 30 50]; % 持续时间
S = [0.5 1.0 1.5 2.0]; % 迁移速度
Py = 5; % y方向速度和持续时间不变

% 定义模拟参数
dt = 0.01; % 时间步长
T = 1000; % 总时间
N = T / dt; % 时间步数
x = zeros(N,length(P)); % x坐标数组
y = zeros(N,length(P)); % y坐标数组

% 模拟各向异性细胞迁移
for i = 1:length(P)
    Px = P(i);
    Sx = S(i);
    for j = 2:N
        x(j,i) = x(j-1,i) + Sx*dt*randn();
        y(j,i) = y(j-1,i) + Sy*dt*randn();
        if rand() < dt/Px
            theta = rand()*2*pi;
            x(j,i) = x(j-1,i) + cos(theta);
            y(j,i) = y(j-1,i) + Py*sin(theta);
        end
    end
end

% 绘制轨迹图
figure
for i = 1:length(P)
    plot(x(:,i),y(:,i))
    hold on
end
xlabel('x')
ylabel('y')
legend('Px=10, Sx=0.5','Px=20, Sx=1.0','Px=30, Sx=1.5','Px=50, Sx=2.0')

代码中使用了 `randn()` 函数来生成均值为0，方差为1的正态分布随机数，模拟各向异性的细胞迁移。在每个时间步长内，细胞都有一定的概率停留在当前位置，概率为 `dt/Px`，同时有一定的概率沿某个方向移动，移动的方向是随机的。在 x 方向上，移动的速度是 `Sx`，在 y 方向上，移动的速度是 `Sy`。