Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ОиАС-1

.pdf
Скачиваний:
29
Добавлен:
26.03.2016
Размер:
1.02 Mб
Скачать

Поэтому реальное движение описывается функциями:

x (t)=x

*(t) + δx (t);

u (t) = u

*(t) + δu (t) (i= 1… n; k=1… m)

(1.2.1)

i

i

i

k

k

k

 

отклонения (возмущения)

отклонения реального

фактического движения от

управления от

программного

программного

 

 

Числа δxi(t0) (i= 1… n) – достаточно малые, но неизвестные числа, являющиеся случайными погрешностями при реализации заданных начальных условий (1.1.2).

Об этих погрешностях обычно известно лишь, что они удовлетворяют неравенству

n

(t0 )≤ ε 2

 

∑δxi2

(1.2.2)

i=1

 

 

известное число

Нетрудно получить уравнения (уравнения возмущенного движения),

описывающие отклонения фактического движения от программного

(невозмущенного) движения.

21

Принимая во внимание, что функции (1.2.1) удовлетворяют (1.1.1), и вычитая из уравнений

x&i* (t )+δx&i (t ) = ϕi [x1* (t )+δx1 (t ),K, xn* (t )+δxn (t ), u1* (t )+δu1 (t ),K, um* (t )+δum (t ),t ]

тождества

x&* (t ) = ϕ

[x* (t ),K, x* (t ),u* (t ),K,u*

(t ),t] (i =

 

)

1, n

i

i

1

n

1

m

 

 

 

получим уравнения возмущенного движения

δx& (t ) = δϕi x1 (t ),Kxn (t ),δu1 (t ),Kum (t ),t ] (i =

 

)

(1.2.3)

1, n

δϕi x1 (t ),Kxn (t ),δu1 (t ),Kum (t ),t]=

i [x1* (t )+δx1 (t ),K, xn* (t )+δxn (t ),u1* (t )+δu1 (t ),K,um* (t )+δum (t ),t]− −ϕi [x1* (t ),K, xn* (t ),u1* (t ),K,um* (t ),t]

22

Если функции δφi(t) (i= 1… n) разложить в ряд Тейлора в окрестности точки

 

 

 

 

 

x *, ..., x

*, u *, ..., u *,

 

 

 

 

 

 

1

n

1

m

 

 

 

 

то уравнения (1.2.3) примут вид

 

 

 

 

 

 

 

n

 

 

m

 

 

 

 

 

 

 

 

(t x j

 

(t uk

+ oi

x1 ,Kxn u1 ,Kum ,t ) (i =1, n)

 

δx&i (t ) = ∑aij

+

bik

(1.2.4)

j=1

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

*

 

 

*

 

 

 

 

 

a

 

(t )= ∂ϕi

b (t ) =

∂ϕi

 

 

 

 

 

 

ij

 

x j

 

ik

uk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

символ |* означает, что частные производные вычисляются в точке xi = xi*, uk = uk*

(i= 1… n, k=1… m);

oix1, ..., δxn, δu1, ..., δum) (i= 1… n) – функции, разложение которых в ряд Тейлора начинается с членов второго порядка малости.

Отбрасывая в (1.2.4) нелинейные члены, получим уравнения первого приближения

n

 

m

 

 

 

 

 

δx&i (t ) = ∑aij

