О компьютерном моделировании случайных величин

Статья: О компьютерном моделировании случайных величин

О компьютерном моделировании случайных величин

Непрерывная случайная величина  имеет равномерное распределение на отрезке , если ее функция распределения задается следующей формулой:

,

Плотность распределения вероятностей при этом имеет вид:

Математическое ожидание и дисперсия случайной величины  соответственно равны [3]:

, .

Обозначим буквой  случайную величину с равномерным распределением на отрезке . Для этой случайной величины функция распределения и плотность распределения вероятностей соответственно имеют вид:

 ,

Если , то вероятность

Моделировать случайную величину  можно многими способами [1].

Мы рассмотрим метод псевдослучайных последовательностей, который наиболее просто реализуется в компьютере. Для получения псевдослучайной последовательности используем алгоритм, который называется методом середины квадратов [4]. Поясним его на примере.

Возьмем некоторое число . Пусть  Возведем его в квадрат:  Выберем четыре средние цифры этого числа и положим  Затем возводим  в квадрат:  и снова выбираем четыре средние цифры. Получаем  Далее находим     и т. д.

Последовательность чисел  принимают за последовательность значений случайной величины  имеющей равномерное распределение на отрезке .

Для оценки степени приближения последовательности  к последовательности случайных чисел с равномерным распределением используют статистические критерии, например, аналогичные критерию, который используется в работе [2].

2. Моделирование последовательности независимых случайных испытаний

Пусть проводится последовательность  независимых испытаний. В результате каждого испытания может произойти одно из  несовместных событий  объединение которых совпадает с пространством элементарных событий . Известна вероятность появления каждого события , , которая не изменяется при переходе от одного испытания к другому. Очевидно, что .

Моделирование последовательности испытаний проводится следующим образом. Разделим отрезок  на  участков  длины которых соответственно равны  Получаем последовательность значений  случайной величины  Если , то считаем, что в -м испытании наступило событие , так как

.

3. Моделирование случайной величины дискретного типа

А. Общий алгоритм моделирования.

Если случайная величина  дискретная, то ее моделирование можно свести к моделированию независимых испытаний. В самом деле, пусть имеет место следующий ряд распределения:

Обозначим через  событие, состоящее в том, что случайная величина  примет значение , при этом . Тогда нахождение значения, принятого случайной величиной  в результате испытания, сводится к определению того, какое из событий  появится.

Так как события  несовместны и вероятность появления каждого из них не изменяется от испытания к испытанию, то для определения последовательности значений, принятых случайной величиной  можно использовать алгоритм моделирования последовательности независимых испытаний.

Б. Моделирование случайной величины с биномиальным распределением.

Случайная величина  считается распределенной по биномиальному закону, если

где ; — вероятность появления некоторого события  в каждом отдельно взятом испытании;  — вероятность появления события  в  независимых испытаниях  раз.

Введем случайную величину  — число появлений события  в -ом испытании,  Для этой величины имеет место:

, . (1)

Тогда случайное число  появлений события  в  испытаниях определяется по формуле

. (2)

Исходя из формул (1) и (2), значения случайной величины  определяются следующим образом:

1) находят последовательность значений  случайной величины

2) для каждого числа ,  проверяют, выполняется ли неравенство  если неравенство выполняется, то полагают  в противном случае считают

3) находят сумму значений  случайных величин  которая совпадает со значением

Повторяя этот алгоритм, получим последовательность значений  случайной величины с биномиальным законом распределения.

В. Моделирование случайной величины, распределенной по закону Пуассона.

Распределением Пуассона называется распределение вероятностей дискретной случайной величины, задаваемое формулой:

, ,

где — число событий простейшего потока, наступающих за некоторый промежуток времени. Распределение Пуассона применяется вместо биномиального распределения тогда, когда число  независимых испытаний велико (порядка нескольких сотен), а вероятность  появления события в каждом отдельно взятом испытании мала, при этом желательно, чтобы имело место .

