Математическое описание системы............................................................................ 18

Результаты моделирования и аналитического решения............................................ 19

Заключение......................................................................................................................... 20

Литература.......................................................................................................................... 21

Приложения........................................................................................................................ 22

Текст программы........................................................................................................... 22

Введение.

В настоящей курсовой работе рассматривается проблема моделирования процессов в Q-схемах – одном из важнейших, с точки зрения применения на практике, классов математических схем, разработанных для формализации процессов функционирования систем массового обслуживания (СМО) в теории массового обслуживания. Предметом изучения в теории массового обслуживания являются системы, в которых появление заявок (требований) на обслуживание и завершение обслуживания происходит в случайные моменты времени, т.е. характер их функционирования носит стохастический характер. Следует отметить, что СМО описывают различные по своей физической природе процессы функционирования экономических, производственных, технических и других систем, например потоки поставок продукции некоторому предприятию, потоки деталей и комплектующих изделий на сборочном конвейере цеха, заявки на обработку информации в ЭВМ от удаленных терминалов и т.д.

САПР состоит из ЭВМ и трех терминалов. Каждый из проектировщиков формирует задания на расчет в интерактивном режиме. Набор строки задания занимает 10±5с. После набора 10 строк задание считается сформированным и поступает на рещение, при этом в течение 10±3с ЭВМ прекращает принимать другие задания. Анализ результата занимает у проектировщиков 30 с, после чего цикл повторяется. Данные по всем проектировщикам одинаковы.

Смоделировать работу системы в течение 6 часов. Определить вероятность простоя проектировщика из-за занятости ЭВМ и коэффициент загрузки ЭВМ.

Для изучения Q-схем используются два подхода: аналитический и имитационный. При аналитическом подходе подлежащая анализу схема описывается с помощью формул, отражающих зависимости между ее различными параметрами. Однако, следует отметить, что разработанные методы аналитического изучения Q-схем подходят далеко не для каждой конкретной системы, они пригодны лишь для систем общего типа. Поэтому при аналитическом изучении систем их необхродимо упрощать до систем основных типов, что в последствии конечно-же сказывается на результатах исследования.  При имитационном подходе ставится эксперимент на машинной модели системы, которая предварительно реализуется на одном из созданных специально для этого языков имитационного моделирования (например, SIMULA, SIMSCRIPT, GPSS и др.) или на языке общего назначения (BASIC, PASCAL, FORTRAN, C++ и др.).

Учитывая условие, построим структурную схему данной системы.

Рис. 1 Структурная блок-схема системы.

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

Рис. 2 Временная диаграмма.

На временной диаграмме:

- оси 1, 2, 3 – возникновение заявки  соответственно у 1-го, 2-го или 3‑го проектировщика;

- ось 4 – обработка заявок проектировщиков на ЭВМ.

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

Так как, по сути, описанные процессы являются процессами массового обслуживания, то для формализации задачи используем символику Q‑схем [2]. В соответствии с построенной концептуальной моделью и символикой Q‑схем структурную схему данной СМО (рис. 1) можно представить в виде, показанном на рис. 3, где И – источник, К – канал.

Рис. 3 Структурная схема системы в символике Q-схем.

Источники И1, И2, И3 имитируют поступление заявок от проектировщиков 1,2 и 3 соответственно. Канал К1 имитирует процесс обработки заявок на центральной ЭВМ. Если канал К1 занят, то клапан 1 закрыт. Источники генерируют заявки, идущие затем на ЭВМ. Если ЭВМ занята, то заявка остается в источнике дожидаться своей очереди на обработку.

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

После формализации задачи можно переходить к построению моделирующего алгоритма.

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

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

- обеспечивать одновременную и независимую работу необходимого числа элементов схемы;

- укладываться в приемлемые затраты ресурсов ЭВМ для реализации машинного эксперимента;

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

- гарантировать выполнение рекуррентного правила – событие, происходящее в момент времени tk может моделироваться только после того, как промоделированы все события, произошедшие в момент времени, произошедшие в момент времени tk-1< tk.

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

