Файл: Шепелев, И. Г. Математические методы планирования и управления в строительстве конспект лекций.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.10.2024
Просмотров: 59
Скачиваний: 0
Для получения (генерирования) чисел | |
РСП |
(0,1) |
суще |
ствует несколько методов. |
|
(см. |
прило |
1.. Разрабатываются таблицы случайных чисел |
|||
жение V). Таблицы составляются с помощью |
какого-либо ге |
||
нератора случайных чисел, например, рулетки (метод |
назван |
«Монте-Карло» по имени города, где особенно распространена игра в рулетку), аппарата жеребьевки, волчка с пронумерован ными гранями, игральной кости и т. д. До появления ЭВМ. таб
лицы были наиболее рас |
&) |
||
пространенными |
источни- |
||
ками |
случайных |
чисел. |
|
При широком использова |
р |
||
нии ЭВМ табличный ме |
|
||
тод получения случайных |
|
||
чисел |
стал применяться |
|
|
реже, так как для таблиц |
|
||
требуется большой объем |
|
||
памяти |
электронных вы |
z p ( • { |
числительных машин.
2. Случайные числа вы рабатываются с помощью генератора случайных чи сел, преобразующего ре зультаты .случайного фи зического процесса в по следовательность двойных разрядов ЭВМ. Обычно в
качестве генератора используют собственные шумы электронных приборов или излучение радиоактивных веществ. Электронный прибор используется в- качестве приставки к ЭВМ. При этом говорят, что машина имеет физический генератор случайных чи сел, получаемых в неограниченном .'количестве, высокого каче ства с большой скоростью. Однако ЭВМ оборудуются датчиками случайных чисел только в том случае, если на ней предусматри вается регулярное решение задач с применением метода МонтеКарло.
3. Наиболее распространенным стал метод получения случайных чисел с помощью программы. Для этого создаются специальные программы, преобразующие первоначально задан ное число в некоторую последовательность чисел, которая Может рассматриваться как случайная. Например, задаются числом £о» из этого числа путем арифметических и логических операций по лучаются следующее число £i и т. д. Таким образом, источником получения -каждого последующего числа является предыдущее.
93
Так как последовательность чисел g получается в результа те программы, эти числа не являются «истинно» случайными, их называют псевдослучайными. Однако эти числа могут ус пешно применяться для моделирования равномерного случай ного процесса, если удовлетворяют требованиям случайности и равномерности.
Существуют специальные методы проверки псевдослучайных чисел на случайность и равномерность [2 , 3].
§ 7.3. Получение случайных величин, заданных законами распределения, с помощью чисел РСП (0,1)
Обычно реальные случайные величины редко задаются рав номерной случайной последовательностью РСП (0,1). Реаль ные случайные величины чаще всего задаются:
а) равномерной случайной последовательностью с предела ми от а до Ь;
б) нормальным законом распределения; в) экспоненциальным законам распределения;
г) пуассоновским законом распределения и другими зако нами.
Очень часто реальные распределения случайных величин не подчиняются вообще никаким законам и могут быть заданы гистограммой или таблицей, полученными в результате обра ботки статистических данных.
В § 8.4 |
рассмотрен |
пример задания |
случайного |
процесса |
|
цепями Маркова. |
перехода от распределения РСП |
(0,1) |
|||
Рассмотрим методы |
|||||
к любому |
другому распределению случайных чисел. |
Способ |
|||
перехода к |
равномерному распределению |
с пределами |
а и b |
рассмотрен в § 7.4. Для перехода к другим законам распреде ления применяют основное соотношение перехода: если случай ная величина г] имеет плотность распределения f (х), то случай
ная величина | равна |
|
E= Jf(jc)dJC |
(7.3.1) |
—СЛ |
|
и имеет равномерный закон распределения в интервале (0 ,1 ). Но обычно задача состоит не в нахождении g g РСП (0,1), а в нахождении гр Зная закон распределения f (х) и выработав g, необходимо взять определенный интеграл (7.3.1) и решить по лученное выражение относительно т). К сожалению этот метод, имеющий универсальное значение, очень редко дцет практиче
94
ские результаты, потому, что не все интегралы берутся, а если
даже интеграл возьмется, то соотношение решить трудно из-за предела—со в выражении (7.3.1).
Однако |
отдельные преобразования случайных |
величин |
|||||||
РСП (0 ,1 ) |
в случайные величины |
(х) |
возможны. |
||||||
Рассмотрим некоторые из них. |
|
|
|
|
|||||
а) Экспоненциальное распределение. |
|
|
|||||||
Если |
случайная величина |
задана |
экспоненциальным зако |
||||||
ном распределения |
|
|
если л > |
|
|
||||
|
|
|
( |
- , п — А х |
0 |
|
|||
|
|
|
J |
Хе |
|
(7.3.2) |
|||
|
|
|
f (х) = |
0 , |
если Л" |
0 , |
|
||
|
|
|
| |
|
|
||||
|
|
|
|
|
|
|
|
А х |
|
ТО |
|
|
J f (х) Ах — |
| |
0 • 6.x + j хе dx. |
(7.3.3) |
|||
Первый |
интеграл |
выражения |
(7.3.3) |
равен |
нулю. |
Взяв второй |
|||
интеграл и сделав |
соответствующие преобразования, получим |
||||||||
|
|
|
&= |
1 |
|
|
|
(7.3.4) |
|
тогда |
|
|
ri = - |
- ^ М |
1 |
|
|
(7.3.5) |
|
|
|
|
|
к |
|
|
|
|
так как число | есть случайное число, принадлежащее РСГ1 (0,1), то (1 — |) можно заменить Окончательная формула для преобразования равномерной случайной последовательно сти в экспоненциальную последовательность, заданную зако ном (7.3.2), имеет вид:
(7.3.6)
б) Распределение типа
|
X |
х. |
— е |
lz‘, при ■х > |
|
f(x) = |
сг2 |
|
0 , |
при х < |
0
(7.3.7)
0
Графики распределений а) и б) приведены на рис. 16.
с = |
0 + |
Г■e~*r— |
d x = - е ~ ^ , |
(7.3.8) |
|
|
] |
а 2 |
|
|
|
|
0 |
|
|
|
|
тогда. |
= |
а ]/"— 2 |
In (1 —- £), |
(7.3.9) |
|
или |
Tj |
= а |
— 2 In %. |
(7.3.10) |
/ 95
Подобно показанному можно аналитически выразить огра ниченное количество распределений случайных величин через случайные числа РСП (0,1). Если такое выражение получено, то для формирования случайной величины rig f (х) необходимо
S)
Рис. 16. Графики распределений случайной величины
выработать случайное |
число |
РСП |
(0,1), подставить это |
число в формулы типа |
(7.3.6) |
и (7.3.10) |
и провести необходи |
мые вычисления. |
|
|
|
§ 7.4 Получение случайных чисел приближенными методами
Если для получения случайных чисел не удается использо вать основное соотношение перехода (7.3.1), пользуются при ближенным методом.
Рассмотрим применение приближенных методов на примере
нормального распределения. Нормальное распределение |
|
Нх) = -----Г д - е 2за“. |
(7.4.1) |
3 1 Л2 к |
|
Это наиболее часто встречающееся распределение случайных величин. Однако применение основного соотношения' (7.3.1) не возможно из-за того, что в выражении
9 6
не удается избавиться от предела —<*>.
Основой получения чисел тр распределенных нормально из чисел £ 6 РСП (0,1), является центральная предельная теорема теории вероятностей, в соответствии с которой, если величина а
1=П
имеет любое распределение, то Pj—2 ai ПРИ достаточно боль- i=i
шом п распределена нормально. Тогда для получения числа 7 , распределенного нормально, с математическим ожиданием 0 и
стандартным |
отклонением единицы, |
пользуются формулой: |
||||
|
|
1=П |
|
|
|
|
|
|
|
1 ), |
|
|
(7.4.3) |
|
|
1=1 |
|
|
|
|
здесь п — количество суммируемых |
случайных |
чисел |
(2 ;* 1 — 1 ), |
|||
имеющих равномёрное распределение с параметрами |
(—1 ; 1 ), |
|||||
с математическим |
ожиданием 0 и стандартным отклонением |
|||||
о = —. Для |
получения удовлетворительных нормальных чисел ^ |
|||||
3 |
(7.4.3) |
необходимо, чтобы п ^ |
10, |
т. е. для получе |
||
по формуле |
||||||
ния каждого числа -( необходимо |
выработать |
и преобразовать |
||||
по формуле (7.4.3) не менее 10 чисел | £ РСП |
(0,1). |
Для со |
||||
кращения затрат времени на выработку | |
применяют формулы |
|||||
коррекции, позволяющие уменьшить п. |
|
|
|
|||
Одна из таких формул: |
|
|
|
|
||
|
|
‘к— Тк-----(ЗТк — |
|
|
(7.4.4) |
|
|
|
20п |
|
|
|
|
Если пользоваться этой формулой коррекции, то можно приме
нять п = 5. Однако числа |
принадлежат нормальной случай |
|||
ной последовательности |
с математическим |
ожиданием |
0 и |
|
стандартным отклонением |
единицы |
НСП |
(0,1). Для |
полу |
чения нормальных чисел с параметрами а и а необходимо дестандартизировать эти числа по формуле:
|
|
% = |
a -f a 7k, |
(7.4.5) |
здесь г]£ |
НСП (а, |
а). |
|
чисел не |
Для |
получения |
последовательности нормальных |
||
обходимо: |
|
РСП (0,1); |
|
|
1) выработать п чисел ££ |
|
97
2) |
вычислить 7 к по формуле (7.4.3); |
(7.4.4); |
||
3) |
откорректировать ^ |
по |
формуле |
|
4) |
получить г| g НСП |
(а, а) по формуле (7.4.5); |
||
5) |
снова выработать п чисел £ и т. д. |
|
||
Для получения чисел РСП |
(а, Ь) (равномерно распределен |
|||
ных |
с пределами а и Ь) |
при |
помощи |
РСП (0,1) пользуются |
несложной формулой: |
|
|
|
|
|
■4, = |
а -И , (6 - а ) . |
(7.4.6) |
|
II |
|
О? |
о) |
чи |
|
Kj |
Рк
в/ аг а"
Рис. 17. Статистический закон распределения случайной величины
а) гистограмма; б) равновероятная гистограмма; в) интегральный график
Часто распределение случайной величины задано гистограм мой (рис. 17) или таблицей.
На рис. 17, а) изображены вероятности Р;, соответствующие каждому отрезку (щ— аг-i) случайной величины тр Для моде лирования необходимо преобразовать гистограмму а) таким
образом, чтобы каждому |
интервалу |
( |
ак — aK-i) |
соответствова |
|
ла равная вероятность Рк, т. е. Pi = |
Р2 |
... Рп. |
|
||
|
i p k = |
i. |
|
|
|
Тогда гистограмма будет |
иметь вид |
(рис. 17,6). |
Для получе |
ния чисел тр соответствующих статистическому закону, изобра женному гистограммой рис. 17, необходимо: 1) выработать чис
98