Алгоритм моделирования случайной величины , распределенной по закону Пуассона при заданном параметре  можно представить следующим образом:

1) выбираем  таким образом, чтобы вероятность  была достаточно малой, например, меньше 0, 01;

2) получаем последовательность значений случайной величины , равномерно распределенной на отрезке ;

3) для каждого числа ,  проверяем, выполняется ли неравенство ; если это неравенство выполняется, то полагают , в противном случае считаем ;

4) вычисляем сумму  которая совпадает со значением случайной величины  распределенной по закону Пуассона.

4. Моделирование случайной величины

абсолютно непрерывного типа

А. Метод обратных функций.

Пусть случайная величина  имеет монотонно возрастающую функцию распределения . Известно, что  значит, случайная величина  с монотонно возрастающей функцией распределения  связана со случайной величиной  соотношением

.

Отсюда следует, что значение  случайной величины является решением уравнения

, (3)

где — значение случайной величины  т. е.

.

Последовательности значений  случайной величины соответствует последовательность  значений случайной величины  с функцией распределения .

Б. Моделирование случайной величины с равномерным распределением на отрезке .

Пусть случайная величина  имеет равномерное распределение на отрезке . Тогда ее функция распределения имеет вид:

.

Составим уравнение (3), получим

,

откуда

.

Последовательности значений  случайной величины  соответствует последовательность значений

, , …

случайной величины  равномерно распределенной на отрезке .

В. Моделирование случайной величины с показательным распределением.

Пусть случайная величина  имеет показательное распределение с параметром . Тогда функция распределения этой случайной величины

, .

Составим уравнение (3). Имеем

. (4)

Решаем уравнение (4) относительно  получаем

. (5)

Так как — случайная величина, равномерно распределенная на , то и  является также случайной величиной, распределенной по равномерному закону на отрезке . Поэтому вместо формулы (5) для моделирования случайной величины  можно использовать формулу

.

Г. Моделирование случайной величины с нормальным распределением.

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

,

где  и  — параметры.

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

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

приближается к нормальному. Требуется найти значения случайной величины  распределенной по нормальному закону с математическим ожиданием  и дисперсией .

Пусть — независимые случайные величины, равномерно распределенные на отрезке . Обозначим

. (6)

Учитывая   , найдем:

 .

При достаточно большом  можно считать, что случайная величина  имеет нормальный закон распределения с математическим ожиданием  и дисперсией .

Пронормируем случайную величину , получим:

. (7)

Для случайной величины  имеет место

, .

Перейдем от случайной величины  к стандартной нормально распределенной случайной величине

.

Тогда

.

Учитывая (6) и (7), получаем:

Например, при

.

Отсюда значение  случайной величины  определится по формуле

, (8)

где  — значения случайной величины , равномерно распределенной на отрезке .

Таким образом, имея 12 значений случайной величины  и подставляя их в формулу (8), получаем значение случайной величины  имея следующие 12 значений величины  и подставив их в формулу (8), получим следующее значение случайной величины  и т. д.

Список литературы

1. Калинина В.Н., Панкин В.Ф. Математическая статистика. М.: Высш. шк., 2001.

2. Кретов М.В. Вероятностные методы оценки прочности строительных материалов // Международная научная конференция «Инновация в науке и образовании—2003». Калининград, 2003. С. 228.

3. Кретов М.В. Теория вероятностей и математическая статистика. Калининград: Янтарный сказ, 2004.

4. Нейман Ю. Вводный курс теории вероятностей и математической статистики. М.: Наука, 1968.

Источник: http://refeteka.ru/r-101286.html

О компьютерном моделировании случайных величин

О компьютерном моделировании случайных величин

М.В.Кретов

1. Моделирование случайнойвеличины,  распределенной по равномерномузакону

Непрерывнаяслучайная величина  имеет равномерное распределение на отрезке , если еефункция распределения задается следующей формулой:

,

Плотность распределениявероятностей при этом имеет вид:

Математическоеожидание и дисперсия случайной величины  соответственно равны [3]:

, .

Обозначим буквой  случайную величину с равномернымраспределением на отрезке . Для этойслучайной величины функция распределения и плотность распределения вероятностейсоответственно имеют вид:

 ,

Если , товероятность

Моделироватьслучайную величину  можно многими способами [1].

Мы рассмотримметод псевдослучайных последовательностей, который наиболее просто реализуетсяв компьютере. Для получения псевдослучайной последовательности используемалгоритм, который называется методом середины квадратов [4]. Поясним его напримере.

Возьмем некоторое число . Пусть  Возведем его в квадрат:  Выберем четыре средние цифры этого числа иположим  Затем возводим  в квадрат:  и снова выбираем четыре средние цифры.Получаем  Далее находим     и т. д.

Последовательность чисел  принимают за последовательность значений случайнойвеличины  имеющей равномерное распределение на отрезке .

Для оценкистепени приближения последовательности  к последовательности случайных чисел с равномернымраспределением используют статистические критерии, например, аналогичныекритерию, который используется в работе [2].

2. Моделированиепоследовательности независимых случайных испытаний

Пусть проводитсяпоследовательность  независимых испытаний. В результате каждогоиспытания может произойти одно из  несовместных событий  объединение которых совпадает с пространствомэлементарных событий . Известнавероятность появления каждого события , , которая неизменяется при переходе от одного испытания к другому. Очевидно, что .

Моделированиепоследовательности испытаний проводится следующим образом. Разделим отрезок  на  участков  длины которых соответственно равны  Получаем последовательность значений  случайной величины  Если , то считаем,что в -м испытаниинаступило событие , так как

.

3. Моделированиеслучайной величины дискретного типа

А. Общий алгоритммоделирования.

Если случайнаявеличина  дискретная, то ее моделирование можно свести кмоделированию независимых испытаний. В самом деле, пусть имеет место следующийряд распределения:

Обозначим через  событие, состоящее в том, что случайнаявеличина  примет значение , при этом . Тогданахождение значения, принятого случайной величиной  в результате испытания, сводится к определениютого, какое из событий  появится.

Так как события  несовместны и вероятность появления каждого изних не изменяется от испытания к испытанию, то для определенияпоследовательности значений, принятых случайной величиной  можно использовать алгоритм моделированияпоследовательности независимых испытаний.

Б. Моделированиеслучайной величины с биномиальным распределением.

Случайная величина считается распределенной по биномиальномузакону, если

где ; — вероятностьпоявления некоторого события  в каждом отдельно взятом испытании;  — вероятность появления события  в  независимых испытаниях  раз.

Введем случайнуювеличину  — число появлений события  в -ом испытании, Для этой величины имеет место:

, . (1)

Тогда случайноечисло  появлений события  в  испытаниях определяется по формуле

. (2)

Исходя из формул(1) и (2), значения случайной величины  определяются следующим образом:

1) находятпоследовательность значений  случайной величины

2) для каждогочисла ,  проверяют, выполняется ли неравенство  если неравенство выполняется, то полагают  в противном случае считают

3) находят суммузначений  случайных величин  которая совпадает со значением

Повторяя этоталгоритм, получим последовательность значений  случайной величины с биномиальным закономраспределения.

В. Моделированиеслучайной величины, распределенной по закону Пуассона.

РаспределениемПуассона называется распределение вероятностей дискретной случайной величины,задаваемое формулой:

, ,

где — числособытий простейшего потока, наступающих за некоторый промежуток времени.Распределение Пуассона применяется вместо биномиального распределения тогда,когда число  независимых испытаний велико (порядканескольких сотен), а вероятность  появления события в каждом отдельно взятомиспытании мала, при этом желательно, чтобы имело место .