Известно, что существуют два основных принципа построения моделирующих алгоритмов: «принцип Dt» и «принцип dz». При построении моделирующего алгоритма Q-схемы по «принципу Dt», т.е. алгоритма с детерминированным шагом, необходимо для построения адекватной модели определить минимальный интервал времени между соседними состояниями Dt’ = {ui} (во входящих потоках и потоках обслуживания) и принять, что шаг моделирования равен Dt’. В моделирующих алгоритмах, построенных по «принципу dz», т.е. в алгоритмах со случайным шагом, элементы Q-схемы просматриваются при моделировании только в моменты особых состояний (в моменты появления из И изменения состояний К). При этом длительность шага Dt = var зависит как от особенностей самой системы, так и от воздействий внешней среды. Моделирующие алгоритмы со случайным шагом могут быть реализованиы синхронным и асинхроным способами. При синхронном способе один из элементов Q-схемы выбирается в качестве ведущего, и по нему «синхронизируется» весь процесс моделирования. При асинхронном способе построения моделирующего алгоритма ведущий (синхронизирующий) элемент не используется, а очередному шагу моделирования (просмотру элементов Q-схемы) может соответствовать любое особое состояние всего множества элементов И, Н, К. при этом просмотр элеменов Q-схемы организовани так, что при каждом особом состоянии либо циклически просматриваются все элементы, либо спорадически, - только те, которые могут в этом случае изменить свое состояние.

Разработку моделирующего алгоритма удобно производить в 2 этапа:

1) разработка укрупненного алгоритма;

2) разработка детального алгоритма.

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

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

На рис. 4 изображена укрупненная схема моделирующего алгоритма.

Рис. 4 Укрупненная схема моделирующего алгоритма.

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

Для описания элементов схемы использовался объектно-ориентированный подход, основными принципами которого являются

- объединение в рамках одной структуры данных полей и методов (процедур), оперирующих над ними;

- наследование – порождение новых классов объектов, при этом наследники получают те же поля и методы, что были объявлены в классе непосредственного предка и его предков;

- полиморфизм – возможность изменения (частичного или полного) методов одноименных с методами предков (т.н. виртуальные методы).

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

Таким образом, в ходе построения моделирующего алгоритма были выстроены следующие классы объектов.

TQSheme – класс схемы. На него возложены основные функции по проведению эксперимента, а именно:

- управлению системным временем

- нахождением порядка опроса элементов в зависимости от  структуры схемы

- опросу элементов в каждом цикле

- учету заявок, находящихся в системе

- учету особых состояний, происходящих в системе

Так, например, именно класс TQSheme реализует блоки 2 и 3 укрупненного алгоритма при помощи своих методов соответственно InitEmulation и Analize, а также блоки 4-7 при помощи метода Emulation. Блок-схемы этих методов приведены ниже.

Рис. 5 Блок-схема процедуры TQSheme.InitEmulation.

Описание:

- блок 1 – происходит инициализация переменных, необходимых для анализа системы;

- блок 2 – создание объектных списков, необходимых для анализа системы: установление конечных элементов, установление порядка очередности опроса элементов схемы;

- блок 3 – инициализация списков заявок и событий, подготовка их к предстоящему прогону схемы.

Рис. 6 Блок-схема процедуры TQSheme.Analize.

Описание:

- блок 1 – нахождение порядка опроса элементов с занесением порядковых номеров элементов в массив порядка опроса FOrder;

- блок 2 – нахождение наименьшего временного интервала, необходимого для анализа схемы по «принципу Dt»

Рис. 7 Блок-схема процедуры TQSheme.Emulation.

Описание:

- блок 1 – процедура инициализации процесса моделирования с установлением начальных значений для всех переменных;

- блок 2 – вызов процедуры Analize;

- блок 3 – проверка окончания моделирования;

- блок 4 – просмотр всех элементов схемы.

Рис. 8 Блок-схема процедуры TQSheme.Step

Описание:

- блок 1 – процедура изъятия заявок из конечных элементов схемы;

- блок 2 – опрос всех элементов схемы в порядке, указанном в массиве FOrder;

- блок 3 – увеличение системного времени на величину Dt;

Каждый элемент схемы, будь то источник, накопитель или канал, также представлен соответствующим классом (TSource, TAccumulator, TChannel). Однако все классы элементов схемы являются наследниками одного общего класса TElement, который обеспечивает общую функциональность для всех типов элементов схемы, а именно:

- содержит информацию о том, из каких элементов схемы заявки поступают на данный элемент;

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

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

- определяет порядок, по которому заявки передаются данному элементу от элементов-источников. (Алгоритм метода AskForParcels показан на рис. 4. Этот метод вызывается объектом класса TQSheme, которому принадлежит этот элемент на этапе моделирования перехода заявок).

Классы TSource, TAccumulator, TChannel доопределяют функции своего общего предка TElement для обеспечения более конкретного поведения, характерного для элементов соответствующего класса.

Рис. 9 Блок-схема процедуры TElement.AskForParcel.

Описание:

