
- •С.А. Фоменков, д.А. Давыдов, в.А. Камаев. Математическое моделирование системных объектов
- •Isbn 5-230-04689-9
- •7.1. Понятие «агрегат» в теории систем ……………………………………….
- •7.3. Моделирование процесса функционирования агрегата ………………..
- •7.4. Кусочно-линейные агрегаты ……………………………………………….
- •1. Понятие о моделях и моделировании
- •1.1. Определение понятия модель. Свойства моделей
- •1.2. Классификация моделей
- •Идеальные модели
- •1.3. Классификация математических моделей по свойствам обобщенного объекта моделирования
- •1.4. Адекватность и эффективность математических моделей
- •1.5. Методы построения моделей
- •1.5.1. Общая логика построения моделей
- •1.5.2. Аналитические модели
- •1.5.3. Идентифицируемые модели
- •1.6. Вопросы для самопроверки
- •2. Построение математических моделей по экспериментальным данным
- •2.1. Постановка задачи идентификации
- •2.2. Идентификация моделей с помощью регрессионного метода
- •2.2.1. Идентификация статических линейных систем с несколькими входами
- •2.2.2. Идентификация нелинейных систем
- •2.2.3. Достоверность (адекватность) регрессионной модели
- •2.3. Построение моделей идентификации поисковыми методами
- •2.4. Вопросы для самопроверки
- •3. Математическое моделирование сложных неоднородных систем
- •3.1. Математические модели элементов
- •3.2. Математические модели взаимодействия элементов сложной системы
- •3.3. Вопросы для самопроверки
- •4. Моделирование по схеме марковских случайных процессов
- •4.1. Классификация марковских процессов
- •4.2. Расчет марковской цепи с дискретным временем
- •4.3. Марковские цепи с непрерывным временем
- •4.3.1. Уравнения Колмогорова
- •4.3.2. Поток событий. Простейший поток и его свойства
- •4.3.3. Пуассоновские потоки событий и непрерывные марковские цепи
- •4.3.4. Предельные вероятности состояний
- •4.3.5. Схема гибели и размножения
- •4.4. Вопросы для самопроверки
- •5. Теория массового обслуживания
- •5.1. Классификация смо и их основные характеристики
- •5.2. Одноканальная смо с отказами
- •5.3. Многоканальная смо с отказами
- •5.4. Одноканальная смо с ожиданием
- •5.5. Вопросы для самопроверки
- •6. Сети Петри
- •6.1. Основные определения
- •6.2. Сети Петри для моделирования
- •6.3. Анализ сетей Петри
- •6.3.1. Задачи анализа сетей Петри
- •6.3.2. Методы анализа
- •6.4. Обобщения сетей Петри
- •6.4.1. Временные сети событий (всс)
- •6.4.3. Комби-сети
- •6.5. Вопросы для самопроверки
- •7. Агрегативное описание систем
- •7.1. Понятие «агрегат» в теории систем
- •7.2. Процесс функционирования агрегата
- •7.2.1. Функционирование агрегата общего вида
- •7.2.2. Функционирование смо как агрегата
- •7.3. Моделирование процесса функционирования агрегата
- •7.4. Кусочно-линейные агрегаты
- •7.4.1. Понятие о кусочно-линейном агрегате.
- •7.4.2. Процесс функционирования кла.
- •7.4.3. Примеры представления систем в виде кла.
- •8.2. Способы организации единичного жребия
- •Появилось или нет событие а?
- •Какое из нескольких возможных событий появилось?
- •3. Какое значение приняла случайная величина?
- •4. Какую совокупность значений примет система случайных величин?
- •8.3. Современное содержание терминов «имитация», «имитационная модель».
- •8.4. Приемы построения и эксплуатации имитационных моделей
- •8.5. Вопросы для самопроверки
- •9. Когнитивные подходы к решению слабоструктурированных и плохо формализованных задач
- •9.1. Когнитивные модели
- •9.2. Вопросы для самопроверки
- •Список использованной литературы
5.2. Одноканальная смо с отказами
Простейшая задача. Пусть СМО состоит только из одного канала (n=1) и на нее поступает пуассоновский поток заявок с интенсивностью , зависящей в общем случае от времени =(t) (5.1). Заявка, заставшая канал занятым, получает отказ и покидает систему. Обслуживание заявки продолжается в течение случайного времени Тоб, распределенного по показательному закону с параметром f(t)= e-t (t0) (5.2).
Из этого следует, что «поток обслуживаний» - простейший, с интенсивностью m. Требуется найти: абсолютную (А) и относительную (q) пропускные способности.
Рассмотрим единственный канал обслуживания как физическую систему S, которая может находиться в одном из двух состояний: S0 – свободен, S1 – занят. Обозначим вероятности состояний p0(t) и p1(t). Очевидно:
t p0(t)+p1(t)=1 (5.3).
Граф состояний системы изображен на рис. 5.1.
Рис. 5.1. Граф состояний одноканальной СМО с отказами
По графу состояний системы составим дифференциальные уравнения Колмогорова:
(5.4)
В соответствии с (5.3) одно уравнение в (5.4) лишнее. Отбросим второе уравнение, а первое перепишем с учетом (5.3):
или
(5.5).
Это
уравнение естественно решать при
начальных условиях p0(0)=1;
p1(0)=0.
Уравнение (5.5) легко может быть решено
не только для простейшего потока заявок
(=const),
но и для случая =(t).
Приведем решение (5.5) только для случая
=const:
. Графически результаты решения показаны
на рис. 5.2.
Рис. 5.2. Результаты решения системы уравнений (5.4)
Для нашего случая вероятность p0 есть не что иное, как q.
Действительно, p0 есть вероятность того, что в момент t канал свободен, иначе вероятность того, что заявка, пришедшая в момент t, будет обслужена. А значит, для данного момента времени t среднее число обслуженных заявок к числу поступивших также равно p0: q= p0.
В
пределе, при t,
когда процесс обслуживания уже
установится, предельное значение q
будет равно
.
Легко
найти и А,
зная q.
Они связаны очевидным соотношением:
.
В пределе, при t,
А тоже установится и будет равна
.
Зная
q
(вероятность того, что пришедшая в момент
t
заявка будет обслужена) легко найти
вероятность отказа: Pотк
=1-q.
Pотк
есть не что иное, как средняя доля
необслуженных заявок среди поданных.
В пределе, при t
.
5.3. Многоканальная смо с отказами
Рассмотрим n-канальную СМО с отказами. Будем нумеровать состояния системы по числу занятых каналов (или, что в данном случае то же, по числу заявок, связанных с системой). Состояния будут:
S0 – все каналы свободны;
S1 – занят ровно один канал, остальные свободны;
……
Sk – заняты ровно k каналов, остальные свободны;
…….
Sn – заняты все n каналов.
Граф состояний имеет следующий вид (рис. 5.3). Слева направо систему переводит один и тот же поток – поток заявок с интенсивностью .
Рис. 5.3. Граф состояний многоканальной СМО с отказами
Очевидно, если обслуживанием занято 2 канала, а не один, поток обслуживаний, переводящий систему по стрелке S2S1, будет вдвое интенсивнее (2), если занято k- каналов – в k раз интенсивнее (k). Процесс такого вида представляет собой частный случай процесса гибели и размножения. Составляем уравнения Колмогорова:
(5.6).
Уравнения (5.6) называются уравнениями Эрланга. Естественными начальными условиями являются:
p0(0)=1; p1(0)=p2(0)=…=pn(0)=0
Интегрировать (6) в аналитическом виде довольно сложно, на практике решают численно с использованием ЭВМ. Такое решение дает нам все вероятности состояний как функции времени: p0(t), p1(t), …, pn(t).
Больше всего интересны предельные вероятности состояний, характеризующие установившийся режим работы СМО (при t). Воспользуемся готовым решением, полученным для схемы гибели и размножения:
(k=1,2,..n)
(5.7).
Обозначим
и будем называть величину
«приведенной интенсивностью» потока
заявок. Физический смысл её таков:
величина
представляет
собой среднее число заявок, приходящих
в СМО за среднее время обслуживания
одной заявки. С учетом этого (5.7) принимает
вид:
(5.8). Формулы
Эрланга.
Теперь можно найти характеристики эффективности СМО: q, А, Ротк.
Заявка получает отказ, если приходит в момент, когда все n каналов заняты. Вероятность этого равна:
.
Вероятность того, что заявка будет принята к обслуживанию (она же q) дополняет Ротк до 1: q = 1-pn. И наконец: А= q=(1- pn).
Одной
из важных характеристик СМО с отказами
является среднее число занятых каналов
(в данном случае оно совпадает со средним
числом заявок, находящихся в системе).
Обозначим это среднее число
.
Величину
можно вычислить непосредственно по
формуле:
как математическое ожидание дискретной случайной величины, принимающей значение 0,1, …n с вероятностями p0, p1…pn.
Однако значительно проще выразить через А . А есть не что иное, как среднее число заявок, обслуживаемых в единицу времени; один занятый канал обслуживает в среднем за единицу времени заявок; следовательно, среднее число занятых каналов
или
.