Опубликован: 16.11.2010 | Уровень: специалист | Доступ: свободно
Лекция 4:

Статистическое моделирование

3.9. Моделирование совместных зависимых событий

Пусть события A и B имеют вероятности свершения P(A) и P(B) соответственно. Условная вероятность P(B/ A) известна.

Покажем способ моделирования совместных зависимых событий на примере.

Пример 3.9. При испытании нового автомата определены вероятности горизонтального и вертикального отклонений пробоин от точки прицеливания P(A) = P(r _{1} \le  10см ) и P(B) = P(r_{2}\le  10см ).

Вероятность отклонения пробоин по высоте относительно тех, которые уложились в пределы допустимого бокового отклонения, равна:

P(B/\overline{A})=\cfrac{P(B) - P(A)P(B/A)}{1-P(A)}

Соответствующий фрагмент модели приведен на рис. 3.18.

Алгоритм моделирования совместных зависимых событий

Рис. 3.18. Алгоритм моделирования совместных зависимых событий

Пример 3.10. В ремонтное подразделение поступают вышедшие из строя средства связи (СС). В каждом СС могут быть неисправными в любом сочетании блоки A, B, C . Вероятности выхода из строя блоков P_{A}, P_{B}, P_{C}  соответственно. Ремонт производится путем замены неисправных блоков исправными блоками. В момент поступления неисправного СС вероятности наличия исправных блоков P_{HA}, P_{HB} , P_{HC} соответственно. При отсутствии хотя бы одного из исправных блоков A, B, C ремонт неисправного СС не производится.

Построить алгоритм имитационной модели с целью определения абсолютного и относительного количества отремонтированных СС с неисправными блоками A, B, C и A, B из общего количества R поступивших в ремонт СС.

Решение

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

Таблица 3.5. Вероятности появления неисправных блоков
Q_i ABC AB\overline{C} С другими блоками
P_i P_{A}P_{B}P_{C} P_{A} P_{B} (1- P_{C} )
l_r l_{1} =P_{A} P_{B} P_{C} l_{2} = l_{1} + P_{A} P_{B} (1- P_{C} ) 1
Таблица 3.6. Вероятности наличия исправных блоков
Q_{Hi} ABC AB\overline{C} С другими блоками
P_{Hi} P_{HA} P_{HB}  P_{HC} P_{HA} P_{HB} (1- P_{HC})
l_{Hr} l_1 = P_{HA} P_{HB}  P_{HC} l_{2} =l_{1} + P_{HA} P_{HB} (1- P_{HC}) 1

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

Алгоритм имитационной модели приведен на рис. 3.19.

В алгоритме приняты следующие обозначения:

N_{0} - заданное количество реализаций модели;

N - счетчик количества реализаций модели; M - счетчик числа отремонтированных СС за N реализаций модели;

A - абсолютное количество отремонтированных СС;

D - относительное количество отремонтированных СС.

Алгоритм модели функционирования системы ремонта

увеличить изображение
Рис. 3.19. Алгоритм модели функционирования системы ремонта

Согласно постановке задачи в блоках 3…7 по данным табл. 3.5 разыгрывается, с какими неисправными блоками поступает СС в ремонт. В результате розыгрыша определяется номер интервала (столбца табл. 3.5) и запоминается в переменной k .

Аналогично в блоках 8…11 разыгрывается по данным табл. 3.6 наличие в ремонтном подразделении необходимых блоков для замены.

Если такие блоки имеются, т. е. выполняется условие k =i в блоке 12, в счетчик M (блок 13) добавляется единица.

3.10. Классификация случайных процессов

Случайная величина X (t ), зависящая от одного неслучайного вещественного аргумента t , называется случайным процессом. X (t ) является случайной величиной при каждом фиксированном

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

Определим два понятия, присущие случайным процессам: сечение и реализация (рис. 3.20).

Сечением случайного процесса X (t ) называется случайная величина x(t_{j} ), являющаяся значением случайного процесса в фиксированный момент времени t_{j} .

Реализацией случайного процесса X (t ) называется функция времени x_{i} (t ), описывающая течение процесса в некотором i -м опыте.

Реализации и сечения случайного процесса

Рис. 3.20. Реализации и сечения случайного процесса

Случайный процесс X (t ) и аргумент t могут быть дискретными или непрерывными.

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

  • разбить интервал исследования на M временных точек t_j (j = \overline{1, M}), которых должно быть столько, чтобы обеспечить необходимую точность воспроизведения исследуемого процесса;
  • выполнить одну реализацию случайного процесса, то есть для каждого момента времени t_j определить сечение, разыграв случайное число, обладающее характеристиками случайного процесса;
  • определить аналогичные сечения для каждой из N реализаций случайного процесса (число N выбирается таким, чтобы обеспечить необходимые точность и достоверность результатов).

Случайные процессы могут быть:

  • стационарные;
  • нестационарные.

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

M[x(t)]=M[x];\,\, D[x(t)] = D[x].

Такие процессы называют стационарными.

Что касается автокорреляционной функции K_{k}(t_{k},t_{3}), то ее значение в стационарном процессе зависит только от разности t_{k}-t_{3} и не зависит от того, в каком месте временной оси находятся точки t_{k} и t_{3}.

Для стационарного процесса нет необходимости определять искомые характеристики для всех M сечений, а достаточно только для одного сечения N реализаций случайного процесса. То есть вместо M\times N измерений выполнить только N измерений. По данным этих измерений рассчитываются оценки x(t_j) и S^2_{x(t_j )}, которые в силу стационарности и являются оценками характеристик всего случайного процесса M[x] и D[x].

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

В работе модели стационарного процесса может присутствовать и нестационарный период. Это разного рода переходные процессы. Например, характеристики начального периода работы модели нестационарные из-за того, что начальные установки характеристик процесса были не равны характеристикам, значения которых установятся в дальнейшем. Естественно, речь идет о средних значениях характеристик.

Важнейшим свойством случайного процесса является свойство эргодичности.

Свойство эргодичности заключается в том, что все реализации случайного процесса имеют одинаковые статистические характеристики. Отсюда следует, что одна реализация случайного процесса характеризует весь случайный процесс X (t ), следовательно, для определения статистических характеристик процесса достаточно выполнить одну реализацию.

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

\cfrac{1}{N}\sum\limits_{i=1}^{N}{x_i(t_j)}=
\cfrac{1}{M}\sum\limits_{j=1}^{M}{x_i(t_j)}

Свойством эргодичностиобладают многие случайные процессы и, в том числе, все стационарные.

Таким образом, можно сформулировать определение эргодиче-ского процесса.

Случайный процесс X (t ) называется эргодическим, если его основные характеристики M [x] и D[x] могут быть получены не только усреднением по множеству реализаций, но и усреднением по времени одной реализации.

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

Владислав Нагорный
Владислав Нагорный
Высшее образование
Лариса Парфенова
Лариса Парфенова
Экстерн
Петр Гончар-Зайкин
Петр Гончар-Зайкин
Россия
Борис Борисов
Борис Борисов
Казахстан, Алматы, Казахский государственный университет, 1983