CC BY
12МОДЕЛИ ИМИТАЦИИ В НЕАРХИМЕДОВОМ ВРЕМЕНИ: ВРЕМЯ, СИСТЕМНЫЕ ДИНАМИКИ
М. И. Нечепуренко*
УДК 510.23:519.681
Предлагается специальная структура ньютонова неархимедова модельного времени, позволяющая рассматривать имитационные модели в рамках общей теории систем (в терминах временных функций) и являющаяся естественной формой синхронизации процессов в структурированных имитационных моделях.
Ключевые слова: имитационные модели, синхронизация процессов.
The special structure of Newton's non-Archimedean model time is proposed. This structure allows consider simulation models in the scope of general theory of systems (in terms of time functions) and is a natural form of synchronization of processes in structural simulation models.
Key words: simulation models, synchronization of processes.
Введение. Широкое применение имитационного моделирования (simulation) к исследованию сложных технических, биологических, экономических и социальных систем привело к развитию теории имитации и, в частности, к формулировке ряда классов систем [1-7] со специфической структурой и сингулярной динамикой поведения. Под сингулярностью динамики здесь понимается патологическая возможность появления в одной реализации эволюции системы последовательности событий, относящихся к одному и тому же моменту модельного времени, но упорядоченных по выполнению (можно было бы говорить о неординарности динамики).
Теория имитации [6, 7] развивается независимо от общей теории систем [8], так как системы с сингулярной динамикой не являются временными (в вещественном времени) и их поведение обычно описывается поеобытийно с помощью конструкций, адекватных механизму упорядочения в программах имитации.
Один из возможных путей сближения теории имитации и общей теории динамических систем состоит в обобщении понятия временного объекта, в частности, в рассмотрении последовательностей вида
(ao,to), (ai,ti), (a2 ,t2),... (1)
как кусочно-постоянных „обобщенных функций" времени. В (1) tj — моменты модельного времени, в которые происходят события ai-1 ^ aj изменения состояний объекта; t0 ^ t1 ^ t2 ^ ... Не более чем счетпость числа элементов поеледовательноети (1) обеспечивает ее имитируемоеть — возможность рекуррентного воспроизведения любого ее начального
(1)
может быть записана в виде слова
* Статья была опубликована в сб. науч. трудов „Эффективность и структурная надежность информационных систем" (Системное моделирование — 7). Новосибирск, 1982. С. 43458.
(¿оК°а?1а?2 ..
(2)
где = ¿¿+1 — интерпретируемого следующим образом,
В начальный момент ¿0 объект находится в состоянии а0. В момент ¿1 = ¿0 + а0 состояние объекта изменилось на а1 и т. д. Ограничившись рассмотрением стационарных объектов, можно в (2) опустить элемент (¿0). Для последовательностей такого типа естественно определяется понятие отрезка последовательности, а также операции: сочленения двух последовательностей, первая из которых конечна; сокращения последовательностей, например, а?ат = а?+т при а,т > 0, но а?ат = а?; внешнего объединения двух последовательностей, например,
Сокращенное сочленение отрезков сокращенной последовательности восстанавливает эту последовательность. Сокращенные последовательности совпадают с соответствующими сокращенными проекциями их внешнего объединения, В терминах последовательно-(2)
стямп в динамике, в частности, кусочно-линейные агрегаты [1, 2],
В данной работе развивается другой подход, предложенный в сообщении [9]: определяется специальная структура модельного времени, позволяющая, в частности, давать описания имитационных моделей в терминах временных функций,
1. Имитационные программы и имитируемые системы.
1,1, Общая модем,ь имитации.
1,1,1, Имитацией называют воспроизведение на цифровых вычислительных машинах отдельных реализаций эволюции в модельном времени исследуемой замкнутой динамической системы (включая среду) и вычисление функционалов, определенных на этих реализациях,
В дальнейшем мы ограничимся рассмотрением детерминированных стационарных динамических систем, эволюционирующих на конечных или бесконечных справа интервалах модельного времени. Обобщение на стохастические системы может быть проведено обычным образом,
С процессом имитации связаны две системы: имитируемая система (точнее, ее модельное представление в имитационной программе) и имитационная программа в целом, рассматриваемая как динамическая система. Обе системы замкнуты (автономны). Имитационная программа рекуррентно воспроизводит фазовую траекторию имитируемой системы (изменение во времени ее состояния). При этом за одну реализацию процесса имитации вычисляются значения состояния имитируемой системы лишь для конечного числа моментов модельного времени. Поведение имитируемой системы между этими моментами фактически не рассчитывается имитационной программой, хотя оно должно быть однозначно определено и заложено в алгоритмы имитации таким образом, чтобы подобный расчет мог быть выполнен для произвольного промежуточного момента. Кроме собственно воспроизведения фазовой траектории имитируемой системы, имитационная программа выполняет ряд вспомогательных функций, связанных непосредственно с организацией имитации (продвижение модельного времени, проверка условий окончания имитации и
др.), а также е документированием эволюции имитируемой системы. Текущее модельное время и состояние системы являются проекциями вектора состояния программы,
1.1.2, В детерминированном стационарном случае имитационная программа может быть описана шестеркой I = (P, C, T, f, g, h), где P х C х T — множество состояний программы, C — множество состояний имитируемой системы, T — множество моментов модельного времени, функции f : P х C м Рид : C м C определяют переходы (события) в программе, функция h : C м T продвигает модельное время. Обычно предполагается вещественность и неотрицательность модельного времени, то есть T = R+, Здесь и далее R+ = [0; то].
Для заданного начального состояния а0 = (p0, с0, t0) G Pх Cх T программа I при своем выполнении порождает последовательность (1), в которой а = (pi, с, t) и Pi+i = f (pi, с), C+i = g(ci) ti+1 = ti + h(ci) для всех i G N = (0; 1; 2;...). Переход ai м ai+1 называется i-м событием в реализации программы G с начальным состоянием а0. Величина ti+1 = Рr3(ai+1) интерпретируется как момент модельного времени i-ro события в программе, Значение ti для каждого i G N обычно называется текущим модельным временем, соотносимым функцией
time = рr3 : P х C х T м T, ai м time(ai) каждому состоянию ai программы. Значение h(ci) определяет текущий шаг имитации:
time(ai+1) = time(ai) + h^r2(ai)).
1.1.3. Последовательность (1), порождаемая программой I = (P, C, T, f, g, h), определяет пособытийную эволюцию имитируемой системы I = (C, T, f, g,5), где семейство функций
5 = : C м C|t G T}
задает динамику системы между событиями в программе. Семейство 5 удовлетворяет для всех с G C и а, т G T условиям 50A = A, 5СТ (5ТA) = 5CT+r A,
Событие а м ai+1 в программе порождает событие Ai м Ai+1 в системе, соотносимое моменту ti+1 модельного времени; состояние ci отнесено моменту t^. Если ti < ti+1, то для каждого t G [ti, ti+1] системе приписывается состояние 5t-tici.
5h
влетворять дополнительному условию
h(c) = h(5CTс) + а для 0 ^ а ^ h(c).
В данную модель имитации, относящуюся к так называемым системам с дискретными событиями, вкладываются более структурированные модели, в частности, кусочно-линейные агрегаты [2], концептуальные модели языков моделирования [4, 5] и математические модели теории имитации [6, 7].
1.2. Временные сингулярности в имитации. Имеются две основные причины появления в имитационных моделях событий, разделенных временными интервалами нулевой длины: квазипараллелизм имитации событий и допустимость нулевых (мгновенных) реакций в компонентах моделируемых систем [7, 9].
1.2.1, Первая из этих причин — чисто технологическая и присуща только имитационному моделированию. Необходимость введения механизма квазипараллельного выполнения событий возникает в случае, когда число процессоров в вычислительной системе, на которой производится имитация, меньше максимально возможного числа одновременных событий, назначенных на один и тот же момент модельного времени. Обычно семантика языков автоматизации программирования имитационных моделей ориентирована на однопроцессорное выполнение программы имитации. Исключение составляет лишь язык моделирования ДИС [5, 10,11], допускающий (с точки зрения наблюдателя) параллельную имитацию одновременных событий,
1.2.2, Вторая причина выходит за рамки имитационного моделирования и присуща математическому моделированию в широком смысле. Нулевые времена реакции систем и временные интервалы нулевой длины появляются при описании явлений, разномасштабных по протяженности во времени. Интервалы нулевой длины являются не менее естественной абстракцией, чем разрывные временные функции. Существенно также, что при математическом моделировании на ЦВМ с использованием вещественного модельного времени разномасштабные интервалы в одной модели практически нереализуемы: если значения текущего модельного времени £ и интервала Д£ очередного продвижения модельного времени отличаются на несколько десятков десятичных порядков, то, хотя они и могут быть одновременно представлены в памяти машины, тем не менее, продвижение времени не будет выполнено, так как результат операции £ + Д£ из-за округления будет равен Дело не спасает использование вычислений с многократной точностью, если в продвижении модельного времени будут участвовать интервалы, значительно различающиеся по порядку,
1,3, Об адекватности модельного времени,. Вопрос о целесообразности (с точки зрения адекватности моделирования) присутствия разномасштабных величин в одной математической модели выходит за рамки данной работы, примыкающей к существенно более узкой проблематике — „миру имитации" как таковому, вне связи его с реальным „имитируемым миром".
Тем не менее, представляется уместным привести здесь следующие высказывания, относящиеся к общей проблеме формализации понятия времени,
1,3,1, Разномасштабные объекты и дискретные модели, „,,,Какого бы взгляда ни придерживаться относительно строения материи, не подлежит сомнению то, что она расчленена на ряд больших, хорошо отграниченных групп с относительно различными размерами масс, так что члены каждой отдельной группы находятся со стороны своей массы в определенных, конечных отношениях друг к другу, а к членам ближайших к ним групп отно-
"
количества, столь колоссально отличные друг от друга, что между ними прекращается
всякое рациональное отношение, всякое сравнение, и что они становятся количественно
"
„Особенности всевозможного рода, разрывы и т.д. появляются, помимо воли исследо-
"
описапии явлений природы постоянно имеют дело, смотря по обстоятельствам, то с непрерывной средою, то с дискретными моментами. Разумеется, речь идет только о приближениях; но в том-то и дело, что такого рода приближения с одинаковым успехом дают нам
"
„Что можно сказать об еще более коротких интервалах времени? Имеет ли смысл вообще говорить о них, если невозможно не только измерить, но даже разумно судить о
процессах, происходящих в течение столь коротких интервалов? Возможно, нет. Это один
"
1.3.2, Понятие времени. Имеется два различных подхода к определению времени как формы координации материальных явлений. При одном из них время предполагается актуально заданным. Этот подход (И, Ньютона) реализуется в математических моделях, „Время — нечто абсолютное, ни от чего не зависящая чистая длительность, как таковая,
равномерно текущая от прошлого к будущему. Оно является пустым „вместилищем со"
на течение времени. Время универсально, одномерно, непрерывно, бесконечно однородно
"
лейбницевекие подходы — прим. авт.).
При другом подходе время рассматривается лишь в связи с эволюцией наблюдаемых
процессов. Этот подход (Г, Лейбница) по сути дела реализуется в имитационном модели"
„Время... есть созерцаемое становление", „...Не во времени все возникает и проходит, а само время есть это становление, есть возникновение и прехождение, сущее абстрагирование, веепорождающий и уничтожающий свои порождения Кронос ". „Если бы все остановилось.., то... времени не'было бы" (Г. Б. Ф. Гегель, [16]).
"
1.3.3. Архимедовость времени. „...Так как говорить о неархимедовом времени не представляется целесообразным ни с какой точки зрения, то отсюда вытекает, что употребление „переменного" ж, развертывающего перед нами все элементы последовательности и в порядке реального времени (не фиктивного, то есть архимедова), сильно ограничивает класс рассматриваемых последовательностей М(Н,Н, Лузин, [17]).
„Что касается аксиомы (имеется в виду аксиома Архимеда — прим. авт.), то она может
"
величин „не слишком малых" (Н. Бурбаки, [18]).
Наконец отметим, что в общей теории стационарных динамических систем [8] о множестве моментов времени предполагается лишь, что оно является правым интервалом некоторой линейно-упорядоченной абелевой группы (вообще говоря, неархимедовой). 2. Модельное время. 2.1. Цепь модельного времени,.
2.1.1. Далее, К =] —; Ж = (0; ±1; ±2;...). Множество К2 с обычной операцией сложения векторов
(£ъ П1) + (6, П2) = (6 + П1 + П2)
(3)
и лексикографическим линейным отношением порядка
(£ъ П1) ^ (6, П2) ^ (6 < 6) или (^1 = ^ Ц1 ^ П2)
образует линейно-упорядоченную, абелеву группу. (3)
n(£ П) = (n£ nn) (5)
В связи с этим R2 можно также рассматривать как линейно-упорядоченное линейное векторное пространство над кольцом целых чисел Z,
Как обычно, будем использовать обозначения O = (0, 0); — (£, п) = (—— n); ti —t2 = ti + ( — ¿2).
Здесь и в (3)-(5) £, п, П е R; n е Z; t е R2.
Таким образом, для всех m, n е Zh t, t, G R2 выполняются следующие определяющие свойства линейно-упорядоченного линейного векторного пространства:
ti + t2 = t2 + ti; ti + (t2 + ta) = (ti + t2) + ta;
t + O = t; t — t = O;
(mn)t = m(nt); 1 * t = t;
m(ti + t2) = mti + mt2; (m + n)t = mt + nt
или^ = t2, ми ti < t2, или t2 < ti;
ti ^ t2, t2 ^ ta ^ ti ^ ta; ti ^ t2, t2 ^ ti ^ ti = t2; ti < t2 ^ t2 — ti > 0.
Отсюда следует, в частности, что
Ot = O, (—m)t = —mt; t > O ^ —t < O;
ti < t2, ta ^ t4 ^ ti + ta < t2 + t4;
ti < t2, ^ mti < mt2 при m > 0; ti < t2, ^ nt2 < mti при n < 0.
2,1,2, Множество T = R x Z с индуцируемыми на нем из R2 операциями векторного
сложения и умножения и умножения на целые числа является подпространством про-R2
t G T
равенетвом t = (Re t, In t).
Очевидно, что не вызывает недоразумений использование записи £ для обозначения элементов вида (£, 0) е T, В этих обозначениях имеет место представление
t = Re t + In t ■ в,
где символ в резервируется за элементом в = (0, 1) е T. Очевидно также, что не вызывает недоразумений использование единого символа 0 для обозначения нулевых элементов 0 е R и O е T,
TT
T
названа сингулярным временем,
TT
значения:
[а; в] = {t|a ^ t ^ в}, [а; в] = {t|a ^ t < в}, ]а; в] = {t|a < t ^ в}, ]а; в[= {t|a < t < в}.
Очевидно, что
[а; в] = 0 при а < в ; [а; в [=]а; в ] =]а; в[= 0 при а ^ в ; [а; в[= [а; в - 0], ]а; в] = [а + в], ]а; в[= [а + 0; в - 0].
Сечения и склеивания интервалов задаются равенствами
[а; y] = [а; в[и[в; y] при а < в ^ Y,
[а; y] = [а; в[и]в; y] пРи а ^ в < Y.
Для любого множества T G T и t G T через T + t будем обозначать сдвиг T на t, то есть T + t = {s|s - t G T}.
При t > 0 будем говорить о сдвиге множества T вправо, а при t < 0 — влево. Очевидно, что [а; в] + t = [а + t; в + t].
2.1.3, Для дальнейшего целесообразно также рассматривать цепь T = T U {то} модельного времени, где бесконечный момент то есть единица решетки T, то ееть то = sup T, Можно положить то = (+то, 0) и (в соответствии с принятой договоренностью) не различать tog T и +то G R , Элеме нт то обладает для л юбых t G T и m G Z следующими свойствами:
t + то = то + t = то; тто = тот = то при m > 0.
T
Для а G T полагаем
[а; то[= {^а ^ t < то}, [а; то] = {^а ^ t ^ то}.
Такие интервалы будем называть бесконечными (справа).
Далее будут рассматриваться только конечные или бесконечные интервалы цепи T, Только такие интервалы будут временными, T-интервалами, Интервалы вида
] — то; а[= {t|t < а}, ] — то; а] = {t|t ^ а}
как области определения функций T-времени не будут рассматриваться. Это не мешает записать:
T =] — то; +то[, T =] — то; +то[.
Отметим, что при имитации момент t = той бесконечные приращения модельного времени Д t = то реализуются большим (по еще предетавимым в памяти машины) числом, причем выход имитационной программы на t = то влечет окончание имитации.
Далее запись tk ^ то используется в обычном смысле: последовательность tk G T не ограничена,
2.1.4, Множество T Ç T будем называть финитным на интервале I Ç T, если на каждом конечном подынтервале J Ç T имеется лишь конечное число элементов множества T П J < +то.
Очевидно, что множество финитно на конечном интервале тогда и только тогда, когда оно конечно.
2,2, Характеристика модельного времени.
2.2.1, Введенная структура модельного Т-временн обладает следующими свойствами. Время актуально задано (ньютоново). Каждое событие в модели имитации привязывается к определенному Т-моменту, События в модели, назначенные на один и тот же Т-момент, с точки зрения наблюдателя одновременны. Поэтому в моделях имитации для большинства систем автоматизации программирования имитационных моделей на каждый Т-момент приходится те более одного события, В ДИС-моделях [5] Т-время обеспечивает синхронизацию параллельных процессов.
Время всюду дискретно, то есть для любого Т-момента t существуют покрывающий его элемент (момент t + в, последующий за t) и элемент, который t покрывает (момент t — в, предшествующий t).
Время всюду однородно в том смысле, что для любых T С Т, a G Т отображение T ^ T + a t ^ t + a является изоморфизмом. Однородность времени обеспечивает возможность имитации когерентных процессов, В [9], кроме данной структуры модельного времени, рассмотрена также иная, не обладающая свойством однородности.
Время неархимедово, так как, например, тв < 1 для любых m G Z; то есть в является слабой единицей [19], Генератор вида t ^ t + в те порождает всю времениую цепь Т и поэтому не может служить целям измерения временных интервалов (то есть не является часами),
2.2.2, Вообще говоря, ничего не мешает рассматривать более тонкие структуры модельного времени, обеспечивающие большие описательные возможности синхронизации параллельных когерентных процессов.
Например, в качестве цепи модельного времени можно взять пространство R2 или даже множество числовых последовательностей t = (£i, ,...), G R с операциями
i
+ ы = &+п), = Е j)
j=i
и отождествлением & = (£, 0, 0,...),
В данной работе эти варианты не рассматриваются,
3. Временные объекты и системные динамики. Далее используются стандартные терминология и обозначения категорий (см., например, [8, 20]), В частности, K(A, B) — множество морфизмов из объекта A в объект B категории K, idA — единичный морфизм в K(A, A), ObK — класс K-объектов, Set — категория множеств с морфизмами-отображениями. Ряд определений и результатов в данном разделе формируется лишь для Set
пых алгебраических систем [21] с морфизмами — гомоморфизмами систем.
При описании временных объектов и системных динамик не обойтись понятиями и терминологией, сложившимися в теории систем [8, 21], Например, из-за неархимедовости цепи модельного времени вместо морфизмов, определяющих системные динамики, пришлось рассматривать однопараметрические полугруппы морфизмов.
Далее вместо временных функций [8] (то есть функций, определенных на подмножествах множества Т) в основном используются конструкции, аналогичные счетным степеням объектов [20], Такой подход естественен при рассмотрении стационарных динамических систем.
Так как в дальнейшем рассматривается только положительная полуось модельного времени, то для упрощения обозначений положено:
Рис. 1.
Рис. 2.
т = т+, т = Т+.
3.1. Степени объектов категорий.
3.1.1 Определение. Будем говорить, что объект А категории К является объектом с бесконечной степенью, если выполняются следующие условия:
1) для каждого Ь € Т, Ь = 0 определены К -объект А4, называемый ¿-й степенью объекта АА
{П : А4 ^ А|0 ^ 5 < ¿},
что для любого семейства К-морфизмов {fs : В ^ А|в < ¿} с общей областью определения В существует единственный К-морфизм f : В ^ А+, для которого коммутативны все диаграммы (рис. 1);
2) определена ассоциативная операция умножений (сочленения) степеней: As•А4 = А^4 для всех 0 < 5 < то, 0 < £ ^ то, причем коммутативны все диаграммы (рис. 2).
Последние диаграммы согласуют семейство морфизмов п4 с порядком записи сомножителей в произведении степеней объектов, а также, очевидно, однозначно определяют морфизмы по морфизмам Далее верхний индекс в обозначениях морфизмов П
будет опускаться, так как их область определения всегда будет ясна из контекста. А4 А
А4
физма, в частности,
Ав ~ А.
Пример. В имеет место представление = {(а^ет а € А для всех 5 € Т} с покоординатными операциями, индуцируемыми из А Морфизм п4 является ¿-й проекцией, то есть
п4(^ет = а4.
Имеют место изоморфизмы
~ А[4;те[ при люб ом Ь € Т
Рис. 3.
$ е Т
1 ос
А А
Рис. 7.
Рис. 4.
С ч" Г~1
А-— А \А°
Рис. 6.
в.
ос
ос
ос
А
П
В —А
Рис. 8.
А2 - А[0;21 ноА2 ф А х А ф А20.
3.1.2. Определение. Пусть существует К-объект для Для 0 < Ь ^ то левым сужением до ¿объекта будем называть единст венный К-морфиз м г4 : ^ А4, для которого коммутативны все диаграммы (рис. 3).
Морфизм г8 продолжается на каждый объект А4 для 5 < Ь ^ то как единственный К-морфизм г8 : А4 ^ А8, для которого коммутативна диаграмма (рис. 4).
3.1.3. Определение. Пусть существует К-объект Для Ь € Т левым сдвигом на Ь объекта будем называть единст венный К-морфиз м : ^ для которого коммутативны все диаграммы (рис. 5).
Отметим коммутативность диаграмм (рис. 6).
Морфизм г8 продолжается на каждый объект А4 для 0 ^ 5 < Ь как единственный К-морфизм г8 : А4 ^ А4-8, для которого коммутативна диаграмма (рис. 7).
Лемма. г0 = ; г80г4 = г8+4 для всех Ь € Т.
3.1.4. Определения и свойства морфизмов сужений, левых сдвигов и умножений степеней приводят к следующим аналогиям между морфизмами в общих категориях и преобразованиями временных функций:
о сужение функций с[0; то на то[; о сужение функций то на [0; Ь[; • о сочленение двух функций с точкой сочленения Ь.
3.1.5. Лемма. Пусть существуют К-объекты и В
Для любого морфизма F : B ^ A существует единственный морфизм F^ : B^ ^ A^, для которого коммутативны все диаграммы (рис, 8),
Доказательство непосредственно следует из определения 3,1,1, если выбрать ft = F•п : B~ ^ A,
3,2, Потоки. В данном подразделе без оговорок рассматривается категория Set.
3.2.1, Определение. Пусть A — множество, |A| ^ 2, и выделен элемент Oa G A, Элемент x G будем называть финитным (A, Оа)-потоком, если множество {t|ntx = Oa} финитно на T,
Каждую пару (ntx, t), такую, что t G T и ntx = Oa, будем называть событием в потоке x в момент t,
x
пом времени: два его события могут относиться к моментам s = t, таким, что Res = Ret.
3.2.2, Определение. Пусть A и B — множества, a x и y — финитные (A, Oa)- и (B, Ob)-потоки соответственно. Внешним объединением потоков x и y будем называть финитный (A х B, (Oa, Ob))-поток x * y, определяемый равенствами
Отметим, что события в потоке х * у могут иметь вид ((а, О в), ¿) и ((Оа , Ь), где а = Оа и Ь = Ов.
3,2,3, Определение. Пусть х и у — финитн ые (2, 0)-потоки, где 2 = {0; 1} (булевы потоки), Слиянием потоков х и у будем называть финитный (2, 0)-поток х Л у, определяемый для Ь € Т равенствами
3,2,4, Определение. Пусть А и В — множества. Отображениие : ^ Вбудем называть морфизмом ^^^^^^^^^щирования, а (В, Ов)-поток у = х ^-дифференциалом элемента х € (функции х : Т ^ А), если выполняются условия:
1) еюръекция F : Ах А ^ В такова, что F-1({Ов}) = Аа, где Аа С Ах А — диагональ;
2) морфизм определен равенствами п4 о х = F (п4х, п4+1х) для всех х € и
Если х : Т ^ А — кусочно-постоянная временная функция с не более чем счетным числом скачков и моменты скачков образуют финитное множество на Т, то условие 1 гарантирует финитноеть потока х , В этом смысле можно говорить, что F-дифференциал отфильтровывает скачки, В общем случае F-дифференциал не является финитным потоком,
3.2.5, Для дальнейшего важен следующий частный случай дифференциала.
Определение. Для каждого множества А существует единственная сюрьекция F :
А х А ^ 2 со свойств ом F -1({О}) = Да, то есть диффере нциал определяется единственным образом, В этом случае для элемента х € будем писать ¿х вместо х и называть морфизм й индикатором скачков,
3.2.6, Теорема. Для любых х € , у € и морфизмов дифференцирования : ^ С и : ^ С выполняются равенства:
1) о г4 = г4 о , Ь € Т;
2) (г4х • г4у) = (г4йрх) • (г4йрх), Ь € Т;
nt(x * y) = (ntx, nty), t G T.
nt(x Л y) =
1, если ntx = ntx = 1; 0, в противном случае.
tT
Рис. 9.
Рис. 10.
3) dFxg(x * y) = dFx * dGy, t G T.
3.3. Автономные динамики.
3.3.1. Определение. Автономной динамикой в категории К будем называть пару (A, S), где A — объект вК, а однопараметрическая полугруппа морфизмов
S = (S4 : A ^ A|t G T}
удовлетворяет условиям:
1) Ss о S4 = Ss+i для всех s G T, t G T;
2) So = idA.
Последнее означает, что So является не только единицей полугруппы S, как это следует из условия 1, но также и единичным морфизмом в K(A, A).
Следующий пример показывает, что из-за неархимедовости T-времени морфизм Sö, вообще говоря, не определяет однозначно динамику (A, S),
Пример. К = Set, A = R, S\a = a + Ret, St,a = a + 2 Ret. Здесь S1 = SÎ, при Ret > 0, хотя = Sg = 1A.
Так как S4 = SRe 4 о (Sö)In то каждый морфизм S4 допускает факторизацию на мор-физмы, описывающие движение в вещественном времени (t = Re t) и микровремени (t = в ■ In t). Тем самым справедлива следующая
Теорема. Произвольная автономная динамика (A, S) допускает факторизацию (представление) (A,SRe,SIn), где
SRe = (SRe : A ^ A|t G R+}, Sm = (SRe : A ^ A|m G N},
так, что S4 = SRe4 о Sjn4 для всех t G T; SRe = S4, S^" = Smd для всех t G R+, m G N.
3.3.2. Теорема. Пусть (A, S) — автономная динамика в категории К, и существует Aœ, Тогда существует единственный К-морфизм а : A ^ Aœ, такой, что коммутативны все диаграммы (рис. 9).
Не вызывает недоразумений использовать символ S вместо а для обозначения такого морфизма.
Пусть (A, S) — автономная динамика в Set. Образ x элемента a G ab кообласти Aœ морфизма S будем называть автономной траекторией в A, Множество автономных траекторий, определяемых динамикой (A, S) для различных начальных состояний a G A, изоморфно множеству ((Sia)teT|a G A}, а значит, и множеству временных функций (xa : m ^ A|a G A}, определяемых равенствами xa(t) = S4a для всех t G T. Каждую из этих реализаций также можно называть автопомпой траекторией.
3.3.3. Определение. Обозначим через DynK класс автономных динамик в категории К. Динаморфизмом g : (A, S) ^ (B, е) будем называть каждый морфизм, индуцируемый на Dyn-К теми К-морфизмами g : A ^ B, для которых коммутативны все диаграммы (рис. 10).
Теорема. DynK является категорией с объектами — автономными динамиками в категории Кис морфизмами — динаморфизмами,
3.3.4, Определение. Автономную динамику (A, 6) в категории Set будем называть кусочно-непрерывной и записывать четверкой (A, f, G, h), если функция h : А ^ Т, однопараметрическое семейство f = {f4|t е Т} морфизмов
f : {а ^ A|h(a) ^ t} ^ А (6)
и морфизм
G : {а е A|h(a) = 0} ^ А (7)
удовлетворяют для всех а е А и s,t еТ следующим условиям:
1) h(a) = то ^ f4a = a;
2) f°a = a;
3) s + t ^ h(a) ^ fs ■ f4a = f s+ta;
4) t ^ h(a) ^ h(a) = h(f4a) +1;
5) h(a)0 ^ h(Ga) > 0, 66a = f6Ga;
6) h(a) > 0 ^ 64a для t < h(a);
7) h(an) = то гДе ao = a an+1 = G о f h(a,n)ara дая n ^ 0,
ra^ü
Замечание 1, Определения (6) и (7) предполагают „измеримость" функции h в категории К в следующем слабом смысле:
h-1 {(0)}, h-1 ([t; то]) е ObK для всех t е T.
Замечание 2, Условие 7 является необходимым и достаточным для рекуррентного определения 6 по f и G, По сути дела, это означает возможность поеобытийной имитации на ЭВМ за конечное время сколь угодно длинных отрезков кусочно-непрерывных автономных динамик.
Пример. Опишем начальный отрезок автономной траектории, определяемой кусочно-непрерывной динамикой (A, f, G, h) с начальным состоянием а е А, таким, что h(a0) > 0, h(a1) = 0 и h(a2) > 0 гДе an определены по (8). Тогда
64a0 = f4a0, h(64a0) = h(a0) — t при 0 ^ t ^ h(a0);
64a0 = f6a1, h(64a0) = 0 при t = h(a0) + 0;
¿4 = f, h(64) = h(a2) — t при 0 ^ t — h(a0) — 0 ^ h(a2).
3.3.5. Определение. Прямым произведением (A,6) x (B, e) автономных динамик (A,6) и (B, e) в категории Set будем называть единственную автономную динамику (A x В), (6 x e), для которой коммутативны все диаграммы (рис. 11).
Лемма. Если D1 = (A1, f1, G1, h1) D2 = (A2, f2, G2, h2) — кусочно-непрерывные автономные динамики в категории Set, то автономная динамика D1 x D2 также кусочно-непрерывна, причем D1 x D2 = (A, f, G, h), где A = A1 x A2, функция h : A ^ T определяется равенствами
5 X £
я
г
Ах В
Ах В
Ах в ^Т
1
Рис. 11.
Л,(а, Ь) = ш1п(^1(а), Л.2(а)) для а € А, Ь € В;
а морфизмы /¿и С — отображениями
/(а, Ь) = (/1а, /2Ь) для 0 ^ £ ^ %, Ь);
Заключение. Введенный в предыдущих разделах формализм позволяет вложить исследование систем с сингулярной динамикой в общую теорию систем |8|, |22|,
Остановимся кратко па последователыюетпом подходе к определению и исследованию стационарных динамических систем, эволюционирующих в Т-времени.
Такие системы могут быть описаны как специальный класс сетей синхронных автоматов. Описание однозначно характеризует отношение входных и выходных функций. Неархимедовость времени исключает возможность говорить при этом о рекуррентной эволюции системы с шагом 9, Подобная эволюция со счетным числом шагов не обеспечивает потенциальную наблюдаемость системы па интервалах ненулевой вещественной
Т
ментами £ + п9, где значение £ € Т фиксирован о, ап € 2 (цепь модельного времени Т не является совершенно упорядоченной).
В связи с этим ноеледователыюетный подход, потенциально реализующий эволюцию Т
телыюсти моментов. В частности, имитируемость предполагает фипитпость потоков и кусочно-непрерывных траекторий эволюции системы.
Рассмотрим следующий пример. Пусть А\{0а} = {2-г|г € N1}, 0а = 0. Поток х : [0; 1] ^ А такой, что х(£) = а1 = 2-1 при £ = £г = 1 — 2-г-1, оде г € N и х(£) = 0 при других £ е [0; 1], является финитным потоком на любом интервале [0; а], 0 < а < 1 и не является таковым на [0; 1]. События в потоке х имитируют поведение наблюдателя в апории Зенона „Ахиллес и черепаха". Значен не аг события (аг,£г) равно расстоянию между Ахиллесом и черепахой в момент ¿г; начальное событие (а0,£0) = (1;0); скорости Ахиллеса и черепахи равны, соответственно, +2 и +1. Интервал ¿г+1 — £г между событиями (аг, £г) и (аг+1, £г+1) равен времени, необходимому Ахиллесу, чтобы пробежать расстояние аг. Поток х не имитируем на интервале [0; 1] (Ахиллес не догонит черепаху за конечное время имитации). Парадокс порожден свойствами модели, а не свойствами реального (моделируемого) движения.
Во второй части данной работы будет сформулирован класс динамических систем, названных дискретными динамическими системами (ДДС), Определение ДДС будет дано
в терминах их специальной схемной реализации. Определение основывается на известных понятиях сети [1, 23-27], автоматов [26, 28-29] и динамических систем [8, 22], Базовая модель ДДС реализует функции асинхронных автоматов,
ДДС при определенных значениях параметров их динамики имитируемы, то есть их эволюция может быть воспроизведена на ЭВМ за конечное время. Если в имитацию эволюции сетей ДДС заложить механизм синхронизации параллельных процессов, идентичный механизму, принятому в языке моделирования ДИС [5, 8], то будет обеспечена независимость результатов имитации сети от порядка последовательной имитации совпадающих (в модельном времени) событий.
Список литературы
1. Бусленко Н. П. К теории сложных систем // Изв. АН СССР: Техническая кибернетика. 1963. № 5. С. 7-18.
2. Коваленко И.Н. О некоторых классах сложных систем // Изв. АН СССР. Техническая кибернетика. 1964. № 6. С. 3-9; 1965. № 1. С. 14-20; № 3. С. 3-11.
3. Бусленко Н. П., Осетинский Н. И. О динамических системах с последействием // Кибернетика. 1972. № 5. С. 123-132.
4. Dahl О.-J., Nygaard К. SIMULA — an ALGOL based simulation of discrete information systems // Comm. of ACM. 1966. V. 9. N 9. P. 671-678.
5. Netchepurenko M.I., Hairutdinov A.H., Shaboov G.V. Simulation of discrete information systems // Simulation programming languages. N.-H. 1968. P. 117-129.
6. KlNDLER E. Dynamic systems and theory of simulation // Kvbernetika. 1979. V. 15. N 2. P. 77-87.
7. Barnden J. A. Nonsequentialytv and concrete activitvphases in discrete-event simulation languages // ACM Trans, on programming languages and systems. 1981. V. 3. N 3. P. 293-317.
8. Месарович \!.. Такахара Я. Общая теория систем: математические основы. М.: Мир, 1978.
9. Нечепуренко М.И. Сингулярное время в имитации. В сб.: Вычислительная техника и системный анализ. Новосибирск, 1982, С. 9-12.
10. Нечепуренко М. И., Хайрутдинов А. X., Шавров Г. В. Язык моделирования ДИС / Докл. V-й Межвуз. конф. по физическому и математическому моделированию. \!.. 1968.
11. Дровышевич С. Г., Хайрутдинов А.Х., Чинин Г. Д. Язык моделирования дискретных информационных систем ДИС-68 / Первая всесоюз. конф. по программированию. Алгоритмические языки. Киев, 1968, С. 89-109.
12. Энгелвс Ф. Диалектика природы. М.: ГНИЛ, 1948.
13. Бэр Р. Теория разрывных функций. М.-Л.: ГТТИ, 1932.
14. Фейнман Р., Лейтон Р., Сэндс М. Фейнмановские лекции по физике. Т. 1. М.: Мир, 1965.
15. Штейнман Р. Я. Пространство и время. БСЭ (3-е изд.). Т. 21. М.: Советская энциклопедия, 1975, С. 117-120.
16. Гегель Г. В. Ф. Энциклопедия философских наук. Т. 2. Философия природы. М.: Мысль, 1975.
17. Лузин Н. Н. Теория функций действительного переменного. М.: Учпедгиз, 1948.
18. Бурваки Н. Общая топология. Топологические группы. Числа и связанные с ними группы и пространства. М.: Наука, 1969.
19. Бирхгоф Г. Теория структур. М.: ИЛ, 1952.
20. Arbib М.А., Manes Е. G. Foundation of svstem theorv: decomposable svstems //Automática. 1974. V. 10. P. 285-302. Перевод в [22]: С. 7-48.
21. Мальцев А. И. Алгебраические системы. М.: Наука, 1970.
22. Математические методы в теории систем: сб. статей. М.: Мир, 1979.
23. Бусленко В. Н., Лукьянов А. Н., Осетинский Н. И. Сложные системы с обобщенным оператором сопряжения. Теория сложных систем и методы их моделирования: Тр. семинара. \!.. 1980, С. 4-11.
24. Явлонский С. В. Основные понятия кибернетики. Проблемы кибернетики. \!.. 1959, вып. 2., С. 7-38.
25. Лупанов О. Б. О возможностях синтеза схем из произвольных элементов. Тр. МИЛИ им. В. А. Стеклова, 1958, Т. 51, С. 158-173.
26. Айзе рман М. А., Гусев Л. А. и др. Логика, автоматы, алгоритмы. \!.: ГИФМЛ, 1963.
27. Нечепуренко М. И. Модели структурного резервирования систем. Прикладные задачи на графах и сетях. Новосибирск, 1981, С. 57-86.
28. Лазар ев В. Г., Пийль В. Г. Синтез асинхронных конечных автоматов. \!.: Наука, 1964.
29. Астрахановский А. Г., Варшавский В. И. и др. Апериодические автоматы. М.: Наука, 1976.
Дата поступления — 5.05.2015