- блок 1 – задается цикл по всем источникам для данного элемента;

- блок 2, 3 – функции определения выражений для клапанов, установленные пользователем;

- блок 4 – проверка возможности принятия данным элементом заявки;

- блок 5 – прием заявки;

- блок 6 – отказ заявке в приеме.

Главное окно программы показано на рисунке 10.

Рис. 10 Главное окно программы.

Данная система представляет собой одноканальную СМО с отказами.

Интенсивность потока заявок для нее будет определяться следующим выражением:

l = l1+l2+l3   (1),

где l1,l2,l3 – интенсивность потока заявок от каждого проектировщика соответственно. С учетом того, что данные по всем проектировщикам одинаковы выражение (1) примет следующий вид:

l=3l1 (2).

Интенсивность потока заявок найдем с помощью выражения

где
– среднее время поступления заявок от одного проектировщика, которое в свою очередь с учетом условий задачи вычисляется как

 (3),

где t1 – среднее время набора проектировщиком одной строки задания,

t2 – время анализа результата проектировщиком

С учетом выражений (2),(3) выражение (1) примет следующий вид:

l=3/(10t1+ t2) (4).

Интенсивность потока обслуживаний для данной системы определяется по формуле

 (5),

где
– среднее время обработке ЭВМ одной задания.

Чтобы найти пропускную способность ЭВМ, воспользуемся формулой для пропускной способности одноканальной СМО с отказами:

Q = m/(m +l)  (6)

Чтобы найти вероятность простоя проектировщика, воспользуемся формулой для нахождения вероятности отказа в одноканальной СМО:

Ротк=l/(m +l)  (7)

Результаты моделирования и аналитического решения.

Воспользовавшись выраженияи (6) и (7), найдем коэффициент загрузки ЭВМ и вероятность простоя проектировщика из-за занятости ЭВМ.

Коэффициент загрузки ЭВМ Q=0,8125

Вероятность простоя проектировщика из-за занятости ЭВМ Р=0,1875.

При моделировании работы системы на ЭВМ были получены следующие результаты:

Коэффициент загрузки ЭВМ Q=0,256

Вероятность простоя проектировщика из-за занятости ЭВМ Р=0,1.

Несовпадение результатов аналитических расчетов и машинного моделирования объясняется тем, что в для расчета аналитическим методом данная система была упрощена и приведена к виду одноканальной СМО с отказами, что не учитывает всех особенностей функционирования исходной системы.

Результате данной работы стало построение программы, моделирующей процесс функционирования заданной ситсемы. Были рассчитаны (аналитически и при помощи построенного моделирующего алгоритма) показатели эффективности данной системы: коэффициент загрузки и вероятность простоя проектировщика из-за занятости ЭВМ. Выявлены основные закономерности и способы взаимодействия элементов Q-схем, а также причины несовпадения рассчетных показателей с результатами прогона моделирующего алгоритма на ЭВМ.

1. Кремер «Исследование операций в экономике» -М.: «Экономика»1997 г.

2. Советов Б.Я., Яковлев С. А. Моделирование систем. - М.:ВШ,1995.

3. Советов Б.Я., Яковлев С.А. Моделирование систем. Практикум. - М.:ВШ,1999.

4. Вентцель Е.С. Теория вероятностей. –М.:Наука, 1969.

5. Вентцель Е.С. Исследование операций. – М.:Сов. Радио, 1972.

Текст программы.

unit QSheme;

interface

uses Classes, SysUtils, Forms, QSObjs, QSTypes, Utils, ExtCtrls;

const

  MaxElementCount = 10000;

