步骤1:
假设我想生成离散的均匀随机数,取值为-1或1.换句话说,我想生成具有以下分布的数字:
P(X = -1) = 0.5
P(X = 1) = 0.5
要生成这些数字中的100个数组,我可以编写以下代码:
n = 100
DV = [-1,1]; % Discrete value
RI = unidrnd(2,n,1); % Random uniform index
DUD = DV(RI); % Discrete uniform distribution
我的DUD数组看起来像:[ – 1,1,1,1,-1,-1,1,-1,…]
第2步:
现在我想生成10个等于和(DUD)的数字,因此10个数字的分布对应于离散均匀分布后的100个数的总和.
我当然可以这样做:
for ii = 1:10
n = 100;
DV = [-1,1]; % Discrete value
RI = unidrnd(2,n,1); % Random index
DUD = DV(RI); % Discrete uniform distribution
SDUD(ii) = sum(DUD);
end
同
SDUD =
2 2 -6 -2 -4 2 4 4 0 2
有没有数学/ matlab技巧呢?不使用for循环.
SDUD的直方图(10000个值,n = 100)如下所示:
奖金:
如果可以修改原始离散值,那将是很好的.因此,代替[-1,1],离散值可以是,例如,[0,1,2],每个都具有出现p = 1 / number_of_discrete_value,因此在该示例中为1/3.