绘制对数正态随机变量的100种不同实现

Question

我输出数字有问题。

pr = pt - plo - 10*gamma*log_base10(d/d0) - x_sigma

鉴于：pt = 0，plo = 55，gamma = 4.7，d0 = 1，sigma = 4.6。 这里，sigma是正态随机变量x_sigma的标准偏差。

对于每个距离值x_sigma，我必须对对数正态随机变量d进行100次不同的实现。

目标是在图表上绘制pr vs d。 这是我的代码，我不知道我在哪里犯了错误：

 Pt = 0;
 Pl0 = 55;
 gamma = 4.7;
 d0 =1;
 d= [1:20:2000]
 x_sigma = lognrnd(0, 4.6, 1, 100)


 Pr = Pt - Pl0 - 10.*gamma.* log10(d/d0) - x_sigma
 figure; plot(d, Pr); hold on;