type

  TQSheme = class

  private

    FParcelsClass : TParcelsClass;

    FStepCount : integer;

    FSysPeriod : TCustTime;

    FStepPeriod : TCustTime;

    FSysTime : TCustTime;

    FElements : TList;

    FFinishElements : TList;

    FLastElements : TList;

    FSourceElements : TList;

    FParcels : TList;

    FOrderElementCount : integer;

    FEventQueue : TList;

    FOrder : array [0..MaxElementCount-1] of integer;

    FDiagram : TPaintBox;

  protected

    function GetElement(Index : integer): TElement;

    function GetElementCount: integer;

    procedure InitAnalize;

    procedure FreeAnalize;

    procedure InitEmulation;

    procedure FreeEmulation;

    procedure FindFinishElements;

    procedure GetRecipientsOrder;

    procedure Step;

    procedure TakeParcelsFromFinishElements;

    function IsFinishElement(Element : TElement): Boolean;

    function FastestStepPeriod : TCustTime;

    procedure ClearParcelList;

    procedure ClearEventQueue;

    function GetCounts(Index : integer): integer;

    function GetParcelCount: integer;

    procedure DrawElementLines;

    procedure DisplayEvents;

  public

    function NewParcel: Pointer;

    procedure NewEvent(AEvent : integer;  ASender, ASource : TObject; AInfo : TInfo);

    procedure RedrawDiagram;

    procedure Emulation;

    procedure Analize;

    constructor Create;

    destructor Destroy; override;

    procedure AddElement(Element : TElement);

    procedure GetOrder;

    procedure DelElement(Element : TElement);

    property SysTime : TCustTime read FSysTime;

    property SysPeriod: TCustTime read FSysPeriod write FSysPeriod;

    property StepPeriod : TCustTime read FStepPeriod write FStepPeriod;

    property Counts[Index : integer]:integer read GetCounts;

    property BornParcelCount : integer index 0 read GetCounts;

    property StoreParcelCount : integer index 1 read GetCounts;

    property WorkParcelCount : integer index 2 read GetCounts;

    property PassedParcelCount : integer index 3 read GetCounts;

    property RefusedParcelCount : integer index 4 read GetCounts;

    property ParcelCount:integer read GetParcelCount;

    property StepCount : integer read FStepCount write FStepCount;

    property ParcelsClass : TParcelsClass read FParcelsClass write FParcelsClass;

  published

    property Diagram : TPaintBox read FDiagram write FDiagram;

    property ElementCount : integer read GetElementCount;

    property Elements[Index : integer] : TElement read GetElement;

  end;

implementation

uses MainFrm;

constructor TQSheme.Create;

begin

  FElements := TList.Create;

  FParcelsClass := TParcel;

  FParcels := TList.Create;

  FEventQueue := TList.Create;

end;

destructor TQSheme.Destroy;

begin

  FElements.Free;

  ClearEventQueue;

  FEventQueue.Free;

  ClearParcelList;

  FParcels.Free;

  inherited;

end;

function TQSheme.GetElement(Index : integer): TElement;

begin

  Result := FElements[Index];

end;

function TQSheme.GetElementCount: integer;

begin

  Result := FElements.Count

end;

procedure TQSheme.AddElement(Element: TElement);

begin

  if Assigned(Element) then begin

    FElements.Add(Element);

    Element.Sheme := Self;

  end;

end;

procedure TQSheme.DelElement(Element: TElement);

var i,j : integer;

begin

  if Assigned(Element) then begin

    for i := 0 to ElementCount - 1 do

      for j:= Elements[i].SourceCount-1 downto 0 do

        if Elements[i].Sources[j] = Element then

          Elements[i].DelSource(Element);

    FElements.Remove(Element);

  end;

end;

function TQSheme.IsFinishElement(Element: TElement):Boolean;

var j,s : integer;

begin

  Result := False;

  for j := 0 to ElementCount-1 do begin

    for s := 0 to Elements[j].SourceCount-1 do begin

      if Element = Elements[j].Sources[s] then Exit;

    end;

  end;

  Result := True;

end;

procedure TQSheme.FindFinishElements;

var i : integer;

begin

  for i := 0 to ElementCount-1 do

    if IsFinishElement(Elements[i]) then begin

      FFinishElements.Add(Elements[i]);

      FLastElements.Add(Elements[i]);

    end;

end;

function TQSheme.FastestStepPeriod: TCustTime;

var i : integer;

    Min : TCustTime;

begin

  Min := FSysPeriod;

  for i := 0 to ElementCount-1 do

    if (Elements[i] is TShop) then

      with TShop(Elements[i]).Generator do

        if Mean-Disp < Min then Min := Mean-Disp;

{$ifndef Precision}

  Result := Min;

{$else}

  Result := Min div 10;

{$endif}

end;

procedure TQSheme.InitAnalize;

begin

  FSysTime := 0;

  FStepCount := 0;

  FOrderElementCount := 0;

  FLastElements := TList.Create;

  FSourceElements := TList.Create;

end;

procedure TQSheme.FreeAnalize;

begin

  FLastElements.Free;

  FSourceElements.Free;

end;

procedure TQSheme.GetRecipientsOrder;

var i,s : integer;

    LastElement : TElement;

