- •Оглавление
- •Предисловие
- •Введение
- •1. Модели и их свойства
- •Основные понятия и определения
- •Вопросы к разделу 1.1
- •Целенаправленность моделей
- •Вопросы к разделу 1.2
- •Свойства моделей
- •Вопросы к разделу 1.3
- •Способы реализации моделей
- •1.4.1. Материальные модели
- •1.4.2. Виды подобия
- •1.4.3. Идеальные модели
- •Вопросы к разделу 1.4
- •1.5. Виды моделей
- •1.5.1. Познавательные и прагматические модели
- •1.5.2. Детерминированные и вероятностные модели
- •1.5.3. Непрерывные и дискретные модели
- •1.5.4. Статические и динамические модели
- •1.5.5. Линейные и нелинейные модели
- •1.5.6. Стационарные и нестационарные модели
- •1.5.7. Сосредоточенные и распределенные модели
- •1.5.8. Классификация видов моделей
- •Вопросы к разделу 1.5
- •Кибернетические модели систем
- •1.6.1. Модель типа «черный ящик»
- •1.6.2. Модель состава системы
- •1.6.3. Модель структуры системы
- •1.6.4. Графы
- •1.6.5. Структурная схема системы
- •1.6.6. Итоги анализа моделей систем.
- •Вопросы к разделу 1.6
- •2. Аналитические математические модели систем
- •2.1. Общая математическая модель динамической системы
- •Вопросы к разделу 2.1
- •2.2. Частные математические модели динамических систем
- •2.2.1. Модели детерминированных линейных непрерывных систем
- •Модели дискретных систем. Конечные автоматы
- •Вопросы к разделу 2.2
- •2.3. Свойства динамических систем
- •2.4. Линейная непрерывная детерминированная модель многомерной динамической системы в переменных состояния
- •Вопросы к разделу 2.4
- •2.5.7. Определитель Грама
- •Вопросы к разделу 2.5
- •Линейное векторное пространство
- •2.6.1. Характеристические числа и характеристические векторы
- •2.6.2. Формула Бохера
- •2.6.3. Модальная матрица
- •2.6.4. Диагонализация квадратной матрицы
- •Вопросы к разделу 2.6
- •Управляемость и наблюдаемость
- •Вопросы к разделу 2.7
- •Компьютерное моделирование. Численное интегрирование дифференциальных уравнений
- •3.1. Ошибки усечения и округления
- •3.2. Метод Эйлера
- •3.3. Методы Рунге–Кутта
- •Сравнение различных методов решения. Контроль величины шага и устойчивость
- •Вопросы к главе 3
- •4. Динамика развития и использования моделей
- •4.1. Сложности алгоритмизации моделирования
- •Вопросы к главе 4
- •5. Аналитические вероятностные математические модели систем
- •5.1. Аналитические модели систем массового обслуживания
- •5.1.1. Важнейшие выходные параметры смо
- •5.1.2. Простейшие модели смо
- •5.1.3. Общая характеристика моделей смо
- •5.1.4. Дисциплины обслуживания
- •5.1.5. Характеристики входного потока заявок
- •5.1.6. Функция распределения Пуассона
- •5.1.7. Характеристики обслуживания
- •5.1.8. Показательный закон распределения времени обслуживания
- •5.1.9. Показатели качества обслуживания
- •5.1.10. Согласование источника заявок с каналом обслуживания
- •5.1.11. Оценка эффективности многоканальной смо
- •Вопросы к разделу 5.1
- •5.2. Сети Петри
- •5.2.1. Маркировка
- •5.2.2. Правила срабатывания переходов
- •5.2.3. Разновидности сетей Петри
- •5.2.4. Конфликтные ситуации в сетях Петри
- •5.2.5. Пример сети Петри для работы группы пользователей на одной рабочей станции
- •5.2.6. Пример сети Петри для системы обнаружения и устранения неисправностей в технической системе
- •5.2.7. Анализ сетей Петри
- •Вопросы к разделу 5.2
- •6. Имитационные модели систем
- •6.1. Имитационный эксперимент
- •Недостатки имитационного моделирования
- •6.2. Развитие имитационного моделирования
- •Основные фазы развития средств им
- •6.3. Этапы имитационного моделирования
- •6.4. Подходы к построению имитационных моделей
- •6.4.1. Событийный подход
- •6.4.2. Подход сканирования активностей
- •6.4.3. Процессно-ориентированный подход
- •6.5. Разработка программ им
- •6.5.1. Использование для им универсальных языков программирования
- •6.5.2. Использование для им специализированных языков моделирования
- •6.6. Имитационное моделирование систем массового обслуживания
- •6.6.1. Событийный метод моделирования
- •6.5.3. Создание и использование проблемно-ориентированных систем моделирования
- •6.6.2. Схема реализации событийного метода имитационного моделирования
- •Вопросы к главе 6
- •7. Метод «ресурсы–действия–операции» (рдо)
- •7.1. Основные положения метода рдо
- •7.1.1. Ресурсы сложной дискретной системы
- •7.1.2. Действия в сдс
- •7.1.3. Операции в сдс
- •7.1.4. Основные положения рдо-метода
- •7.2. Представление сдс в рдо-методе
- •7.3. Базовая структура инструментальной среды интеллектуальной системы
- •7.4. Продукционный имитатор
- •7.5. Моделирование в среде рдо
- •7.5.1. Основные понятия
- •7.5.2. Объекты исходных данных и объекты, создаваемые рдо-имитатором при выполнении прогона
- •7.5.3. Состав объектов модели
- •7.5.4. Назначение объектов модели
- •7.6. Интегрированная среда моделирования рдо
- •7.6.1. Состав функций исм
- •7.6.2. Главное окно исм рдо
- •7.6.3. Инструментальная панель
- •7.6.4. Работа с рдо-имитатором
- •Описание кадра анимации
- •Пример описания кадра анимации
- •Вопросы к главе 7
- •8. Краткое описание языка gpss
- •8.1. Оператор generate
- •8.2. Оператор function
- •8.3. Операторы split и assemble
- •8.4. Операторы seize и release
- •8.5. Оператор advance
- •8.6. Операторы enter и leave
- •8.7. Операторы queue и depart
- •8.8. Оператор test
- •8.9. Операторы start и terminate
- •8.10. Оператор transfer
- •8.11. Оператор assigne
- •8.12. Операторы управления движением заявок
- •8.13. Вычислительный оператор variable
- •8.14. Оператор синхронизации матсн
- •8.15. Пример программы на языке gpss для смо
- •Программа к примеру смо
- •Вопросы к главе 8
- •9. Планирование компьютерных экспериментов с моделями систем
- •9.1. Основные понятия теории планирования экспериментов
- •9.2. Модели планирования эксперимента
- •9.3. Виды планов экспериментов
- •Вопросы к главе 9
- •10. Обработка и анализ результатов компьютерного моделирования
- •10.1. Методы оценки
- •10.2. Статистические методы обработки
- •10.3. Задачи обработки результатов моделирования
- •10.3.1. Критерий согласия Колмогорова
- •10.3.2. Критерий согласия Пирсона
- •10.3.3. Критерий согласия Смирнова
- •10.3.4. Критерий согласия Стьюдента
- •10.3.5. Критерий согласия Фишера
- •10.4. Анализ и интерпретация результатов компьютерного моделирования
- •10.4.1. Корреляционный анализ результатов моделирования
- •10.4.2. Регрессионный анализ результатов моделирования
- •10.4.3. Дисперсионный анализ результатов моделирования
- •Вопросы к главе 10
- •Заключение
- •Список литературы
- •Приложение 1 Некоторые сведения из теории матриц
- •Основные типы матриц
- •Специальные типы матриц
- •Операции над матрицами Сложение матриц
- •Умножение матриц
- •Дифференцирование матриц
- •Интегрирование матриц
- •Определители
- •Свойства определителей
- •Нуль-граф и полный граф
- •Изоморфные графы
- •Плоские графы
- •Число ребер графа
- •Формула Эйлера для числа вершин, ребер и граней плоского графа
- •Распределение Лапласа
- •Вырожденное (причинное) распределение
- •Приложение 4 Краткие сведения о специализированных языках и проблемно-ориентированных системах имитационного моделирования
- •Предметный указатель
- •Список сокращений
Формула Эйлера для числа вершин, ребер и граней плоского графа
Формула Эйлера связывает число вершин и ребер плоского графа с числом его граней. Гранью называется область плоскости, ограниченная ребрами плоского графа, не содержащая внутри себя ни ребер, ни вершин.
Итак, формула Эйлера:
где n – число вершин, m – число ребер графа, f – число граней графа.
Исходя из этой формулы, был сформулирован ряд следствий:
Следствие 1. В любом простом планарном графе существует вершина, степень которой не больше пяти.
Следствие 2. Каждый планарный граф G с n ≥ 4 вершинами имеет, по крайней мере, четыре вершины со степенями, не превышающими 5.
Следствие 3. Если G – связный простой планарный граф с n ≥ 3 вершинами и m ребрами, то m ≤ 3n – 6.
Приведенные следствия определяют зависимость планарности графа от числа его вершин и ребер и задают границы интервала по числу ребер, при попадании в который необходимо проводить дополнительные исследования, чтобы получить достоверный ответ на вопрос, является ли планарным исследуемый граф.
ПРИЛОЖЕНИЕ 3
Непрерывные одномерные распределения вероятностей
Нормальное распределение (Гаусса)
Непрерывная случайная величина x распределена нормально с математическим ожиданием (центром) μ и дисперсией σ2, если
Здесь F(X) – функция распределения (рис. 1), а w(X) – плотность распределения (рис. 2).
Рис. 1. Вид функции нормального распределения вероятностей F(X)
Рис. 2. Вид плотности нормального распределения вероятностей w(X)
Равномерное (прямоугольное) распределение
Функция распределения на интервале (μ-a, μ+a) равна:
и имеет вид, показанный на рис. 3.
Рис.3. Вид функции равномерного распределения вероятностей F(X)
Плотность равномерного распределения на интервале (μ-a, μ+a) равна:
и имеет вид, показанный на рис. 4.
Рис.4. Вид плотности равномерного распределения вероятностей w(X)
Дисперсия для равномерного распределения имеет значение:
D(X)=a2/3.
Распределение Лапласа
Функция распределения равна:
и имеет вид, показанный на рис.5.
Рис.5. Вид функции распределения Лапласа для μ=0,5 и β=0,1
Плотность распределения Лапласа равна:
и имеет вид, показанный на рис.6.
Рис.6. Вид плотности распределения Лапласа для μ=0,5 и β=0,1
Значение дисперсии составляет: D(X)=2β2.
Вырожденное (причинное) распределение
Функция распределения равна единичной ступенчатой функции (рис.7):
Рис.7. Вид функции вырожденного распределения F(X)
Плотность распределения равна дельта-функции Дирака:
и имеет вид, показанный на рис.8.
Рис.8. Вид плотности вырожденного распределения w(X)
Приложение 4 Краткие сведения о специализированных языках и проблемно-ориентированных системах имитационного моделирования
Реализация требований к имитационным моделям в рамках универсального языка программирования приводит к громоздким и неудобным для практического использования программам. В большинстве таких программ могут разобраться только их авторы, а любое изменение в постановке задачи требует переработки значительной части текста программы.
Более того, даже при создании специализированных языков имитационного моделирования в зависимости от концептуальной направленности их авторов можно выделить четыре основных класса, которые перечислим в соответствии с классообразующим признаком:
1 – доработка универсального языка группой операторов, реализующих необходимые функции. Этот путь наиболее простой и в литературе приводится много языков, незначительно отличающихся друг от друга.
Наибольшее распространение в этом классе получил язык GASP, который может быть реализован на любой вычислительной системе, имеющей в программной среде компилятор с FORTRAN-IV. Существенным недостатком класса является отсутствие концептуальной выразительности и средств для проверки логики имитационной модели;
2 – ориентация на системы дифференциальных уравнений, удобная для реализации группы исследовательских моделей, породила класс языков, ярким представителем которого является DYNAMO и распространенный в России его диалект ИМИТАК.
Подробно идеология языка ИМИТАК изложена в монографии с большим количеством примеров. Но концепция, опирающаяся на понятия математической физики, оказалась не очень удобной для экономических расчётов.
Оставшиеся языки многие исследователи делят на две школы: школу фирмы IBM, где применяется язык GPSS, использующий в качестве дескрипторов схемные символы, и школу, использующую операторные языки;
3 – применение операторов приводит к построению моделей на основе аналитического расчленения изучаемого процесса на условия, необходимые для выполнения действий, результаты этих действий, а также временные взаимозависимости между элементами системы и действиями, в которых они участвуют.
Примером такого языка может служить SIMSCRIPT, пока не нашедший широкого применения в исследовательской практике;
4 – в языке GPSS, ориентированном на процессы, моделирование системы производится путем движения транзактов через блоки действий. Решения принимаются в результате воспроизведения логики функционирования моделируемой системы, которая представляется в виде блоков.
При моделировании исследователем, не владеющим программированием, необходимо отдать предпочтение наиболее наглядному и простому языку GPSS. Не менее важным достоинством GPSS является его распространенность и наличие версии языка для персональных компьютеров.
По причине значительной трудоемкости имитационное моделирование применялось только как метод исследований в научной деятельности и крайне редко – в практической работе. Однако, задачи практики требовали создания языков моделирования, более применимых в работе по управлению объектами более широкого класса, нежели уникальные научные установки. И к середине 80-х годов было создано в основном американскими фирмами, по данным Емельянова А.А., более 20 различных систем. К числу них относились уже упомянутые выше GASP, SIMULA, GPSS и SLAM.
Process Charter-1.0.2
Process Charter-1.0.2 (Процесс Чартер), разработчик – компания Scitor, Калифорния, США.
Процесс Чартер ориентирован на дискретное моделирование, удобный и простой в построении модели, самый дешевый из представленных продуктов, но ориентирован в основном на распределение ресурсов, слабая поддержка моделирования непрерывных компонент и недостаточно средств для построения диаграмм.
Powersim-2.01
Powersim-2.01 (Пауерсим), разработчик компания Modell Data, Берген, Норвегия.
Пакет Пауерсим, наоборот, является лучшим средством для построения непрерывных моделей, имеет множество встроенных функций, многопользовательский режим для коллективной работы группы исследователей с моделью, средства обработки массивов для упрощения создания аналогичных моделей. К недостаткам можно отнести сложную специальную систему обозначений и ограниченную поддержку дискретного моделирования.
Lthink-3.0.61
Lthink-3.0.61 (Лтинк), разработчик компания High Performans Systems, Ганновер, Германия.
Пакет Лтинк комбинирует достоинства первых двух программных средств и эффективно работает как с дискретными, так и с непрерывными моделями посредством встроенных блоков, имеется подробная обучающая программа, поддерживаются значительное количество форматов входных данных, развиты средства анализа чувствительности, имеется поддержка моделирования слабо подготовленными пользователями. К недостаткам можно отнести сложную систему обозначений и поддержку меньшего числа функций по сравнению с Пауерсим.
Extend+BPR-3.1
Extend+BPR-3.1 (Экстенд+БПР), разработчик – Imagine That!, Сан-Хосе, США.
Пакет, Экстенд+БПР, специально предназначен для анализа бизнес-проектов и также поддерживает дискретное и непрерывное моделирование. К числу достоинств нужно отнести понятную среду построения моделей с помощью блоков, множество встроенных блоков и функций для облегчения конструирования моделей, средства построения дополнительных функций посредством встроенного языка. К сожалению, данный пакет используется в достаточном объеме только на компьютерах Макинтош и имеет высокую стоимость порядка 1200 долларов.
ReThink
ReThink (Ретинк), разработчик – Gensym, Кембридж, Массачусетс, США.
Пакет Ретинк обладает основными свойствами Экстенд и имеет хороший графический транслятор для создания моделей, но работает под управлением экспертной системы G2, работающей в режиме реального времени. Недостатками являются слабая поддержка создания непрерывных компонент, необходимость экспертной системы, которая далеко не всегда необходима при моделировании, и высокая стоимость.
Piligrim-2.1
Piligrim-2.1 (Пилигрим) совместной разработки МегаТрон, Москва, Keisy, Гаага, Нидерланды и Энит, Таллин, Эстония.
Попытки создать универсальную систему моделирования привели к созданию системы Пилигрим, которая популярна в Европе и России. Пакет обладает широким спектром возможностей имитации временной, пространственной и экономической динамики моделируемых объектов. Разрабатываемые модели имеют качество коллективного управления процессом моделирования. Стоимость пакета ниже, чем Ретинк и экспертная система G2. К недостаткам можно отнести лишь отсутствие графического транслятора.
Arena-BE 3.6.1
Arena-BE 3.6.1, разработчик – Systems Modeling (Rockwell Software),
Одним из наиболее эффективных инструментов имитационного моделирования является система Arena компании Systems Modeling Arena позволяет строить имитационные модели, проигрывать их и анализировать результаты такого проигрывания. Имитационное моделирование – это универсальное средство для оптимизации процессов, поэтому модели с помощью Arena могут быть построены для самых разных сфер деятельности – производственных технологических операций, складского учета, банковской деятельности, обслуживания клиентов в ресторане и т.д. и т.п.
Имитационная модель компании Systems Modeling включает следующие основные элементы: источники и стоки (Create и Dispose), процессы (Process) и очереди (Queue). Источники – это элементы, от которых в модель поступает информация или объекты. Скорость поступления данных или объектов от источника обычно задается статистической функцией. Сток – это устройство для приема информации или объектов. Понятие очереди близко к понятию хранилища данных – это место, где объекты ожидают обработки. Времена обработки объектов (производительность) в разных процессах могут быть разными. В результате перед некоторыми процессами могут накапливаться объекты, ожидающие своей очереди. Часто целью имитационного моделирования является минимизация количества объектов в очередях. Тип очереди в имитационной модели может быть конкретизирован. Очередь может быть похожа на стек: пришедшие последними в очередь объекты первыми отправляются на дальнейшую обработку (LIFO: last-in-first-out). Альтернативой стеку может быть последовательная обработка, когда первыми на дальнейшую обработку отправляются объекты, пришедшие первыми (FIFO: first-in-first-out). Могут быть заданы и более сложные алгоритмы обработки очереди. Процессы – это аналог работ в функциональной модели. В имитационной модели может быть задана производительность процессов.
