- •Оглавление
- •Предисловие
- •Введение
- •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.5.3. Создание и использование проблемно-ориентированных систем моделирования
- •6.6. Имитационное моделирование систем массового обслуживания
- •6.6.1. Событийный метод моделирования
- •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 Краткие сведения о специализированных языках и проблемно-ориентированных системах имитационного моделирования
- •Предметный указатель
- •Список сокращений
Умножение матриц
Число столбцов матрицы A должно быть равно числу строк матрицы B, тогда говорят, что эти две матрицы согласуются по форме, и произведение AB существует. Если матрица A имеет размерность (m×n), а матрица B имеет размерность (n×k), то матрица C, являющаяся результатом произведения AB=C, будет иметь размерность (m×k). Условно это обозначим так:
.
Для матриц A (m×n) и B (n×m) существует как произведение AB, так и произведение BA. Произведение AB имеет размерность (m×m), а произведение BA - размерность (n×n). Естественно, что они в общем случае не равны. Даже в случае m=n, а значит, при одинаковой размерности (m×m) произведений AB и BA, эти произведения не обязательно равны. Если же оказывается, что они равны, т.е. AB=BA, то в этом случае говорят, что матрицы коммутативны.
Пример 3. Вычислить произведения указанных матриц:
;
.
Свойства умножения матриц. Умножение в общем случае не коммутативно, ассоциативно и дистрибутивно.
1. Некоммутативность:
AB≠BA.
2. Ассоциативность:
(AB)C=A(BC).
3. Дистрибутивность:
(A+B)C=AC+BC.
Умножение на скаляр.
При умножении на скалярную величину каждый элемент матрицы умножается на него.
Умножение на диагональную матрицу.
Умножение слева матрицы A на диагональную матрицу D эквивалентно операции эквивалентную операции со строками A. При умножении справа матрицы A на диагональную матрицу D операции производятся со столбцами матрицы A.
Умножение транспонированных матриц (транспонирование произведения матриц):
(A∙B)T = BT∙AT.
Умножение на единичную матрицу.
Умножение как слева, так и справа на единичную матрицу не изменяет исходную матрицу, т.е.
Дифференцирование матриц
Пусть
A(t)
– матрица (m x
n), элементы которой
aij
есть дифференцируемые функции
скалярной переменной t.
Производная от матрицы A(t)
по переменной t есть
матрица, элементами которой являются
:
.
Производная от суммы двух матриц равна сумме производных от этих матриц:
.
Производная от произведения матриц:
.
При этом должен сохраняться первоначальный порядок следования сомножителей произведения.
.
Интегрирование матриц
Интеграл от матрицы определяется как матрица, образованная из интегралов от элементов исходной матрицы. Следовательно,
Для обозначения интеграла от матрицы обычно используется символ Q=∫ ( )dt. Если оператор Q снабжен индексами (сверху t, а снизу t0), то они указывают нижний и верхний пределы интегрирования:
Пример
4. Найти
Определители
Определители существуют только для квадратных матриц.
В общем случае используется разложение Лапласа определителя n порядка по элементам строки (столбца) на сумму n определителей (n–1) порядка.
Например, для n = 3:
Свойства определителей
1. Определитель равен единице, если матрица А – единичная.
2. Определитель равен нулю, либо если все элементы матрицы равны нулю, либо все элементы строки (или столбца) равны нулю, или равны между собой или пропорциональны элементы произвольных двух строк (или двух столбцов).
3. Величина определителя остается неизменной по модулю при перестановке местами его строк (или столбцов).
4. Знак определителя изменяется на противоположный при замени местами его двух строк (или столбцов).
5. Значение определителя умножается на постоянную k, если все элементы какой-либо его строки (столбца) умножаются на k.
6. Значение определителя не изменяется, если к какой-либо его строке (или столбцу) прибавить умноженные на k соответствующие элементы другой строки (или столбца).
Миноры и алгебраические дополнения
Если в определителе |А| вычеркнуть i-ю строку и j-й столбец, то оставшиеся элементы образуют определитель, порядок которого на единицу меньше, чем у |А|, называемый минором элемента aij.
.
Минор Мij, взятый со знаком (–1)i+j называется алгебраическим дополнением элемента aij.
Cij = (–1)i+jMij .
Разложение определителя по Лапласу можно представить:
Присоединенная матрица
Присоединенная матрица образуется из алгебраических дополнений исходной матрицы А с последующим ее транспонированием:
Обратная матрица
Обратная матрица находится как частное от деления присоединенной матрицы Adj(A) на определитель |A|:
.
Произведение определителей
Определитель произведения двух квадратных матриц равен произведению их определителей:
.
Производная от определителя
Производная от определителя по одному из элементов равна алгебраическому дополнению этого элемента:
Произведение обратных матриц
Матрица, обратная произведению матриц, равна произведению матриц, обратных матрицам-сомножителям, взятому в обратном порядке.
Рассмотрим произведение двух матриц:
Для того, чтобы найти С–1, умножим обе части равенства слева на В–1А–1, а справа на С–1:
Производная от обратной матрицы
Для значений t, при которых A(t) дифференцируема и существует обратная матрица, производная от A–1(t) имеет вид:
.
Это выражение можно получить, рассматривая
Раскрывая, получим:
откуда уже окончательно получаем выражение для производной от A–1(t).
Некоторые специальные обратные матрицы
1. Инволютивная матрица. Матрица А называется инволютивной, если она совпадает со своей обратной матрицей:
2. Ортогональная матрица. Матрица А называется ортогональной, если ее транспонирование равносильно нахождению обратной матрицы:
A–1 = AT.
3. Унитарная матрица. Матрица А называется унитарной, если матрица, обратная А, равна матрице, сопряженной с А:
ПРИЛОЖЕНИЕ 2
Некоторые сведения из теории графов
Описание спортивных соревнований с помощью графов
Предположим, что футбольная команда вашего вуза участвует в соревнованиях и играет с командами других вузов. Пусть общее число команд равно шести. Вашу команду обозначим буквой А, а другие команды – буквами В, С, D, E и F. Через несколько недель после начала соревнований окажется, что некоторые из команд уже сыграли друг с другом, например:
А с С,D,F,
В с С, Е, F,
С с А, В,
D с А, Е,F,
Е с В,D,F,
F с А,В,D,Е.
Это можно изобразить при помощи такой геометрической схемы. Каждую команду представим точкой или маленьким кружочком и соединим отрезком те пары точек, которые соответствуют командам, уже игравшим друг с другом. Тогда для данного списка проведенных игр получим схему, изображенную на рис. 1.
Схема такого вида называется графом. Она состоит из нескольких точек А, В, С, D, Е, F, называемых вершинами, и нескольких соединяющих эти точки отрезков, таких как АС или ЕВ, называемых ребрами графа.
Из рис. 1 видно, что точки пересечения некоторых ребер графа могут не являться его вершинами; это происходит потому, что наш граф изображен на плоскости.
Рис. 1. Граф G сыгранных футбольных матчей в ходе соревнований
Возможно, удобнее было бы представлять себе его ребра нитями, проходящими друг над другом в пространстве. При изображении на плоскости вершины графа во избежание путаницы должны отмечаться достаточно отчетливо (например, кружочками).
Рис. 2. Граф состязания восьми команд
Каждую совокупность игр любого турнира можно представить соответствующим графом. Наоборот, если задан некоторый граф, т. е. фигура, состоящая из точек – вершин, соединенных прямолинейными отрезками – ребрами, то его можно рассматривать как схему такого состязания. В качестве примера рассмотрим граф, изображенный на рис. 2. Его можно представлять себе как граф, соответствующий состязанию восьми команд, где команда А уже играла с командами В, Е, D, команда В играла с А, F, G, С и т. д.
Упражнения
1. Начертите граф игр, сыгранных к середине сезона вашими футбольными или волейбольными командами.
2. Дайте полный список проведенных игр, соответствующий графу рис. 2.
3. Сколько ребер и сколько вершин имеют графы рис. 1 и 2?