begin

  if FLastElements.Count = 0 then Exit;

  for i := 0 to FLastElements.Count-1 do begin

    LastElement := TElement(FLastElements[i]);

    FOrder[FOrderElementCount] := FElements.IndexOf(LastElement);

    Inc(FOrderElementCount);

    for s := 0 to LastElement.SourceCount - 1 do

      if FSourceElements.IndexOf(LastElement.Sources[s])<0 then

        FSourceElements.Add(LastElement.Sources[s]);

  end;

  SwapPointers(Pointer(FSourceElements),Pointer(FLastElements));

  FSourceElements.Clear;

  GetRecipientsOrder;

end;

procedure TQSheme.GetOrder;

begin

  FindFinishElements;

  GetRecipientsOrder;

end;

procedure TQSheme.TakeParcelsFromFinishElements;

var i : integer;

    Parcel : TParcel;

begin

  for i := 0 to FFinishElements.Count-1 do

    with TElement(FFinishElements[i]) do

      if CanDrop then begin

        Parcel := Container;

        NewEvent(EV_PASS,nil,FFinishElements[i],Parcel.Info);

        DoBeforeDrop(FFinishElements[i]);

        DropParcel;

        DoAfterDrop(FFinishElements[i]);

        Parcel.State := psPassed;

      end;

end;

procedure TQSheme.Step;

var i : integer;

begin

  TakeParcelsFromFinishElements;

  for i := 0 to FOrderElementCount-1 do Elements[FOrder[i]].AskForParcel;

  Form1.Gauge1.Progress := Round(FSysTime/FSysPeriod*100);

  Inc(FSysTime,FStepPeriod);

  Inc(FStepCount);

end;

procedure TQSheme.Analize;

begin

  try

    try

      InitAnalize;

      GetOrder;

      FStepPeriod := FastestStepPeriod;

    finally

      FreeAnalize;

    end;

  except

    on EInvalidPointer do raise;

  end;

end;

procedure TQSheme.ClearEventQueue;

var i : integer;

begin

  if Assigned(FEventQueue) then begin

    for i := 0 to FEventQueue.Count - 1 do FreeMem(FEventQueue[i],SizeOf(TEventRec));

    FEventQueue.Clear;

  end;

end;

procedure TQSheme.ClearParcelList;

var i : integer;

begin

  if Assigned(FParcels) then begin

    for i := 0 to FParcels.Count - 1 do TParcel(FParcels[i]).Free;

    FParcels.Clear;

  end;

end;

procedure TQSheme.InitEmulation;

var i : integer;

begin

  ClearParcelList;

  ClearEventQueue;

  for i := 0 to ElementCount - 1 do

    Elements[i].ClearContainer;

  FFinishElements := TList.Create;

end;

procedure TQSheme.FreeEmulation;

begin

  FFinishElements.Free;

end;

procedure TQSheme.Emulation;

begin

  try

    InitEmulation;

    Analize;

    while FSysTime < FSysPeriod do Step;

    Form1.Gauge1.Progress := 100;

    //RedrawDiagram;

  finally

    FreeEmulation;

  end;

end;

function TQSheme.NewParcel: Pointer;     

var P : Pointer;

begin

  P := FParcelsClass.Create;

  FParcels.Add(P);

  Result := P;

end;

procedure TQSheme.NewEvent(AEvent : Integer; ASender, ASource: TObject; AInfo : TInfo);

var P : PEventRec;

begin

  GetMem(P,SizeOf(TEventRec));

  with P^ do begin

    Event := AEvent;

    Sender := ASender;

    Source := ASource;

    Info := AInfo;

    SysTime := FSysTime;

  end;

  FEventQueue.Add(P);

end;

function TQSheme.GetCounts(Index : integer): integer;

var i : integer;

begin

  Result := 0;

  for i := 0 to FParcels.Count-1 do

    if Ord(TParcel(FParcels[i]).State) = Index then Inc(Result);

end;

function TQSheme.GetParcelCount: integer;

begin

  Result := FParcels.Count;

end;

const //DrawConstants

  Top = 20;

  Left = 20;

  Interval = 20;

procedure TQSheme.DrawElementLines;

var i : integer;

    Y : integer;

begin

  for i := 0 to ElementCount-1 do begin

    Y :=Top + interval *i;

    with Diagram.Canvas do begin

      TextOut(0,Y + Font.Height,Elements[i].Name);

      MoveTo(0,Y);

      LineTo(Diagram.ClientWidth,Y)

    end;

  end;

end;

procedure TQSheme.DisplayEvents;

{var i : integer;

    s : string;}

