Create a multinomial distribution object using the default parameter values.

pd = makedist('Multinomial')

pd = MultinomialDistribution Probabilities: 0.5000 0.5000

创建一个多项分布对象说tribution with three possible outcomes. Outcome 1 has a probability of 1/2, outcome 2 has a probability of 1/3, and outcome 3 has a probability of 1/6.

pd = makedist('Multinomial','probabilities',[1/2 1/3 1/6])

pd = MultinomialDistribution Probabilities: 0.5000 0.3333 0.1667

Generate a random outcome from the distribution.

rng('default');% for reproducibilityr = random(pd)

r = 2

The result of this trial is outcome 2. By default, the number of trials in each experiment, $$n$$, equals 1.

Generate random outcomes from the distribution when the number of trials in each experiment, $$n$$, equals 1, and the experiment is repeated ten times.

rng('default');% for reproducibilityr = random(pd,10,1)

r =10×12 3 1 3 2 1 1 2 3 3

Each element in the array is the outcome of an individual experiment that contains one trial.

Generate random outcomes from the distribution when the number of trials in each experiment, $$n$$, equals 5, and the experiment is repeated ten times.

rng('default');% for reproducibilityr = random(pd,10,5)

r =10×52 1 2 2 1 3 3 1 1 1 1 3 3 1 2 3 1 3 1 2 2 2 2 1 1 1 1 2 2 1 1 1 2 2 1 2 3 1 1 2 3 2 2 3 2 3 3 1 1 2

Each element in the resulting matrix is the outcome of one trial. The columns correspond to the five trials in each experiment, and the rows correspond to the ten experiments. For example, in the first experiment (corresponding to the first row), 2 of the 5 trials resulted in outcome 1, and 3 of the 5 trials resulted in outcome 2.