close

Вход

Забыли?

вход по аккаунту

?

Безударное приведение механических систем за конечное время в голономное программное многообразие в условиях неопределённости.

код для вставкиСкачать
УДК 531.31:62-56
Безударное приведение механических систем за конечное
время в голономное программное многообразие в условиях
неопределённости
И. А. Мухаметзянов
Кафедра теоретической механики
Российский университет дружбы народов
улица Миклухо-Маклая, 6, Москва, Россия, 117198
Строится алгоритм управления процессом безударного приведения в голономное программное многообразие фазового состояния механических систем любой конфигурации
за конечное время при произвольно действующих на них не управляющих активных сил
и ограниченных случайных возмущений.
Ключевые слова: управление, программное многообразие, алгоритм управления,
безударное приведение в голономное программное многообразие, конечное время.
1.
Постановка задачи
Рассмотрим механическую систему, движения которой описываются следующими уравнениями Лагранжа второго рода:
(︂
)︂
d 

−
=  + ′ ,
(1)
d  ˙

где  — кинетическая энергия системы вида
 =
1 
˙ (, )˙ + ˙ (, ) + 0 (, ),
2
(2)
 — -мерный вектор обобщённых координат, (, ) — ( × ) матрица, (, ) —
-мерный вектор, 0 (, ) — скалярная функция,  — -мерный вектор управляющих сил, ′ — -мерный вектор неуправляющих активных сил и случайных
возмущающих сил, ограниченных по величине.
Элементы вектора (, ) и матрицы (, ), а также функцию 0 (, ) и их
производные по  и  в области (, ) функционирования системы (1) будем считать ограниченными и непрерывными. Заметим также, что (, ), (, ), 0 (, ),
определяющие конфигурацию системы, кроме этих условий, не стеснены другими
ограничениями.
Пусть невозмущённое состояние системы (1) задано в виде ( − )-мерного
многообразия, образованного голономными программными связями:
(, ) = 0,
(3)
где  — -мерный вектор с непрерывными и линейно независимыми в области
(, ) элементами, непрерывно дифференцируемыми по  и  в этой области.
Заметим, что  6 .
Задача заключается в построении управляющей обобщённой силы  в виде
комбинации непрерывных и ступенчатых функций от  и ,
˙ обеспечивающей
приведение за конечное время фазового состояния системы (1) в многообразие
(3) при любых начальных условиях 0 , ˙0 , 0 , независимо от конкретного вида
(, ), (, ), 0 (, ), ′ .
Статья поступила в редакцию 20 марта 2012 г.
Работа выполнена при финансовой поддержке РФФИ, проект 10-01-00381-а.
Мухаметзянов И. А. Безударное приведение механических систем за ко‌ . . .
105
Отыскание вектора  в виде функции от  и ˙ объясняется тем, что их значения в невозмущённом состоянии (3) равны нулю, а при отклонениях от него
становятся отличными от нуля. Следовательно, вектор  может быть принят в
качестве меры отклонения от многообразия (3).
Заметим, что при размерности вектора , равной количеству степеней свободы
системы, решение поставленной задачи возможно по принципу декомпозиции [1],
развитому в работах [2–5]. Важной особенностью цели данной работы является
решение поставленной задачи при минимальной размерности вектора , равной
размерности  вектора  или при любой размерности  вектора , удовлетворяющей условию  6  6 .
2.
Алгоритм управления укороченной системой
Для решения задачи переходим от обобщённых координат 1 , 2 , . . . ,  к другим обобщённым координатам 1 , 2 , . . . ,  , являющимся элементами вектора ,
и координатам 1 , 2 , . . . , − , ортогональным к ним.
В силу ортогональности этих групп координат члены  матрицы (, , )
кинетической энергии в новых координатах будут равны нулю. Следовательно,
кинетическая энергия системы будет иметь следующую структуру:
где
 =  +  + 0 (, , ),
(4)
1 
˙  (, , )˙ + ˙   (, , ),
2
1
 = ˙  (, , )˙ + ˙  (, , ).
2
(5)
 =
Заметим, что матрица  является определённо положительной. При этом
уравнение (1) в новых координатах разбивается на две части:
(︂
)︂

d 
−
(6)
=  + ′ ,
d  ˙

(︂
)︂
d 

