从截断正态分布绘制伪随机

Question

x = 0.5 + 0.1*randn() 

从平均值0.5和标准差0.1的正态分布中得出伪随机数。

既然如此，下面的Matlab代码是否等价于从一个在0处截断的正态分布在1处进行采样？

    while x <=0 || x > 1

    x = 0.5 + 0.1*randn();

    end

Answer 1

使用MATLAB的概率分布对象使得从截断分布中进行采样变得非常容易。

您可以使用makedist()和truncate()函数定义对象，然后修改（截断它）为random()函数准备对象，该函数允许从它生成随机变量。

% MATLAB R2017a
pd = makedist('Normal',0.5,0.1)     % Normal(mu,sigma)
pdt = truncate(pd,0,1)              % truncated to interval (0,1)
sample = random(pdt,numRows,numCols)  % Sample from distribution `pdt`

一旦创建了对象（这里是pdt，即pd的截断版本），您就可以在各种函数调用中使用它。

figure, hold on
h = histogram(cr,'Normalization','pdf','DisplayName','@Luis Mendo samples');
X = 0:.01:1;
p = plot(X,pdf(pdt,X),'b-','DisplayName','Theoretical (w/ truncation)');