24.07.2014 Views

Полная версия - Самарский государственный аэрокосмический ...

Полная версия - Самарский государственный аэрокосмический ...

Полная версия - Самарский государственный аэрокосмический ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!