ÐÐ¾Ð»Ð½Ð°Ñ Ð²ÐµÑÑÐ¸Ñ - СамаÑÑкий гоÑÑдаÑÑÑвеннÑй аÑÑокоÑмиÑеÑкий ...
ÐÐ¾Ð»Ð½Ð°Ñ Ð²ÐµÑÑÐ¸Ñ - СамаÑÑкий гоÑÑдаÑÑÑвеннÑй аÑÑокоÑмиÑеÑкий ...
ÐÐ¾Ð»Ð½Ð°Ñ Ð²ÐµÑÑÐ¸Ñ - СамаÑÑкий гоÑÑдаÑÑÑвеннÑй аÑÑокоÑмиÑеÑкий ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Авиационная и ракетно-космическая техника<br />
быстрой переменной (выбор структуры управления<br />
на витке) и определение законов<br />
изменения параметров этой программы от<br />
витка к витку.<br />
Критерий оптимальности представим в<br />
следующем виде:<br />
T<br />
0<br />
0<br />
( , )<br />
I = a∫ F xu dt . (1.13)<br />
Перейдем в системе (1.12) от времени t<br />
к быстрой переменной ϕ :<br />
dx<br />
dϕ<br />
dτ<br />
dϕ<br />
= aXx<br />
= a⋅ω<br />
−1<br />
( , ϕ,<br />
u()<br />
t) ⋅ω<br />
( x,<br />
ϕ)<br />
−1<br />
( x,<br />
ϕ).<br />
,<br />
(1.14)<br />
Здесь τ - так называемое «медленное» время,<br />
τ = a⋅ t. В соответствии с принципом<br />
максимума Понтрягина введем вектор сопряженных<br />
переменных Ψ и запишем Гамильтониан<br />
системы (1.14):<br />
H = ψ<br />
+ ψ<br />
Г<br />
T<br />
⋅aω<br />
−1<br />
( aω<br />
( x, ϕ ) ⋅ X( x, ϕ ,u( t)<br />
))<br />
+<br />
−1<br />
−1<br />
( x, ϕ ) − aω<br />
( x, ϕ ) ⋅ F ( x, ϕ ,u) =<br />
( x, ϕ,<br />
ψ u)<br />
= aF , . (1.15)<br />
Определим локально-оптимальное управление<br />
u% из условия максимума Гамильтониана<br />
на отрезке ϕ∈<br />
[ 0,2π]<br />
0<br />
. Проведем затем<br />
процедуру усреднения исходной и сопряженной<br />
систем по быстро меняющейся переменной<br />
ϕ. Усредненная система уравнений<br />
будет иметь вид<br />
∧<br />
2π<br />
dx<br />
a<br />
∧ ∧ ∧<br />
−1<br />
⎛<br />
⎞<br />
= ω ( , ϕ)<br />
⋅ , ϕ, ⎛<br />
x, ψϕ ,<br />
⎞<br />
⎜ ⎟ dϕ,<br />
dϕ<br />
2π<br />
∫ x X⎜x u%<br />
⎟<br />
⎝ ⎝ ⎠⎠<br />
∧<br />
dψ<br />
dϕ<br />
dτ<br />
dϕ<br />
0<br />
2π<br />
∧ ∧ ∧<br />
a ⎛<br />
⎞<br />
= Fx<br />
, yu ,<br />
⎛<br />
x, ψϕ ,<br />
⎞<br />
⎜ ⎟ dϕ,<br />
2π<br />
∫ ⎜x<br />
% ⎟<br />
⎝ ⎝ ⎠⎠<br />
0<br />
2π<br />
a<br />
∧<br />
−1<br />
= ω<br />
⎛<br />
x, ϕ<br />
⎞<br />
⎜ ⎟dϕ,<br />
2π<br />
∫<br />
⎝ ⎠<br />
0<br />
а усредненный критерий оптимальности<br />
(1.16)<br />
T 2π<br />
a<br />
∧ ∧ ∧<br />
⎛<br />
=<br />
0<br />
,<br />
⎛<br />
, ψϕ ,<br />
⎞<br />
2π<br />
⎜ ⎜ ⎟<br />
0 0 ⎝ ⎝ ⎠<br />
J ∫∫ F xu% ⎞<br />
x ⎟dϕ<br />
dt<br />
⎠<br />
. (1.17)<br />
Интегралы в правых частях системы<br />
(1.16) образуют совокупность так называемых<br />
«усредняющих интегралов»:<br />
2π<br />
0<br />
()<br />
∫ , ( ) T<br />
j<br />
Φ= ⋅ dϕ<br />
Φ= Φ , j = 1,2n+ 1. (1.18)<br />
После усреднения правые части системы<br />
(1.16) не содержат циклической переменной<br />
ϕ. Поэтому модель «медленной» эволюции<br />
вектора состояния и вектора сопряженных<br />
переменных может быть представлена в<br />
виде системы интегро-дифференциальных<br />
уравнений с «медленным» временем в качестве<br />
независимой переменной:<br />
∧<br />
dq<br />
dτ<br />
j<br />
=<br />
Φ<br />
∧<br />
Φ<br />
j⎜z<br />
⎟<br />
2n+<br />
1<br />
j = 1,2n+ 1,<br />
⎛ ⎞<br />
∧ ∧ ∧<br />
T<br />
⎝ ⎠<br />
∧<br />
, где q =<br />
⎛<br />
xψ ,<br />
⎞<br />
⎜ ⎟<br />
⎛<br />
z<br />
⎞ ⎝ ⎠ ,<br />
⎜ ⎟<br />
⎝ ⎠<br />
2π<br />
∧<br />
−1<br />
⎛<br />
2n 1<br />
x,<br />
⎞<br />
+<br />
d<br />
0 ⎝ ⎠<br />
Φ = ∫ ω ⎜ ϕ⎟<br />
ϕ . (1.19)<br />
Исходная оптимизационная задача сводится,<br />
таким образом, к решению краевой<br />
задачи для системы (1.16, 1.17). Однако вычисление<br />
усредняющих интегралов (1.18)<br />
представляет самостоятельную проблему и<br />
требует разработки специальных процедур.<br />
1.4. Приближенный метод решения<br />
динамической задачи (метод разбиения)<br />
Пусть движение КА описывается системой<br />
дифференциальных уравнений (1.2).<br />
Откажемся от получения универсального решения<br />
для всего пространства переменных<br />
и поставим цель определить ряд упрощенных<br />
решений для каждой отдельной выделенной<br />
области X пространства состояний. Разобьем<br />
допустимую область фазового пространства<br />
переменных на т подобластей таких, что<br />
X<br />
⊆ X1 ∪ X<br />
2<br />
∪ X<br />
3<br />
∪ ... ∪ X m<br />
. Заменим исходную<br />
двухточечную краевую задачу на многоэтапную<br />
последовательность переходов:<br />
x<br />
0<br />
→ x → x<br />
1<br />
2<br />
→ ... → x<br />
m−1<br />
→ x<br />
К<br />
,<br />
41