我正在尝试使用 MATLAB 中的 Levy 分布生成随机数。我正在使用逆变换方法生成随机数，然后绘制其直方图以进行 pdf 估计。然后我分析地绘制了 Levy 分布的 PDF。

问题：估计的 PDF 和分析的 PDF 不匹配。我将我的代码放在这里以供参考，其中我采用了位置参数 = 5 和比例参数 = 5。

clear all;
close all;
clc;
for loopr=1:1:10^8
X11(loopr)=5+5/(2*(erfcinv(rand))^2);
end
bin=[0:0.1:100];
l=0;
for loop1=5:0.1:100
l=l+1;
pdf12(l)=sqrt(5/(2*pi))*exp(-5/(2*(loop1-5)))/(loop1-5)^1.5;
end
loop1=5:0.1:100;
plot(loop1, pdf12)
hold on;
histogram(X11, bin, 'normalization', 'pdf');
hold off;

我得到的数字是：