步骤1.定义分配的参数。

创建矢量p包含每个结果的概率。结果1具有1/2的概率，结果2具有1/3的概率，和结果3具有的1/6的概率。试验在每个实验数ñ是5，并且实验的重复次数代表8。

P = [1/2 1/3 1/6];N = 5;代表= 8;

第2步：生成一个随机数。

从生成的多项分布，这是一个单次试验的结果一个随机数。

RNG（'默认'）％用于重现R = mnrnd（1，P，1）

R =1×30 1 0

返回的向量[R包含三个元素，显示每个可能的结果计数。这种单一的试验导致的结果2。

第3步：生成随机数的矩阵。

您还可以生成随机数从多项分布，其中报道了多个实验，每个包含多个试验结果的矩阵。生成一个包含实验的结果矩阵N = 5试验和代表= 8重复。

R = mnrnd（N，P，代表）

R =8×31 1 3 3 2 0 1 1 3 0 4 1 5 0 0 1 2 2 3 1 1 3 1 1

在产生的矩阵中的每一行包含数为每个 $$ķ$$多项箱。例如，在第一个实验中（对应于第一行），这五个试验之一导致的结果1，五个试验之一导致结果2，和三个五个试验导致结果3。

第四步：计算的PDF文件。

由于多项功能与库计数工作，创造一切可能的结果组合的多维数组，并利用计算PDFmnpdf。

count1 = 1: n;是从= 1:n;(x1, x2) = meshgrid (count1,是从);x3 = n - (x1 + x2);y = mnpdf ([x1 (:), x2 (:), x3 (:)), repmat (p (n) ^ 2, 1));

第5步绘制的PDF文件。

创建3-d柱状图形象化为结果的频率的每个组合的PDF。

Y =重塑（Y，N，N）;BAR3（y）的组（GCA，'XTickLabel'，1：N）;集（GCA，'YTickLabel'，1：N）;xlabel（将“x_1频率”）ylabel（“X_2频率”）zlabel（“概率质量”）

该图显示了结果的每个可能的组合的概率质量。它不显示 $$X_3$$，这是由约束来确定 $$X_1+X_2+X_3=ñ$$。