Моделирование систем
Министерство общего и профессионального образования
Северо-Кавказский
Государственный Технический Университет
Факультет
информационных систем и технологий
Кафедра
автоматизированных систем обработки информации и управления
ПОЯСНИТЕЛЬНАЯ ЗАПИСКА
К курсовой работе по дисциплине
«Моделирование систем»
На
тему_______________________________________________________
(наименование темы в
соответствии с заданием)
Студент группы________________
_______________________________
______________________
(подпись)
(Ф. И. О.)
Руководитель___________________
______________________
(подпись)
(Ф. И. О.)
Оценка__________________________
2000 г.
Министерство
общего и профессионального образования
Северо-Кавказский
Государственный Технический Университет
Факультет
информационных систем и технологий
Кафедра
автоматизированных систем обработки информации и управления
«УТВЕРЖДАЮ»
заведующий кафедрой
_____________________
«_____»_________2000 г.
ЗАДАНИЕ
на курсовую
работу по дисциплине «Моделирование систем»
Студент____________________________группа_____________________
(Ф. И. О.)
Руководитель__________________________________________________
(ученое звание,
степень, Ф. И. О.)
Тема курсовой
работы__________________________________________
Задание выдано «_____»_____________2000 г.
Руководитель___________________
(подпись)
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++ и др.).
Учитывая условие, построим структурную схему
данной системы.
Рис. 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;
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;
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.