lognrnd

Логарифмически нормальные случайные числа

Синтаксис

r = lognrnd(mu,sigma)
r = lognrnd(mu,sigma,sz1,...,szN)
r = lognrnd(mu,sigma,[sz1,...,szN])

Описание

пример

r = lognrnd(mu,sigma) генерирует случайное число от логарифмически нормального распределения с параметрами распределения mu (среднее значение логарифмических значений) и sigma (стандартное отклонение логарифмических значений).

пример

r = lognrnd(mu,sigma,sz1,...,szN) или r = lognrnd(mu,sigma,[sz1,...,szN]) генерирует sz1 -by-szN массивом логарифмически нормальных случайных чисел.

Примеры

свернуть все

Найдите параметры распределения от среднего значения и отклонения логарифмически нормального распределения и сгенерируйте логарифмически нормальное случайное значение от распределения.

Найдите параметры распределения mu и sigma от среднего значения и отклонения.

m = 1; % mean
v = 2; % variance
mu = log((m^2)/sqrt(v+m^2))
mu = -0.5493
sigma = sqrt(log(v/(m^2)+1))
sigma = 1.0481

Сгенерируйте логарифмически нормальное случайное значение.

rng('default') % For reproducibility
r = lognrnd(mu,sigma)
r = 1.0144

Сохраните текущее состояние генератора случайных чисел. Затем создайте вектор 1 на 5 логарифмически нормальных случайных чисел от логарифмически нормального распределения с параметрами 3 и 10.

s = rng;
r = lognrnd(3,10,[1,5])
r = 1×5
109 ×

    0.0000    1.8507    0.0000    0.0001    0.0000

Восстановите состояние генератора случайных чисел к s, и затем создайте новый вектор 1 на 5 случайных чисел. Значения те же, что и прежде.

rng(s);
r1 = lognrnd(3,10,[1,5])
r1 = 1×5
109 ×

    0.0000    1.8507    0.0000    0.0001    0.0000

Создайте матрицу логарифмически нормально распределенных случайных чисел с тем же размером как существующий массив.

A = [3 2; -2 1];
sz = size(A);
R = lognrnd(0,1,sz)
R = 2×2

    1.7120    0.1045
    6.2582    2.3683

Можно объединить предыдущие две строки кода в одну строку.

R = lognrnd(1,0,size(A));

Входные параметры

свернуть все

Среднее значение логарифмических значений для логарифмически нормального распределения, заданного как скалярное значение или массив скалярных значений.

Чтобы сгенерировать случайные числа от нескольких дистрибутивов, задайте mu и sigma с помощью массивов. Если и mu и sigma являются массивами, то размеры массивов должны быть тем же самым. Если или mu или sigma являются скаляром, то lognrnd расширяет скалярный аргумент в постоянный массив, одного размера в качестве другого аргумента. Каждый элемент в r является случайным числом, сгенерированным от распределения, заданного соответствующими элементами в mu и sigma.

Пример: [0 1 2; 0 1 2]

Типы данных: single | double

Стандартное отклонение логарифмических значений для логарифмически нормального распределения, заданного как неотрицательное скалярное значение или массив неотрицательных скалярных значений.

Если sigma является нулем, то вывод r всегда равен exp(mu).

Чтобы сгенерировать случайные числа от нескольких дистрибутивов, задайте mu и sigma с помощью массивов. Если и mu и sigma являются массивами, то размеры массивов должны быть тем же самым. Если или mu или sigma являются скаляром, то lognrnd расширяет скалярный аргумент в постоянный массив, одного размера в качестве другого аргумента. Каждый элемент в r является случайным числом, сгенерированным от распределения, заданного соответствующими элементами в mu и sigma.

Пример: [1 1 1; 2 2 2]

Типы данных: single | double

Размер каждой размерности, заданной как целые числа или вектор - строка из целых чисел. Например, определение 5,3,2 или [5,3,2] генерирует 5 3 2 массивами случайных чисел от распределения вероятностей.

Если или mu или sigma являются массивом, то заданные измерения sz1,...,szN должны совпадать с общими размерностями mu и sigma после любого необходимого скалярного расширения. Значения по умолчанию sz1,...,szN являются общими размерностями.

  • Если вы задаете одно значение sz1, то r является квадратной матрицей размера sz1.

  • Если размером какой-либо размерности является 0 или отрицательный, то r является пустым массивом.

  • После второго измерения lognrnd игнорирует последующие измерения с размером 1. Например, lognrnd (3,1,1,1) производит вектор 3 на 1 случайных чисел.

Типы данных: single | double

Выходные аргументы

свернуть все

Логарифмически нормальные случайные числа, возвращенные как скалярное значение или массив скалярных значений с размерностями, заданными sz1,...,szN. Каждый элемент в r является случайным числом, сгенерированным от распределения, заданного соответствующими элементами в mu и sigma.

Больше о

свернуть все

Логарифмически нормальное распределение

Логарифмически нормальное распределение является распределением вероятностей, логарифм которого имеет нормальное распределение.

Средний m и отклонение v логарифмически нормальной случайной переменной являются функциями логарифмически нормальных параметров распределения µ и σ:

m=exp(μ+σ2/2)v=exp(2μ+σ2)(exp(σ2)1)

Кроме того, можно вычислить логарифмически нормальные параметры распределения µ и σ от среднего m и отклонения v:

μ=журнал(m2/v+m2)σ=журнал(v/m2+1)

Альтернативная функциональность

  • lognrnd является функционально-специализированным к логарифмически нормальному распределению. Statistics and Machine Learning Toolbox™ также предлагает родовой функции random, который поддерживает различные распределения вероятностей. Чтобы использовать random, создайте распределение вероятностей LognormalDistribution, возражают и передают объект как входной параметр или задают имя распределения вероятностей и его параметры. Обратите внимание на то, что специфичный для распределения функциональный lognrnd быстрее, чем родовая функция random.

  • Чтобы сгенерировать случайные числа в интерактивном режиме, используйте randtool, пользовательский интерфейс генерации случайных чисел.

Ссылки

[1] Marsaglia, G. и В. В. Цанг. “Быстрый, Легко Реализованный метод для Выборки от Уменьшения или Симметричных Одномодовых Функций плотности”. SIAM Journal на Научном и Статистическом Вычислении. Издание 5, Номер 2, 1984, стр 349–359.

[2] Эванс, M., Н. Гастингс и Б. Пикок. Статистические Дистрибутивы. 2-й редактор, Хобокен, NJ: John Wiley & Sons, Inc., 1993.

Расширенные возможности

Представлено до R2006a