67313

Методы определения характеристик моделируемых систем

Лекция

Экономическая теория и математическое моделирование

По каждой такой характеристике y определяется N значений по которым строится гистограмма относительных частот вычисляется математическое ожидание дисперсия и моменты более высокого порядка определяются средние по времени и максимальные значения. В случае стационарного эргодического процесса...

Русский

2014-09-07

1.13 MB

2 чел.

Лекция № 14

Методы определения характеристик моделируемых систем.

Измеряемые характеристики моделируемых систем.

При имитационном моделировании можно измерять значения любых характеристик, интересующих исследователя. Обычно по результатам вычислений определяются характеристики всей системы, каждого потока и устройства.

Для всей системы производится подсчёт поступивших в систему заявок, полностью обслуженных и покинувших систему заявок без обслуживания по тем или иным причинам. Соотношения этих величин характеризует производительность системы при определённой рабочей нагрузке.

По каждому потоку заявок могут вычисляться времена реакций и ожидания, количества обслуженных и потерянных заявок. По каждому устройству определяется время загрузки при обслуживании одной заявки м число обслуженным устройством заявок, время простоя устройства в результате отказов и количество отказов, возникших в процессе моделирования, дины очередей и занимаемые ёмкости памяти.

При статистическом моделировании большая часть характеристик — это случайные величины. По каждой такой характеристике y определяется N значений, по которым строится гистограмма относительных частот, вычисляется математическое ожидание, дисперсия и моменты более высокого порядка, определяются средние по времени и максимальные значения. Коэффициенты загрузки устройств вычисляются по формуле:

k=Vk*Nok/Tm   (1)

Vk- среднее время обслуживания одной заявки к-ым устройством;

 Nok - количество обслуженных заявок устройством за время моделирования Tm.

Определение условий удовлетворения стохастических ограничений при имитационном моделировании производится путём простого подсчёта количества измерений,  вышедших и  не вышедших за допустимые пределы.

Расчёт математического ожидания и дисперсии выходной характеристики.

В случае стационарного эргодического процесса функционирования системы вычисление М(у) и Д(у) выходной характеристики у производится усреднением не по времени, а по множеству Nзнач., измеренных по одной реализации достаточной длительности. В целях экономия ОЗУ ЭВМ М(у) и Д(у) вычисляются по рекуррентным формулам:

   mn=mn-1*(n-1)/n + y/n;    (2)

где mn-1 - математическое ожидание, вычисленное на предыдущем шаге.

   dn=dn-1*(n-2)/(n-1) + 1/n*(yn-mn-1)2  (3)

здесь dn-1 - дисперсия, вычисленная на предыдущем шаге.

При большом количестве измерений эти оценки являются состоятельными и несмещёнными.

Расчёт среднего по времени значения выходной характеристики.

 Например, средняя длина очереди к каждому устройству вычисляется по формуле:

       (4)

где i - номер очередного изменения состояния очереди (занесение заявки в очередь или исключение из очереди); N - количество изменений состояния очереди; - интервал времени между двумя последними изменениями очереди.

Ёмкость накопитель:     (5)

где - ёмкость накопителя, занятая в интервале между двумя последними обращениями к накопителю для ввода-вывода заявки.

Построение гистограммы для стационарной системы.

 Г - эмпирическая плотность распределения вероятностей. Задаются границы изменения интересующей характеристики. уi[yнв], числом интервалов Ng. Определяется ширина интервала =( yн - ув)/Ng.

Затем в процессе моделирования по мере появления значений уi определяется число попаданий этой случайной величины в каждый из интервалов Ri гистограммы. По этим данным вычисляется относительная частота по каждому интервалу: Gi=Ri/(N*), где N - общее число измерений у. Площадь гистограммы равна единице, равна сумме площадей:

, т.к.

При необходимости выдвигается гипотеза о том, что эмпирическое распределение согласуется с некоторым теоретическим распределением. Эта гипотеза проверяется по тому или иному критерию. Например, при использовании критерия 2 в качестве меры расхождения используется выражение   

 (6);

где -  определяется из выбранного теоретического распределения вероятность попадания случайной величины в i-ый интервал.

  (7).

Из теоремы Пирсона следует, что для любой функции распределения F(y) случайной величины у при N распределения величины 2 имеет вид:

, где z - значение случайной величины 2 ,

k=Ng-(r +1) - число степеней свободы распределения 2 . r - количество параметров теоретического распределения, Г(к/2) - гамма функция.

Функция распределения 2 табулирована. По вычисленному значению 2 и числу степеней свободы с помощью таблиц определяется вероятность Р(2<Z). Если она превышает заданный уровень значимости С, то выдвинутая гипотеза принимается.

