- •Построение математических моделей методом регрессионного анализа
- •Идентификация линейной функции.
- •Идентификация внутренне линейных функций.
- •Моделирование по схеме марковских случайных процессов
- •Марковские цепи с дискретным временем
- •Марковские цепи с непрерывным временем
- •Предельные вероятности состояний
- •Системы массового обслуживания
- •Системы массового обслуживания с отказами
- •Системы массового обслуживания с ожиданием
- •Агрегативное описание систем
- •Моделирование стохастических процессов методом статистических испытаний
Агрегативное описание систем
Лекции. Разделы 11.1-11.3.
Под агрегатом будем понимать объект <T,X,U,Y,Z,H,G>, где T – фиксированное подмножество действительных чисел, X,U,Y,Z – множества любой природы, H,G – операторы (вообще говоря, случайные) переходов и выходов. Элементы указанных множеств назовем: tT – моментом времени; xX - входным сигналом; uU- управляющим сигналом; yY – выходным сигналом; zZ - состоянием. Состояния, входные, управляющие и выходные сигналы, рассматриваемые как функции времени, обозначим z(t), x(t), u(t) и y(t). Операторы переходов и выходов H и G реализуют функции z(t) и y(t).
Для агрегатов справедливы следующие два предположения.
Предположение 1. За конечный интервал времени в агрегат поступает конечное число входных и управляющих сигналов и вырабатывается конечное число выходных сигналов.
Предположение 2. Из особых состояний агрегат может переходить в новое состояние скачком.
Под особыми состояниями будем понимать его состояния в момент получения входного либо управляющего сигналов или выдачи выходного сигнала. Все остальные состояния агрегата будем называть неособыми.
Рассмотрим оператор переходов агрегата.
Пусть z(t*) – некоторое особое состояние агрегата, а us – последний управляющий сигнал usU. Примем следующие обозначения для операторов, являющихся частными видами оператора H и определяющих состояние агрегата в момент t*+0. Если t* - момент поступления входного сигнала x, то
z(t*+0)=V[z(t*),x,us] (5.1)
Аналогично, если t* - момент поступления управляющего сигнала u, то
z(t*+0)=V[z(t*),u] (5.2)
При одновременном поступлении x и u
z(t*+0)=V [z(t*),x,u] (5.3)
Наконец, если t* - момент выдачи выходного сигнала y, то
z(t*+0)=W[z(t*),us] (5.4)
В интервале между особыми состояниями, значение z(t) определяется при помощи операторов U, вид которых в общем случае зависит от особого состояния, являющегося для данного интервала времени начальным состоянием:
z(t*+0)=
[t,z(t*+0),us]
(5.5)
Здесь t* - момент особого состояния, являющегося исходным для данного интервала времени. В общем случае операторы U,V,V,V и W являются случайными операторами.
Теперь рассмотрим оператор выходов агрегата. Во множестве Z состояний z(t) агрегата выделим класс подмножеств {Zy}, обладающих следующими свойствами. Выходной сигнал y выдается в момент t в тех случаях, когда: 1) z(t)Zy; z(t-0)Zy и 2) z(t+0)Zy, но z(t’)Zy. Тогда, оператор выходов G можно представить в виде совокупности двух операторов: G, вырабатывающего выходной сигнал
y=G[z(t),us] (5.6)
и G, проверяющего для каждого t принадлежность z(t) к одному из подмножеств Zy. В общем случае, оператор G является случайным оператором. Это значит, что данным t,z(t),u ставится в соответствие не одно определенное значение выходного сигнала, а некоторое множество значений y с соответствующим распределением вероятностей, задаваемых оператором G.
1.
Постановка задачи.
Описать в общем виде процесс функционирования агрегата в интервале [t0,T].
t0 – начальный момент времени, состояние агрегата z0, значение управляющего сигнала u0.
t1>t0 – поступает управляющий сигнал u1.
t2>t1 – поступает входной сигнал x1 и выдаётся выходной сигнал y1
t2<t3<T – выдаётся выходной сигнал y2 (во всех остальных случаях оснований для выдачи выходного сигнала нет).
Решение.
1. В полуинтервале (t0,t1] изменение состояния агрегата описываются оператором (5.5) (промежуток времени между особыми состояниями, начальное состояние z0 условно считаем начальным особым состоянием этого полуинтервала).
z(t)=
[t,z0,u0].
(t0<tt1)
2. В момент t1 поступает управляющий сигнал u1:
z(t1+0)=V [z(t1), u1].
3. В полуинтервале (t1,t2] снова действует оператор (5.5):
z(t)=
[
t,z(t1+0),u1].
(t1<tt2)
4. В момент t2 поступает входной сигнал x1 и выдаётся выходной сигнал y1. В первую очередь выдается выходной сигнал:
y1=G [z(t2), u1].
Состояние агрегата меняется скачком, которое описывается оператором (5.4):
z(t2+0)=W[z(t2), u1].
После этого происходит изменение состояния агрегата при принятии входного сигнала x1, которое описывается оператором (5.1):
z(t2+0+0)=V [z(t1+0), x1,u1].
5. В полуинтервале (t2,t3] действует оператор (5.5):
z(t)=
[
t,z(t2+0+0),
u1].
(t2<tt3)
6. В момент t3 выдаётся выходной сигнал y2:
y2=G [z(t3), u1]
Cостояние z меняется:
z(t3+0)=W[z(t3), u1].
7. В полуинтервале (t3,T] снова действует оператор (5.5):
z(t)=
[
t,z(t3+0),
u1]. (t3<tT)
2.
Постановка задачи.
Описать процесс функционирования следующей однолинейной системы массового обслуживания как агрегата. В моменты времени tj, представляющие собой случайный поток однородных событий с заданным законом распределения, в систему поступают заявки, каждая из которых характеризуется параметром aj. Если линия обслуживания занята, заявка попадает в очередь и может находиться там (ждать) не более чем τj(ж)=f(aj,ui), где ui – значение управляющего параметра. Величины τj образуют случайный поток однородных событий с заданным законом распределения. Заявка должна быть принята к обслуживанию не позднее, чем в момент времени tj +τj(ж), в противном случае она получает отказ. Время обслуживания заявки или время занятости линии τj(З)=s(aj,ui). Исследуемый интервал времени [t0,T].
t0 – система занята обслуживанием заявки x1, в очереди ещё три заявки (заявки обслуживаются в порядке очерёдности), значение управляющего параметра u1.
t1>t0 – закончилось время ожидания заявки x3 (заявка получила отказ).
t1<t2<T поступает заявка x5 с параметром a5.
Решение.
Представим систему массового обслуживания в виде агрегата. Состояние агрегата будем характеризовать вектором z(t)Z со следующими компонентами:
z1(t) – время, оставшееся до конца обслуживания заявки (для заявки, которая находится на обслуживании);
z2(t) – значение величины ui;
z3(t) – количеств заявок в очереди;
z4(t) – оставшееся время ожидания первой заявки из очереди (время до момента отказа);
z5(t) – оставшееся время ожидания второй заявки из очереди (время до момента отказа);
z6(t) – оставшееся время ожидания третьей заявки из очереди (время до момента отказа).
1. В полуинтервале (t0,t1] изменение состояния агрегата описывается оператором (5.5). Для нашего случая он имеет следующий вид:
2. В момент времени t1 выдаётся выходной сигнал y:
,
где
- представляет собой признак – заявка
получила отказ,
- выходной параметр заявки.
Состояние агрегата в этот момент изменяется скачком, который описывается оператором (5.4). В нашем случае оператор W имеет следующий вид:
3. В полуинтервале (t1,t2] изменение состояния агрегата описывает оператор U (5.5):
4. В момент времени t2 поступает заявка x5, состояние агрегата изменяется скачком под действие оператора (5.1). Оператор V' в нашем случае имеет следующий вид :
5. В полуинтервале (t2,T] изменение состояния агрегата описывает оператор U (5.5):
Практическое занятие 6