begin

  {Form1.mResults.Items.Clear;

  for i := 0 to FEventQueue.Count - 1 do begin

    with TEventRec(FEventQueue[i]^) do begin

      case Event of

        EV_TAKE: S := '+++:';

        EV_REFUSE: S := '------:';

        EV_PASS: S := 'PASS:';

      end;

      S := S + IntToStr(Info);

      S := S + '['+IntToStr(SysTime)+ '] ';

      if Assigned(Source) then S := S + TElement(Source).Name

      else S := S+'nil';

      S := S + '->';

      if Assigned(Sender) then S := S + TElement(Sender).Name

      else S := S+'nil';

    end;

    Form1.mResults.Items.Add(S);

  end;}

end;

procedure TQSheme.RedrawDiagram;

//var i : integer;

begin

  //Diagram.Canvas.FillRect(Rect(0,0,Diagram.Width,Diagram.Height));

  //DrawElementLines;

  DisplayEvents;

end;

initialization

  Randomize;

end.

unit QSObjs;

interface

uses Classes,QSTypes,SysUtils, Utils;

type

  TElement = class;

  TIsRightElement = function(Element : TElement): Boolean of object;//far;

  TBeforeAfterAction = procedure (Sender : TElement) of object;

  TElement = class

  private

    FId : integer;

    FName : string;

    FSources : TList;

    FSheme : TObject;

    FContainer : TParcel;

    FOnSourceValidate : TIsRightElement;

    FOnDestinationValidate : TIsRightElement;

    FBeforeTake: TBeforeAfterAction;

    FAfterTake: TBeforeAfterAction;

    FBeforeDrop: TBeforeAfterAction;

    FAfterDrop: TBeforeAfterAction;

    procedure SetSheme(ASheme : TObject);

    function GetSourceCount: integer;

    function GetSource(Index : integer): TElement;

    function GetParcelPresent: Boolean;

    function GetCanDropParcelFor(Destination : TElement): Boolean;

    function GetCanTakeParcelFrom(Source: TElement): Boolean;

    procedure Pass(SourceIndex : integer); virtual;

  protected

    function GetCanTake: Boolean; virtual; abstract;

    function GetCanDrop : Boolean; virtual; abstract;

  public

    constructor Create;virtual;

    destructor Destroy; override;

    procedure AddSource(Element : TElement);

    procedure DelSource(Element : TElement);

    procedure AskForParcel; virtual;

    procedure ClearContainer; virtual;

    procedure RefuseParcel(SourceIndex : integer);

    procedure DropParcel;virtual;

    procedure TakeParcel(SourceIndex : integer); virtual;

    procedure DoBeforeDrop(Sender : TElement);

    procedure DoBeforeTake(Sender : TElement);

    procedure DoAfterDrop(Sender : TElement);

    procedure DoAfterTake(Sender : TElement);

    property CanDropParcelFor[Destination : TElement]: Boolean read GetCanDropParcelFor;

    property CanTakeParcelFrom[Source : TElement]: Boolean read GetCanTakeParcelFrom;

    property Container : TParcel read FContainer write FContainer;

    property ParcelPresent : Boolean read GetParcelPresent;

    property CanTake : Boolean read GetCanTake;

    property CanDrop : Boolean read GetCanDrop;

    property Id: integer read FId write FId;

  published

    property Name : string read FName write FName;

    property Sheme: TObject read FSheme write SetSheme;

    property SourceCount : integer read GetSourceCount;

    property Sources[Index : integer]: TElement read GetSource;

    property OnSourceValidate : TIsRightElement read FOnSourceValidate write FOnSourceValidate;

    property OnDestinationValidate : TIsRightElement read FOnDestinationValidate write FOnDestinationValidate;

    property BeforeTake: TBeforeAfterAction read FBeforeTake write FBeforeTake;

    property AfterTake: TBeforeAfterAction read FAfterTake write FAfterTake;

    property BeforeDrop: TBeforeAfterAction read FBeforeDrop write FBeforeDrop;

    property AfterDrop: TBeforeAfterAction read FAfterDrop write FAfterDrop;

  end;

  TElementClass = class of TElement;

  TGenerator = class

  private

    FMean : TCustTime;

    FDisp : TCustTime;

    FRandomType : TRandomType;

    function GetRandom: TCustTime;

  public

    constructor Create;

    property Mean : TCustTime read FMean write FMean;

    property Disp : TCustTime read FDisp write FDisp;

    property RandomType : TRandomType read FRandomType write FRandomType;

    property Time : TCustTime read GetRandom;

  end;

  TShop = class(TElement)

  private

    FGenerator : TGenerator;

    FEndWorkTime : TCustTime;

    procedure Pass(SourceIndex : integer); override;

    function GetProcessed: Boolean;

  protected

    function GetCanTake : Boolean; override;

    function GetCanDrop : Boolean; override;

    property EndWorkTime : TCustTime read FEndWorkTime write FEndWorkTime;

  public

    constructor Create; override;

    destructor Destroy;override;

    procedure DropParcel; override;

    property Generator : TGenerator read FGenerator;

    property Processed : Boolean read GetProcessed;

    procedure Work; virtual;

  end;

  TChannel = class(TShop)

  public

    procedure Pass(SourceIndex : integer); override;

  end;

  TSource = class(TShop)

  private

    procedure TakeParcel(SourceIndex: integer);override;

  public

    procedure Pass(SourceIndex : integer); override;

    procedure AskForParcel; override;

  end;

  TAccumulator = class(TElement)

  private

    FParcels : TList;

    FLimited : Boolean;

    FCapacity : integer;

    function GetParcel(Index : integer): TParcel;

    function GetFreeSpacePresent : Boolean;

    function GetEmpty: Boolean;

    procedure SetCapacity(Value : integer);

    function GetCapacity : integer;

    function GetParcelCount : integer;

    procedure Pass(SourceIndex : integer); override;

    function GetCanTake : Boolean; override;

    function GetCanDrop : Boolean; override;

  public

    constructor Create; override;

    destructor Destroy; override;

    procedure ClearContainer; override;

    procedure DropParcel; override;

    property ParcelCount : integer read GetParcelCount;

    property Parcels[Index : integer]: TParcel read GetParcel;

    property FreeSpacePresent: Boolean read GetFreeSpacePresent;

    property Empty : Boolean read GetEmpty;

    procedure TakeParcel(Index : integer); override;

  published

    property Capacity : integer read GetCapacity write SetCapacity;

    property Limited : Boolean read FLimited write FLimited;

  end;

  TAccumulatorClass = class of TAccumulator;

