Санкт-Петербургский государственный университет
Кафедра теории управления
Медведев Кирилл Сергеевич
Выпускная квалификационная работа бакалавра
Оценка допустимых погрешностей в одной
задаче стабилизации
Направление 010400
Прикладная математика, фундаментальная информатика
и основы программирования
Научный руководитель,
доктор физ.-мат. наук,
профессор
Харитонов В. Л.
Санкт-Петербург
2016
Содержание
Введение . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3
Глава 1. Ранее полученные результаты . . . . . . . . . . . . . . . . . . .
4
1.1. Дополнительные утверждения . . . . . . . . . . . . . . . . . . . .
7
Глава 2. Задача робастного анализа . . . . . . . . . . . . . . . . . . . . .
8
2.1. Возмущение запаздывания в состоянии . . . . . . . . . . . . . . .
8
2.2. Возмущение запаздывания в управлении . . . . . . . . . . . . . . 13
Пример решения задачи . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
Заключение . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
Список литературы . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
Приложение . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2
Введение
Основной целью данной работы является робастный анализ систем с запаздыванием в управлении и состоянии. Подобные задачи возникают в ситуациях, когда неизвестно точное значение параметров системы. Известно,
что при отклонении параметров замкнутой системы от расчетных значений
новая система может утратить свойство устойчивости. Именно поэтому важно знать, в каких интервалах могут меняться параметры системы так, чтобы система сохраняла устойчивость. Оценка будет производиться на основе
рассмотрения характеристических матриц возмущенной и невозмущенной систем.
3
Глава 1. Ранее полученные результаты
Рассмотрим систему с запаздыванием
ẋ(t) = Ax(t − h) + Bu(t − τ ).
(1)
Здесь матрицы A ∈ Rn×n и B ∈ Rn×m , пара (A, B) — полностью управляема,
запаздывание h > 0, τ ∈ [h, 2h].
По формуле Коши [1] решение данной системы можно представить как
R0
x(t + τ − h) = K(τ − h)x(t) +
K(τ − 2h − θ)Ax(t + θ)dθ +
−h
+
t+τ
R −h
K(t + τ − h − ξ)Bu(ξ − τ )dξ,
t>0
t
где K(t) — фундаментальная матрица системы. Известно, что она должна
удовлетворять уравнению
K̇(t) = AK(t − h),
t > 0,
и начальным условиям
K(t) =
0n×n ,
t < 0,
E ,
n×n
t = 0.
Так как (τ − h) ∈ [0, h], то фундаментальная матрица K(τ − h) = En×n .
Кроме того, на интервале t ∈ [0, h] справедливо
K(t − h) = 0n×n ,
K̇(t) = 0n×n ,
откуда с учетом K(0) = En×n следует
K(t) = En×n ,
4
t ∈ [0, h].
Таким образом решение может быть записано в виде
τZ−2h
Z−h
Ax(t + θ)dθ + Bu(t + ξ)dξ.
x(t + τ − h) = x(t) +
(2)
−τ
−h
Управление для системы выберем следующим образом
u(t) = F x(t + τ − h).
(3)
Используя формулу (2), получим
τZ−2h
Z−h
Ax(t + θ)dθ + F Bu(t + ξ)dξ.
u(t) = F x(t) + F
−τ
−h
Характеристическая функция замкнутой системы (1), (3) имеет вид
f0 (s) = det[sEn×n − e−sh (A + BF )].
Лемма [1]: Корни функции f0 (s) имеют отрицательную вещественную
часть тогда и только тогда, когда собственные числа матрицы A + BF
лежат в области
n
π π o
Γ = s = α + iβ | α = −ω sin ωh, |β| < |ω cos ωh|, ω ∈ − ,
2h 2h
комплексной плоскости.
Далле будем полагать, что выбор матрицы F обеспечивает выполнение
данной леммы.
Замена интегралов в формуле (3) конечными суммами может приводить к
неустойчивости замкнутой системы [2]. Для того чтобы не допустить этого в
работе [4] было предложено использовать дополнительный фильтр. Выберем
матрицу Гурвица [5] G ∈ Rm×m , и вычислим разность
5
du(t)
− Gu(t) = F ẋ(t + τ − h) − GF x(t + τ − h) =
dt
= F (Ax(t + τ − 2h) + Bu(t − h)) − GF x(t)−
τZ−2h
Z−h
Ax(t + θ)dθ − GF Bu(t + ξ)dξ.
− GF
−h
−τ
В результате получим замкнутую систему вида
ẋ(t) = Ax(t − h) + Bu(t − τ ),
τZ−2h
du(t)
x(t + θ)dθ +
= Gu(t) + F Bu(t − h) − GF (x(t) + A
dt
−h
−h
Z
+ B u(t + ξ)dξ) + F Ax(t + τ − 2h).
−τ
6
(4)
1.1.
Дополнительные утверждения
Утверждение 1: Для того чтобы E+D была невырождена достаточно,
чтобы выполнялось условие ||D|| < 1
Докажем от противного: пусть выполнено условие ||D|| < 1, но при этом
у матрицы E + D есть нулевое собственное число. Пусть x 6= 0 — собственый
вектор, соответствующий нулевому собственному числу. Тогда справедливы
следующие выражения
(E + D)x = 0x,
Dx = −x.
Из этих равенств с учетом max ||Dx|| = ||D|| получим
||x||=1
||x|| = ||Dx|| 6 ||D|| < 1,
что противоречит условию выбора x. Доказано
7
Глава 2. Задача робастного анализа
Рассмотрим возмущенную систему
dx(t)
= Ax(t − h1 ) + Bu(t − τ1 ),
dt
τZ−2h
du(t)
= Gu(t) + F Bu(t − h) − GF (x(t) + A
x(t + θ)dθ+
dt
−h
Z−h
+B u(t + ξ)dξ) + F Ax(t + τ − 2h).
(5)
−τ
где h1 = h + δ1 и τ1 = τ + δ2 . Здесь δ1 ∈ R и δ2 ∈ R — значения возмущений,
причем
|δ1 | < ρ1 ,
|δ2 | < ρ2 .
Задача робастного анализа заключается в поиске величин ρ1 и ρ2 таких, что
возмущенная система (5) сохраняет устойчивость.
2.1.
Возмущение запаздывания в состоянии
В данном разделе рассматривается случай, когда возмущение содержится
только в запаздывании состояния, то есть |δ2 | = 0. Тогда система (5) примет вид
dx(t)
= Ax(t − h1 ) + Bu(t − τ ),
dt
τZ−2h
du(t)
= Gu(t) + F Bu(t − h) − GF (x(t) + A
x(t + θ)dθ+
dt
−h
−h
Z
+B u(t + ξ)dξ) + F Ax(t + τ − 2h).
−τ
8
(6)
Характеристическая матрица системы (6) имеет вид
R(s) = R0 (s) + ∆(s) =
=
sEn×n − Ae−sh
GF +
GF A(es(τ −2h) −e−sh )
s
−Be−sτ
− F Aes(τ −2h)
sEm×m − G + Ω(s)
(7)
A(e−sh − e−sh1 ) 0n×m
,
+
0m×n
0m×m
где Ω(s) = −F Be−sh +
GF B(e−sh −e−sτ )
,
s
R0 (s) — характеристическая матрица
системы (4) без возмущений, ∆(s) — матрица, возникшая в силу добавления
возмущения.
Характеристической функцией возмущенной системы (6) является
q(s) = det(R(s)).
Для получения оценки ρ1 будем рассматривать ситуацию, когда нарушается устойчивость системы (6), то есть у характеристической функции появляется корень на мнимой оси
q(s) = det(R(iω)) = 0,
ω∈R
Это утверждение эквивалентно тому, что матрица R(iω) будет вырожденной. Воспользуемся тем фактом, что R(s) удовлетворяет соотношению (7) и
вынесем за скобки множитель R0 (iω)
R0 (iω) E(m+n)×(m+n) + R0−1 (iω)∆(iω) .
9
Известно, что R0 (s) — характеристическая матрица невозмущенной системы
(4), следовательно она не может иметь нулевого собственного числа на мнимой
оси. Это означает что такое собственное число должна иметь матрица
E(m+n)×(m+n) + R0−1 (iω)∆(iω)
Воспользовавшись вспомогательным утверждением 1, получим
||R0−1 (iω)|| ||∆(iω)|| > 1.
Это означает, что, если ∀ ω ∈ R выполнено
||R0−1 (iω)|| ||∆(iω)|| < 1,
то ни одно из собственных чисел системы (6) не попадает на мнимую ось.
Рассмотрим более подробно норму ||∆(iω)||. Известно, что
−sh
Ae
− Ae−sh1 0n×m
,
∆(s) =
0m×n
0m×m
откуда непосредственно получаем, что
||∆(s)|| = ||A|| |e−sh − e−sh1 |,
||∆(iω)|| = ||A|| |e−iωh − e−iωh1 |.
Дальнейшую оценку величины |e−iωh − e−iωh1 | можно провести двумя способами. Первый способ является более грубым, поскольку мы воспользуемся
свойством |eiφ | = 1, а значит получим
||∆(iω)|| 6 2||A||.
10
Основываясь на этой оценке можно сформулировать следующее утверждение.
Теорема 1: Если для системы (5) выполнены условия δ2 = 0, и ∀ ω ∈ R
||R0−1 (iω)|| <
1
,
2||A||
то эта система сохраняет устойчивость при любых допустимых значениях h1 .
Второй способ оценки заключается в том, что
|e−iωh − e−iωh1 | = |e−iω
Далее, учитывая, что |e−iω
h+h1
2
h+h1
2
(e−iω
| = 1 и |e−iω
h−h1
2
h−h1
2
− eiω
− eiω
h−h1
2
h−h1
2
)|.
1
| = |2sin(ω h−h
2 )|, по-
лучим
||∆(iω)|| 6 ||A|| |2sin(ω
h − h1
)|
2
Перед тем, как сформулировать еще одну теорему, рассмотрим матрицу
R0 (iω). Представим ее в виде
R0 (iω) = iω E − O(iω) ,
где
O11
O(iω) =
O21
O12
,
O22
e−iωh
O11 =
A
iω
e−iωτ
O12 =
B
iω
O21 = −
1
eiω(τ −2h) − e−iωh
eiω(τ −2h)
GF +
GF
A
−
FA
iω
ω2
iω
11
O22
1
e−iωh
e−iωh − e−iωτ
= G+
FB +
GF B
iω
iω
ω2
Нетрудно заметить, что ∃ r такое что ∀ |ω| > r ||O(iω)|| < 1.
Рссмотрим теперь норму ||R0−1 (iω)||. Для нее справедливо
||R0−1 (iω)||
∞
∞
1 X
1
1
1 X k
||O (iω)|| 6
||O(iω)||k =
<
=
|ω| ||E − O(iω)|| |ω|
|ω|
k=0
=
1
1
→ 0,
|ω| 1 − ||O(iω)||
k=0
ω → ∞.
Это означает, что ∀ M ∈ R ∃ Ω такое что ∀ |ω| > Ω ||R0−1 (iω)|| > M . Выбрав M =
1
2||A|| ,
теперь можно сформулировать вторую теорему.
Теорема 2: Если для системы (5) выполнены условия δ2 = 0, и ∃ Ω ∈ R
такое что ∀ |ω| > Ω
||R0−1 (iω)|| 6
1
,
2||A||
то для того, чтобы данная возмущенна система сохраняла устойчивость,
необходимо, чтобы ∀ |ω| < Ω выполнялось неравенство
h−h
1
1
|<
.
|sin ω
2
2||A||
12
2.2.
Возмущение запаздывания в управлении
В данном разделе рассматривается случай, когда возмущение содержится
только в запаздывании управления, то есть |δ1 | = 0. Тогда система (5) примет
вид
dx(t)
= Ax(t − h) + Bu(t − τ1 ),
dt
τZ−2h
du(t)
= Gu(t) + F Bu(t − h) − GF (x(t) + A
x(t + θ)dθ+
dt
−h
Z−h
+B u(t + ξ)dξ) + F Ax(t + τ − 2h),
(8)
−τ
где τ1 = τ + δ2 , δ2 ∈ R — значение возмущения, причем
|δ2 | < ρ2 .
Задача робастного анализа в данном случае заключается в поиске величины
ρ2 такой, что возмущенная система (8) сохраняет устойчивость.
Характеристическая матрица данной системы имеет вид
R(s) = R0 (s) + ∆(s) =
=
sEn×n − Ae−sh
GF +
GF A(es(τ −2h) −e−sh )
s
−Be−sτ
− F Aes(τ −2h)
sEm×m − G + Ω(s)
0n×n B(e−sτ − e−sτ1 )
,
+
0m×n
0m×m
13
(9)
где Ω(s) = −F Be−sh +
GF B(e−sh −e−sτ )
,
s
R0 (s) — характеристическая матрица
системы (3) без возмущений, ∆(s) — матрица, возникшая в силу добавления
возмущения.
Характеристической функцией возмущенной системы (8) является
q(s) = det(R(s)).
Для получения оценки ρ2 будем рассматривать ситуацию, когда нарушается устойчивость системы (8), то есть у характеристической функции появляется корень на мнимой оси
q(iω) = det(R(iω)) = 0,
ω∈R
Это утверждение эквивалентно тому, что матрица R(iω) будет вырожденной. Воспользуемся тем фактом, что R(s) удовлетворяет соотношению (9) и
вынесем за скобки множитель R0 (iω)
−1
R0 (iω) E + R0 (iω)∆(iω) .
Известно, что R0 (s) — характеристическая матрица невозмущенной системы (4), следовательно она не может иметь нулевого собственного числа на
мнимой оси. Это означает что такое собственное число должна иметь матрица
E + R0−1 (iω)∆(iω).
Воспользовавшись вспомогательным утверждением 1, получим
||R0−1 (iω)|| ||∆(iω)|| > 1.
Это означает, что если ∀ ω ∈ R выполнено
||R0−1 (iω)|| ||∆(iω)|| < 1,
14
то ни одно из собственных чисел системы (8) не попадает на мнимую ось.
Рассмотрим более подробно норму ||∆(iω)||. Известно, что
0n×n B(e−sτ − e−sτ1 )
,
∆(s) =
0m×n
0m×m
откуда непосредственно получаем, что
||∆(s)|| = ||B|| |e−sτ − e−sτ1 |,
||∆(iω)|| = ||B|| |e−iωτ − e−iωτ1 |.
Дальнейшую оценку величины |e−iωτ − e−iωτ1 | можно провести двумя способами. Первый способ является более грубым, поскольку мы воспользуемся
свойством |eiφ | = 1, а значит получим
||∆(iω)|| 6 2||B||.
Основываясь на этой оценке можно сформулировать следующее утверждение.
Теорема 3: Если для системы (5) выполнены условия δ1 = 0, и ∀ ω ∈ R
||R0−1 (iω)|| <
1
,
2||B||
то эта система сохраняет устойчивость при любых допустимых значениях τ1 .
Второй способ оценки заключается в том, что
|e−iωτ − e−iωτ1 | = |e−iω
Далее, учитывая, что |e−iω
τ +τ1
2
τ +τ1
2
(e−iω
τ −τ1
2
| = 1 и |e−iω
τ −τ1
2
− eiω
− eiω
получим
||∆(iω)|| 6 ||B|| |2sin(ω
15
τ −τ1
2
τ − τ1
)|.
2
)|.
τ −τ1
2
1
| = |2sin(ω τ −τ
2 )|,
Перед тем, как сформулировать еще одну теорему, рассмотрим матрицуR0 (iω).
Очевидно, что
R0 (iω) = iω E − O(iω) ,
где
O12
,
O22
O11
O(iω) =
O21
O11
e−iωh
=
A
iω
O12 =
O21
e−iωτ
B
iω
1
eiω(τ −2h) − e−iωh
eiω(τ −2h)
= − GF +
GF A −
FA
iω
ω2
iω
O22
e−iωh
e−iωh − e−iωτ
1
GF B
FB +
= G+
iω
iω
ω2
Нетрудно заметить, что ∃ r такое что ∀ |ω| > r ||O(iω)|| < 1.
Рссмотрим теперь норму ||R0−1 (iω)||. Для нее справедливо
||R0−1 (iω)||
∞
∞
1
1
1 X
1 X
<
=
||O(iω)|| 6
||O(iω)||k =
|ω| ||E − O(iω)|| |ω|
|ω|
k=0
=
k=0
1
1
→ 0 ω → ∞.
|ω| 1 − ||O(iω)||
Это означает, что ∀ M ∈ R ∃ Ω такое что ∀ |ω| > Ω ||R0−1 (iω)|| > M . Выбрав M =
1
2||B|| ,
теперь можно сформулировать вторую теорему.
Теорема 4: Если для системы (5) выполнены условия δ1 = 0, и ∃ Ω ∈ R
такое что ∀ |ω| > Ω
||R0−1 (iω)|| 6
16
1
,
2||A||
(10)
то для того, чтобы данная возмущенная система сохраняла устойчивость
необходимо, чтобы ∀ |ω| < Ω выполнялось неравенство
τ −τ
1
1
|sin ω
|<
.
2
2||A||
17
Пример решения задачи
Рассмотрим пример решения конкретной задачи. Пусть
−2
4 5
, B = , h = 1, τ = 1.5
A =
1
3 2
(11)
Выбрем так же
F = 8 8 ,
G = −1
благодаря чему собственные числа матрицы A + BF будут соответственно
λ1 = λ2 = −1
В данном случае из рисунков 1 и 2 следует
max||R0−1 (iω)|| < 835
ω∈R
||A|| = 7.2854
откуда, выбрав m = 835, получим
835 >
1
.
14.5708
Это означает, что условие первой теоремы не выполнено и нужно воспользоваться второй теоремой. На интервале ω ∈ [−39.4, 39.4]
||R0−1 (iω)|| >
1
,
2||A||
а значит на нем должно выполнятся неравенство
ω
1
|sin [h − h1 ] | <
2
2||A||
отсюда, выбрав ω = 39.4 следует оценка
|δ1 | < 0.0035.
Аналогично можно получить оценку для δ2
|δ2 | < 0.0249.
18
Заключение
В данной работе найдены достаточные условия робастной стабилизации
для класса дифференциальных уравнений с запаздываниями в управлении и
состоянии. Полученные теоремы позволяют найти оценку допустимых отклонений для запаздываний в управлении и состоянии, сохраняющих устойчивость возмущенной системы.
В качестве направлений дальнейших исследований отметим распространение полученных результатов на случаи, когда возмущение присутствует как
в управлении, так и в состоянии.
19
Список литературы
1. Bellman R., Cooke K. L. Differential Difference Equations. New York:
Academic Press, 1963. P. 462
2. Kharitonov V. L., 2015. Prediction Based Control: Implementation Issue.
Differential Equations and Control Processes, No 4. P. 52–65.
3. Manitius A. Z., Olbrot A. W. (1979). Finite spectrum assignment for systems
with delay, IEEE Trans. on Automatic Control, 24, 541-553.
4. Mondie S., Michiels W.,2003. Finite spectrum assignment of unstable timedelay sytems with a safe implementation, IEEE Trans. on Automatic Control,
No 12. P. 2207-2212.
5. Емельянов С. В., Коровин С. К., Ильин А. В., Фомичев В. В., Фурсов
А. С. Математические методы теории управления. Проблемы устойчивости, управляемости и наблюдаемости. М.: ФИЗМАТЛИТ, 2013. 200 с.
20
Приложение
Рис. 1: График функции f (ω) = ||R0−1 (iω)||.
21
Рис. 2: Графики функций f (ω) = ||R0−1 (iω)|| (синий) и g(ω) =
22
1
2||A||
(красный).
Отзывы:
Авторизуйтесь, чтобы оставить отзыв