УДК 519.62+681.3.06
Сибирский журнал науки и технологий. 2017. Т. 18, № 4. С. 810-819
АЛГОРИТМЫ СИМВОЛЬНЫХ ВЫЧИСЛЕНИЙ НА ОСНОВЕ КОРНЕВЫХ ДЕРЕВЬЕВ ДЛЯ ОЦЕНКИ ВОЗМОЖНОСТЕЙ УПРАВЛЕНИЯ
А. А. Рогалев
Сибирский федеральный университет, Институт космических и информационных технологий Российская Федерация, 660074, г. Красноярск, ул. Академика Киренского, 26 E-mail: a.rogalyov@yahoo.com
Представлены методы и алгоритмы для оценки возможностей управления, в которые входит определение возможности приведения управляемой системы в заданную точку y* фазового пространства. Для этого вычисляются границы множества достижимости управляемой системы - множества точек фазового пространства, в которые возможно сдвигаться вдоль траектории управляемой системы, и проверяется включение точки y во множество достижимости. Подобные задачи часто встречаются на практике, например, при оценке маневренных возможностей летательных аппаратов. Оценки маневренных возможностей должны быть достаточно точными и реализуемыми за малое время. Описанные алгоритмы основаны на символьных формах представления решений. При реализации этих алгоритмов подвергаются обработке, а также хранятся числовые данные и символьные данные. Для оценки областей значений выражений, построенных формулами, исследуются структуры символьных данных, включающие формальные линейные комбинации корневых помеченных деревьев. Операции над корневыми помеченными поддеревьями выполняются при вложении структуры данных в ассоциативный набор данных. Возможны сокращения членов символьных формул, если неперестановочные операторы записываются через перестановочные операторы. Такие алгоритмы уменьшают время вычисления операторов, в которых присутствуют производные, часто экспоненциально. Приводятся примеры применения этих методов. Среди этих примеров можно отметить оценки для модели управления «машина Дуббинса», представляющей систему обыкновенных дифференциальных уравнений 3-го порядка и описывающей движение беспилотного летательного аппарата (БПЛА). Поскольку скорости полета находятся в интервале 30-65 км/ч, который типичен для БПЛА, а скорость ветра на высоте 30-200 м над уровнем земли почти всегда превышает 18 км/ч, то БПЛА должны эффективно маневрировать в воздушном потоке. Вводя параметризацию кривой длинной дуги, постановка этой задачи формулируется геометрически. Это позволяет исследовать границы множеств достижимости в каждый момент времени. Анализ структуры множеств достижимости дает полезную информацию для оценки возможностей управления.
Ключевые слова: возможности управления, множества достижимости, области значений, символьные данные, корневые деревья.
Siberian Journal of Science and Technology. 2017, Vol. 18, No. 4, P. 810-819
ALGORITHMS OF SYMBOLIC COMPUTATIONS BASED ON ROOT TREES FOR EVALUATING THE CAPABILITY OF THE CONTROL PROCESS
A. A. Rogalev
Siberian Federal University, Institute of Space and Information Technologies 26, Academica Kirenskogo Str., Krasnoyarsk, 660074, Russian Federation E-mail: a.rogalyov@yahoo.com
The article presents the methods and algorithms for evaluating management capabilities. These methods also determine the possibility of bringing the controlled system to a given point y* in the phase space. A set of reachability is a set of points in the phase space into which it is possible to move along the trajectory of the controlled system. It is necessary to check whether the point y* is included in the reachable set. Such problems are often encountered in practice, for example, when assessing the maneuverability of aircraft. Estimates of maneuverability should be accurate enough and implemented in a short time. The algorithms described in the article are based on the symbolic forms of solutions. When these algorithms are implemented, numerical data and character data are processed, and also are stored. To evaluate the ranges of values of expressions constructed by formulas, the article explores the structures of character data, including formal linear combinations of root labeled trees. Operations on the root tagged subtrees are
performed by nesting the data structure in an associative data set. Reductions of the members of symbolic formulas are possible, if non-commutative operators are written through permutation operators. Such algorithms reduce the time of computation of operators in which derivatives are present, often exponentially. The article gives the examples of the application of these methods. Among these examples, we can mention the "Dobbins machine" control model, which represents the system of ordinary differential equations of the third order and describes the motion of an unmanned aerial vehicle (UAV). Since the flight speeds are in the range of 30-65 km/h, which is typical for UAVs, and the wind speed at an altitude of30-200 m above the Earth level almost always exceeds 18 km/h, UAVs must effectively maneuver in the air flow. This allows us to examine the boundaries of reachable sets at each instant of time. An analysis of the reachable sets provides useful information for evaluating management capabilities.
Keywords: control capabilities, reachability sets, ranges of values, symbolic data, root trees.
Введение. Практическое построение управления объектом выполняется как в детерминированном случае, так и при априорной неопределенности внешних возмущений, а также неопределенности состояния текущего объекта. Появление неопределенности зависит от влияния различных неопределенных факторов: возмущающих воздействий, не поддающихся контролю погрешностей параметров в определении начальных условий, и некоторых других причин. Вероятностные (стохастические) характеристики во многих случаях не могут обеспечить эффективную оценку поведения управляемой системы с неопределенными факторами [1-3]. Причины этого заключаются в том, что их использование принципиально не гарантирует результат одного конкретного опыта (испытания) и, кроме того, требует знания статистических характеристик, которые не всегда доступны в достаточной мере [4-6].
В данной статье за основу выбран подход, в котором для оценки возможностей управления строятся формулы решения задачи (в общем случае приближенные с контролем их точности), что помогает оценить области значений этих формул. Области значений этих формул определяют гарантированные границы всех возможных фазовых состояний системы управления с учетом влияния всех факторов неопределенности. Эти границы являются верхней и нижней границами множеств достижимости. В статье описываются методы, повышающие эффективность работы с символьными формулами, на которых основан расчет включений множеств достижимых множеств. Символьные формулы решений лежат в основе этих алгоритмов, их объем может достигать значительной величины. Символьная формула (аналитическое выражение) определяется как последовательность имен переменных и знаков операций, которые нужно проделать в определенном порядке над значениями этих переменных, чтобы получить значение выражения. Поэтому символьный метод (аналитический метод) представляет метод преобразования символьной информации (символьных формул) на языке анализа математических выражений. В эти формулы могут включаться числовые константы с отложенным выполнением арифметических действий над ними. Можно провести некоторую аналогию с определением из [7], где предлагается разработать аппарат для преобразования уравнений с помощью простых математических приемов, большинство из которых подобно обычным алгебраическим алгоритмам. В нашем случае этот аппарат рассчитан на использование ком-
пьютеров, поэтому используются структуры данных, составленные из формальных линейных комбинаций корневых помеченных деревьев для эффективных символьных вычислений композиций операторов. Можно напомнить, что дерево - широко распространённая программная единица для хранения и обработки однотипных данных, имитирующая древовидную структуру в виде набора связанных узлов [8], являющаяся связным графом, не содержащим циклы. В статье будут использоваться корневые конечные деревья, в которых выделена одна вершина (корень дерева).
Формально корневое дерево определяется как конечное множество Т одного или более узлов со следующими свойствами:
а) существует один корень дерева Т;
б) остальные узлы (за исключением корня) распределены среди т > 0 непересекающихся множеств Т1, ...,Тт, и каждое из множеств является корневым деревом; деревья Т1, ...,Тт называются поддеревьями данного корня Т.
Постановка задачи оценки множеств достижимости. Структура управляемого объекта описывается системой дифференциальных уравнений
di = f (t, y, u) dt
(1)
с заданным классом допустимых управлений и начальным и конечным состоянием управляемого объекта при условиях, налагаемых на правую часть системы дифференциальных уравнений:
и С) £ и с Яг, t £ Т = [/0, гР ], (2)
где У0,и - выпуклые компактные множества. Полагаем, что функция /(/, у, и) вместе со своими производными по у, и является непрерывной относительно аргументов ($, у, и) на произведении Т х К" х и и удовлетворяет условию Липшица по у £ К" с константой Ь > 0 V у1, у2 £ Я" на Т х Я" х и . Кусочно-непрерывная функция и(Г), t £ Т, - это управляющая функция, удовлетворяющая (2) и такая, что при этом управлении существует непрерывная кусочно-дифференцируемая траектория у^) системы (1) при t £ Т и заданном начальном состоянии у0. Множество допустимых управлений обозначается через и .
Для описания управляемых систем, описываемых обыкновенными дифференциальными уравнениями
[1-3], во многих задачах требуются множества достижимости Y(t) - множества точек в фазовом пространстве Rn, в которые возможно сдвинуться вдоль траектории решения задачи (1), (2) за интервал времени [t0, t ].
Из определения множества достижимости ясно, что практически необходимо численно приближать множества достижимости и вычислять границы областей, содержащих эти множества, т. е. строить внешние оценки множеств достижимости. Более детально изучены множества достижимости линейных систем дифференциальных уравнений с выпуклыми ограничениями на управления и состояния системы. В работах [4-6] представлены численные методы, позволяющие вычислять границы множеств достижимости для систем такого типа, кроме того, методы нахождения внешних оценок множеств достижимости для линейных систем описаны в работах [3; 9].
Для нелинейных систем эта задача решается значительно более сложно, возможно оценивать множества достижимости, причем предлагаемые методы обычно требуют больших вычислительных затрат. Множества достижимости нелинейных систем в общем случае невыпуклые, несвязные. Созданы методы аппроксимации множеств достижимости, использующие уравнение Беллмана, а также принцип сравнения с реализацией метода эллипсоидов, известны схемы [2; 5] аппроксимации множеств достижимости за счет вычисления характерных точек множеств непосредственно через решение серии задач оптимального управления [10; 11].
Основная задача статьи - применить методы символьных вычислений на основе корневых деревьев для улучшения реализации гарантированных методов решения дифференциальных систем с учетом управляющих воздействий и основанных на символьном представлении решений. Более полное описание этого подхода можно найти в [12-18].
Пусть справедлива равномерная оценка | y(t) |< b
для всех решений (1) на интервале t е [t0,T], где b = const > 0; множество достижимости Y (t) = f (t, y,^) для всех у, t e[s,T], является компактным и выпуклым множеством.
Как было отмечено выше, символьная формула решений определяется сдвигом вдоль траектории решения и получается в процессе последовательного преобразования формул вида
Yn = pn (tо fn yо Y1 Yn) = = Sn (Y 0) 0...0 S2 (Yi_1) о S1 (Y'). (3)
Пусть ф(y0) обозначает однозначное отображение единичного интервала из R1 на прямоугольный параллелепипед из Rn, которое каждой точке t е R ставит в соответствие некоторую точку y = ф^).
Такое отображение определяет процедуру исполнения, в которой для каждой точки t е R определяется
формула отображения F (y/
0 у0 Y 2 ,
, Y
1) и процесс ее
сборки по адресам. Для этого возможно использовать ф(у01 - непрерывное, однозначное отображение
единичного интервала на «-мерный куб, называемое кривой Пеано, заполняющей пространство и представляющей непрерывную, нигде не дифференцируемую кривую, которая проходит через все точки единичного гиперкуба в пространстве К". Это соответствие задано отображением элементов конечного множества отрезков из единичного интервала и элементами конечного множества гиперкубов, входящих в К". Формула будет представлять рекурсивную структуру, размер которой изменяется. Для записи такой формулы в компьютере используются линейные динамические структуры.
В рассматриваемых задачах управление является измеримой функцией, что не допускает применения методов рядов Тейлора для построения опорных траекторий.
В символьный метод оценки множеств достижимости входит запись компонент символьных формул решений как векторных функций, состоящих из символьных компонент 5 (1:к ), зависящих от символьных
форм начальных данных у°, . .., у° и описывающих сдвиг вдоль кривой, аппроксимирующей траекторию решения системы, на каждом шаге. Сдвиг вдоль траектории решений определяется на основе построенных символьных формул решений. Символьные формулы не преобразуются на каждом шаге алгоритма, они хранятся в памяти компьютера, для этого полезным являются символьные формулы кусочно-полиномиальных функций. Используя символьную формулу У, аппроксимирующую оператор сдвига вдоль траектории, в методе определяются прообразы экстремальных значений (верхних и нижних границ для множеств решений) как решения экстремальной задачи. Эти точки определяют границы множеств достижимости.
Выполнение преобразований символьных формул -это первый этап гарантированного метода включения решений, за которым следует этап числовых вычислений. Символьные преобразования реализуются в системах компьютерной алгебры. Однако такое применение компьютерной алгебры, как правило, является крайне затратным по времени исполнения, особенно если вычисления производятся в цикле. Найденные формулы часто являются громоздкими, что и приводит к неэкономным вычислениям, особенно если речь идет о вычислениях в цикле. Поэтому в статье предлагается хранить и обрабатывать формулы, хранящиеся в виде помеченных деревьев.
Преобразования символьных формул. Вычисления, рассматриваемые в статье, представлены как алгебраические выражения над входными данными, которые сохраняют важнейшие соотношения между входными данными и результирующей аналитической информацией. В предлагаемых методах символьных вычислений моделируются операторы присваивания, условные операторы, циклы и процедуры. Это помогает повысить эффективность и точность при настойчивом применении процедур упрощения.
Определим умножение на укоренившихся помеченных деревьях, тем самым создавая набор этих структур данных как ассоциативную алгебру. Алгебра гомоморфизмов задана на исходной алгебре операторов, отображаемой в данную алгебру деревьев. Если операторы представить на этой структуре данных, то появятся сокращения, связанные с тем, что неперестановочные операторы выражаются через коммутирующие операторы. Это приводит к алгоритму, применение которого для определения формул производных ускорит вычисления оператора экспоненциально. Многие вопросы, описанные в этом разделе, опираются на работы [19; 20].
Кратко напомним основные определения, используемые в этой работе.
Алгебра над полем - это векторное пространство, снабженное билинейным произведением. Отсюда следует, что алгебра над полем является одновременно векторным пространством и кольцом, причём эти структуры согласованы.
Алгебра Ли - это векторное пространство Ь над полем к; на пространстве введено билинейное отображение [-, -]: Ь х Ь ^ Ь , удовлетворяющее условиям:
[ х, у] + [у, х] = 0 для всех х, у е Ь ; [ х,[ у, г ]]+[ у, [ 2, х ]] + [ 2, [ х, у ]] = 0 для всех х, у, г е Ь .
Так, пример алгебры Ли выглядит следующим образом: билинейное отображение [-, -] определено
как [ х, у ] = ху - ух для х, у е А . Дифференцированием в алгебре и называется линейное отображение 5 : и ^ и , удовлетворяющее правилу Лейбница дифференцирования произведения
5{аЬ ) = 5(а )Ь + а5(Ь ).
Пусть Я - это коммутативная алгебра с единичным элементом над полем к характеристики О. Пусть Д, ...,Бн - это N перестановочных производных Я, т. е. для I,у = 1, ...,N , ООуа = ОуОа для всех а е Я .
Предположим, что мы также получили М производных О//1, ...,Б1/М е Я , которые могут быть выражены как Я линейных комбинаций производных Ц, т. е. для у = 1, ..., М,
Difj ' a1 е R
ц=1
(4)
Diff (Dj, ...,Dn;R) является алгеброй всех отображений R ^ R , порожденных отображениями D ^ L(a), a е R (отображение L(a), a е R определяется зависимостью L(a)(b) = ab для b е R). Пусть
X : k(Ej, ..., EM) ^ Diff (Dj, ..., Dn ; R) (5) обозначает отображение, которое отображает p е k{Ej, . ., EM} в линейный дифференциальный
оператор %(p), полученный подстановкой (1) и упрощением, с использованием того факта, что D^ являются
производными R. Пусть p е k(E1, ...,EM} имеет вид
i
p = ^pi, где каждый член pi имеет степень m. Для
i=1
расчета значения %(p) будем вычислять %(pi). Это
даст значение формулы lm! Nm членов. Предположим, что Expenditure(p) - это затраты применения
алгоритма A для упрощения p е k (E1,..., Em ). Затраты пропорциональны количеству дифференцирований и умножений.
Тогда Expenditure(p) = O (imn! Nm ) .
Simple
Проанализируем алгоритм, выполняющий предварительную обработку выражения p, целью которой является устранение любых членов, сокращаемых после подстановок (1), из вычислений.
Возможно сформулировать следующую теорему. Теорема 1. Предположим:
1) p - это сумма l = 2m_1 членов, каждый из которых является однородным относительно степени m;
2) X(p) - это линейный дифференциальный оператор степени 1;
3) m, N таким образом, что 2m << Nm .
Тогда
ExpenditureBetter (p)
= O
1
(6)
Для создания эффективного символьного алгоритма необходимо записать производные высокого порядка, порожденные О//1, ..., Б1/М, с помощью членов перестановочных производных О1, ...,DN .
Обозначим свободную ассоциативную алгебру символом О/1, ..., О1/М, и пусть (О1, ..., ; Я) -это пространство формальных линейных дифференциальных операторов с коэффициентами из Я, т. е.
ЕхрепйиигеШтр1г (р) ^ т • 2
Напомним, что скобка Ли [ , ]: g х g ^ g - это билинейная, кососимметрическая операция, удовлетворяющая тождеству Якоби: [[х, у], г] + [[г, х], у] + + [[у, 2], х] = 0.
Значение теоремы 1 состоит в том, что алгоритм, представленный для упрощения вычисления производных, входящих в формулу х(Р), является методом, который часто экспоненциально улучшает эффективность вычислений. Конечно, если мы априорно знаем определенные свойства %(р), то мы можем использовать эти свойства, чтобы вычисления были эффективнее. Тем не менее, наш алгоритм не требует знания свойств %(р).
Сокращения, которые появляются, когда неперестановочные операторы выражены через коммутирующие операторы, возникают естественно. Это приводит к алгоритму, применение которого для операторов, таких как производные, ускоряет вычисления экспоненциально.
Пусть {е^ ...,Ем | - множество символов. Будем
называть дерево помеченным метками |e1; ..., EM j, если каждому узлу дерева, отличному от корня, соответствует элемент |е1, ...,EM j, присвоенный ему. Обозначим множество всех деревьев, помеченных метками |е1, ...,EM j, через Forest{,...,EM j.
Пусть k |Forest{, ..., EM jj обозначает векторное
пространство над полем k с базисом Forest {, ..., Em j.
Определим умножение в k[Forest[е1, ...,EM jj следующим образом. Поскольку набор помеченных деревьев образует базис для k[Forest[е1, ...,EM jj
векторного пространства, достаточно описать произведение двух помеченных деревьев.
Пусть ti, t2 - два помеченных дерева. Пусть si,s2, ...,sr - дети корня ti. Если t2 имеет n +1 узлов (считая корень), существуют (n + l)r способов присоединения r поддеревьев ti, для которых si, ...,sr используются как корни помеченного дерева t2. Присоединим каждый si как потомок некоторого узла t2, сохраняющий исходные метки. Произведение деревьев tit2 определяется как сумма этих (n + l)r помеченных деревьев. Можно показать, что это произведение ассоциативно и что дерево, состоящее только из корня, является мультипликативным тождеством (элементом, являющимся единичным для мультипликативной операции) (см. [19; 20] для доказательства). Определим ко-умножение над k {Forest {Ei, ..., ЕМ jj
следующим образом.
Пусть t - помеченное дерево, и пусть si, ..., sr -
дети корня t. Если P - подмножество Ct = {si, ...,sr j,
полагаем, что tp - помеченное дерево, образованное
путем создания (размещения) элементов из P как потомков нового корня. При этом сохраняются исходные метки. Определим ко-умножение как операцию A(t) = ^ t ®tc xP , где X \ Y обозначает операцию
PcQ t
теоретико-множественного относительного дополнения Y в X. Определим ко-единицу измерения, полагая е(/) равным 1, если t имеет только один узел (его корень), и равным 0 в противном случае. Можно показать, что это делает k{Forest(Ei, ...,EM )j ко-коммутативной ко-алгеброй. Возможно определить шкалу градации на k{Forest(Ei, ...,EM)j, полагая k{Forest(Ei, ...,EM)j
подпространством k{Forest(Ei, ...,EM )j, натянутым на деревья с n +1 узлами и обозначенным через
k {Forestn (Ei, ..., Ем )j.
Помеченное упорядоченное дерево является помеченным деревом, для которого существует линейное упорядочение потомков каждого узла. Обозначим векторное пространство над k с базисом из множества помеченных упорядоченных деревьев через k{LOF(E1, ...,EM)} . Определения операций умножения и ко-умножения для k{LOF(E1, ...,EM)} аналогичны определениям, приведенным выше.
Следующая теорема доказана в работе [20].
Теорема 2. Набор помеченных деревьев t, корень которых имеет ровно одного потомка, является базисом для P(k {Forest {Ej, ..., Eu }}) .
Доказательство. Непосредственно получим, что любое дерево, корень которого имеет только одного потомка k, охватывает примитивные элементы. Пусть O = k{Forest{Ej, ...,Eu}}, и определим линейное отображение л : O ® O ^ O следующим образом: если tj и t2 - помеченные деревья, то л ( ® t2 ) - это помеченное дерево, образованное путем идентификации корней tj и t2. Другими словами, л ( ® t2 ) является помеченным деревом, для которого поддеревья корня образуют все поддеревья корней tj и t2. Легко видеть, что если t - помеченное дерево, корень которого имеет r потомков, то л о Д (t) = 2r t. С другой стороны, если a = ^ ait е P(O), мы получаем, что л о Д(а) = 2а . Так как деревья t являются линейно независимыми, следует, что ai = 0 , если корень t имеет более одного потомка. Это завершает доказательство теоремы.
Далее опишем, как помеченные деревья могут быть использованы для упрощения вычислений дифференциальных операторов. Начнем с определения отображения
у : k{Forest{Ej, ...,EM }} ^ Diff (Dj, ...,DN;R) следующим образом.
Шаг 1. Задав помеченное дерево t е LFm (Ej, ...,EM ), ставим в соответствие корню число 0 и назначаем остальным узлам числа j, ...,m . Мы в дальнейшем идентифицируем узел с числом, назначенным узлу. Чтобы определить отображение, мы просуммируем индексы ...,цт , связанные с каждым узлом дерева, кроме корня. Закрепим узел k дерева t, и пусть l, ...,l' обозначают детей этого узла.
Пусть
ID.....D,, аЦ1, если k - не корень,
R(k; ц, ..., Цг ) = | w Цг Yk
[Dw -D, если k - к°рень.
Обозначим эту величину R(k). Заметим, что R(k ) е R при k > 0.
Шаг 2. Определим
y(t) = £ R(m)...R(l)R(0)
ЦЬ -,Цт =j
и расширим у на все k{Forest{Е1, ...,ЕМ }} в силу его линейности.
Следующие три предложения описывают основные свойства отображения у. Первое предложение является примером упрощения факторизации х с помощью множества помеченных деревьев. Легко получить, что часто менее затратно вычислять у и ф совместно, чем вычислять %.
Предложение 1
а) отображение у является гомоморфизмом алгебры;
б) Х = V ° ф.
Доказательство. Доказательство (а) получено непосредственным вычислением по правилу Лейбница и содержится в работе [19]. Так как х и у°ф согласованы (располагаются) на генерируемом множестве Е1, ..., EM , то часть (б) предложения 1 следует из части (а). Это завершает доказательство предложения.
На самом деле верно более сильное утверждение: отображение у соблюдает взаимодействие ко-умножения на k {Forest{Е1, ..., ЕМ }} -ассоциативной алгебре символов над помеченными деревьями и умножение в R в следующем смысле.
Предложение 2
Для всех a,b е R и для всех t е k{LF(Е1, ...,EM)}
((у ® у) ° A(t)) (a ® b ) = у (t) (ab ).
Для многих приложений важно знать, ограничиваются ли действия этих гомоморфизмов областью алгебры Ли, порожденной производными Е1, ...,Ем .
Следующее ниже утверждение показывает, что если известны вычисленные элементы алгебры Ли, то сле-
дует рассматривать только деревья, корень которых содержит только один узел - потомка.
Предложение 3
а) пусть Ь (Е1, . .., ЕМ) обозначает алгебру Ли, порожденную производными Е1, ...,ЕМ, тогда
ф(Ь ((, ...,Ем )) Р(к{((, ...,Ем )}),
Р(к{ЬЕ((, ...,Ем )}))с Оег Я;
б) элементы Ь (Е1, . . ., Ем) отображаются под действием ф в суммы деревьев, корень которых имеет только одного потомка.
Завершим этот раздел простым примером вычисления производной третьего порядка с использованием деревьев.
Фиксируем производные
Еу =£аУО, 1< ] < 3
Ц=1
и рассмотрим производные высшего порядка
Р = Е3Е2Е1 - Е3Е1Е2 - Е2Е1Е3 + Е1Е2Е3.
Вычисление %( Р) «в лоб» включает в себя вычисление 24N3 членов: 8N3 - члены-сомножители О, , ,
2 N3 из них сокращаются, 12N3 членов, содержащих множители О, , все из них сокращаются, и 4N3
членов, включающих множители О,..... , все из них
сокращаются.
Вычисляя ф( р), получим
ф(р) = о
Е-
Е,
Е3 Е1
О + о
О Е1 О Е О Е3 о
Ез
О Е, о Е1 о Е1 о Е
Ö Ез Ö Ез О Е, О Е
+
Все сокращения членов при вычислениях %(р) = У(ф(р)) происходят при вычислении ф(р).
Сделаем следующие предположения: р ек(^Е1, ...,Ет} (элемент ассоциативной алгебры символов) есть сум-
ма
p=У Р'
где каждый член pt имеет степень m;
'=1
затраты на умножение составляют одну единицу; затраты на дифференцирование - одна единица; затраты на сложение равны нулю; затраты добавления узла к дереву - одна единица, так что стоимость построения дерева : е ЬЕт (Е1, ...,Ет) равна т единиц. Предложение 1
Вычислительные затраты на вычисление %(р)
равны величине 21тт! Ыт.
Доказательство. Предположим, что р^ имеет вид
E,(m... E для некоторых индексов 1 < уь Тогда x(p) равно
, У m < M .
( N
У
гЦтД,
4m
А ( N
У О11 D
L^ у/ Ц/ VVi =1
Л
членов,
x(p ) = Ух(Рг) составляют 2lmm! N"
Cost,
BETTER
Используя теорему 3, можно построить еще более эффективные алгоритмы для вычисления результата воздействия элементов k(E1, ...,EM}, которые, как
известно, являются дифференциальными операторами. В этом случае деревья, для которых корень имеет более одного потомка, даже не нужно строить.
Приложения
1. Рассмотрим движение управляемого объекта на плоскости, описываемого следующей системой дифференциальных уравнений:
dyi I/ —1 = V cos ф, dt
^ = V sin ф,
dt
d Ф k ¡I, • = — u,u < 1,
(7)
После разложения появляются т! Ыт каждый из которых включает в себя т дифференцирований и т умножений (включая умножения и дифференцирования, связанные с применением оператора
Ф)).
Предложение 2
Затраты на вычисление ф( р) равны 1тт!.
Доказательство. Моному степени т отображение ф ставит в соответствие сумму (объединение) т! помеченных деревьев. Этот факт легко доказать по индукции (доказательство в [19]). По предположениям, представленным выше, затраты на построение помеченного дерева с т узлами (в дополнение к корню) составляют т единиц. Поэтому общие затраты равны 1тт!.
Предложение 3
Пусть ст = ф(р) обозначим через |ст| число помеченных деревьев с ненулевыми коэффициентами с. Тогда затраты на вычисление равны 2т|ст| Ыт .
Доказательство. Фиксируем помеченное дерево : е ЬЕт (Е1, ...,Ем). Из определения отображения у мы видим, что затраты на вычисление у(:) равны 2тЫт, и, следовательно, полные затраты равны 2т |ст| Ыт.
Объединяя эти три утверждения, получим теорему.
Теорема 3. В соответствии с принятыми выше предположениями затраты Соз:$ 1е (р) вычисления
dt V V = const, k = const > 0.
Значения угла ф рассматриваются на интервале (-л, л). Здесь y1, y2 - координаты объекта, отождествляемого с точкой на плоскости; ф - угол между вектором скорости объекта и осью x; u - управляющий параметр, удовлетворяющий указанному ограничению и характеризующий скорость изменения угла ф; k - максимальное боковое ускорение; V - величина k
скорости; а = ^ ■
Неравенство в (7) ограничивает радиус кривизны траектории объекта, а именно, радиус кривизны не может больше единицы. Рассматриваемая система использовалась Р. Айзексом при постановке задачи «шофер-убийца» (рис. 1). Состояние z0 = (y1(t0), y2(t0), ф(^)) объекта в начальный момент времени предполагается заданным.
Множество достижимости G(T) в момент времени T > 0 есть совокупность всех точек фазового пространства, в каждую из которых возможен перевод системы (7) в момент времени T > 0 при помощи
некоторого допустимого управления на промежутке
*
T = [0, t ] из начальной точки z0.
2. Рассмотрим систему обыкновенных дифференциальных уравнений, в правую часть которой входит управляющее воздействие:
dyL dt
= 2 y2 + y4,
dyi dt
= -2 yi + y3 + u (t),
dy^ dt
(8)
= y2 + 10 y4,
тогда как
(p) вычисления L = у о ф(p) равны
lmm!+ 2m ст N"
аУГ = 2 У1 " 2 Уз + и2(:). а:
Система (8) - это гамильтонова система, имеющая гамильтониан
н = У12 + У2У4 + 5У4 + У12 - У1 Уз + Уз2 - и1(:)У1 - и2(:^з. В этой формуле У1, у2 - обобщенные координаты, у2 , у4 - обобщенные импульсы.
'=1
На возмущающие силы системы и1 (/), и2 (/) наложены ограничения
{{ < 1,Ы < 1}.
В [1] для этой системы были выбраны начальные данные в виде параллелепипедов размерности 4. Динамические оценки множеств достижимости были построены в виде эллипсоидов [1]. Поэтому полезно
предложить гарантированные оценки множеств достижимости, качественно совпадающие с результатами [1], показывающими отсутствие экспоненциального роста границ эллипсоидов. В [1] графики построены в логарифмической шкале с несколькими числовыми метками. Других числовых данных для этой системы в [1] не приводится (рис. 2, 3).
Рис. 1. Гарантированные границы множества достижимости задачи (7) - проекция на оси t-y1-y2 задачи «шофер-убийца»
Fig. 1. Guaranteed set boundaries of tasks (7) -projection on an axis t-y1-y2 "driver-killer" tasks
Рис. 2. Гарантированные границы множества достижимости задачи (8) - проекция на оси Г-у1: 1 - верхняя граница проекции множества достижимости; 2 - нижняя граница проекции множеств достижимости
Fig. 2. Guaranteed set boundaries of tasks (8) - projection on an axis t-y 1. 1 - upper limit of set projection of reachability; 2 - lower limit of set projection of reachability
6
-6 J
Рис. 3. Гарантированные границы множества достижимости задачи (8) - проекция на оси t-y2: 1 - верхняя граница проекции множества достижимости; 2 - нижняя граница проекции множеств достижимости
Fig. 3. Guaranteed set boundaries of tasks (8) - projection on an axis t-y2: 1 - upper limit of set projection of reachability; 2 - lower limit of set projection of reachability
Заключение. В статье описаны символьные алгоритмы, улучшающие применение гарантированных методов оценки решений дифференциальных уравнений для вычисления включений множеств достижимости управляемых систем. В статье применяется операция умножения на укоренившихся помеченных деревьях, что создает набор этих структур данных как ассоциативную алгебру. Затем определяется алгебра гомоморфизма из исходной алгебры операторов в данную алгебре деревьев. Сокращения, которые появляются, когда неперестановочные операторы выражены через коммутирующие операторы, возникают естественно, когда операторы представлены с использованием этой структуры данных. Это приводит к алгоритму, применение которого для операторов, таких как производные, ускоряет вычисления оператора, что повышает эффективность применения гарантированных методов.
Библиографические ссылки
1. Куржанский А. Б. Управление и наблюдение в условиях неопределенности. M. : Наука, 1977. 392 с.
2. Черноусько Ф. Л. Оценивание фазового состояния динамических систем. М. : Наука, 1988. 319 с.
3. Chernousko F. L. State Estimation for Dynamic Systems. Boca Raton : CRC Press, 1994. 304 p.
4. Овсеевич А. И., Шматков А. М. К вопросу о сопоставлении вероятностного и гарантированного подходов к прогнозу фазового состояния динамических систем // Известия Академии наук. Теория и системы управления. 2007. № 4. С. 11-16.
5. Черноусько Ф. Л. Эллипсоидальные аппроксимации множеств достижимости управляемых линейных систем с неопределенной матрицей // Прикладная математика и механика. 1996. Т. 60, № 6. С. 940-950.
6. Куржанский А. Б., Фурасов Б. Д. Задачи гарантированной идентификации билинейных систем с дискретным временем // Известия Академии наук. Теория и системы управления. 2000. № 4. С. 5-12.
7. Кнут Д. Э. Искусство программирования. Т. 2. Основные алгоритмы. Т. 2. Основные алгоритмы. 3-е изд. М. : Вильямс, 2000. 832 с.
8. Шеннон К. Работы по теории информации и кибернетике. М. : Иностранная литература, 1963. 832 с.
9. Пацко Б. В., Пятко С. Г., Федотов А. А. Трехмерные множества достижимости нелинейных управляемых систем // Известия Академии наук. Теория и системы управления. 2003. № 3. С. 8-16.
10. Тятюшкин А. И., Моржин О. В. Численное исследование множеств достижимости нелинейных управляемых дифференциальных систем // Автоматика и телемеханика. 2011. Вып. 6. С. 160-170.
11. Финкельштейн Е. А., Горнов А. Ю. Алгоритм квазиравномерного заполнения множества достижимости нелинейной управляемой системы // Известия Иркутского государственного университета. Сер. «Математика». 2017. Т. 19. С. 217-223.
12. Новиков В. А., Рогалев А. Н. Построение сходящихся верхних и нижних оценок решения систем обыкновенных дифференциальных уравнений // Жур-
нал вычислительной математики и математической физики. 1993. Т. 33, № 2. С. 219-231.
13. Рогалев А. Н. Включение множеств решений дифференциальных уравнений и гарантированные оценки глобальной ошибки // Вычислительные технологии. 2003. Т. 8, № 6. С. 80-94.
14. Rogalyov A. N. Computation of reachable sets guaranteed bounds // Proceedings of the IASTED International Conferences on Automation, Control, and Information Technology - Control, Diagnostics, and Automation (ACIT - CDA 2010). Calgary, Canada : ACTA Press, 2010. B6. Р. 132-139.
15. Рогалев А. Н. Вычисление гарантированных границ множеств достижимости управляемых систем // Автометрия. 2011. Т. 47, № 3. С. 100-112.
16. Rogalev A. N. Calculation of Guaranteed Boundaries of Reachable Sets of Controlled Systems // Optoelectronics, Instrumentation and Data Processing. 2011. Vol. 47. № 3. P. 287-296.
17. Рогалев А. Н., Рогалев А. А. Численный расчет включений фазовых состояний в задачах наблюдения за движением самолета // Вестник СибГАУ. 2012. № 1(41). С. 53-57.
18. Рогалев А. Н., Рогалев А. А. Численные оценки предельных отклонений траекторий летательных аппаратов в атмосфере // Вестник СибГАУ. 2015. Т. 16, № 1. С. 97-104.
19. Grossman R. Evaluation of expressions involving higher order derivations // J. Math. Systems, Estimation, and Control. 1991. № 1. Р. 91-106.
20. Grossman R., Larson R. G. Hopf-algebraic structures of families of trees // J. Algebra. 1989. Vol. 126. Р. 184-210.
References
1. Kurzhanskiy A. B. Upravlenie i nablyudenie v us-loviyakh neopredelennosti [Control and observation under uncertainty]. Moscow, Nauka Publ., 1977, 392 p.
2. Chernous'ko F. L. Otsenivanie fazovogo sostoy-aniya dinamicheskikh sistem [Dynamical systems phase state evaluating]. Moscow, Nauka Publ., 1988, 319 p.
3. Chernousko F. L. State Estimation for Dynamic Systems. Boca Raton, CRC Press, 1994, 304 p.
4. Ovseevich A. I., Shmatkov A. M. [Probability and guaranteed approaches for dynamical systems phase state forecasting]. Izvestiya Akademii Nauk. Teoriya i sistemy upravleniya. 2007, No. 4, P. 11-16 (In Russ.).
5. Chernous'ko F. L. [Ellipsoidal approximation of reachable sets of controlled systems with matrix of uncertainty]. Prikladnaya matematika i mekhanika. 1996, Vol. 60, No. 6, P. 940-950 (In Russ).
6. Kurzhanskij A. B., Furasov B. D. Zadachi garanti-rovannoj identifikatsii bilinejnykh sistem s diskretnym vremenem, Izvestiya Akademii Nauk. Teoriya i sistemy upravleniya.2000, Vol. 4, P. 5-12 (in Russ).
7. Donal'd Knut. Iskusstvo programmirovaniya. Osnov-nye algoritmy [Art of programming. Basic algorithms]. Moscow, Vil'yams Publ., 2000, 832 p.
8. Shennon K. Raboty po teorii informatsii i kiber-netike [Notes on information theory and cybernetics]. Moscow, Inostrannaya literature Publ., 1963, 832 p.
9. Packo B. V., Pyatko S. G., Fedotov A. A. [Three-dimensional reachable sets of nonlinear controlled systems]. Izvestiya Akademii Nauk. Teoriya i sistemy upravleniya. 2003, Vol. 3, P. 8-16 (In Russ.).
10. Tyatyushkin A. I., Morzhin O. V. [Nonlinear controlled differential systems reachable sets numerical estimation]. Avtomatika i telemekhanika. 2011, No. 6, P. 160-170 (In Russ).
11. Finkel'shtejn E. A., Gornov A. Yu. [Algorithm of quasiuniform fill of Nonlinear controlled systems reachable sets]. Izvestiya Irkutskogo gosudarstvennogo univer-siteta. Seriya Matematika. 2017, Vol. 19, P. 217-223 (In Russ.).
12. Novikov V. A., Rogalev A. N. [Upper and lower converging bounds estimation constructing for ODE systems solutions]. Zhurnal vychislitel'noy matematiki i matematicheskoy fiziki. 1993, Vol. 33, No. 2, P. 219-231 (In Russ.).
13. Rogalev A. N. [Differential equations solutions sets inclusions and guaranteed estimation of global error].
Vychislitel'nye tekhnologii. 2003, Vol. 8, No. 6, P. 80-94 (In Russ.).
14. Rogalyov A. N. [Computation of reachable sets guaranteed bounds]. Proceedings of the IASTED International Conferences on Automation, Control, and Information Technology - Control, Diagnostics, and Automation. ACTA Press, B6, Calgary, Canada. 2010, P. 132-139.
15. Rogalev A. N. [Estimation of guaranteed bounds of controlled systems reachable sets]. Avtometriya. 2011, Vol. 47, No. 3, P. 100-112 (In Russ.).
16. Rogalev A. N. [Calculation of Guaranteed Boundaries of Reachable Sets of Controlled Systems ]. Optoelectronics, Instrumentation and Data Processing. Allerton Press, 2011, Vol. 47, No. 3, P. 287-296.
17. Rogalev A. N., Rogalev A. A. [Numerical computations of phase states inclusions for problems of aircraft displacement inspection]. VestnikSibGAU. 2012, No. 1(41), P. 53-57 (In Russ.).
18. Rogalev A. N., Rogalev A. A. [Numerical estimations of maximum deviations for aircraft trajectories in the atmosphere]. Vestnik SibGAU. 2015, Vol. 16, No 1, P. 97-104 (In Russ.).
19. Grossman R. [Evaluation of expressions involving higher order derivations]. J. Math. Systems, Estimation, and Control. 1991, No. 1, P. 91-106.
20. Grossman R., Larson R. G. Hopf-algebraic structures of families of trees. J. Algebra. 1989, Vol. 126, P. 184-210.
© Рогалев А. А., 2017