
- •С.А. Фоменков, д.А. Давыдов, в.А. Камаев. Математическое моделирование системных объектов
- •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.4. Одноканальная смо с ожиданием
На СМО поступает поток заявок с интенсивностью , интенсивность обслуживания (т.е. в среднем непрерывно занятый канал будет выдавать обслуженных заявок в единицу времени), n=1. Заявка, поступившая в момент, когда канал занят, становится в очередь и ожидает обслуживания. Предположим, что количество мест в очереди ограничено числом m (в дальнейшем, при m можно получить характеристики одноканальной СМО без ограничений по длине очереди). Будем нумеровать состояния СМО по числу заявок, находящихся в системе (как обслуживаемых, так и ожидающих обслуживания):
S0 – канал свободен; S1 – канал занят, очереди нет; S2 – канал занят, одна заявка стоит в очереди; Sk – канал занят, k -1 стоят в очереди; Sm+1 – канал занят, m заявок стоят в очереди. Граф состояний (размеченный) представлен на рис. 5.4.
Рис. 5.4. Размеченный граф состояний одноканальной СМО с ожиданием
Снова схема гибели и размножения. Пользуясь общим решением, напишем выражения предельных вероятностей состояния:
или
(5.9).
В знаменателе выражения для р0 стоит геометрическая прогрессия, суммируя её находим
(5.10).
(5.10)
справедливо только при
(иначе неопределенность вида
).
Но сумму геометрической прогрессии со
знаменателем
найти ещё проще, чем по (5.10); она равна
(m+2)
и в этом случае p0=1/(m+2)
{то же самое получится, если раскрыть
неопределенность по правилу Лапиталя}.
Определим
характеристики СМО: Ротк
q,
А, среднюю
длину очереди
,
среднее число заявок, связанных с
системой
.
Очевидно, заявка получает отказ только в том случае, когда канал занят и все m мест в очереди – тоже:
(5.11).
(5.12).
А=q.
Найдем среднее число
,
находящихся в очереди; определим эту
величину как математическое ожидание
дискретной случайной величины R
– числа заявок , находящихся в очереди:
.
.
Выведем
формулу для суммы. Эта сумма не что иное,
как производная по
суммы
,
а для этого выражения мы можем
воспользоваться формулой суммы
геометрической прогрессии:
Продифференцируем её по и проведя преобразования, найдем
(5.13).
Тогда
.
Подставляем p0 из (5.10) и получаем
(5.14).
Выведем
теперь формулу для
.
Рассмотрим общее число заявок К,
связанных с системой, как сумму двух
случайных величин: числа заявок, стоящих
в очереди и числа заявок, находящихся
под обслуживанием:
.
По теореме сложения математических ожиданий:
,
где
- среднее число заявок в очереди;
- среднее число заявок под обслуживанием.
Найдем
.
Т.к. канал у нас один, то случайная
величина
может принимать только два значения: 0
или 1. Значение 0 она принимает, если
канал свободен; вероятность этого равна
.
Значение 1 она принимает, если канал
занят; вероятность этого равна
.
Отсюда находим:
.
,
где
находим из (5.14).
Выведем
выражение еще для одной существенной
характеристики СМО с ожиданием: среднего
времени ожидания заявки в очереди.
Обозначим его
.
Пусть заявка приходит в систему в
какой-то момент времени. С вероятностью
p0
канал обслуживания не будет занят и ей
не придется стоять в очереди (tож=0).
С вероятностью p1
она придет
в систему во время обслуживания какой-то
заявки, но перед ней не будет очереди,
и заявка будет ждать начала своего
обслуживания в течение времени 1/
(среднее время обслуживания одной
заявки). С вероятностью p2
в очереди перед ней будет стоять еще
одна заявка и время ожидания в среднем
будет 2/
и т.д. Вообще, с вероятностью pk
пришедшая
заявка застанет в системе k
заявок и будет ждать в среднем k/
единиц времени (1km).
При k=m+1
(в очереди m
заявок, вероятность этого pm+1)
tож=0
(заявка не обслуживается).
.
Подставим сюда выражения для p1,p2,…pm из (5.9).
.
Преобразуем сумму в скобках, используя (5.13)
Или, выражая p0 через
.
Сравнивая это выражение с (5.14) , замечаем, что
,
(5.15)
т.е. среднее время ожидания равно среднему числу заявок в очереди, деленному на интенсивность потока заявок.
Выведем ещё одну формулу для среднего времени пребывания заявки в системе. Обозначим случайную величину – время пребывания заявки в СМО через Тсист.. Она складывается из двух слагаемых (тоже случайных):
Тсист.=Тож
+
,
где Тож
- время ожидания заявки в очереди,
случайная величина, равная времени
обслуживания Тоб,
если заявка обслуживается и 0, если она
не обслуживается (получает отказ). По
теореме сложения математических
ожиданий:
,
но в наших обозначениях
,
а
.
Отсюда находим:
или с учетом (5.15)
.
Формула Литтла (первая): для любой СМО, при любом характере потока заявок, при любом распределении времени обслуживания, при любой дисциплине обслуживания среднее время пребывания заявки в системе равно среднему числу заявок в системе, деленному на интенсивность потока заявок.
Формула Литтла (вторая): для любой СМО, при любом характере потока заявок, при любом распределении времени обслуживания, при любой дисциплине обслуживания среднее время пребывания заявки в очереди равно среднему числу заявок в очереди, деленному на интенсивность потока заявок.
Теперь можно рассмотреть работу одноканальной СМО с ожиданием при m (неограниченная очередь). Совершить предельный переход m. Можно рассмотреть работу многоканальной СМО с ожиданием. Состояние системы будем нумеровать по числу заявок, связанных с системой:
S0 – все каналы свободны;
S1 – занят один канал, остальные свободны;
… … …
Sk – заняты k каналов, остальные свободны;
… … …
Sn – заняты все n каналов;
Sn+1 – заняты все n каналов, одна заявка стоит в очереди;
… … …
Sn+r – заняты все n каналов, r заявок в очереди;
… … …
Sn+m – заняты все n каналов, m заявок в очереди.
Размеченный граф состояний представлен на рис. 5.5.
Рис. 5.5. Размеченный граф состояний многоканальной СМО с ожиданием
Теперь необходимо написать уравнения Колмогорова, найти вероятности состояний и с их помощью рассчитать все интересующие величины. Затем опять можно рассмотреть и случай m.
Можно
рассмотреть СМО с ограниченным
временем ожидания
(на каждую заявку, стоящую в очереди
действует как бы «поток уходов» с
интенсивностью
(
-
среднее время пребывания в очереди)).
Фрагмент размеченного графа состояний
для такой системы приведен на рис. 5.6.
Рис. 5.6. Фрагмент размеченного графа состояний многоканальной СМО с ограниченным временем ожидания
Существуют и другие разновидности СМО: замкнутые СМО (интенсивность потока поступающих заявок зависит от состояния самой СМО), СМО с «взаимопомощью» между каналами (незанятые каналы «помогают» занятому в обслуживании) и т.д. Результаты расчетов такого рода систем можно найти в соответствующей литературе, указанной в списке источников.