implementation

uses QSheme;

constructor TElement.Create;

begin

  FSources := TList.Create;

end;

destructor TElement.Destroy;

begin

  FSources.Free;

  inherited;

end;

procedure TElement.SetSheme(ASheme : TObject);

begin

  if Assigned(ASheme) then FSheme := ASheme;

end;

procedure TElement.AddSource(Element : TElement);

begin

  if Assigned(Element) then FSources.Add(Element);

end;

procedure TElement.DelSource(Element: TELement);

begin

  if Assigned(Element) then FSources.Remove(Element);

end;

function TElement.GetSourceCount: integer;

begin

  Result := FSources.Count;

end;

function TElement.GetSource(Index: integer): TElement;

begin

  Result := FSources[Index];

end;

procedure TElement.TakeParcel(SourceIndex : integer);

begin

  FContainer := Sources[SourceIndex].FContainer;

  TQSheme(Sheme).NewEvent(EV_TAKE,Self,Sources[SourceIndex],FContainer.Info);

  Sources[SourceIndex].DropParcel;

end;

procedure TElement.Pass(SourceIndex: integer);

var Source : TElement;

begin

  if SourceIndex <> -1 then Source := Sources[SourceIndex];

  DoBeforeTake(Self);

  if SourceIndex <> -1 then Source.DoBeforeDrop(Source);

  TakeParcel(SourceIndex);

  DoAfterTake(Self);

  if SourceIndex <> -1 then Source.DoAfterDrop(Source);

end;

function TElement.GetCanDropParcelFor(Destination: TElement): Boolean;

begin

  Result := CanDrop;

  if Assigned(OnDestinationValidate) then

    Result := Result and OnDestinationValidate(Destination)

end;

function TElement.GetCanTakeParcelFrom(Source : TElement) : Boolean;

begin

  if Assigned(OnSourceValidate) then

    Result := OnSourceValidate(Source)

  else Result := True;

end;

procedure TElement.AskForParcel;

var i : integer;

    Source : TElement;

begin

  for i := 0 to SourceCount - 1 do begin

    Source := Sources[i];

    if Source.CanDropParcelFor[Self] and CanTakeParcelFrom[Source]  then

      if CanTake then begin

        Pass(i);

        if Self is TShop then Exit;

      end

      else

        if not (Source is TAccumulator) then RefuseParcel(i);

  end;//for

end;

function TElement.GetParcelPresent: Boolean;

begin

  Result := FContainer <> nil;

end;

procedure TElement.ClearContainer;

begin

  DropParcel;

end;

procedure TElement.RefuseParcel(SourceIndex: integer);