Моделирование систем с использованием типовых математических схем

Блочные иерархические модели процессов функционирования систем

Рассмотрим машинную модель Mm, системы S как совокупность блоков {mi}, i=1,2…n. Каждый блок модели можно охарактеризовать конечным набором возможных состояний {Z0}, в которых он может находиться. Пусть в течение рассматриваемого интервала времени (0,Т) блок i  изменяет состояние в моменты времени tijТ , где j - номер момента времени. Момент времени можно разделить на три группы:

случайные, связанные с внутренними свойствами блока;

случайные, связанные с изменением состоянием других блоков, имитирующая воздействие среды Е;

детерминированные моменты, связанные с заданным расписанием функционирования блока.

Моментами смены состояний модели Мм в целом t(k) Т будем считать все моменты изменения блоков {mi}, рис. 8.1.  см. ниже.

Рис. 8.1. Смена состояний модели для случаев 3-х блоков

При этом моменты ti(j) и tk являются моментами системного времени, т.е. времени, в котором функционирует система S. При машинной реализации модели Мм её блки представляются соответствующими программными модулями.

Особенности реализации процессов с использованием Q-схем

При моделировании Q-схем следует адекватно учитывать как связи, отражающие движения заявок (сплошные линии) так и управляющие связи (пунктирные линии).

Рассмотрим фрагмент Q-схемы (Рис. 8.2.):

Рис. 8.2. Фрагмент Q-схемы.

Примерами управляющих связей являются различные блокировки обслуживающих каналов (по входу и по выходу): "клапаны" изображены в виде треугольников, а управляющие связи пунктирными линиями. Блокировка канала по входу означает, что этот канал отключается от входящего потока заявок, а блокировка канала по выходу указывает, что заявка обслуженная блокированным каналом, остаётся в этом канале до момента снятия блокировки. В этом случае, если перед накопителем нет "клапана", то при его переполнении будут иметь место потери заявок.

Моделирующий алгоритм должен отвечать следующим требованиям:

обладать универсальностью относительно структуры, алгоритмов функционирования и параметров системы S;

обеспечивать одновременную и независимую работу системы S;

укладываться в приемлемые затраты ресурсов ЭВМ. (памяти, времени расчёта для реализации машинного эксперимента);

проводить разбиение на достаточно автономные логические части (блоки);

гарантировать выполнение рекуррентного правила расчётов;

При этом необходимо иметь виду, что появление одной заявки входящего потока в некоторый момент времени ti может вызвать изменение состояния не более чем одного из элементов Q-схемы, а окончание обслуживания заявки в момент ti в некотором канале К может привести в этот момент времени к последовательному изменению состояний нескольких элементов (Н,К), т.е. будет иметь место процесс распространения смены состояний в направлении противоположном движению заявки в системе S. Поэтому просмотр элементов Q-схемы должен быть противоположным движению заявок.

Все виды моделирующих алгоритмов Q-схемы можно классифицировать следующим образом (см. Рис. 8.3.):

Рис. 8.3. Виды моделирующих алгоритмов Q-схемы.

Алгоритмы моделирующие Q-схему по принципу "t" являются детерминированными (по шагу), а по принципу особых состояний – стохастические. Последние могут быть реализованы синхронным и асинхронным способами.

При синхронном способе один из элементов Q-схемы (И, Н или К) выбирается в качестве ведущего и по нему "синхронизируется" весь процесс моделирования.

При асинхронном способе — ведущий (синхронизирующий) элемент не используется, а очередному шагу моделирования (просмотру элементов Q-схемы) может соответствовать любое особое состояние всего множества элементов И, Н и К. При этом просмотр элементов Q-схемы организован так, что при каждом особом состоянии либо циклически просматриваются все элементы, спорадически - только те элементы, которые в этом случае могут изменить своё состояние. (просмотр с прогнозированием)

Построение и реализация моделирующих алгоритмов Q-схем

Прежде чем использовать какой либо язык для моделирования Q-схемы, необходимо глубже вникнуть в суть процесса построения и реализации М.А.

Пример. Рассмотрим Q-схему (Рис. 8.4.):

Рис. 8.4. Трехфазная Q-схема.

Примем обозначения:

Р - вероятность потери заявки (Р=N1/(N1+N3));

tm - время появления очередной заявки из источника;

tk,j - время окончания обслуживания заявки каналом Кк,j, k=1,2,3…; j=1,2…;

zi, zk,j - состояния накопителей  и каналов обслуживания;

tn - текущее время моделирования;

Li - ёмкость i-ого накопителя;

Lkm - число каналов в к-ой фазе;

N1, N2 - число выходных заявок;

Т - интервал моделирования;

