Файл: Шепелев, И. Г. Математические методы планирования и управления в строительстве конспект лекций.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.10.2024

Просмотров: 59

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
Рис. 15. Закон распределения равномер­ ной случайной последовательности:
а) дифференциальный закон; б) интегральный закон

Для получения (генерирования) чисел |

РСП

(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