
- •Основные понятия теории моделирования Моделирование как метод научного познания
- •Понятие модели и моделирования
- •Понятие информационной модели
- •Виды моделирования
- •Математическое и компьютерное моделирование
- •Пример моделирования простейшей системы
- •Формализация систем и классификация моделей Некоторые термины системного анализа
- •Формальная модель объекта.
- •Классификация моделей
- •Основные требования к математическим моделям
- •Этапы процесса моделирования
- •Статические регрессионные модели Понятие черного ящика и регрессии
- •Линейная одномерная регрессионная модель
- •Линейная множественная регрессионная модель
- •Нелинейные регрессионные модели
- •Динамические модели Динамические системы
- •Динамическая система первого порядка
- •Звено второго порядка (колебательное звено)
- •Динамические регрессионные модели
- •Модель в виде фильтра Каллмана
- •Статистическое моделирование Понятие статистического моделирования
- •Метод Монте-Карло.
- •Генераторы случайных чисел
- •Проверка качества работы генератора
- •Моделирование случайной величины с заданным законом распределения
- •Моделирование системы случайных величин
- •Потоки случайных событий
- •Системы массового обслуживания (смо)
- •Принципы компьютерного моделирования смо
- •Задания
- •1. Моделирование простейшей системы
- •2. Статическая регрессионная модель
- •3. Динамическая регрессионная модель
- •4. Метод Монте-Карло и генераторы случайных чисел
- •5. Моделирование смо
- •Приложения
- •Нормальный закон распределения (закон Гаусса)
- •Распределение Пуассона
- •Список источников
- •Моделирование процессов и систем
- •153000 Г. Иваново, пр. Ф. Энгельса, 21 Оглавление
Нормальный закон распределения (закон Гаусса)
Наиболее часто на практике встречается нормальный закон распределения, часто называемый законом Гаусса. Он играет исключительно важную роль в теории вероятностей и занимает среди других законов распределения особое положение. Главная его особенность состоит в том, что он является предельным законом, к которому приближаются другие законы распределения при весьма часто встречающихся условиях. Большинство случайных величин (таких как ошибки измерений, стрельбы и т.д.) могут быть представлены как суммы большого числа сравнительно малых слагаемых – элементарных ошибок, каждая из которых вызвана действием одной какой-то причины. Каким бы законам распределения ни были подчинены отдельные элементарные ошибки, в сумме большого числа слагаемых их особенности нивелируются, и сумма оказывается подчиненной закону распределения, близкому к нормальному.
Одна из формулировок центральной предельной теоремы: «Для большого числа N случайных величин X с любым законом распределения их сумма есть случайное число с нормальным законом распределения при условии, что влияние на сумму отдельных слагаемых является равномерно малым1».
Формула распределения вероятности значений случайной величины x по нормальному закону имеет вид (рис. 16):
где x – случайная величина;
f
(x)
– вероятность принятия случайной
величиной значения x;
mx – математическое ожидание;
σx – среднее квадратичное отклонение.
Максимальная
ордината кривой соответствует точке
и равна
По
мере удаления от точки
кривая
асимптотически приближается к оси
абсцисс. Среднее квадратичное отклонение
σx
характеризует форму кривой: при уменьшении
σx
кривая вытягивается вверх, сужаясь с
боков.
Н
ормализованным
нормальным распределением (рис.
17) называется
такое нормальное распределение, у
которого mx
= 0 и σx
= 1. Из
нормализованного распределения можно
получить любое другое нормальное
распределение с заданными mx
и σx
по формуле:
Распределение Пуассона
Говорят, что случайная дискретная величина X распределена по закону Пуассона, если вероятность того, что она примет определенное значение m, выражается формулой
Здесь a – параметр Пуассона или математическое ожидание случайной величины X. Кроме того, для распределения Пуассона
Приложение 4. Пример моделирования СМО
Условие задачи.
Составы с песком приходят на разгрузочный пункт с интенсивностью 6 шт. в сутки (простейший поток). В каждом составе от 5 до 14 вагонов, количество вагонов распределено равномерно. В каждом вагоне количество песка распределено нормально, m=20, σ=3. На пункте несколько погрузчиков, но с одним вагоном работает только один. Время, необходимое для разгрузки 1 тонны песка 6 мин. Время на подъезд погрузчика к вагону 3 мин. Время подачи состава с места ожидания на площадку 6 мин. Поезд не встает в очередь (уходит), если перед ним 3 состава. Определить простой погрузчиков и поездов, потерянную выручку из-за отчаявшихся машинистов и т.д.
Некоторый анализ системы.
Заявки – это, очевидно, подходящие составы. Поданные под разгрузку вагоны не образуют ординарный поток, но ординарным можно считать поток составов. Обслуживанием можно назвать процесс разгрузки состава. Можно считать, что заявке отказано в обслуживании, если в момент ее поступления в очереди 3 заявки. Время на подготовку к обслуживанию не требуется, следовательно, каналы горячие. Дисциплина обслуживания – First In, First Out.
Обсуждение алгоритма
Предлагается
реализовать принцип последовательной
проводки заявок, т.к. процесс ожидания
заявки и процесс обслуживания не
представляют интереса и являются
довольно длительными, а канал обслуживания
один. Считаем, что первый состав приходит
на пункт в момент времени
.
К этому моменту все погрузчики свободны.
Определяем
случайным образом время
прибытия
всех остальных поездов, считая поток
поездов простейшим:
где
случайное
равномерно распределенное число. Получим
массив моментов появления заявок.
Следуя принципу последовательной проводки заявок, рассматриваем каждую заявку отдельно, сохраняя в специально созданном массиве время отправления разгруженного состава, т.к. оно определяет количество поездов в очереди на момент поступления очередной заявки и, следовательно, определяет, будет ли в обслуживании отказано.
Итак, для каждой i-й заявки определим сначала величину очереди: перебирая поезда с 1-го по i–1-й, подсчитаем количество составов, время отправления которых оказалось больше, чем время прибытия рассматриваемого i-го. Если в очереди больше трех составов, заявка отклоняется, а в качестве времени отправления ушедшему без обслуживания составу записывается ноль. Для обслуживаемой заявки определяем при помощи датчика случайных чисел количество вагонов и количество песка в каждом вагоне, назначаем исполнителя – один из погрузчиков. В массиве для каждого погрузчика указывается время его освобождения. Исполнителем становится либо свободный погрузчик (его время освобождения меньше времени начала разгрузки состава), либо тот, время освобождения которого минимально. Назначение погрузчику текущего вагона как задания есть переопределение времени освобождения.
Псевдокод алгоритма
ввод величин из условия задачи
n=1;
;
;
Пока не закончится время эксперимента
{n=n+1; определяем как время наступления события в простейшем потоке}
Для каждого
погрузчика (
количество
погрузчиков):
Для каждого i-го поезда (i=1..n)
{ определение длины очереди – подсчет количества еще не отправленных среди ранее пришедших составов;
если длина очереди >3, то
;
составу отказано в обслуживании}
иначе
{если длина очереди =0, то
иначе
};
количество вагонов в составе = равномерно распределенному от 5 до 14 случайному числу;
для каждого j-го вагона (j=1...количество вагонов)
{количество песка = нормально распределенному случайному числу;
определение свободного или освобождающегося первым погрузчика с номером k;
кол-во песка*скорость
разгрузки + время для подхода к вагону}
количество
погрузчиков
} //конец обработки обслуженной заявки
} //конец обработки i-й заявки
время
эксперимента =
конец
Интересующие исследователя величины (такие как количество отклоненных заявок, среднее время простоя составов и погрузчиков) можно рассчитывать, вводя в алгоритм дополнительные переменные. Для сбора статистических данных лучше всего 1-й поезд не учитывать.