莱维飞行的matlab代码

由于莱维飞行是一种随机游走过程，其代码实现方式可以有很多种，下面是其中一种实现方式的MATLAB代码：

% 设定初始位置和步长
x0 = 0;
y0 = 0;
step = 0.1;

% 设定迭代次数和初始方向
N = 5000;
theta = 2*pi*rand;

% 初始化位置数组
x = zeros(1, N+1);
y = zeros(1, N+1);
x(1) = x0;
y(1) = y0;

% 生成随机步长和方向
for i = 2:N+1
    r = step * sqrt(-2*log(rand));
    phi = 2*pi*rand;
    x(i) = x(i-1) + r*cos(theta+phi);
    y(i) = y(i-1) + r*sin(theta+phi);
end

% 绘制轨迹图
plot(x, y, 'LineWidth', 1.5);
axis equal;
grid on;
xlabel('X');
ylabel('Y');
title('Levy Flight');

运行以上代码，即可得到莱维飞行的轨迹图。需要注意的是，由于莱维飞行是一种长程依赖的随机过程，因此其迭代次数要比普通的随机游走要多一些，才能得到比较典型的结果。