МАТЕМАТИКА
J
УДК 57.519
ПРЕОБРАЗОВАНИЕ АФФИННЫХ СИСТЕМ СО СКАЛЯРНЫМ УПРАВЛЕНИЕМ К КВАЗИКАНОНИЧЕСКОМУ ВИДУ
С.Б. Ткачев, А.А. Шевляков
МГТУ им. Н.Э. Баумана, Москва e-mail: aash29@gmail.com
Для аффинных систем со скалярным управлением получены локальные условия, при которых они преобразуются в систему квазиканонического вида. Разработан алгоритм поиска функций, определяющих такое преобразование, и рассмотрены особенности его реализации с использованием систем аналитических вычислений. Приведен пример приложения теории.
Ключевые слова: теория управления, аффинные системы, нелинейные системы, квазиканонический вид.
TRANSFORMATION OF AFFINE SYSTEMS WITH SINGLE INPUT TO QUASI-CANONICAL FORM
S.B. Tkachev, A.A. Shevlyakov
Bauman Moscow State Technical University, Moscow e-mail: aash29@gmail.com
For affine systems with single input, local conditions are obtained, under which they are transformed to a system of quasi-canonical form. An algorithm is developed for searching the functions that define this transformation, and peculiarities of its implementation using the systems of analytical calculations are considered. An example of theory application is given.
Keywords: control theory, affine systems, nonlinear systems, quasi-canonical form.
Введение. Один из подходов к решению задачи управления нелинейной динамической системой основывается на преобразовании системы к специальному виду, для которого метод решения соответствующей задачи управления известен. Примером являются системы канонического вида [1-3], которые с помощью линеаризации обратной связью можно преобразовать в линейную систему, записанную в канонической форме Бруновского [4].
Условия приводимости к каноническому виду хорошо известны [3], однако не всякую аффинную систему можно к этому виду преобразовать. Поэтому среди аффинных систем выделяют системы, которые преобразуются к квазиканоническому виду [5]. Такие системы содержат подсистему, которая линеаризацией обратной связью преобразуется в каноническую форму Бруновского, и подсистему общего вида.
Основные теоретические положения о преобразовании к квазиканоническому виду в некоторой открытой области аффинных систем со скалярным управлением приведены в работе [5], а с векторным управлением — в [6].
Представляет интерес получение различных локальных условий существования требуемых преобразований, а также условий, при выполнении которых квазиканонический вид имеет специальные свойства.
Проверка условий, при выполнении которых аффинная система преобразуется к квазиканоническому виду, нахождение соответствующей замены переменных, если она существует, а также запись системы в новых переменных требуют выполнения значительных объемов аналитических вычислений.
Компьютерную алгебру начали применять для решения задач нелинейной теории управления в 1990-х гг. [7]. В настоящее время существуют пакеты программ, предназначенные для решения задачи стабилизации положения равновесия путем преобразования к каноническому виду [8]. Однако актуальными остаются разработка алгоритма и создание пакета программ, позволяющих в автоматическом или полуавтоматическом режимах проверить условия существования и выполнить преобразование аффинной системы к квазиканическому виду.
В рамках проводимого исследования использованы системы Maple и MATLAB, выбор которых обусловлен распространенностью системы MATLAB и обширными возможностями Maple для проведения символьных вычислений.
Квазиканонический вид аффинной системы. Рассмотрим аффинную систему со скалярным управлением
x = A(x) + B (x)u, (1)
где x E Rn; u E R1; A(x) = (ai(x),..., an(x))T, A(0) = 0; B(x) =
= (b1(x),... ,bn(x))T; a,t(x), ht(x) E С^(П); П — открытое множество, содержащее положение равновесия x = 0.
Пусть система (1) в области П преобразуется к виду
Zi = Z2, . . . , Zr-1 = Zr,
zr = f (z,n) + g(z,n)u, (2)
n = q(z, n) + p(z, n)u,
где z = (z1,..., zr)T, z E Rr; n = (n1,..., nn-r)T, n E Rn-r; q(z, n) =
= (q1(z,n),...)qn-r (z,n))T; p(z,n) = (P1(z,n),...,Pn-r (z,n))T.
Указанный вид при г < п называют квазиканоническим видом [5], а при г = п — каноническим видом [1] системы (1). Число г называют индексом приводимости.
Система (2) квазиканонического вида регулярна в точке (г°,п°), если коэффициент при управлении д(г,п) не обращается в точке (г°, п°) в нуль.
Говоря о том, что в области П стационарная аффинная система преобразуется к стационарной системе того или иного вида, предполагают [3], что существует диффеоморфизм Ф : П ^ Ф(П), который переводит траектории одной системы в траектории другой системы, соответствующие тем же управлениям.
Пусть замена переменных (г, п) = Ф(х), задающая преобразование системы (1) к квазиканоническому виду (2), выбрана так, что Ф(0) = (0, 0). Тогда точка (г, п) = (0, 0) будет положением равновесия системы (2). Если система (2) является минимально фазовой [9, 10], а ее квазиканонический вид регулярен в этой точке, то управление
и = п) + Е ад^ |п) (3)
обеспечивает локальную стабилизацию положения равновесия (г, п) =
= (0, 0).
Удобным для анализа свойств аффинных систем является аппарат дифференциальной геометрии. В рамках дифференциально-геометрического подхода системе (1) на области П можно взаимно-однозначно сопоставить гладкие векторные поля
п Я п Я
А = Е *(*)^ в = ЕВД^ (4)
г=1 г г=1 г
Напомним, что для гладкой функции <^(ж) определена производная этой функции по векторному полю [3]. Например, производная функции по векторному полю А, записанная в координатах х, имеет вид
А<(ж) = аг(ж)
г=1
Если в замене переменных (г, п) = Ф(х) координатные функции П1 = п1(х),... ,Пп-г = Пп-Г(х) выбраны таким образом, что
В^г(ж) = 0, i = 1, n — r, (5)
то система (2) в переменных z, п примет вид
Zi = Z2, . . . , Zr-1 = Zr,
zr = / (z,n) + g(z,nK (6)
п = q(z, п),
который будем называть специальным квазиканоническим видом. Такой вид более удобен для анализа свойств системы.
Условия существования квазиканонического вида. Следуя работам [5, 6], приведем результаты, относящиеся к преобразованию аффинных систем (1) к квазиканоническому виду (2).
Напомним, что коммутатором гладких векторных полей X и У называют [11] векторное поле [X, У], координаты которого можно вычислить через координаты X (ж) и У (ж) векторных полей X и У по формуле
КУК.)- ^ X (ж) - ^ У (ж).
Зададим последовательность коммутаторов векторных полей А и
В: а^В - В, а^В - [А, В], а^+1В - [А, а^В], к > 1.
Теорема 1. [5] Для того чтобы в некоторой области П для аффинной системы (1) существовали переменные, в которых она имеет квазиканонический вид (2), необходимо и достаточно, чтобы существовала функция <^(ж) € Сте(П), удовлетворяющая в П системе уравнений
а^В^(ж) - 0, к - 0, г - 2, (7)
и существовали функции п - П (ж) € Сте(П), 3 - 1,п — г, которые вместе с функциями
г - Аг-1^(ж), г -1~Т, (8)
задают в П гладкую невырожденную замену переменных (г, п) - <^(ж).
Следствие 1. Пусть выполнены условия теоремы 1 и дополнительные координатные функции п - П(ж), 3 - 1,п — г, удовлетворяют в области П соотношениям (5). Тогда в П аффинная система (1) заменой переменных (г, п) - <^(ж) преобразуется к специальному квазиканоническому виду (6).
Локальные условия существования квазиканонического вида. Локальные условия, при которых аффинная система преобразуется к квазиканоническому виду, задает следующая теорема.
Теорема 2. Для того чтобы в некоторой окрестности точки ж0 для аффинной системы (1) существовали переменные, в которых она имеет квазиканонический вид (2), необходимо и достаточно, чтобы существовала гладкая функция <^(ж), удовлетворяющая в окрестности точки ж0 системе уравнений (7), и ранг матрицы Якоби отображения
г - Аг-1^(ж), г -1~г, (9)
в точке ж0 был равен г.
Если система в некоторой окрестности точки х° преобразуется к квазиканоническому виду (2), то согласно теореме 1 в этой окрестности определена функция <^(х), удовлетворяющая системе уравнений (7), причем соотношения (9) задают в этой окрестности невырожденную замену по части переменных г. Следовательно, ранг матрицы Якоби отображения (9) в точке х° был равен г.
Пусть в некоторой окрестности точки х° существует гладкая функция <^(х), удовлетворяющая в этой окрестности системе уравнений (7), а ранг матрицы Якоби отображения (9) в точке х° был равен г. Дополним множество функций (9) п — г гладкими в окрестности точки х° функциями пj = П (х), 3 = 1,п — г, так, чтобы ранг матрицы Якоби построенного отображения (¿,п) = Ф(х) в точке х° был равен п. Такие функции всегда существуют. Тогда в некоторой окрестности точки х° отображение (г, п) = Ф(х) будет задавать гладкую невырожденную замену переменных и в силу теоремы 1 в некоторой окрестности точки х° в этих переменных система (1) имеет квазиканонический вид (2).
Следствие 2. Если В(х°) = 0, выполнены условия теоремы 2 и в окрестности точки х° гладкие функции пj = П (х), 3 = 1,п — г, выбраны так, что выполнены условия (5) и ранг матрицы Якоби построенного отображения (г, п) = Ф(х) в точке х° равен п, то аффинная система (1) в окрестности этой точки преобразуется к специальному квазиканоническому виду (6).
Обозначим через ик (х) столбец координат векторного поля (-1)*а^В. Матрицу и(х) = (и°(х),..., ип-1(х)) называют [1] матрицей управляемости системы (1).
Пусть [5, 6]
7 (х) = ВАг-1^(х) = (-1)г-1ааА-1В^(ж), (10)
где функция <^(х) есть решение системы (7). Поскольку
¿г = АГ ^(ж)|я=^-1(г,п) + иВАг-1^(х)|ж=^-1(г,п), (11)
то функция 7(х) есть коэффициент при управлении в системе квазиканонического вида, записанный в переменных х.
Теорема 3. [9] Пусть иг (х) — матрица, состоящая из первых г столбцов матрицы управляемости и(х), <^'(х) — матрица Якоби замены по части переменных = Аг^(х), г = 1,г, где <^(х) —решение системы (7) линейных однородных дифференциальных уравнений первого порядка. Тогда
Ие%'(х)иг (х))| = |7 (х)|т. (12)
Следующая теорема содержит локальные условия существования преобразования аффинной системы (1) к квазиканоническому виду (2),
который обобщает известные результаты для систем, эквивалентных каноническому виду [3].
Теорема 4. Если при некотором г, 2 < г < п, для аффинной системы (1) существует такое решение <^(ж) € С ^(П) системы линейных однородных уравнений в частных производных (7), для которых соответствующая функция 7(ж) - ВАг-1^(ж) в точке ж0 пространства состояний системы удовлетворяет условию
т(ж°) - 0, (13)
то в некоторой окрестности точки ж0 :
— подматрица иг (ж) матрицы управляемости системы (1) имеет ранг г;
— матрица Якоби <^'(ж) функций г, - Аг-1 <^(ж), г - 1,г, имеет ранг г;
— аффинная система (1) преобразуется к регулярному квазиканоническому виду (2).
Пусть 7(ж0) - 0. Рассмотрим произведение матрицы Якоби <^'(ж) на подматрицу иг (ж) матрицы управляемости. По теореме 3 имеем |^'(ж)иг(ж)| - |7Г(ж)|. Поскольку 7(ж0) - 0, то определитель матрицы Ф' (ж)иг (ж) отличен от нуля в точке ж0 ив силу непрерывности его элементов также отличен от нуля в некоторой окрестности точки ж0. Следовательно, в этой окрестности матрица ^'(ж)иг(ж) имеет ранг, равный г, и ранг каждого сомножителя не меньше г. Поскольку матрица иг содержит ровно г столбцов, а матрица Якоби Ф'(ж) г строк, то в рассматриваемой окрестности ранги этих матриц постоянны и равны г.
Поскольку ранг матрицы Якоби <^'(ж) равен г в точке ж0, то по теореме об обратной функции [11] функциональные соотношения (8) в некоторой окрестности О точки г0 - ^(ж0) разрешимы относительно г переменных из числа ж1,..., жп. Не нарушая общности, примем, что
это ж1,..., жг:
ж^ - ж^ (г1,...,гг, жг+1,..., жп), 3 - 1,г. (14)
Положим
Пк - жг+(%, к - 1,п — г. (15)
Соотношения (14) и (15) в некоторой окрестности точки ж0 задают гладкую обратимую невырожденную замену переменных, и в новых переменных г, п аффинная система (1) запишется в квазиканоническом виде (2).
Поскольку согласно (11) функция 7(ж) определяет коэффициент при управлении в полученной системе квазиканонического вида, то в
некоторой окрестности точки x0 квазиканонический вид будет регулярным.
Замечание 1. Если для системы (1) в окрестности точки x0 существует регулярный квазиканонический вид, то B(x0) = 0.
Замечание 2. Если условия теоремы (4) выполнены при r = n, то аффинная система (1) преобразуется в некоторой окрестности точки x0 к регулярному каноническому виду [1, 3].
Замечание 3. Любая аффинная система может рассматриваться как система, записанная в квазиканоническом виде с индексом приводимости r = 1, если в качестве координатной функции zi выбрать любую координатную функцию x^, такую что соответствующее ей уравнение аффинной системы (1) содержит в правой части управление. При этом такой квазиканонический вид не обязательно будет регулярным в фиксированной точке x0.
Индекс приводимости не может быть больше n, поэтому для каждой аффинной системы (1) в окрестности точки x0 существует максимальный индекс приводимости к квазиканоническому виду rmax, 1 < rmax < n, и для преобразования системы к квазиканоническому виду можно выбрать любое r, 1 < r < rmax.
При фиксированном r для системы (1) эквивалентная ей система квазиканонического вида (2) не единственна. По крайней мере, имеется некоторая свобода выбора замены по части переменных
xr+1 , . . . , xn.
Если B(x0) = 0, то в окрестности точки x0 для системы (1) существует эквивалентная система регулярного квазиканонического вида с индексом приводимости r = 1 . Перейдем к определению максимального значения индекса r приводимости системы (1) к регулярному квазиканоническому виду, 1 < f < rmax.
Локальные условия существования регулярного квазиканонического вида задает следующая теорема.
Теорема 5. Пусть в некоторой окрестности точки x0 размерность инволютивного замыкания Fr-1 распределения
Fr-1 = span(B, adAB,..., ad^-2B),
r > 2, постоянна и равна m,r — 1 < m<n, и для векторного поля adA-1B имеет место adA-1B(x0) £ Fr-1(x0).
Тогда в некоторой окрестности точки x0 аффинная система (1) преобразуется к регулярному квазиканоническому виду (2) с индексом приводимости r.
Поскольку инволютивное замыкание Fr-1 распределения Fr-1 по условию теоремы регулярно (dim Fr-1 = m) и инволютивно в некоторой окрестности точки x0, то у этого распределения существует
к - п — т функционально независимых в окрестности точки ж0 первых интегралов [3, 11].
Среди них существует интеграл <^(ж), такой что аа^-1 В^(ж0) - 0. Действительно, если для всех первых интегралов распределения ^г-1 имеет место ааА-1В^(ж0) - 0, то в точке ж0 аа^-1В(ж0) € ^г-1(ж0), что противоречит последнему условию теоремы.
Указанная функция <^(ж) является решением системы уравнений (7). Поскольку при этом
7(ж0) - (—1)г-1аа^-1В^(ж0) - ВАг-1^(ж0) - 0,
то в силу теоремы 4 в некоторой окрестности точки ж0 система (1) преобразуется к квазиканоническому виду с индексом приводимости г.
Следствие 3. Пусть Г < п — максимальное число, для которого в окрестности точки ж0 выполнены условия теоремы 5. Тогда г — максимальный индекс приводимости системы (1) к регулярному квазиканоническому виду (2).
Алгоритм преобразования к регулярному квазиканоническому виду. Приведем основные этапы исследования проблемы преобразования системы (1) в окрестности заданной точки ж0 к регулярному квазиканоническому виду (2) с использованием системы компьютерной алгебры.
Проверяем условие В (ж0) - 0. Если это условие выполнено, то аффинную систему (1), как отмечалось выше, можно рассматривать как систему, записанную в регулярном квазиканоническом виде с индексом приводимости г - 1. Если условие не выполнено, прекращаем поиск.
Переходим к поиску индекса приводимости г > 1. Если указанное неравенство выполнено, то распределение - эрап(В) регулярно в некоторой окрестности точки ж0 и инволютивно, так как [В, В] - 0.
Далее для к - 2,3,... последовательно вычисляем векторное поле ааА-1В, проверяем регулярность распределения
^к - эрап(В, аа^В,..., аа^-1В)
в окрестности точки ж0, определяем размерность его инволютивного замыкания и регулярность. Процесс останавливается, как только либо распределение ^к не будет регулярным в окрестности точки ж0, либо инволютивное замыкание ^к не будет регулярным, либо размерность станет равной п — размерности системы (1).
Пусть остановка произошла при к - ктах. Возвращаемся к анализу распределения ^ктах-1. Это распределение инволютивно и регулярно.
Проверяем условие ааАтах-1В(ж0) € ^ктах-1. Если это условие выполнено, то найден максимальный индекс гг - ктах приводимости к
регулярному квазиканоническому виду. Если это условие не выполнено, переходим к рассмотрению распределения ^ктах-2.
Проверка будет продолжаться до тех пор, пока не найдется такой первый индекс ктах — 1 > т > 2, что будет выполнено условие а^т-1В(ж0) Ф ^т-1. Тогда Г = т. Если же и при т = 2 указанное неравенство не выполняется, то система приводится к регулярному квазиканоническому виду с г = 1.
Проверка регулярности или инволютивности заданного распределения сводится к исследованию ранга соответствующей ему функциональной матрицы, по столбцам которой записаны координаты векторных полей, порождающих это распределение, и может выполняться с использованием алгоритма, предложенного в [2].
Решая методами компьютерной алгебры соответствующую систему уравнений в частных производных, для найденного индекса приводимости г находим такой первый интеграл <^(х) распределения ^ который удовлетворяет условию а^4_ 1^(ж°) = 0.
Если в аналитическом виде его получить не удается, то аналитически преобразовать систему (1) к регулярному квазиканоническому виду с указанным индексом невозможно.
Пусть удалось найти <^(х) в аналитическом виде. Тогда для выбранной функции <^(х) с использованием системы компьютерной алгебры вычисляем функции Аг-1^(х), г = 1,г.
В режиме диалога задаем функции пг = Пг(х), г = 1,п — г, такие что ранг матрицы Якоби множества функций гг = Аг-1^(х), г = 1,г, П] = П](х), ] = 1, п — г, равен п в точке х0.
Для получения квазиканонического вида (2), в котором р(г, п) = 0, функции Пг следует выбрать из числа первых интегралов векторного поля В, но это осуществимо не всегда.
Для сформированной замены переменных (г,п) = Ф(х) ищем обратную замену х = Ф-1(г,п) (если она выражается аналитически) с использованием системы компьютерной алгебры.
В режиме диалога определяем (если это возможно) область, в которой построенная замена переменных (г, п) = Ф(х) является диффеоморфизмом, на основе анализа гладкости прямого и обратного отображений и других их свойств.
Оценкой сверху для искомой области может быть область, в которой не вырождена матрица Якоби отображения (г, п) = Ф(х).
Если прямая и обратная замены получены в аналитическом виде, то с использованием системы компьютерной алгебры получаем квазиканонический вид в переменных г, п.
Замечание 4. Для преобразования системы (1) к регулярному квазиканоническому виду с индексом приводимости г < г следует найти
первые интегралы распределения ^г-1 и среди них выбрать такой первый интеграл <^(х), который удовлетворяет условию а^—1^(х0) — 0. Если такую функцию <^(х) удалось найти, далее действуем согласно приведенному выше алгоритму.
Замечание 5. Неоднозначность регулярного квазиканонического вида при фиксированном индексе приводимости г < Г связана в том числе с тем, что могут существовать несколько первых интегралов распределения ^г-1, для которых аёА-1<£(х0) — 0. Пример. Рассмотрим динамическую систему
х1 — Х2 х 1, х 2 — 2х1х2 — 2x1 + и, А з — х1 — 2х1х2 + 2х^, хс 4 — х2 — х? + х1и
в окрестности точки х0 — (1, 0, 0, 0). Векторное поле В имеет вид
(16)
д
+ жг
д
дх2 дх4
и не обращается в нуль в точке х0. Следовательно, распределение регулярно и инволютивно в окрестности точки х0.
Перейдем к анализу распределения ^2 — эрап(В, ааАВ). Найдем векторное поле
adA B =
dxi
Поскольку ранг матрицы
д д ^ д , 2 ^ д — 2xi —— 2xi —--+ (x2 — Xj — 1)
дх2
дх3
дх4
(17)
(B(x), adAB(x)) =
/ 0 —1 \
1 —2xi
0 2xi
V Xi 2 X2 — Xi — 1
(18)
в окрестности точки х0 равен двум, распределение ^2 регулярно. Для анализа инволютивности распределения вычислим вектор-
д
ное поле Х1 — [В, ааАВ] — 2—— и сформируем соответствующую
дх4
функциональную матрицу
/0 —1 0 \
1 2х 0
(В(х), аа^В(х),Х1(х))— ^ 01х1 0 (19)
0 —1 0
1 —2xi 0
0 2xi 0
^ Xi X2 — xi — 1 2
Поскольку ранг матрицы (19) равен трем, распределение — — эрап(В, а^В) не является инволютивным.
Найдем инволютивное замыкание распределения ^2. Для этого вычислим коммутаторы векторных полей, входящих в распределение
эрап(В, аёАВ, Х^.
Получим [В, XI] = 0 и [аё^В, Х1] = 0. Таким образом, инволютивное замыкание распределения ^2 порождается векторными полями В, аёАВ и Х1. Ранг матрицы (19), а следовательно, и размерность равны трем.
Вычислим векторное поле
д
ааА В = —. (20)
дхз
Можно видеть, что ранг матрицы (В(ж), аёАВ(ж), аё^В(ж)) равен трем и поэтому распределение ^3 = эрап(В, В, аё^В(ж)) регулярно.
Перейдем к построению инволютивного замыкания распределения ^3. Поскольку ранг функциональной матрицы
/0 -1 0 0 \ 1 -2ж1 0 0 0 2ж1 0 1 ж1 ж - ж1 - 1 2 0 равен четырем, т.е. размерности системы, процесс поиска инволю-тивного и регулярного распределения, имеющего размерность г < п, завершен и искомым распределением является ^ 2.
Поскольку ранг матрицы (21) равен четырем и в самой точке ж0,
то
а^В(ж0) / ^(ж°).
Следовательно, у системы (16) существует регулярный квазиканонический вид с максимальным индексом приводимости Г = 3. В данном случае Г = Гтах.
Составим систему уравнений для нахождения первых интегралов распределения ^2:
+ ж^ = 0,
дж2 дж4
д^ о д^ ^ д<£ 1 , 2 1 ^ д^
(B(x), adAB(x),X1(x), adAB(x)) =
(21)
- 2ж^ + 2ж^ + (ж2 - ж2 - 1^ = 0, (22) дж1 дж2 дж3 дж4
2 ^ =°-
дж4
Поскольку размерность распределения ^2 равна трем, то система (22) [11] имеет 4 - 3 = 1 функционально независимое решение. Найдя его с использованием системы компьютерной алгебры, получим искомый
первый интеграл <^(х) — х1 + х3. Множество решений системы (22) имеет вид Л,(^(х)), где — произвольная гладкая функция скалярного аргумента.
Положим г1 — <^(х). Тогда ¿2 — ¿1 — х1, ¿3 — ¿2 — х2 — х^. Поскольку существование регулярного квазиканонического вида уже установлено, проверять ранг матрицы Якоби замены по части переменных г не нужно.
В качестве п выберем х4. Получим набор функций
2
¿1 — х2 + хз,
¿2 — X1, 2 (23)
¿3 — х2 — х2,
П — х4,
который задает замену переменных в окрестности точки х0. Обратная к (23) замена имеет следующий вид:
х1 — ¿2,
х2 — ¿3 + ^ (24)
хз — ¿1 — ¿2 ,
х4 — п.
Поскольку прямая и обратная замены задаются гладкими в М4 функциями, а отображение (24) взаимно-однозначно отображает М4 на К4, то это отображение является диффеоморфизмом Л4.
Записав систему в переменных (23), получим квазиканонический
вид
¿1 — ¿2,
¿2 — ^ (25)
¿3 — и, П — ¿3 + ¿2«.
Он определен всюду в Л4 и является всюду регулярным.
Выбор п, сделанный ранее, не является единственно возможным. Чтобы в последнее уравнение системы (25) не входило управление, требуется выбрать п из числа первых интегралов векторного поля В. Таких функционально независимых первых интегралов три: р — х1, д — х3 и 8 — —х1 х2 + х4.
Выбрав п — —х1х2 + х4, получим невырожденную замену переменных
¿1 — х2 + х3,
Z2 = Xi,
2
— X2 — Xi П = —x1X2 + X4.
Обратная к (26) замена имеет вид
ж1 = ¿2,
2
ж2 = ¿3 + ¿2, 2
ж3 = ¿1 - ¿2 ,
2
ж 4 = П + ¿2^3 + ¿2
Соответствующее отображение является диффеоморфизмом Л4. В переменных (26) система (16) имеет специальный квазиканонический вид
¿1 = ¿2,
¿2 = ¿^ (26) ¿3 = и,
П = -¿2 - 3¿з¿2 + ¿3,
поскольку в последнее уравнение управление не входит.
Заключение. Получены новые локальные условия существования квазиканонического (в том числе регулярного) вида для аффинных систем со скалярным управлением, а также приведен алгоритм преобразования к регулярному квазиканоническому виду, основанный на использовании систем компьютерной алгебры.
Работа выполнена при финансовой поддержке РФФИ (проекты 12-07-00329 и 11-01-733) и Министерства образования и науки Российской Федерации (соглашение № 14.B37.21.0370).
СПИСОК ЛИТЕРАТУРЫ
1. Жевнин А. А., Крищенко А. П. Управляемость нелинейных систем и синтез алгоритмов управления // Докл. АН СССР. - 1981. - T. 258, № 4. -С. 805-809.
2. К р и щ е н к о А. П. Преобразование многомерных аффинных управляемых систем // Управляемые нелинейные системы. - 1991. - № 2. - C. 5-14.
3. КраснощеченкоВ. И., Крищенко А. П. Нелинейные системы: геометрические методы анализа и синтеза. М.: Изд-во МГТУ им. Н.Э. Баумана, 2005. - 520 с.
4. BrunovskyP.A. A classification of linear controllable systems // Kybernetika.
- 1970. - Vol. 6. - P. 176-188.
5. Крищенко А. П. Преобразование нелинейных систем и стабилизация программных движений // Труды МВТУ им. Н.Э. Баумана. - 1988. - № 512. - С. 6987.
6. Крищенко А. П., Клинковский М. Г. Преобразование аффинных систем с управлением и задача стабилизации. // Дифференциальные уравнения.
- 1992. - № 1. - Т. 28. - С. 1945-1952.
7. B i r k J., Z e i t z M. Program for symbolic and rule-based analysis and design of nonlinear systems // Lecture Notes in Control and Information Sciences, SpringerVerlag, Berlin. - 1991. - No. 165. P. 115-128.
8.Fischer G. NONLINCON: symbolic analysis and design package for nonlinear control systems // Master's thesis. Eindhoven University of Technology. - 1994. - 176 p. [Электронный ресурс]. URL: http://alexandria.tue.nl/repository/books/639746.pdf (дата обращения: 20.05.2012).
9. I s i d o r i A. Nonlinear control systems. London: Springer-Verlag, 1995. - 587 p.
10. K h a l i l H. K. Nonlinear systems. 3rd edition. Upper Saddle River: Prentice-Hall, 2002. - 750 p.
11. Канатников А.Н.,КрищенкоА. П., Четвериков В. Н. Дифференциальное исчисление функций многих переменных. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2000. - 456 с.
12. Шевляков А. А. Вычисление ранга функциональной матрицы // http://technomag.edu.ru Наука и образование: электронное научное издание. -2011. № 10. [Электронныйресурс]. URL: http://technomag.edu.ru/doc/243762.html (дата обращения: 20.05.2012).
Статья поступила в редакцию 21.06.2012
Сергей Борисович Ткачев — д-р физ.-мат. наук, профессор кафедры "Математическое моделирование" МГТУ им. Н.Э. Баумана. Автор более 60 научных работ в области математической теории управления.
S.B. Tkachev — D. Sc. (Phys.-Math.), professor of "Mathematical Simulation" department of the Bauman Moscow State Technical University. Author of more than 60 publications in the field of mathematical control theory.
Андрей Анатольевич Шевляков — аспирант кафедры "Математическое моделирование" МГТУ им. Н.Э. Баумана. Автор двух научных работ в области теории управления.
A.A. Shevlyakov — post-graduate of "Mathematical Simulation" department of the Bauman Moscow State Technical University. Author of two publications in the field of control theory.