Алгоритммоделирования случайной величины ,распределенной по закону Пуассона при заданном параметре  можно представить следующим образом:

1) выбираем  таким образом, чтобы вероятность  была достаточно малой, например, меньше 0, 01;

2) получаемпоследовательность значений случайнойвеличины , равномернораспределенной на отрезке ;

3) для каждогочисла ,  проверяем, выполняется ли неравенство ; если этонеравенство выполняется, то полагают , в противномслучае считаем ;

4) вычисляем сумму которая совпадает со значением случайнойвеличины  распределенной по закону Пуассона.

4. Моделированиеслучайной величины

абсолютнонепрерывного типа

А. Метод обратныхфункций.

Пусть случайнаявеличина  имеет монотонно возрастающую функциюраспределения . Известно,что  значит, случайная величина  с монотонно возрастающей функциейраспределения  связана со случайной величиной  соотношением

.

Отсюда следует,что значение  случайной величины являетсярешением уравнения

, (3)

где — значениеслучайной величины  т. е.

.

Последовательностизначений  случайной величины соответствуетпоследовательность  значений случайной величины  с функцией распределения .

Б. Моделированиеслучайной величины с равномерным распределением на отрезке .

Пусть случайнаявеличина  имеет равномерное распределение на отрезке . Тогда еефункция распределения имеет вид:

.

Составим уравнение(3), получим

,

откуда

.

Последовательностизначений  случайной величины  соответствует последовательность значений

, , …

случайной величины равномерно распределенной на отрезке .

В. Моделированиеслучайной величины с показательным распределением.

Пусть случайнаявеличина  имеет показательное распределение с параметром. Тогдафункция распределения этой случайной величины

, .

Составим уравнение(3). Имеем

. (4)

Решаем уравнение(4) относительно  получаем

. (5)

Так как — случайнаявеличина, равномерно распределенная на , то и  является также случайной величиной,распределенной по равномерному закону на отрезке . Поэтомувместо формулы (5) для моделирования случайной величины  можно использовать формулу

.

Г. Моделированиеслучайной величины с нормальным распределением.

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

,

где  и  — параметры.

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

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

приближается кнормальному. Требуется найти значения случайной величины  распределенной по нормальному закону сматематическим ожиданием  и дисперсией .

Пусть — независимыеслучайные величины, равномерно распределенные на отрезке . Обозначим

. (6)

Учитывая   , найдем:

 .

При достаточнобольшом  можно считать, что случайная величина  имеет нормальный закон распределения сматематическим ожиданием  и дисперсией .

Пронормируемслучайную величину , получим:

. (7)

Для случайнойвеличины  имеет место

, .

Перейдем отслучайной величины  к стандартной нормально распределенной случайнойвеличине

.

Тогда

.

Учитывая (6) и(7), получаем:

Например, при

.

Отсюда значение  случайной величины  определится по формуле

, (8)

где  — значения случайной величины , равномернораспределенной на отрезке .

Таким образом,имея 12 значений случайной величины  и подставляя их в формулу (8), получаемзначение случайной величины  имея следующие 12 значений величины  и подставив их в формулу (8), получимследующее значение случайной величины  и т. д.

Список литературы

1. Калинина В.Н.,Панкин В.Ф. Математическая статистика. М.: Высш. шк., 2001.

2. Кретов М.В.Вероятностные методы оценки прочности строительных материалов // Международнаянаучная конференция «Инновация в науке и образовании—2003». Калининград, 2003.С. 228.

3. Кретов М.В.Теория вероятностей и математическая статистика. Калининград: Янтарный сказ,2004.

4. Нейман Ю.Вводный курс теории вероятностей и математической статистики. М.: Наука, 1968.

Для подготовкиданной работы были использованы материалы с сайта http://old.albertina.ru/

Дата добавления: 07.04.2010

Источник: https://www.km.ru/referats/1C83D21863854622A431D35B16FC76C9

Refy-free
Добавить комментарий