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

Некоторые примеры неустойчивых задач управления. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
670 Кб
Скачать

вое условие на координату материальной точки (3.53):

x1(1) = a,

(4.2)

где a [0, 1]. Кроме того, мы предполагаем, что выполнено условие (3.55), которое в рассматриваемом случае неотрицательных управлений сводится к требованию

1

u(t)dt ≤ 1.

(4.3)

0

Таким образом, рассматривается материальная точка с импульсным ограничением (4.3) и краевым условием (4.2). Согласно (4.1), наша система в начальный момент t = 0 покоится в нуле. Неотрицательность программного управления в системе (3.53) означает содержательно, что используется нереверсируемый двигатель.

Будем конструировать ОД по скорости (т. е. по второй координате в (3.53)) при наличии вышеупомянутых ограничений; кроме того, будем рассматривать ослабление краевого условия (4.2), т. е. замену (4.2) неравенством

|x1(1) − a| < ε,

(4.4)

где ε ]0, ∞[. Для того чтобы исключить субъективизм в выборе параметра ε в (4.4), будем конструировать множество притяжения, которое в настоящем разделе будет введено достаточно строго (впрочем, данное множество в идейном отношении подобно (3.52)).

Полагаем в дальнейшем, что F есть по определению множество всех неотрицательных в/з функций u = u(·) = (u(t), 0 ≤ t < 1) на промежутке [0, 1[ = R | (0 ≤ ξ) & (ξ < 1)}, для каждой из которых выполнено (4.3). Элементы множества F можно рассматривать в качестве программных управлений.

51

Если u F, то, с учетом (4.1), компоненты траектории системы (3.53) имеют вид (3.57), (3.58). Для всех наших целей будет существенным только (терминальное) состояние систе-

мы в последний момент времени: если u F, то

1

φu,1(1) =

0

(1 − t)u(t)dt,

(4.5)

φu,2(1) =

0

1 u(t)dt.

(4.6)

В виде φu(1) = (φu,1(1), φu,2(1)) имеем состояние управляемой системы (3.53) в последний момент времени. Условие

(4.2) преобразуется в силу (4.5) к требованию

 

0

1(1 − t)u(t)dt = a.

(4.7)

Поэтому множество всех допустимых программных управлений (имеется в виду допустимость в смысле (4.2), (4.3)) имеет вид 1

F[a] = {u F | (1 − t)u(t)dt = a}. (4.8)

0

С учетом (4.6) и (4.7) мы можем рассматривать множество

{∫ 1 }

G[a] = u(t)dt : u F[a] (4.9)

0

как область достижимости по скорости при точном соблюдении условий (4.2), (4.3). Условие (4.3) возмущаться не будет.

Что же касается условия (4.2), то мы будем рассматривать его ослабленные аналоги в виде условия (4.4), где параметр ε, ε > 0, может быть различным. Тогда полагаем

 

 

1

 

(4.10)

F(ε)[a] = {u F|

 

0

 

(1 − t)u(t)dt − a

< ε} ε ]0, ∞[.

 

 

 

 

 

 

 

 

 

 

 

 

52

Считаем,

что

при ε

 

]0,

[ программные управления

(ε)

 

 

 

 

 

 

u F

[a],

и только они допустимы в смысле соблюдения

условия (4.4). Эти управления формируют новую ОД:

G(ε)[a] = {∫0

1

u F(ε)[a]}

ε ]0, ∞[. (4.11)

u(t)dt :

Каждое множество (4.11) содержит ОД (4.9). Для целей «исключения» параметра ε из конструкции (4.11) полагаем,

что

 

 

 

 

 

 

 

 

 

 

 

 

(ε)

 

 

(4.12)

(AS)[a] =

 

G

[a]

ε ]0,

[

 

 

 

 

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

всмысле обычной | · |–топологии вещественной прямой R). Ясно, что

G[a] (AS)[a].

(4.13)

Вложение (4.13) следует из (4.12) и ранее упомянутых свойств ОД (4.11). Тогда [9]

(AS)[a] = [a, 1] = G[a] a ]0, 1[;

(4.14)

свойство (4.14) можно содержательно интерпретировать как устойчивость ОД при замене (4.2) на (4.3) в случае малых ε, ε > 0. Кроме того, из [9] следует, что

(G[0] = {0}) & ((AS)[0] = [0, 1]).

(4.15)

Мы видим, что случай a = 0 соответствует ощутимой неустойчивости в вопросах построения ОД (в левой части (4.12) имеем одноэлементное множество, а в правой – континуум [0, 1]). Итак, при переходе от a ≈ 0 к случаю a = 0 мы имеем скачок ОД. При этом, увеличивая a сколь угодно мало в сравнении с нулем, мы «отстраиваемся» от неустойчивости. Эту возможность можно переформулировать в несколько иных терминах.

53

Рассмотрим систему (3.53) в условиях, когда начальные координаты (первая компонента фазового вектора в начальный момент времени и скорость) могут быть малыми положительными числами.

Зафиксируем a0 ]0, 1[. Кроме того, заменим (4.1) условиями

x1(0) = x10, x2(0) = x20,

(4.16)

где x01 ]0, ∞[ и x02 ]0, ∞[. Тогда рассмотрим значение

x01 + x02 ]0, ∞[,

с помощью которого определяем новое краевое условие

x1(1) = a0,

(4.17)

где, как уже отмечалось, a0 ]0, ∞[; при этом значение x1(1) всякий раз соответствует траектории с начальным условием (4.16). Иными словами, сохраняя ограничение (4.3) на выбор (неотрицательного) программного управления, будем рассматривать ОД по скорости при начальных условиях (4.16) и краевом условии (4.17). В этом случае

1

0 0

F

0

0

 

(1 − t)u(t)dt = a0} (4.18)

