莱维飞行源代码（Levy）

clc
clear
close all

n = 200;
m = 3; 

[x,y,z] = deal(zeros(n,3));

d = levy(n,m,1.4); 

for i = 2:length(d) 
    x(i)= x(i-1)+d(i,1);
    y(i)= y(i-1)+d(i,2);
    z(i)= z(i-1)+d(i,3);
end

line(x,y,z)
    
function [z] = levy(n,m,beta)

    num = gamma(1+beta)*sin(pi*beta/2); % 分子
    
    den = gamma((1+beta)/2)*beta*2^((beta-1)/2); % 分母

    sigma_u = (num/den)^(1/beta);% 标准偏差

    u = random('Normal',0,sigma_u^2,n,m); 
    
    v = random('Normal',0,1,n,m);

    z = u./(abs(v).^(1/beta));
     
end