begin

  Sources[SourceIndex].Container.State := psRefused;

  TQSheme(Sheme).NewEvent(EV_REFUSE,Self,Sources[SourceIndex],Sources[SourceIndex].Container.Info);

  Sources[SourceIndex].DropParcel;

end;

procedure TElement.DropParcel;

begin

  Container := nil;

end;

procedure TElement.DoBeforeDrop(Sender : TElement);

begin

  if Assigned(FBeforeDrop) then FBeforeDrop(Sender);

end;

procedure TElement.DoAfterDrop(Sender : TElement);

begin

  if Assigned(FAfterDrop) then FAfterDrop(Sender);

end;

procedure TElement.DoBeforeTake(Sender : TElement);

begin

  if Assigned(FBeforeTake) then FBeforeTake(Sender);

end;

procedure TElement.DoAfterTake(Sender : TElement);

begin

  if Assigned(FAfterTake) then FAfterTake(Sender);

end;

constructor TGenerator.Create;

begin

  inherited;

  FRandomType := rtPlane;

end;

function TGenerator.GetRandom: TCustTime;

var R : single;

begin

  case FRandomType of

    rtPlane: R := PlaneRND;

    rtNormal: R := NormRND;

    rtExponent: R := ExpRND

  else

    R := Random;

  end;

  Result := FMean - FDisp + Round(R * 2 * FDisp);

end;

constructor TShop.Create;

begin

  inherited;

  FGenerator := TGenerator.Create;

end;

destructor TShop.Destroy;

begin

  FGenerator.Free;

  inherited;

end;

procedure TShop.DropParcel;

begin

  inherited;

  FEndWorkTime := 0;

end;

procedure TShop.Pass(SourceIndex : integer);

begin

  inherited;

  Work;

end;

function TShop.GetProcessed: Boolean;

begin

  Result := (TQSheme(Sheme).SysTime >= FEndWorkTime);

end;

function TShop.GetCanTake: Boolean;

begin

  Result := not ParcelPresent and Processed;

end;

function TShop.GetCanDrop: Boolean;

begin

  Result := ParcelPresent and Processed;

end;

procedure TShop.Work;

begin

  FEndWorkTime := TQSheme(Sheme).SysTime + FGenerator.GetRandom;

end;

procedure TChannel.Pass(SourceIndex: integer);

begin

  inherited;

  Container.State := psWork;

end;

procedure TSource.TakeParcel(SourceIndex: integer);

begin

  Container := TQSheme(Sheme).NewParcel;

end;

procedure TSource.Pass(SourceIndex : integer);

begin

  inherited;

  Container.State := psBorn;

end;

procedure TSource.AskForParcel;

begin

  if CanTake then Pass(-1);

end;

constructor TAccumulator.Create;

begin

  FLimited := False;

  FParcels := TList.Create;

  inherited;

end;

destructor TAccumulator.Destroy;

begin

  FParcels.Free;

end;

function TAccumulator.GetParcel(Index : integer): TParcel;

begin

  Result := FParcels[Index];

end;

function TAccumulator.GetCanDrop: Boolean;

begin

  if Empty then AskForParcel;

  if not Empty then Container := FParcels.First;

  Result := not Empty;

end;

function TAccumulator.GetCanTake: Boolean;

begin

  Result := FreeSpacePresent;

end;

function TAccumulator.GetFreeSpacePresent: Boolean;

begin

  Result := (Capacity <> FParcels.Count) or (not Limited);

end;

function TAccumulator.GetEmpty: Boolean;

begin

  Result := FParcels.Count = 0;

  //if not Result then Container := FParcels.First;

end;

procedure TAccumulator.DropParcel;

begin

  if not Empty then FParcels.Delete(0);

  inherited;

end;

function TAccumulator.GetCapacity : integer;

begin

  Result := FCapacity;

end;

function TAccumulator.GetParcelCount: integer;

begin

  Result := FParcels.Count;

end;

procedure TAccumulator.SetCapacity(Value : integer);

begin

  FLimited := True;

  FCapacity := Value;

end;

procedure TAccumulator.ClearContainer;

begin

  FParcels.Clear;

  inherited;

end;

procedure TAccumulator.Pass(SourceIndex : integer);

begin

  inherited;

  Container.State := psStore;

end;

procedure TAccumulator.TakeParcel(Index : integer);

begin

  FParcels.Add(Sources[Index].Container);

  TQSheme(Sheme).NewEvent(EV_TAKE,Self,Sources[Index],Sources[Index].Container.Info);

  Container := FParcels.Last;

  Sources[Index].DropParcel;

end;

end.