F(x1, x2) = {u

 

| (x1

+ x2) +

0

есть множество всех допустимых программных управлений. Значения x01 и x02 можно будет, в принципе, менять, но толь-

ко так, чтобы

x10 + x20 ≤ a0.

(4.19)

Заметим при этом, что множество всех неотрицательных покомпонентно векторов (x01, x02), удовлетворяющих (4.19), имеет мощность континуума, т. к. a0 > 0; это множество содержит нулевой вектор.

54

Согласно (4.8), определено множество F[a0 (x01 + x02)]; более того, из (4.18) вытекает, что

F(x10, x20) = F[a0 (x10 + x20)].

(4.20)

Соответственно, мы имеем в виде множества

1

0 0

0

+

0 0

G(x1, x2) = {x2

u(t)dt : u F(x1, x2)} =

0

1

={x02 + u(t)dt : u F[a0 (x01 + x02)]} =

0

 

= {x20 + y : y G[a0 (x10 + x20)]}.

(4.21)

Теперь, действуя по аналогии с (4.10), мы

введем при

ε ]0, ∞[ множество всех ε–допустимых (в смысле (4.17)) программных управлений:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(x1, x2, ε) =

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

= {u F|

 

(x10

+ x20) +

 

 

0

(1 − t)u(t)dt − a0

< ε} =

 

 

 

 

 

 

 

 

 

 

 

 

(ε)

 

 

 

 

 

 

 

 

0

0

 

 

 

 

=

 

u

 

F

 

(1

 

 

 

 

 

 

 

(a0

 

(x1

 

 

< ε

 

=

{

 

|

 

 

t)u(t)dt

+ x2))

}

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= F

 

[a0

 

 

(x0

+ x0)];

 

 

 

 

(4.22)

Как следствие (см. (4.10)) при ε ]0, ∞[ мы в виде множества

G(x10, x20, ε) = {x20 +

1

u(t)dt : u F(x10, x20, ε)} =

0

1

={x02 + u(t)dt : u F(ε)[a0 (x01 + x02)]}

0

(см. (4.22)) получаем сдвиг множества, подобного (4.11), откуда, с учетом (4.11), вытекает равенство

G(x01, x02, ε) = {x02 + y : y G(ε)[a0 (x01 + x02)]}. (4.23)

55

Для того чтобы исключить субъективность в выборе ε, ε > 0, мы, как и в (4.12), введем множество притяжения. А именно: полагаем, что

 

 

0 0

 

 

 

(4.24)

(

AS

 

 

0

0

 

 

)[x1; x2] =

 

 

G(x1, x2, ε).

 

 

ε

 

]0, [

 

 

 

Фиксируем произвольные

значения

x10

]0, ∞[ и

x20 ]0, ∞[, для которых

 

 

 

 

 

 

 

 

a0 (x10 + x20) [0, 1[.

 

 

(4.25)

В связи с (4.24) отметим одно простое представление замыканий множеств (4.23):

G(x01, x02, ε) =

= {x02 + y : y G(ε)[a0 (x01 + x02)]} ε ]0, ∞[. (4.26)

Свойство (4.26) является простейшим вариантом весьма общего положения теории топологических векторных пространств: замыкание множества, полученного сдвигом на x0 другого множества, есть сдвиг замыкания последнего. В целях полноты изложения приведем соответствующее рассуждение в следующем замечании.

Замечание 4.1

Фиксируем ε ]0, ∞[. Для краткости полагаем, что

 

0

+ y :

(ε)

0

0

(4.27)

Ω = {x2

y G

[a0 (x1

+ x2)]}.

Кроме того, используем (4.23). Нам требуется установить равенство

G(x10, x20, ε) = Ω.

(4.28)

56

Пусть x G(x01, x02, ε). Тогда можно указать последовательность

 

 

(x(k))k N : N G(x10, x20, ε)

(4.29)

(здесь и ниже

N

 

 

 

= {1; 2; ...}), для которой

 

 

 

(x(k))k N → x .

(4.30)

С учетом (4.23) имеем очевидное свойство

 

x(k) − x20 G(ε)[a0 (x10 + x20)] k N.

(4.31)

Из (4.30) вытекает свойство сходимости (x(k) −x02)k N → y ,

 

0

где y = x − x2. С учетом (4.31) имеем из этого последнего свойства включение

y G(ε)[a0 (x01 + x02)],

а тогда x = x02 + y Ω, согласно (4.27), чем и завершается обоснование вложения:

 

G(x10, x20, ε) Ω.

(4.32)

 

 

0

 

Пусть w . С учетом (4.27) получаем, что w˜ = w − x2

 

 

 

 

G(ε)[a0 (x10 + x20)], и тогда можно указать

последователь-

ность

 

(zk)k N : N → G(ε)[a0 (x10 + x20)],

для которой имеет место сходимость

 

(zk)k N → w˜.

(4.33)

В этом случае, согласно (4.23), справедлива система включений

x20 + zk G(x10, x20, ε) k N.

(4.34)

57

При этом x02 + w˜ = w, а потому из (4.33) следует сходимость (x02 + zk)k N → w. С учетом (4.34) получаем следующее включение

w G(x01, x02, ε).

Поскольку выбор w был произвольным, установлено вложение

G(x01, x02, ε),

откуда, с учетом (4.32), вытекает (4.28), и следовательно (см. (4.27)), требуемое равенство:

G(x01, x02, ε) = {x02 + y : y G(ε)[a0 (x01 + x02)]}.

Поскольку выбор ε был произвольным, свойство (4.26)

установлено.

2

Предложение 4.1 Справедливо равенство

 

(AS)[x10; x20] = {x20 + y : y (AS)[a0 (x10 + x20)]}.

(4.35)

Доказательство. Обозначим для краткости множество в правой части (4.35) через Ξ. Требуется установить равенство

(AS)[x10; x20] = Ξ.

(4.36)

Пусть x (AS)[x01; x02]. Тогда x R удовлетворяет, в силу (4.24), системе включений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x G(x10, x20, ε)

ε ]0, ∞[.

 

(4.37)

Из (4.26) и (4.37) вытекает, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x0

G(ε)[a

0

(x0

+ x0)]

ε

]0,

[.

 

2

 

1

2

 

 

 

 

Поэтому, согласно

(4.12),

имеем включение

x − x20

(AS)[a0 (x10 + x20)], и тогда

 

 

 

 

 

 

x = x02 + (x − x02) Ξ,

58

чем завершается проверка вложения

 

(AS)[x10; x20] Ξ.

(4.38)

Выберем произвольное число ξ Ξ. Тогда, по определению

множества Ξ, мы имеем

 

 

 

 

 

 

ξ − x20 (AS)[a0 (x10 + x20)].

(4.39)

Из (4.12) и (4.39) вытекает система включений

 

 

 

 

 

 

 

 

ξ − x20 G(ε)[a0 (x10 + x20)] ε ]0, ∞[.

 

Тогда, с учетом (4.26), мы получаем

 

 

 

 

 

 

ξ G(x10, x20, ε) ε ]0, ∞[.

 

Иными словами, установлено следующее включение

 

 

 

 

 

 

 

 

 

 

ξ

G(x10, x20, ε).

 

 

 

 

 

 

ε

]0, [

 

 

 

 

С учетом (4.24) получаем включение ξ (AS)[x01; x02]. Поскольку выбор ξ был произвольным, установлено вложение

Ξ (AS)[x01; x02],

из которого, с учетом (4.38), вытекает требуемое неравенство

(AS)[x01; x02] = Ξ.

Используя определение множества Ξ, мы получаем равенство (4.35). 2 Введем в рассмотрение множество X всех упорядоченных

пар (x1, x2) ]0, ∞[ × ]0, ∞[, для каждой из которых

a0 (x1 + x2) [0, 1[.

(4.40)

Поскольку в предыдущих построениях выбор x01 и x02 был произвольным, с соблюдением условия (4.25), мы можем использовать конструкцию (4.24) при замене x01 → x1 и x02 → x2

59

в случае (x1, x2) X: определено множество (AS)[x1; x2] (при упомянутой замене определены также множества G(x1, x2, ε) при ε ]0, ∞[). Кроме того, с учетом (4.40), определено (см. (4.12)) множество (AS)[a0 (x1 + x2)], где (x1, x2) X. Наконец, согласно предложению (4.1), имеем систему

(AS)[x1; x2] =

= {x2 + y : y (AS)[a0 (x1 + x2)]} (x1, x2) X. (4.41) Напомним теперь, что при (x1, x2) X мы, согласно (4.9) и (4.40), можем определить множество G[a0(x1+x2)], причем,

используя (4.21), мы, с одной стороны, определяем

1

G(x1, x2) = {x2 + u(t)dt : u F(x1, x2)}

0

(«точная» ОД для случая начальных условий, определяемых вектором (x1, x2)), а с другой — (также см. (4.21)) имеем равенство

G(x1, x2) = {x2 + y : y G[a0 (x1 + x2)]}. (4.42)

Согласно (4.14), (4.41) и (4.42), мы имеем следующее положение: если (x1, x2) X и при этом

0 < a0 (x1 + x2),

(4.43)

то справедливо равенство

(AS)[x1; x2] = {x2 + y : y G[a0 (x1 + x2)].

(4.44)

Теперь, используя свойство, подобное (4.21), мы имеем (в случае (4.43)) представление

(AS)[x1; x2] = G(x1, x2),

(4.45)

имеющее смысл устойчивости задачи о построении ОД по скорости.

60

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]