−
=  + ′ .
(7)
d  ˙

Систему (6) назовём укороченной системой.
Теперь переходим к преобразованиям, связанным лишь с укороченной системой, считая влияние системы (7) на систему (6) через элементы  и ,
˙ входящими
в неё через  (, , ), (, , ) и их производные по , , , возмущающимися факторами системы (6).
Введём вместо ˙ квазискорости 
˜˙

˜˙ = ˙ + ,
где  — вектор с элементами   ,  > 0 ( = 1, 2, . . . , ).
Введём обозначение 
˜ = .
Подставляя (8) в (5), получим
 = ˜ +  + ˜0 ,
(8)
106
Вестник РУДН. Серия Математика. Информатика. Физика. № 3. 2012. С. 104–114
где
˜ =
 =
1 ˙

˜  (, , )
˜˙ + 
˜˙  ˜ (, , ),
2
1 

˜  
˜ −   ,
2
1 ˙
˜  
˜˙ .
= 
2
1 
˙  (, , )˙ + ˙  (, , ),
2
˜ =  −  
˜,
(2)
˜
˜ = 0 + ˜0 ,
˜0 =
Уравнение (6) представим в квазикоординатах
(︂
)︂
d ˜
˜


−
=  + ′ +
−
˙
d  

˜


˜
˜
(9)
в силу того, что

˜
˜
=
=
,
 ˙
 ˙

˜˙
так как

˜˙
= 1.
 ˙
˜
Имеет место
=  
˜˙ + ˜ . Следовательно, уравнение (9) можно предста
˜˙
вить в виде
(2)
)︀
˜
d (︀
˜,
˜˙ −
˜˙ + ˜ 
 
=  + 
(10)
d

˜
где
(︂
)︂
0
˜

˜  = ′ + ′ +  − ˜ − ˜ ˙ − ˜ + ˜ .
˜ =
− ˜ , 

˜

˜


˜



˜
Умножая скалярно на 
˜˙ уравнение (10), левую часть полученного уравнения
представим в виде
(︃
)︃
(2)
(2)
(2)
)︀

