Моделирование систем

Министерство общего и профессионального образования
Северо-Кавказский Государственный Технический Университет
Факультет информационных систем и технологий
Кафедра автоматизированных систем обработки информации и управления

ПОЯСНИТЕЛЬНАЯ ЗАПИСКА
К курсовой работе по дисциплине «Моделирование систем»
На тему_______________________________________________________
(наименование темы в соответствии с заданием)
Студент группы________________
_______________________________ ______________________
(подпись) (Ф. И. О.)
Руководитель___________________ ______________________
(подпись) (Ф. И. О.)
Оценка__________________________
2000 г.
Аннотация.
В данном курсовом проекте рассматривается задача о машинном моделировании процесса взаимодействия между тремя проектировщиками и центральной ЭВМ. Проводится анализ этой системы и формализация ее в терминах Q-схем, а также приведена программа моделирующего алгоритма и математическое описание системы. Проводится сравнительный анализ аналитического и имитационных методов исследования.
Содержание.
Введение………………………………………………………………………………………………………………. 5
Основная часть……………………………………………………………………………………………………… 6
Постановка задачи…………………………………………………………………………………………….. 6
Возможные пути исследования…………………………………………………………………………. 6
Этап моделирования…………………………………………………………………………………………. 7
Разработка Q-схемы системы………………………………………………………………………… 7
Разработка моделирующего алгоритма и машинная реализация……………………. 10
Математическое описание системы…………………………………………………………………. 18
Результаты моделирования и аналитического решения…………………………………….. 19
Заключение…………………………………………………………………………………………………………. 20
Литература………………………………………………………………………………………………………….. 21
Приложения………………………………………………………………………………………………………… 22
Текст программы…………………………………………………………………………………………….. 22
Введение.
В настоящей курсовой работе рассматривается проблема моделирования процессов в Q-схемах – одном из важнейших, с точки зрения применения на практике, классов математических схем, разработанных для формализации процессов функционирования систем массового обслуживания (СМО) в теории массового обслуживания. Предметом изучения в теории массового обслуживания являются системы, в которых появление заявок (требований) на обслуживание и завершение обслуживания происходит в случайные моменты времени, т.е. характер их функционирования носит стохастический характер. Следует отметить, что СМО описывают различные по своей физической природе процессы функционирования экономических, производственных, технических и других систем, например потоки поставок продукции некоторому предприятию, потоки деталей и комплектующих изделий на сборочном конвейере цеха, заявки на обработку информации в ЭВМ от удаленных терминалов и т.д.
Основная часть. Постановка задачи.
САПР состоит из ЭВМ и трех терминалов. Каждый из проектировщиков формирует задания на расчет в интерактивном режиме. Набор строки задания занимает 10±5с. После набора 10 строк задание считается сформированным и поступает на рещение, при этом в течение 10±3с ЭВМ прекращает принимать другие задания. Анализ результата занимает у проектировщиков 30 с, после чего цикл повторяется. Данные по всем проектировщикам одинаковы.
Смоделировать работу системы в течение 6 часов. Определить вероятность простоя проектировщика из-за занятости ЭВМ и коэффициент загрузки ЭВМ. Возможные пути исследования.
Для изучения Q-схем используются два подхода аналитический и имитационный. При аналитическом подходе подлежащая анализу схема описывается с помощью формул, отражающих зависимости между ее различными параметрами. Однако, следует отметить, что разработанные методы аналитического изучения Q-схем подходят далеко не для каждой конкретной системы, они пригодны лишь для систем общего типа. Поэтому при аналитическом изучении систем их необхродимо упрощать до систем основных типов, что в последствии конечно-же сказывается на результатах исследования. При имитационном подходе ставится эксперимент на машинной модели системы, которая предварительно реализуется на одном из созданных специально для этого языков имитационного моделирования (например, SIMULA, SIMSCRIPT, GPSS и др.) или на языке общего назначения (BASIC, PASCAL, FORTRAN, C++ и др.). Этап моделирования. Разработка Q-схемы системы.
Учитывая условие, построим структурную схему данной системы.

Рис. 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.