- •Введение
- •1 Обзор методов управления в условиях запаздывания и возмущений. Постановка задачи
- •1.1 Обзор методов управления в условиях запаздывания
- •1.1.1 Предиктор Смита
- •1.1.2 Предиктор на основе метода М. Крстича
- •1.2 Обзор методов управления в условиях возмущающих воздействий
- •1.3 Обобщенная постановка задачи
- •2 Методы построения адаптивных наблюдателей мультигармонических сигналов
- •2.1 Алгоритм адаптивной идентификации параметров смещенного синусоидального сигнала
- •2.1.1 Постановка задачи
- •2.1.2 Алгоритм идентификации частоты смещенного гармонического сигнала
- •2.1.3 Алгоритм идентификации смещения, амплитуды и фазы
- •2.1.4 Числовой пример
- •2.2 Алгоритм адаптивной идентификации частот и наблюдатель гармоник мультигармонического сигнала
- •2.2.1 Постановка задачи
- •2.2.2 Алгоритм идентификации частот смещенного мультигармонического сигнала
- •2.2.3 Алгоритм идентификации смещения, амплитуд и фаз гармоник
- •2.2.4 Числовой пример
- •2.3 Заключительные выводы по главе
- •3 Компенсация мультигармонических возмущений для устойчивых объектов с запаздыванием в управлении
- •3.1 Алгоритм компенсации мультигармонического возмущения, действующего на устойчивый линейный объект управления с запаздыванием
- •3.1.1 Постановка задачи
- •3.1.2 Алгоритм адаптивной идентификации частот
- •3.1.3 Синтез закона управления
- •3.1.4 Числовой пример
- •3.2 Алгоритм компенсации мультигармонического возмущения, действующего на устойчивый нелинейный объект управления с запаздыванием
- •3.2.1 Постановка задачи
- •3.2.2 Преобразование нелинейной системы
- •3.2.3 Алгоритм адаптивной идентификации частот
- •3.2.4 Синтез закона управления
- •3.2.5 Числовой пример
- •3.3 Экспериментальные исследования алгоритма управления
- •3.4 Заключительные выводы по главе
- •4 Компенсация мультигармонических возмущений для неустойчивых объектов с запаздыванием в управлении
- •4.1 Постановка задачи
- •4.2 Алгоритм стабилизации неустойчивого объекта управления с запаздыванием
- •4.3 Алгоритм адаптивной идентификации частот и наблюдатель гармоник мультигармонического возмущения
- •4.4 Алгоритм компенсации мультигармонического возмущения, действующего на неустойчивый объект управления с запаздыванием
- •4.5 Числовой пример
- •4.6 Заключительные выводы по главе
- •Заключение
- •Литература
3.1. Компенсация возмущения для линейного объекта 77
Подставляя (3.60) в (3.7) с учетом леммы 3.1, получим
( ) = |
[ ( )](~( ) − 0 0 |
|
( + − |
|
|
|
)) |
|||||||||
( + ) − =1 |
||||||||||||||||
|
|
( ) |
|
1 |
|
|
∑ |
1 |
|
|
+ |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∑ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ 0( ) + |
( ) + ( ) |
|
|
|
|
|
|
|
|||||||
|
=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
= ( ) + 2( ), |
|
[ ( ) ] |
|
|
|
|
|
|
|
|
||||||
где функция ( ) |
= |
~( ) ограничена и стремится к нулю, |
||||||||||||||
|
|
|
|
|
|
( ) |
|
|
|
|
|
|
|
|
|
2( ) — экспоненциально затухающая функция. Теорема 3.1 доказана.
Таким образом, получен реализуемый закон управления (3.43)– (3.47), обеспечивающий основную цель (3.4) lim →∞ ( ) = 0.
З а м е ч а н и е 3.3 Следует отметить, что в случае неточного задания параметров объекта рассматриваемый алгоритм управления может оказаться неудовлетворительным. Разработка робастных алгоритмов, обеспечивающих достижение сформулированной цели управления в случае неточно заданных параметров, является предметом дальнейших исследований.
3.1.4Числовой пример
Для демонстрации эффективности предлагаемого алгоритма управления представим результаты математического моделирования неминемально-фазового объекта (3.1)-(3.3) с относительной степенью − = 2 и параметрами
= |
−35 0 1 |
0 |
, = |
|
0 , |
= |
0 |
, 0 |
= |
1 |
, |
|||||||||||||
|
|
−10 |
1 |
0 |
0 |
|
|
|
0 |
|
|
|
1 |
|
|
|
|
1 |
|
|
|
|||
|
|
24 0 0 |
0 |
|
13 |
|
0 |
|
|
0 |
|
|||||||||||||
|
|
−50 |
0 |
0 |
1 |
|
|
−5 |
|
|
|
|
|
|
|
|||||||||
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, 1 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
= |
10 |
, 1 |
= |
4 |
= |
6 |
. (3.61) |
|||||||||
|
|
|
|
|
|
|
|
|
5 |
|
|
2 |
|
|
|
|
|
3 |
|
|
||||
|
|
|
|
|
|
|
|
|
4 |
|
1 |
|
|
3 |
|
|||||||||
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
78 Глава 3. Устойчивые объекты с запаздыванием
На рис. 3.1 представлены переходные процессы в замкнутой системе для частоты возмущающего воздействия = 1, 5, параметров
закона управления = 10, = 2 и запаздывания = 6. На рис. 3.2а,
3.2б представлены переходные процессы в замкнутой системе для запаздывания = 5, частоты возмущающего воздействия = 3 и
различных значениях параметров алгоритма управления. На рис. 3.2в, 3.2г представлены переходные процессы в замкнутой системе для различных значений запаздывания , включая его отсутствие
= 0, частоты возмущающего воздействия = 3 и параметров алгоритма управления = 10, = 3. На рис. 3.2в, 3.2г представле-
ны переходные процессы в замкнутой системе для возмущающего воздействия с нерегулярной составляющей 1.
Для иллюстрации эффективности разработанного алгоритма
компенсации возмущений рассмотрим случай с двумя гармониками. Выберем параметры возмущения (3.3) 2 = [1 2 3 − 1] и
2 = [2 4 − 1 5].
На рис. 3.3 представлены результаты моделирования замкнутой системы для частот 1 = 2 и 2 = 3, запаздывания = 4, пара-
метров закона управления 1 = 20, 2 = 20, = 2. На рис. 3.4
представлен результат моделирования мультигармонического возмущения, частоты которого изменились скачкообразно с течением времени.
Из рис. 3.1-3.4 видно, что цель управления (3.4) достигается при различных значениях запаздывания, даже при его значительном увеличении. Отсутствие временной задержки является частным случаем рассматриваемой задачи, где = 0. С увеличением
значений параметров закона управления и идентификация ча-
стоты и, как следствие, компенсация возмущения происходит быстрее. Показаны адаптивные и робастные свойства алгоритма компенсации возмущения.
1Нерегулярная составляющая моделировалась с помощью формирующего фильтра ( ) = 1
10 +1 , на вход которого подавалось белошумное воздействие (с ограниченной частотой) мощностью = 0, 1. Матрица входов для этого
воздействия была принята равной вектору .
3.1. Компенсация возмущения для линейного объекта 79
ωˆ |
|
|
|
|
1 |
|
|
|
|
|
|
2,5 |
|
|
|
|
|
|
|
|
|
|
0,5 |
2 |
|
|
|
|
|
1,5 |
|
|
|
|
0 |
1 |
|
|
|
|
|
|
|
|
|
|
−0,5 |
0,5 |
|
|
|
|
|
|
|
|
|
|
t, c |
0 |
10 |
20 |
30 |
40 |
−1 |
0 |
50 |
|
|
|
|
|
a |
|
|
|
|
|
b |
|
|
|
|
|
t, c |
0 |
10 |
20 |
30 |
40 |
50 |
(а) Временная диаграмма оценки ча- |
(б) Выходная переменная ( ) без |
||||||||||||||||
стоты ^( ) |
|
|
|
|
|
управления (a) и с управлением (b) |
|||||||||||
1 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆ |
−0,5 |
|
|
|
|
|
|
|
|
|
ϕˆ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|
ϕˆξ |
0,8 |
|
|
|
|
|
ˆ |
−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
Lξ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−1,5 |
|
|
|
|
|
|
|
|
|
|
0,6 |
|
|
|
|
|
|
−2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,4 |
|
|
|
|
|
|
−2,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−3 |
|
|
|
|
|
|
|
|
|
|
0,2 |
|
|
|
|
|
|
−3,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t, c |
−4 |
|
|
|
|
|
|
|
|
|
t, c |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
10 |
|
20 |
30 |
40 |
50 |
−4,5 |
10 |
|
20 |
|
|
|
|
30 |
40 |
50 |
|
0 |
|
0 |
|
|
|
|
|
||||||||||
(в) |
Временная |
диаграмма |
оценки |
(г) |
Временная |
диаграмма |
оценки |
||||||||||
функций ^ |
|
и |
^ |
|
|
функций ^( ) и ^ |
|
( ) |
|
|
|||||||
|
( ) |
|
( ) |
|
|
|
|
|
|
|
|
|
|
|
|||
ˆ |
|
|
|
|
|
|
25 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ξ0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
kp |
|
|
|
|
|
|
|
20 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
kd |
0,5 |
|
|
|
|
|
|
15 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
−0,5 |
|
|
|
|
|
|
−5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t, c |
−10 |
|
|
|
|
|
|
|
|
|
|
−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t, c |
|
10 |
|
20 |
30 |
40 |
50 |
−15 |
10 |
|
20 |
|
|
|
|
30 |
40 |
50 |
|
0 |
|
0 |
|
|
|
|
|
||||||||||
(д) |
Временная |
диаграмма |
оценки |
(е) |
Временная |
диаграмма |
оценки |
||||||||||
функции ^ |
( ) |
|
|
|
|
функций |
|
( ) и |
|
( ) |
|
|
|||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
Рис. 3.1. Временные диаграммы оценки частоты возмущения = |
|||||||||||||||||
1, 5, выходной переменной, параметров закона управления при за- |
|||||||||||||||||
паздывании = 6 и коэффициентах идентификатора = 2, = 10 |
80 |
|
|
Глава 3. Устойчивые объекты с запаздыванием |
|||||||||
ωˆ |
|
|
|
|
|
|
2 |
|
|
|
|
|
4 |
|
|
|
k = 1, λ = 1 |
|
1 |
|
|
|
|
|
|
|
|
|
k = 10, λ = 5 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||
3 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
U = 0 |
|
|
|
|
|
|
|
−2 |
|
|
|
k = 1, λ = 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
t, c |
|
|
|
|
k = 10, λ = 5 |
t, c |
0 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||
5 |
10 |
15 |
20 |
25 |
30 |
−3 |
|
|
|
|
|
|
0 |
10 |
20 |
30 |
40 |
50 |
|||||||
|
|
|
|
|
|
|
0 |
|||||
(а) Временная диаграмма оценки ча- |
(б) Выходная переменная ( ) при |
|||||||||||
стоты возмущения = 3 |
|
|
запаздывании = 4 |
|
|
|
||||||
ωˆ |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
0,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−0,5 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
h = 0 |
|
|
|
|
|
|
|
−1,5 |
|
|
|
h = 6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
t, c |
|
|
|
|
h = 15 |
t, c |
0 |
|
|
|
|
|
−2 |
|
|
|
|
||
|
10 |
20 |
|
30 |
40 |
10 |
20 |
|
30 |
40 |
||
0 |
|
|
0 |
|
||||||||
(в) Оценки частоты ^( ) при различ- |
(г) Выходная переменная ( ) при |
|||||||||||
ном запаздывании |
|
|
|
различном запаздывании |
|
|||||||
ωˆ |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
0,5 |
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
|
|
||
1,5 |
|
|
|
|
|
|
0 |
|
|
|
|
|
1 |
|
|
|
|
|
|
−0,5 |
|
|
|
|
|
0,5 |
|
|
|
|
|
|
−1 |
|
|
|
|
|
|
|
|
|
|
|
t, c |
|
|
|
|
|
t, c |
0 |
|
|
|
|
|
−1,5 |
|
|
|
|
|
|
|
20 |
40 |
60 |
80 |
100 |
20 |
40 |
60 |
80 |
100 |
||
0 |
|
0 |
(д) Временная диаграмма оценки ча- |
(е) Выходная переменная ( ) без |
стоты возмущения ^( ) с аддитивной |
управления (a) и с управлением (b) |
нерегулярной составляющей |
при запаздывании = 4 и коэффици- |
|
ентах идентификатора = 1 и = 10 |
Рис. 3.2. Временные диаграммы оценки частоты и выходной переменной ( )
3.1. Компенсация возмущения для линейного объекта 81
4 |
|
|
|
|
3,5 |
|
|
|
ωˆ1 |
|
|
|
ωˆ2 |
|
|
|
|
|
|
3 |
|
|
|
|
2,5 |
|
|
|
|
2 |
|
|
|
|
1,5 |
|
|
|
|
1 |
|
|
|
|
0,5 |
|
|
|
|
|
|
|
|
t, c |
0 |
50 |
100 |
150 |
200 |
0 |
(а) Временная диаграмма оценки частот ^ ( )
1,2
|
L1 |
1 |
L2 |
|
Lξ1 |
0,8 |
Lξ2 |
|
|
0,6 |
|
0,4
0,2
t, c
0
0 |
50 |
100 |
150 |
200 |
(в) Временная диаграмма оценки
функций ^ ^
( ) и ( )
ˆ
ξ0
0,5
0
−0,5
t, c
−1
0 |
50 |
100 |
150 |
200 |
(д) Временная диаграмма оценки
функции ^
0( )
1,5
a
1 |
b |
0,5
0
−0,5
−1
t, c
−1,5
0 |
50 |
100 |
150 |
200 |
(б) Выходная переменная ( ) без управления (a) и с управлением (b)
1 |
|
|
ϕˆ1 |
0 |
ϕˆ2 |
|
ϕˆξ1 |
−1 |
ϕˆξ2 |
|
|
−2 |
|
−3
−4
t, c
−5
0 |
50 |
100 |
150 |
200 |
(г) Временная диаграмма оценки функций ^ ( ) и ^ ( )
10
0 |
|
|
|
|
−10 |
|
|
|
|
−20 |
|
|
kp1 |
|
|
|
|
|
|
|
|
|
kp2 |
|
−30 |
|
|
kd1 |
|
|
|
|
|
|
|
|
|
kd2 |
t, c |
−40 |
50 |
100 |
150 |
200 |
0 |
(е) Временная диаграмма оценки функций ( ) и ( )
Рис. 3.3. Временные диаграммы оценки частот возмущения 1 = 1,2 = 2, выходной переменной ( ) и параметров закона управления при запаздывании = 6 и коэффициентах идентификатора = 2,
= 10