ОиАС-1
.pdfПоэтому реальное движение описывается функциями:
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 ),K,δxn (t ),δu1 (t ),K,δum (t ),t ] (i = |
|
) |
(1.2.3) |
1, n |
δϕi [δx1 (t ),K,δxn (t ),δu1 (t ),K,δum (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 ,K,δxn ,δu1 ,K,δum ,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);
oi(δx1, ..., δ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)=rk[δx1(t),…, δxn(t), t] (k=1… m). |
(1.2.10) |
Заметим, что вид этих функций не зависит от начальных условий из множества (1.2.2). Поясним это обстоятельство подробнее.
27
Допустим, что найдено управление
δuk0(t)=rk0[δx1(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 ,K,δxn ,δu1 ,K,δum ,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
устройства реализации алгоритма управления.