При имитации Q-схемы на ЭВМ требуется организовать массив состояний:

zk,j, tk,j, j=1, Lkm; zi - число заявок в накопителе Hi; i=1,2; ti - i-ая заявка из источника.

zk,j = {1- канал занят; 0 - канал свободен; 2 - заблокирован};

Укрупнённая схема детерминированного МА Q-схемы, построенного по "принципу t" представлена на рисунке 8.5.

Рис. 8.5. Блок схема моделирования Q-схемы по принципу "t".

А далее более подробно рассмотрены алгоритмы блоков 4-9.

 

 


 

А также другие работы, которые могут Вас заинтересовать

39519. Здание 70 квартирного жилого дома в городе Слуцке 1.86 MB
  Удельный вес монтажных работ в строительстве увеличивается с каждым годом. Наряду со снижением массы отдельных конструкций они укрупняются и доводятся до максимальной заводской готовности. Процесс монтажа базируется также на дальнейшем росте уровня комплексной механизации, развитие автоматизации, применение прогрессивных технологий.
39520. Управление строительством многоквартирного жилого дома №7 в программной среде Rillsoft Project 3.08 MB
  Куцепалова Объем дипломного проекта: дипломная работа 82 страницы графическая часть . Куцепалова Объем дипломного проекта: дипломная работа 82страницы графическая часть . Целью проекта является исследование использования в современном строительном бизнесе информационных технологий и специализированного программное обеспечения. Областью возможного практического применения ПО Rillsoft Project 2007 являются оценка строительного проекта с точки зрения объемов работ стоимости общей потребности в ресурсах календарного плана работ графика...
39522. Восьмиэтажное административное здание – здание филиала «Приорбанка» на проспекте Победителей в городе Минске 1007 KB
  Вопросы экономики железобетонных конструкций следует решать совместно с вопросами прочности на протяжении всего процесса проектирования: при выборе объемнопланировочной и конструктивной схемы здания; членении конструкций на сборные элементы; выборе формы и размеров сечения элементов; назначении класса бетона класса стальной арматуры.пр  Nпр  RbtАb  pАпр  м из условия анкеровки арматуры колонны hoф20dпрк2000320640 м Принимаем hoф1 м Полная высота фундамента равна: hф hoфa1007107 м Принимаем hф1.4 Расчёт...
39523. Экономическое обоснование цены реализации строительной продукции 829.5 KB
  Составление локальной сметы на общестроительные работы. Составление акта приёмки выполненных работ. Составление локальной сметы на общестроительные работы Сметная документация в строительстве составляется на основе Методических указаний по определению стоимости строительства зданий и сооружений и составлению сметной документации с применением РСН РДС 8.
39524. Процесс проектирования металлических конструкций 1.84 MB
  Подбор сечения затяжки и проверка ее на прочность Наибольшее усилие в затяжке при полной расчетной нагрузке Nз=2137 кН; lз=52 м длина затяжки; Ез=195 ГПа модуль упругости затяжки; Rз=1200 МПа расчетное сопротивление материала затяжки в качестве затяжки принимаем высокопрочные канаты французской фирмы Freyssinet.2 Подбор сечений стержней фермы В качестве материала фермы принята сталь С375 с расчетным сопротивлением стали Ry=365 МПа при толщине в интервале от...
39525. Промышленное и гражданское строительство 369.5 KB
  Методические указания содержат рекомендации по выполнению экономической части дипломных проектов и учитывают особенности определения техникоэкономических показателей строительных конструкций а также сметной стоимости строительства на основе ресурсносметных норм. БНТУ 2006 ВВЕДЕНИЕ Методические указания разработаны для студентов специальности 170 02 01 Промышленное и гражданское строительство и могут быть использованы при расчетах экономической части дипломных проектов а также при курсовом проектировании по дисциплинам...
39526. Разработка стратегии развития и совершенствование нынешней системы стратегического управления СУ-16 2.78 MB
  Если в прошлом многие компании могли весьма успешно функционировать обращая внимание в основном на внутренние проблемы связанные с повышением эффективности использования ресурсов в текущей деятельности то сегодняшнее развитие рыночных отношений делает необходимым изменение сложившихся стереотипов хозяйствования характера управления.Есть ли изменения существенные для фирмы Так туристические и транспортные компании постоянно оценивают динамику цен на топливо. Например реклама коттеджей это тактика строительной компании выбравшей...
39527. Инвестиционный проект предприятия малого бизнеса 7.76 MB
  В настоящее время основными видами деятельности акционерного общества являются производство строительно-монтажных работ, проектно-изыскательские работы, торгово-закупочная деятельность в области жилищного и производственного строительства. Производственная деятельность акционерного общества осуществляется в сотрудничестве с организациями