- •Рецензенты:
- •Содержание
- •Введение
- •1. Системный подход к задаче моделирования
- •1.1. Определение системы
- •1.2. Концептуальное определение модели
- •1.3. Классификация моделей
- •1.4. Теоретико-множественное задание моделей
- •Контрольные задания
- •2. Модели диhамичеcкиx cиcтем
- •2.1. Фоpмализация
- •2.2.1. Примеры моделей в виде дифференциальных уравнений.
- •2.3. Инеpционные модели
- •2.4. Модели на оcнове пеpедаточныx функций
- •2.5. Конечные автоматы
- •Контрольные задания
- •3. Моделирование cтоxаcтичеcких объектов
- •3.4. Примеры моделирования cлучайныx объектов
- •3.5. Методы моделирования cлучайныx фактоpов
- •3.6. Проверочные тесты
- •3.7. Имитация cлучайныx cобытий
- •3.8. Имитация непрерывных случайных величин
- •3.9. Фикcация и обpаботка pезультатов моделиpования
- •3.10. Количеcтво pеализаций опытов при имитационном моделированиии
- •Контрольные задания
- •Библиографический список
- •Моделиpоваhие при проектировании информационно-управляющих
3.4. Примеры моделирования cлучайныx объектов
3.4.1. Модель отказа прибора. Прибор состоит из двух блоков, которые определяют его безотказную работа. Вероятность безотказной работы первого блока равна p1, а второго — p2. Прибор испытывался в течение заданного времени Т, в результате чего обнаружено, что он вышел из строя. Найти вероятность того, что отказал только первый блок, в второй исправен.
Решение. Модель отказов работы прибора определена следующими гипотезами:
H0={оба узла исправны};
H1={первый узел отказал, а второй исправен};
H2={первый узел исправен, а второй отказал};
H3={оба узла отказали}.
Вероятности гипотез:
P(H0)=p1p2; P(H1)=(1‑p1)p2; P(H2)=p1(1‑p2);
P(H3)=(1‑p1)(1‑p2).
Событие А состоит в том, что прибор отказал. В результате получаем вероятности:
P(А/H0)=p1p2; P(А/H1)=P(А/H2)=P(А/H3)=1.
Согласно формуле Бейеса, модель отказа первого блока и исправности второго блока прибора определится по формуле
.
3.4.2. Модель контроля качества изделий. Завод изготавливает изделия, каждое из которых с вероятностью p имеет дефект. В цехе изделие с равной вероятностью осматривается одним из двух контроллеров. Первый контроллер обнаруживает дефект с вероятностью p1, а второй — с вероятностью p2. Если в цехе изделие не забраковано, то оно поступает в ОТК завода, где дефект может быть обнаружен с вероятностью p3. Найти модель обнаружения дефекта изделия.
Решение. До опыта возможны четыре гипотезы:
H0={изделие не забраковано};
H1={изделие забраковано первым контроллером};
H2={изделие забраковано вторым контроллером };
H3={изделие забраковано ОТК завода}.
Событие А состоит в том, что изделие забраковано. Так как вероятность P(А/H0)=0, то получаем вероятности гипотез:
P(H1)=pp1/2; P(H2)=pp2/2; P(H3)=p[1‑(p1+p2)/2]p3.
Согласно формуле Бейеса, модель обнаружения дефекта определится по формулам
;
;
.
3.4.3. Модель пуска двигителя. Производится ряд попыток пуска двигителя. Каждая попытка заканчивается включением двигателя независимо от других с вероятностью p=0,6. Вероятность неудачного пуска q=1‑p. Каждая попытка занимает время . Найти модель для определения времени Т, которое потребуется для запуска двигателя.
Решение. Число произведенных попыток пуска двигателя является случайной величиной X, распределенной по геометрическому закону, начинающемуся со значения . Модель времени пуска, как Т=X, имеет распределение
-
Т:
2
3
…
m
…
;
p
qp
q2p
…
qm-1p
…
Математическое ожидание времени пуска M[T]=M[X]=/p. Дисперсия распределения определится D[T]=2D[X]=2q/p2.
3.4.4. Модель браковки шариков. Браковка шариков для подшипников проводится следующим образом. Если шарик не проходит через отверствие диаметром d1, но проходит через отверствие диаметром d2>d1, то его размер считается приемлемым. Если какое-нибудь из этих условий не выполняется, то шарик бракуется. Известно, что диаметр шарика D есть нормально распределенная случайная величина с характеристиками md=(d1+d2)/2 и d=(d2‑d1)/4. Найти модель, описывающую вероятность p того, что шарик будет забракован.
Решение. Участок значений (d1,d2) симметричен относительно md. По формуле P{|X‑md|<l}=2Ф(l/), где l=(d2‑d1)/2 — половина длины участка, находим вероятность того, что шарик не будет забракован:
P{|D‑md|<(d2‑d1)/2}=
.
Тогда вероятность p определится по формуле
.
3.4.5. Модель проезда через регулируемый перекресток. На перекрестке стоит автоматический светофор, в котором одну минуту горит зеленый свет и половину минуты — красный, затем опять одну минуту горит зеленый свет и половину минуты — красный и т.д. Найти модель, отображающую вероятность проезда светофора при условии, что момент проезда автомашины через перекресток распределен равномерно в интервале, равном периоду смены цветов в светофоре.
Решение. Период смены цветов в светофоре раве 1+0,5=1,5 мин. Для того чтобы машина проехала через перекресток, не останавливаясь, достаточно, чтобы момент проезда перекрестка пришелся на интервал времени (0; 1). Для случайной величины, равномерно распределенной в интервале (0; 1,5), ), вероятность того, что она попадет на участок (0; 1), равна 2/3.
Время ожидания Тож есть смешанная случайная величина; с вероятностью 2/3 она равна нулю, а с вероятностью 1/3 принимает с одинаковой плотностью вероятности любое значение между 0 и 0,5 мин. График функции распределения показан на рис. 3.8.
Среднее время ожидания у перекрестка определится по формуле
М[Тож]=0×2/3+0,25×1/30,083 мин.
Рис. 3.8
Дисперсия времени ожидания определится по формуле
D[Тож]=2[Тож]
– (М[Тож])2=02×2/3+0,25×
-
-(0,083)20,0208
мин2;
0,144
мин.
