Logo GenDocs.ru

Поиск по сайту:  


Загрузка...

Лекции по моделированию систем - файл 10_Псевдослучайные числа, процедуры их получения.doc


Лекции по моделированию систем
скачать (1799.4 kb.)

Доступные файлы (18):

10_Псевдослучайные числа, процедуры их получения.doc127kb.31.03.2005 20:30скачать
11_Моделирование случайных воздействий.doc302kb.21.04.2005 22:59скачать
12_Приближенные способы преобразования.doc279kb.22.04.2005 01:52скачать
13_Имитационное моделирование.doc505kb.28.04.2005 15:43скачать
14_Характеристики мод-х систем и типовые схемы.doc1157kb.04.05.2005 23:18скачать
15_Планирование экспериментов.doc236kb.12.05.2005 16:07скачать
1_введение.doc207kb.07.01.2005 19:19скачать
1_общ_вопр_мод.DOC105kb.26.01.2005 11:17скачать
2_матем_мет_мод.doc89kb.22.02.2005 11:16скачать
3_Сетевые модели.doc21kb.08.02.2005 18:48скачать
6_Системы массового обслуживания.doc234kb.02.03.2005 23:51скачать
7_Сетевые модели Сети Петри.doc264kb.11.03.2005 10:17скачать
8_Обощенные модели А-схемы.doc206kb.18.03.2005 01:16скачать
9_Концептуальные, алгоритмические, статические модели.doc90kb.25.03.2005 13:09скачать
P-схемы.doc137kb.24.02.2005 22:48скачать
Модели данных.doc26kb.08.02.2005 14:10скачать
Непрерывно детерминированные модели.doc58kb.22.02.2005 17:07скачать
Сетевые модели.doc379kb.08.02.2005 18:42скачать

содержание
Загрузка...

10_Псевдослучайные числа, процедуры их получения.doc

Реклама MarketGid:
Загрузка...
Лекция № 10

ПСЕВДОСЛУЧАЙНЫЕ ПОСЛЕДОВАТЕЛЬНОСТИ И ПРОЦЕДУРЫ ИХ МАШИННОЙ ГЕНЕРАЦИИ

При статистическом моделировании систем одним из основных вопросов является учет стохастических воздействий. Количество случайных чисел, используемых для получения статистически устой­чивой оценки характеристики процесса функционирования системы S при реализации моделирующего алгоритма на ЭВМ. Количество случайных чисел колеблется в достаточно широких пределах в зависимости от:

1 класса объекта моделирования;

2.вида оцениваемых характеристик;

3необходимой то­чности и достоверности результатов моделирования. Результаты статистического моделирования существенно за­висят от качества исходных (базовых) последовательностей случай­ных чисел.

На практике используются три основных способа генерации случайных чисел:

  • аппаратный (физический);

  • табличный (файловый);

  • алгоритмический (программный).

Аппаратный способ. Генерация случайных чис­ел вырабатываются специальной электронной приставкой — гене­ратором (датчиком) случайных чисел,— служащей в качестве одно­го из внешних устройств ЭВМ. Реализация этого способа генерации не требует дополнительных вычислительных опе­раций ЭВМ по выработке случайных чисел, а необходима только операция обращения к внешнему устройству (датчику). В основе лежит физический эффект, лежащего в основе таких генераторов чисел, чаще всего используются шумы в электронных и полупроводнико­вых приборах, явления распада радиоактивных элементов и т. д.

Достоинства:

Запас чисел не ограничен;

Расходуется мало операций;

He занимается место в памяти .

Недостатки:

Требуется периодическая проверка;

Нельзя воспроизводить последовательности;

Используется специальное устройство;

Необходимы меры по обеспечению стабильности.

Табличный способ. Случайные числа, представленные в виде таблицы, помещаются в память ЭВМ. Этот способ получения случайных чисел обычно используют при сравнительно небольшом объеме таблицы и файла чисел.

Достоинства:

Требуется однократная проверка;

Можно воспроизводить последовательности.

Недостатки:

Запас чисел ограничен;

Много места в ОЗУ;

Необходимо время для обращения к памяти.

^ Алгоритмический способ. Способ получения последовательности случайных чисел основанный на формировании случайных чисел в ЭВМ с помощью специальных алгоритмов и реализующих их программ. Каждое случайное число вычисляется с помощью соответствующей программы по мере возникновения потребностей при моделирова­нии системы на ЭВМ.

Достоинства:

Требуется однократная проверка;

Многократная воспроизводимость последовательности чисел;

Мало места в памяти и нет внешних устройств.

Недостатки:

Запас чисел ограничен периодом последовательности;

Затраты машинного времени.

Программная имитация случайных воздействий сводится к генерированию некоторых стандартных процессов и их последующего функционального преоб­разования. В качестве базового может быть принят любой удобный для моделирования конкретной системы S про­цесс (например, пуассоновский поток при моделировании Q-схемы). При дискретном моделирований базовым процессом явля­ется последовательность чисел , которые представляют реализации независимых, равномерно распределенных на ин­тервале (0, 1) случайных величин . В стати­стических терминах - повторная выборка из равномерно распре­деленной на интервале (0, 1) генеральной совокупности значений величины .

