我需要根据确定的概率密度函数(截断拉普拉斯分布)生成随机值:
Sigma_Phi = 1;
B = 1/(1-exp(-sqrt(2)*pi/Sigma_Phi));
Phi = -60:0.001:60;
for iter=1:length(Phi)
if Phi(iter) < pi && Phi(iter)>=0
P(iter) = ((B*exp(-abs(sqrt(2)*Phi(iter)/Sigma_Phi)))/(sqrt(2)*Sigma_Phi));
elseif Phi(iter) >= -pi && Phi(iter)<=0
P(iter) = ((B*exp(-abs(sqrt(2)*Phi(iter)/Sigma_Phi)))/(sqrt(2)*Sigma_Phi));
else
P(iter)=0;
end
end
我需要按照该概率密度函数生成随机值,但我不知道如何操作。