(t x j

+ ∑bik

(t uk

(i =1, n)

(1.2.5)

j=1

 

k =1

 

 

 

 

 

23

Понятие об оптимальном стабилизирующем управлении

24

Решения уравнения (1.2.4)

n

 

m

 

 

 

 

 

(t x j

(t uk

+ oi x1 ,K,δxn ,δu1 ,K,δum ,t )

(i = 1, n)

δx&i (t ) = ∑ aij

+ ∑ bik

j=1

 

k =1

 

 

 

 

 

при начальных условиях из множества (1.2.2)

n

(t

) ≤ ε 2

∑δx2

i

0

 

i=1

 

 

описывают отклонения реального движения от программного в каждый момент времени. Для количественной характеристики этих отклонений часто используют значение интеграла

t1

n

δxi2 )dt

 

J = ∫ (∑ qii

(1.2.6)

t0

i=1

 

 

в котором qii (i= 1… n) – положительные числа.

Интеграл (1.2.6) представляет собой взвешенную с помощью коэффициентов qii (i= 1… n) сумму площадей, ограниченных квадратом отклонений истинного движения от программного по каждой переменной состояния. Он характеризует «расстояние» между реальным движением и программным и является «мерой» близости этих движений.

Используем δuk(t) (k=1… m) для сближения этих движений, тогда δuk(t) (k=1… m)

называются стабилизирующими управлениями.

Таким образом, результирующие управления

uk(t) = uk*(t) + δuk(t) (k=1… m)

состоят из программных и стабилизирующих управлений. Подставляя это выражение в (1.1.5), получим ограничения на стабилизирующее управление:

 

u

* u

*(t) ≤ δu (t) ≤ u *

u *(t).

(1.2.7)

 

k

k

k

k

k

 

Обычно |u

*(t)| ≥ |δu (t)| (k=1…

m).

 

 

 

 

k

k

 

 

 

 

 

Это объясняется тем, что программное управление обеспечивает основное (программное) движение системы, а стабилизирующее управление лишь «парирует» малые отклонения от программного движения, обеспечивая, если t1 →∞, устойчивость (отсюда термин «стабилизирующее управление») и требуемую точность осуществления программного движения.

25

В связи с этим часто вместо ограничений (1.2.7), определяющих допустимый «расход» стабилизирующего управления в каждый момент времени, накладывают на стабилизирующие управления интегральные ограничения (ограничения на «энергию»)

t1

 

 

 

 

 

 

 

 

 

 

 

(t )dt

*

(k =1, m)

 

 

 

2

(1.2.8)

∫δuk

Juk

t0

 

 

 

 

 

 

 

 

 

 

 

Для учета ограничений (1.2.8) будем вместо (1.2.6) рассматривать критерии

качества стабилизации

 

 

 

 

 

 

 

 

 

 

 

J = t1

(n

q

δx2 + m γ

kk

δu 2 )dt

(1.2.9)

∫ ∑

ii

i

 

k

t0

i=1

 

 

k =1

 

 

 

 

 

 

 

где числа γkk (k=1… m) определяются значениями

 

 

 

 

 

*

 

(k =1, m)

Juk

Стабилизирующее управление предназначено для минимизации интеграла (1.2.9). Кроме того, если t1 →∞, то для существования этого интеграла стабилизирующее управление должно обеспечивать асимптотическую устойчивость нулевого решения системы (1.2.4).

26

Если отыскивать стабилизирующее управление как явную функцию времени (по аналогии с программным управлением), то для каждого начального условия из множества (1.2.2) получим управления

δuk[t, δx1(t0),…, δxn(t0)] (k=l… m),

для реализации которых необходимо измерять переменные состояния в момент времени t=t0, так как числа δxi(t0) (i=l… n) неизвестны. Кроме того, функции

δuk[t, δx1(t0),…, δxn(t0)]

будут различными для каждого набора δxi(t0) (i=l… n) из множества (1.2.2).

В связи с этим естественно отыскивать стабилизирующее управление не как

явную функцию времени, а как функцию переменных состояния

 

δuk(t)=rkx1(t),…, δxn(t), t] (k=1… m).

(1.2.10)

Заметим, что вид этих функций не зависит от начальных условий из множества (1.2.2). Поясним это обстоятельство подробнее.

27

Допустим, что найдено управление

δuk0(t)=rk0x1(t),…, δxn(t), t] (k=1… m),

при котором (1.2.9)

J = t1

(n

q

δx2

+ m

γ

kk

δu 2 )dt

∫ ∑

ii

i

 

k

t0

i=1

 

 

k =1

 

 

 

принимает наименьшее значение на движениях системы (1.2.4),

n

 

m

 

 

 

 

 

(t x j

(t uk

+ oi x1 ,Kxn u1 ,Kum ,t )

(i =1, n)

δx&i (t ) = ∑aij

+ ∑bik

j=1

 

k =1

 

 

 

 

 

и пусть в начале движения оптимальной системы реализовались начальные условия δxi*(t0) (i=l… n).

Подставляя в (1.2.4) δuk0(t)=rk0 (k=l… m), решим это уравнение и, подставляя его решения в (1.2.9), вычислим значение интеграла.

Получим число J(t0, t1, δx1*(t0), …, δxn*(t0)). Оно является наименьшим из всех значений интеграла (1.2.9) при δuk0(t) ≠ rk0 (k=l… m).

28

Допустим теперь, что в начале движения системы (1.2.4) реализовались начальные условия

δxi**(t0) ≠ δxi*(t0) (i=l… n),

тогда при δuk0(t)=rk0 (k=l… m) получим другое значение (1.2.9), определяемое как

J(t0, t1, δx1**(t0), …, δxn**(t0)).

Это число опять должно быть наименьшим по сравнению со значениями интеграла на траекториях системы (1.2.4) при uk(t) ≠ rk0 (k=l… m) и начальных условиях δxi**(t0) (i=l… n).

Оптимальное стабилизирующее управление − функция переменных состояния и времени, при которых на движениях системы (1.2.4), возбужденных произвольными начальными отклонениями из множества (1.2.2), показатель качества, например (1.2.9), принимает наименьшее

. Если в (1.2.9) верхний предел t1 не ограничен, то стабилизирующее управление должно также обеспечивать асимптотическую

системы.

n δxi2 (t0 )≤ ε 2 i=1

Примечание 1.2.1.

Стабилизирующее управление реализуется регулятором, который является сложным динамическим устройством, состоящим обычно из трех компонент:

измерительных органов,

устройств реализации алгоритма управления (корректирующих контуров),

исполнительных органов.

Здесь и далее известные дифференциальные уравнения, описывающие

.

измерительные и исполнительные органы, включаются в уравнения (1.2.4) (δxi(t)).

Другими словами, уравнения (1.2.4) – это уравнения физического объекта вместе с

измерительными и исполнительными устройствами регулятора. Тогда δxi(t)

(i=l… n) – выходы, измерительных устройств, a δuk(t) (k=l… m) – входы исполнительных органов.

Уравнения (1.2.10) (δuk(t)) описывают устройство реализации алгоритмов управления.

Для упрощения терминологии будем по-прежнему называть уравнениями объекта уравнения (1.2.4) известной (неизменяемой) части системы, состоящей из объекта и элементов регулятора, а уравнениями регулятора – (1.2.10) – называть неизвестную (подлежащую определению) часть системы, состоящую лишь из

30

устройства реализации алгоритма управления.