Непрерывная случайная величина  имеет равномерное рас­пределение в интервале (а, Ь), если ее функции плотности (а) и функция распределения (б) примет вид (Рис. 1):







Рис.1




Числовые характеристики случайной величины , принимающей значения х— это математическое ожидание, дисперсия и среднее квадратическое отклонение соответственно:

При моделировании систем на с случайными числами интервала (0, 1), где границы интервала соответственно а=0 и б = 1. Частным случаем равномерного распределения является функция плотности и функция распределения, соответственно имеющие вид:


Такое распределение имеет математическое ожидание М [] = 1/2

и дисперсию D[] = 1/12.

Это распределение требуется получить на ЭВМ. Но получить его на цифровой ЭВМ невозможно, так как машина оперирует с п-разрядными числами. По­этому на ЭВМ вместо непре­рывной совокупности равно­мерных случайных чисел интервала (0, 1) используют дискретную последователь­ность 2" случайных чисел то­го же интервала. Закон рас­пределения такой дискрет­ной последовательности на­зывают квазиравномерным распределением.

Случайная величина , имеющая квазиравномерное распределе­ние в интервале (0, 1), принимает значения с вероят­ностями , .

Математическое ожидание и дисперсия квазиравномерной слу­чайной величины соответственно имеют вид







На ЭВМ невозможно получить идеальную последовательность случайных чисел хотя бы потому, что на ней можно оперировать только с конечным множеством чисел. Кроме того, для получения значений х случайной величины  используются формулы (алгоритмы). Поэтому такие последовательности, являющиеся по своей сути детерминированными, называются псевдослучайными.

^ Требования к генератору случайных чисел.

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

  • состоять из квазиравномерно рас­пределенных чисел;

  • содержать статистически независимые числа;

  • быть воспроизводимыми;

  • иметь неповторяющиеся числа;

  • получаться с минимальными затратами машинного времени;

  • занимать ми­нимальный объем машинной памяти.

В практике моделирования применяются генерации последовательностей псевдослучайных чисел находят алгоритмы вида (1)

(1)
Данные алгоритмы представляют рекуррентные соотношения первого порядка, для которых начальное число х0 и постоянные параметры уже заданы.

Метод середин­ных квадратов

Пусть имеется 2n-разрядное число, меньшее 1:



1.Возведем его в квадрат:



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

Пример, если начальное число х0=0,2152, то (х0)2=0,04631104,

т. е. Xj=0,6311, затем (х1)2=0,39828721, т. е. х2=0,8287, и т. д.

Недостаток метода:

Наличие корреляции между числа­ми последовательности, в некоторых случаях может отсутствовать.

^ Конгруэнтные процедуры генерации.

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

Два целых числа  и  конгруэнтны или сравнимы по модулю m, m — целое число, тогда и толь­ко тогда, когда существует такое целое число k, что т. е. разность делится на m и числа  и  дают одинаковые остатки от деления на абсолютную величину числа m.

Например,









К



онгруэнтные процедуры являются чисто детерминированными, так как описываются в виде рекуррентного соотношения (1), и имеют вид.

где — неотрицательные целые числа.

Раскроем рекуррентное соотношение (2):





(3)



Если заданы начальные числа (3) последовательность целых чисел {Xi}, составленную из остатков от деления на М членов

последовательности



Таким образом, для любого i>=1 справедливо неравенство Xt<M, получится последовательность рациональных чисел из единичного интервала (0,1)

^ Мультипликативный метод.

Задается последовательность неотри­цательных целых чисел {Xt}, не превосходящих М, рассчитанных по формуле


(4)


т. е. это частный случай соотношения (2) при =0.

В силу детерминированности метода получаются воспроизводи­мые последовательности..

В машинной реализации наиболее удобна версия M=pg, где р — число цифр в системе счисления в ЭВМ; g -— число битов в ма­шинном слове. Тогда вычисление остатка от деления на М сводится к выделению g младших разрядов делимого. Преобразование целого числа Xt в рациональную дробь из интервала осуществляется подстановкой слева от Xi двоичной или десятичной запятой.

Алгоритм построения последовательности для двоичной маши­ны M=pg сводится к выполнению таких операций:

1. Выбрать в качестве X0 произвольное нечетное число.

2. Вычислить коэффициент где t — любое целое поло­жительное число.

3. Найти произведение , содержащее не более 2g значащих
разрядов.

4. Взять g младших разрядов в качестве первого члена последовательности X1 а остальные отбросить.

5. Определить дробь из интервала (0, 1).

6. Присвоить .

7. Вернуться к п. 3.

Смешанный метод.

Позволяет вычислить последовательность неотрицательных целых чисел {Xi}, не превосходящих М, по формуле



Отличием от мультипликативного метода является .

С вычислительной точки зрения смешанный метод генерации сложнее мультип­ликативного на одну операцию сложения. При этом возможность выбора дополнительного параметра позволяет уменьшить возмож­ную корреляцию получаемых чисел.


Скачать файл (1799.4 kb.)

Поиск по сайту:  

© gendocs.ru
При копировании укажите ссылку.
обратиться к администрации
Рейтинг@Mail.ru