d (︀ ˙ 






˜

˜

˜
¨˜  
˜˙ + 
˜˙
˜˙ − 

˜  
+ ˙
+
.
(11)
d

˜


Так как
(2)
(︀
)︀
(2)

˜˙   
˜˙ = 2˜ ,
(2)
(2)


˜
¨˜   

˜˙ + 
˜˙  ˜ + ˙ ˜ +

˜


=
d (︁ (2) )︁

,
d ˜
то из (11) следует
(2)
d˜
d
=
˜˙ 
(︂
)︂
 ˙
1 ˙  d ˙
˜ + 1 
 + 
˜˙ 

˜ − 
˜

˜.
2

˜
2
d
Вектор обобщённых сил управления выберем в виде
⃒
⃒
⃒ d ⃒
1
⃒
⃒,
 =  − 
˜˙ ,  > ⃒
2 d ⃒
где  — постоянная определённо положительная матрица,
(︀
)︀ (︀
)︀
 = − sign 
˜˙ 0 + 
˜˙  0 
˜˙ ,
(12)
(13)
(14)
Мухаметзянов И. А. Безударное приведение механических систем за ко‌ . . .
˙
107
˙
где 0 — постоянный вектор, 
˜ 0 
˜˙ — вектор с элементами 
˜ 0 
˜˙ , 0 — опре˙
делённо положительные матрицы, sign 
˜ — диагональная матрица с элементами
sign 
˜˙  ,  = 1, 2, . . . , . Если 0 , 0 , являющиеся элементами 0 и 0 , удовлетворяют условиям
⃒
⃒
⃒ ⃒
⃒  ⃒
1
⃒ ˜ ⃒
⃒
⃒,
(15)
0 > ⃒ ⃒ , 0 > max ⃒
2

˜ ⃒
то правая часть (12) становится
определённо отрицательной функцией по 
˜˙ . Сле⃒ ⃒
⃒
⃒
˙
˜ со временем будет убывать.
довательно, значение 
Покажем, что при выборе  в виде (14) при любых
конечных начальных
⃒ ⃒
значениях 
˜˙ (0) время обращения в нуль величины ⃒
˜˙ ⃒ будет конечным. Доказательству этого утверждения посвящается следующий раздел.
3.
Оценка времени приведения системы в терминальное
состояние
Так как матрица  положительно определена, то имеет место
1

∑︁
(2)

˜˙ 2 6 ˜ 6 2
=1

∑︁

˜˙ 2 ,
 = const > 0,
 = 1, 2.
(16)
=1
Поэтому из (12) следует неравенство
(2)
d˜
d
(︂ √︁
)︂
(2)
2
6 − 1 ˜ + 2 ˜ ,
(17)
так как правая часть (12), являющаяся определённо отрицательной функцией,
наибольшее значение которой может быть оценено правой частью (17), где 1 и
(2)
2 — положительные постоянные. Обозначая для простоты ˜ через  , из (17)
получим
d
(18)
√ )︁ 6 −d.
√ (︁
 1 + 2 
Имеет место
1
2
1
√ (︁
√ )︁ = √ − (︁
√ )︁ .
1 
 1 + 2 
1 1 + 2 
(19)
Интегрируя левую часть (18) с учётом (19) от 0 до 0, а правую часть от 0 до 1 ,
получим оценку времени приведения системы в состояние 
˜˙ = 0:
(︂
)︂
2
2 √︀
1 6
ln 1 +
0 ,
(20)
2
1
(2)
где 0 — начальное значение ˜ .
√ Заметим, что при 2 = 0, 1 ̸= 0 правая часть (20) имеет значение, равное
2 0
.
1
Начиная с момента времени 1 , переведём систему в «режим торможения»,
учитывая противоположность знаков  (1 ) и ˙  (1 ). Для этого элементы вектора
108
Вестник РУДН. Серия Математика. Информатика. Физика. № 3. 2012. С. 104–114
(8) заменим на
˙ 2 (1 )
 sign ˙  ( )

˜˙ 1 = ˙  ( ) − (
˙ 1) + 
2 ()
( = 1, 2, . . . , ),
(21)
где  =  − 1 . При этом значения  и ˙  обращаются в нуль одновременно в
моменты времени
2 | (1 )|
2 = 1 +
( = 1, 2, . . . , ).
(22)
|˙  (1 )|
Элементы вектора  в (13) можно представить в виде
(︀
)︀
(︀
)︀
 = − 0 + 
˜˙ 0 0 
˜˙ 0 sign 
˜˙ 0 − 
˜˙ 0  ,
где
(23)
(︀
)︀

˜˙ 0 = ˙  sign 
˜˙ 2 + 
˜˙ 1 sign 1 − sign 
˜˙ 2 .
Такое управление обеспечивает приведение системы без удара в терминальное
состояние (3) за конечный промежуток времени 2 = max 2 .
4.
Алгоритм управления исходной системой
Теперь необходимо определить вектор обобщённых сил управления  исходной
системой (1). С этой целью определим зависимость между  и построенной в п.
2 функцией  . Для этого определим сумму элементарных работ всех активных
сил управления
 =  ,
(24)
где  — вектор изохронных вариаций элементов .

Выделим из (24) элементарную работу 
 =  , совершаемую лишь при
вариациях
 = Ω,
(25)
⃦ ⃦
⃦  ⃦
⃦
вытекающих из (3), где Ω = ⃦
⃦  ⃦ — прямоугольная ( × ) матрица.
Из системы  уравнений (25) определим элементы вектора  в количестве 
через  элементов вектора . Для этого вектор  разложим на две составляющие: () — вектор, нормальный к многообразию (3), и () — вектор, касательный к (3). Первый из них ищем в виде () = Ω , где  — -мерный
искомый вектор.
Подставляя
 = () + ()
(26)
в (25), получим ΩΩ +Ω() = . Следовательно, имеем () = Ω (ΩΩ )−1 .
Подставляя в (24) значение (26), получим
 =  Ω (ΩΩ )−1  +  () .
Второй член в правой части этого выражения не зависит от . Следовательно,
частью суммы элементарных работ управляющих сил, совершаемых на элементарных перемещениях , вносящих вклад в вариацию , является
 
 −1

,
 =  Ω (ΩΩ )
откуда
 =  Ω (ΩΩ )−1 .
(27)
Мухаметзянов И. А. Безударное приведение механических систем за ко‌ . . .
109
Если вектор обобщённых сил управления исходной системой (1) задавать в виде  = 0 , где  — -мерный вектор управления,
˙ ) — матрица ( × ),
⃦  ⃦0 (, ,
⃦
⃦
удовлетворяющая в области  условию det 0 0 ̸= 0 и Ω0 ̸= 0, то при
подстановке  = 0  в (27) получим следующую систему  уравнений для определения  элементов вектора :
(ΩΩ )−1 Ω0  =  .
(28)
Заметим, что правая часть этого уравнения была определена в виде (13), где 
имеет вид (14) или (16).
Решение уравнения (28) относительно  можно представить в виде [6]:
(︀
)︀−1
 = Ω̄ Ω̄Ω̄
 +  ,
(29)
где Ω̄ = (ΩΩ )−1 Ω0 , det ‖Ω̄ Ω̄ ‖ ̸= 0,  — -мерный произвольно задаваемый
вектор, удовлетворяющий условию Ω̄ = 0, который можно представить в виде [6]:
[︁
(︀
)︀−1 ]︁
 =  − Ω̄ Ω̄Ω̄
Ω̄ 
˜,
где  — единичная матрица, 
˜ — произвольный вектор. Заметим, что при  = 
(︀
)︀−1
матрица Ω̄ является квадратной, причём  − Ω̄ Ω̄Ω̄
Ω̄ = 0. Следовательно, имеет место  ≡ 0. Отсюда следует, что минимальная размерность вектора
управления  может быть равна размерности  вектора  при  < . Как отмечалось в п.1, в этом заключается принципиальное преимущество предлагаемого
здесь метода управления от принципа декомпозиции [1] при задании невозмущённого состояния системы в виде ( − )-мерного многообразия (3).
Следует отметить также то, что в случае  > , полагая  = 0, в силу произвольности вектора  , получим вектор управления , имеющий минимальную
евклидову норму, в виде
(︀
)︀−1 (︀
)︀
 = Ω̄ Ω̄Ω̄
 − 
˜˙ ,
(30)
где  — ступенчатая функция (14).
В частном случае  =  = 1 матрицы 0 и Ω становятся -мерными векторамистолбцами, а , ,  , 0 ,  — скалярными величинами. При этом из (30) получим скалярное управление
=
)︀
Ω2 (︀
 −  
˜˙ ,

(31)
где  — скалярное произведение векторов 0 и Ω ,  — выражается в виде (14).
5. Управление процессом приведения преследующего тела
в заданную ориентацию при пропорциональной навигации
5.1.
Постановка задачи
Рассмотрим твёрдое тело, жёстко связанное с подвижной системой координат
˜1 построим так, чтобы центр масс
 . Главный вектор управляющих сил 
тела двигался по принципу пропорциональной навигации [7] при погоне за преследуемой точкой  при её произвольном движении в пространстве. При этом
˜2 должен быть таким, чтобы ось тела  из
вектор управляющих моментов 
¯ 3 ) > 0, где 
¯3 —
любого начального положения, удовлетворяющего условию (¯ · 
орт оси , была приведена в положение, совпадающее с вектором скорости ¯
110
Вестник РУДН. Серия Математика. Информатика. Физика. № 3. 2012. С. 104–114
центра масс тела, за конечное время. Следовательно, программная ориентация
тела в этом случае может быть задана выражениями
¯  ¯ = 0,


¯  = ¯
ℓ
( = 1, 2),
(32)
¯ 1, 
¯ 2 — орты осей  и  ; 
где 
¯  — вектор угловой скорости вращения
вектора ¯; 
¯ ℓ — вектор угловой скорости вращения линии визирования ;  —
заданный положительный коэффициент пропорциональности.
Заметим, что при выполнении условий (32) вектор ¯ будет направлен по оси
¯ 3 . Следовательно, при этом вектор 
 тела с ортом 
¯  будет равен сумме ком¯1 и 
¯ 2 . При таком
понентов вектора угловой скорости тела на оси с ортами 
принципе навигации угловая скорость 
¯ ℓ линии визирования  считается доступной измерению в каждый момент времени .
5.2.
Построение главного вектора управляющих сил
¯1 в правой части
Определим выражение главного вектора управляющих сил 
дифференциального уравнения движения центра масс  тела
¯1 ,
¯˙ = ¯ + 
(33)
так, чтобы движение центра масс тела происходило по принципу пропорциональной навигации

¯  = ¯
ℓ ,
(34)
где  — масса тела,  — главный вектор внешних неуправляющих сил, действующих на тело.
Заметим, что вращение вектора скорости ¯ осуществляется благодаря действию силы ¯ , направленной по главной нормали траектории движения центра
масс тела. Определим эту силу ¯ , исходя из условий (34) и ¯ ⊥¯ . Вектор 
¯
2
¯ × ¯



имеет величину  =
и направлен по орту ¯ =
, где  =
,  —

 · 

радиус кривизны траектории точки . Следовательно, имеет место
(︀
)︀
¯ × ¯

¯ =
.
 2
Скалярно умножая (34) на 
¯ ℓ , получим
)︀
1 (︀ ¯
 × ¯ 
¯ ℓ = ℓ2 .
2

Отсюда
(︀
)︀
¯ · 
¯ ℓ × ¯ =  2 ℓ2 .
(35)
¯ ℓ × ¯ , обеспечивающем
Если искать решение этого уравнения в виде ¯ =  
коллинеарность векторов ¯ и (¯
 × ¯ ), то величина силы  будет минимальна.
¯
Подставляя  в (35), определим
(︀
)︀
 2 ℓ2
 = (︀
)︀2 ,

¯ ℓ × ¯
а затем
)︀
 2 ℓ2 (︀
¯ = (︀
¯ ℓ × ¯ .
)︀2 

¯ ℓ × ¯
(36)
Мухаметзянов И. А. Безударное приведение механических систем за ко‌ . . .
111
Покажем, что при движении преследующего тела со скоростью, намного превышающей скорость цели, величина 
¯ ℓ × ¯ при ℓ =
̸ 0, не может быть равна 0.
Для этого заметим, что принцип пропорциональной навигации призван осуществлять равенство компонентов скорости ¯ и скорости цели ¯ , перпендикулярных
к линии визирования , с целью обеспечения поступательности движения .
При этом максимальное значение синуса угла  отклонения вектора ¯ от линии

 не превышает значения sin max =
. А угол между векторами 
¯ ℓ и  яв
ляется прямым. Следовательно, случай max = 2 возможен лишь тогда, когда
 =  , при котором сближение с целью не возможно. Таким образом, в силу
того, что  ≫  при ℓ ̸= 0, случай (¯
ℓ × ¯ ) исключается.
¯1 = ¯ − ¯ , где
Итак, главный вектор управляющих сил выражается в виде 
¯ имеет вид (36), ¯ — компонента ¯ на главную нормаль траектории точки .
Отсюда напрашивается вопрос о том, что случится, если среди составляющих силы  имеются случайные недоступные измерению возмущающие силы. На
этот вопрос возможен такой ответ. Составляющие этих возмущающих сил, направленных по главной нормали траектории точки , возмущают силу ¯ , что
равносильно возмущению коэффициента пропорциональности . Следовательно,
при достаточно больших значениях коэффициента  эти возмущения на процесс
преследующего движения не оказывают существенного влияния. Что касается составляющих возмущающих сил, направленных по касательной к траектории точки , то они приводят лишь к изменению скорости ¯ точки , что не оказывает
влияния на угловую скорость 
¯  , т.е. на процесс преследующего движения. Таким образом, эти возмущения могут оказывать некоторое влияние лишь на время
встречи тела с целью .
5.3.
Построение главного момента управляющих сил
Известно, что вращательное движение тела вокруг центра масс в осях подвижной системы координат   описывается дифференциальным уравнением

˜˙ 0 = ( 
¯ × 
¯0 +  +  ) ,
(37)
где  — тензор инерции тела в точке ; 
¯ 0 (, , ) — мгновенная угловая скорость
тела; , ,  — проекции 
¯ 0 на оси ,  , ;  — главный момент относительно точки  внешних неуправляющих сил, действующих на тело;  — главный
момент управляющих сил относительно . Оси системы координат   будем
считать главными центральными осями инерции тела.
Вектор  определим из первого условия (32), обозначая
¯ 1 · ¯ ,
1 = 
¯ 2 · ¯ .
2 = 
(38)
Введём квазискорости
(︀
)︀
¯  · ¯ ,

˜˙  = ˙  +  
 > 0
( = 1, 2),
(39)
где
)︁
(︀
)︀ (︁
¯  × ¯ + 
¯  · ¯˙ .
˙  = 
¯0 · 
(40)
Из (39) следует
(︀
)︀
¯  × ¯ = 

¯0 · 
˜˙  + 
(︁
)︁
(︀
)︀
¯  · ¯ ( = 1, 2).
¯  · ¯˙ −  
где  = − 
( = 1, 2),
(41)
112
Вестник РУДН. Серия Математика. Информатика. Физика. № 3. 2012. С. 104–114
Решение системы (41) представим в виде
Ω¯
0 = 
˜˙ + ,
(42)
где матрица Ω определяется из
)︀
(︀
)︀
(︀
)︀
(︀
¯3 = 
¯ 1 × ¯ · 
¯1 +  
¯ 1 × ¯ · 
¯2 +  
¯ 1 × ¯ · 
 
˜˙ 1 + 1 ,
(︀
)︀
(︀
)︀
(︀
)︀
¯ 2 × ¯ · 
¯1 +  
¯ 2 × ¯ · 
¯2 +  
¯ 2 × ¯ · 
¯3 = 
 
˜˙ 2 + 2
в виде
)︀
(︀
)︀
(︀
)︀
⃦(︀ ¯
¯2
¯ 1 × ¯ · 
¯ 3⃦
¯1
¯ 1 × ¯ · 
⃦ 1 × ¯ · 
⃦


⃦
⃦.
)︀
(︀
)︀
(︀
)︀
Ω = ⃦(︀ ¯
(43)
¯
¯
¯
¯
¯
¯
¯
¯
2 ×  · 1
2 ×  · 2
2 ×  · 3 ⃦
)︀
(︀
)︀
(︀
)︀
(︀
)︀ (︀
(︀
)︀
¯ 1 × ¯ · 
¯ 1 = 0, 
¯ 1 × ¯ · 
¯ 2 = − ¯ · 
¯3 , 
¯ 1 × ¯ · 
¯ 3 = ¯ · 
¯2 ,
Так
как

(︀
)︀
(︀
)︀ (︀
)︀
(︀
)︀
(︀
)︀
¯ 1 , то (43)
¯ 2 × ¯ · 
¯ 1 = ¯ · 
¯3 , 
¯ 2 × ¯ · 
¯2 = 0 
¯ 2 × ¯ · 
¯ 3 = − ¯ · 

представляется в более компактной форме
(︀
)︀ (︀
)︀
⃦
¯ ·
¯3
¯ ·
¯2 ⃦
⃦
⃦
0
−


)︀
(︀
)︀⃦
Ω=⃦
(44)
⃦(︀¯ · 
¯3
¯1 ⃦ .
0
− ¯ · 
Решение (42) ищем в виде

¯ 0 = Ω ,
(45)
где  — вектор  = (1 , 2 ).
Подставляя (45) в (42), получим
ΩΩ  = 
˜˙ + .
Отсюда
(︀
)︀
 = (ΩΩ )−1 
˜˙ +  .
(46)
Подставляя (46) в (45), получим
(︀
)︀

¯ 0 = Ω (ΩΩ )−1 
˜˙ +  .
Это выражение представим в виде
(︀
)︀

¯0 =  
˜˙ +  ,
(47)
где  = Ω (ΩΩ )−1 .
Теперь определим кинетическую энергию  =
(47). Получим
1 

¯ 
¯ 0 , подставляя 
¯ 0 вида
2 0
)︀
(︀
)︀
1 (︀ ˙ 

˜ +      
˜˙ +  .
2
Следовательно,  имеет следующую структуру:
 =
 = ˜2 + ˜1 + ˜0 ,
(48)
(49)
1 ˙
1

˜  
˜˙ , ˜1 = ˜ 
˜˙ , ˜0 =    ,     , ˜ =    . Теперь
2 (︁
2
)︁
(1)
(2)
вектор  =  ,  в (13) может быть построен с элементами (22).
где ˜2 =
Для выражения вектора управления ( ,  ,  ) в правой части (37) через
 вариации  
˜  ( = 1, 2) представим через вариации  
¯ 0 (, , ). Из (47)
Мухаметзянов И. А. Безударное приведение механических систем за ко‌ . . .
113
следует  
¯0 =  
˜ . Сумму элементарных работ управляющих сил  ,  , 
можно представить в виде
 = 
˜  
¯0 = 
˜   
˜,
где 
˜ =  −1  .
Ту же сумму можно записать в виде  =   
˜ . Следовательно, имеет место



˜  =  . Отсюда
 =   
˜.
(50)
Вектор 
˜ ищем в виде
˜

˜ =  .
(51)
Подставляя (51) в (50), получим
˜ =  .
 
˜ = (   )−1  . Следовательно, из (51) определим выражение 
˜ через
Отсюда 
 в виде 
˜ =  (   )−1  .
Заметим, что в случае равенства нулю одного из компонентов вектора  матрица (44) становится квадратной. Например, при  = 0 имеет место
(︀
)︀
⃦
¯ ·
¯3 ⃦
⃦
⃦
0
−

⃦.
)︀
Ω=⃦
⃦(︀¯ · 
⃦
¯3
0
При этом управление  будет двумерным вектором  = ( ,  ).
В заключение отметим, что время приведения оси  тела в положение, направленное по ¯ , не превышает максимального значения (21).
Литература
1. Пятницкий Е. С. Принцип декомпозиции в управлении механическими системами // Доклады АН СССР. — 1988. — Т. 300, № 2. — С. 300–303.
[Pyatnickiyj E. S. Princip dekompozicii v upravlenii mekhanicheskimi sistemami //
Dokladih AN SSSR. — 1988. — T. 300, No 2. — S. 300–303. ]
2. Матюхин В. И. Универсальные законы управления механическими системами. — М.: МАКС Пресс, 2001. — 249 с. [Matyukhin V. I. Universaljnihe zakonih
upravleniya mekhanicheskimi sistemami. — M.: MAKS Press, 2001. — 249 s. ]
3. Матюхин В. И. // ДАН. — 2009. — Т. 427, № 1. — С. 44–47. [Matyukhin V. I. //
DAN. — 2009. — T. 427, No 1. — S. 44–47. ]
4. Ананьевский И. М. Непрерывное управление по обратной связи возмущёнными
механичскими системами // ПММ. — 2003. — Т. 67, вып. 2. — С. 163–178.
[Ananjevskiyj I. M. Neprerihvnoe upravlenie po obratnoyj svyazi vozmuthyonnihmi
mekhanichskimi sistemami // PMM. — 2003. — T. 67, вып. 2. — S. 163–178. ]
5. Ананьевский И. М. Синтез непрерывного управления механичской системой с
неизвестной матрицей инерции // Изв. РАН. Теория и системы управления. —
2006. — № 3. — С. 24–35. [Ananjevskiyj I. M. Sintez neprerihvnogo upravleniya
mekhanichskoyj sistemoyj s neizvestnoyj matriceyj inercii // Izv. RAN. Teoriya i
sistemih upravleniya. — 2006. — No 3. — S. 24–35. ]
6. Мухаметзянов И. А. Построение уравнений программных движений // Автоматика и телемеханика. — 1972. — № 10. — С. 16–23. [Mukhametzyanov I. A.
Postroenie uravneniyj programmnihkh dvizheniyj // Avtomatika i telemekhanika. —
1972. — No 10. — S. 16–23. ]
114
Вестник РУДН. Серия Математика. Информатика. Физика. № 3. 2012. С. 104–114
7. Кан В. Л., Кельзон А. С. Теория пропорциональной навигации. Л. — Судостроение, 1965. — 423 с. [Kan V. L., Keljzon A. S. Teoriya proporcionaljnoyj navigacii.
L. — Sudostroenie, 1965. — 423 s. ]
UDC 531.31:62-56
Nonimpact Reducing Mechanical Systems to Holonomic
Programmed Set During Finite Time Under Indeterminancy
Conditions
I. A. Mukhametzyanov
Department of Theoretical Mechanics
Peoples’ Friendship University of Russia
6, Miklukho–Maklaya str., Moscow, Russia, 117198
The procedure of constructing the controls algorithm of the nonimpact reduction of the
mechanics systems to holonomic programmed set in finite time under indeterminancy is
proposed.
Key words and phrases: control, programmed set, controls algorithm, nonimpact
reduction of the holonomic set, finite time.
1/--страниц
Пожаловаться на содержимое документа