Герасимов Адаптивное управление
.pdf
где =δ 2 / 2λ +θ2 / 4γ постоянная величина. Решая полученное дифференциальное неравенство, получаем:
V (ε(t)) ≤ e−λtV (0) + λ − λ e−λt ,
откуда с учетом (2.5) следует, что
1 |
|
2 |
|
−λt |
|
|
|
|
|
|
−λt |
|
|
|
|
|
|
|
− λ e |
|
|
|
|||||||
2 |
ε |
|
≤ e |
V (0) + λ |
|
|
|
|
||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ε(t) ≤ |
|
−λt |
+ |
|
|
− |
|
e |
−λt |
(2.6) |
||||
2 e |
V (0) |
λ |
λ |
. |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
Из последнего неравенства следует экспоненциальная сходимость ошибки управления ε
к ограниченному |
множеству с |
границей |
= |
2 |
/ λ . При этом величину |
можно |
уменьшить путем |
увеличения |
коэффициентов |
λ |
и γ . Как следствие, величина θˆ |
||
становится ограниченной.
Таким образом, алгоритм управления (1.6), (2.4) обеспечивает устойчивость в замкнутой системе и является робастным по отношению к внешнему возмущению. В то же время этот алгоритм имеет следующие недостатки, а именно:
при отсутсвии возмущения установившаяся ошибка ε(t) отлична от нуля, что видно из неравенства (2.6);
если подставить (2.4) в (1.6), то нетрудно видеть, что управление
пропорционально величине x2 . Следовательно при росте x амплитуда управления возрастает квадратично, в связи с чем практическая применимость такого закона (1.6) имеет существенные ограничения.
Рассмотрим решение, лишенное недостатков алгоритмов (1.6), (1.9) и (1.6), (2.4) за счет наделения нового алгоритма управления адаптивными и робастными свойствами.
Решение №3
Рассмотрим совместно с настраиваемым регулятором (1.6) алгоритм адаптации, параметрический дрейф в котором ограничивается линейной стационарной обратной связью:
& |
(2.7) |
θˆ = −σθˆ −γxε , |
где σ постоянная положительная величина.
Проведем анализ устойчивости замкнутой системы, представленной объектом (2.1), регулятором (1.6) и алгоритмов адаптации (2.7) с помощью функции Ляпунова (1.8). Возьмем производную от функции и проведем ряд преобразований:
11
|
|
|
|
1 |
|
|
|
|
|
1 ~~& |
|
|
|
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
~ & |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
~ |
|
|
|
|
|
|
|
|
|
|
|
θ~ |
|
|
|
ˆ |
|||||||||||||||||
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
V = 2 |
2εε |
+ |
2γ |
2θθ |
= ε(−λε |
−θ x −δ )− |
γ θθ |
|
= −λε − |
θ xε −δε − γ (− |
σθ −γxε)= |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= −λε2 −δε + |
|
σ |
θ~θˆ |
= −λε2 −δε |
+ |
|
σ |
θ~(− |
θ~ +θ)= −λε2 −δε − σ |
θ~2 + |
σ |
θ~θ = |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
γ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
γ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
γ |
|
|
|
|
γ |
|
|
||||||
|
|
|
|
λ |
2 |
|
λ |
|
2 |
|
|
|
|
|
1 |
|
2 |
|
σ |
~2 |
|
|
σ |
|
~2 |
|
|
|
σ |
|
~ |
|
|
|
|
|
|
|
σ |
θ |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
= − |
2 |
ε |
|
− |
2 ε |
|
|
|
−δε ± 2λ |
δ |
|
− |
2γ |
θ |
|
|
|
− |
2γ θ |
|
|
+ |
|
γ |
θθ ± |
|
2γ |
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
= − |
λ |
ε |
2 |
− |
σ |
|
|
~2 |
|
λ |
ε + |
|
|
1 |
|
|
|
|
|
2 |
|
− |
σ |
|
|
1 ~ |
+ |
|
|
|
1 |
θ |
|
2 |
|
|
|
|
|
1 |
|
δ |
|
2 |
+ |
σ |
θ |
2 |
≤ |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
θ |
− |
|
|
|
|
|
|
δ |
|
|
|
|
|
|
|
|
|
|
θ |
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
2 |
|
|
|
2γ |
|
|
|
|
|
|
2 |
|
|
|
2λ |
|
|
|
|
|
|
|
|
γ |
|
|
2 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2λ |
|
|
|
|
|
|
2γ |
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
λ |
2 |
|
|
|
σ |
~2 |
|
|
1 |
|
2 |
|
|
|
|
σ |
|
|
2 |
|
|
|
λ |
|
2 |
|
|
|
|
σ |
~2 |
|
|
|
|
|
1 |
|
|
|
|
2 |
|
|
σ |
|
|
|
2 |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
≤ − 2 ε |
|
|
|
− |
|
θ |
+ |
|
δ |
|
+ |
|
|
|
|
|
θ |
|
|
≤ − |
|
2 ε |
|
− |
|
|
|
θ |
|
+ |
|
|
δ |
|
|
+ |
|
|
|
θ |
|
. |
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
2γ |
2λ |
|
|
2γ |
|
|
|
|
2γ |
|
|
|
2λ |
|
|
|
2γ |
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
Введем обозначение κ = min{λ, σ}. Тогда, считая λ, σ положительными, имеем: |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|
|
|
1 |
|
|
~ |
2 |
|
|
|
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
σ |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V& ≤ −κ |
|
ε |
|
|
|
|
− |
|
|
|
|
θ |
|
+ |
|
|
|
δ |
|
|
|
+ |
|
|
|
|
|
θ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
2γ |
2λ |
|
|
|
|
2γ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 / 2λ +σθ2 / 2γ |
V& ≤ −κV + |
|
|
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
где |
|
|
|
|
постоянная |
|
величина. Далее, применяя результат, |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
=δ |
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
полученный в ходе второго решения задачи, окончательно получаем: |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
V (ε(t)) ≤ e−κtV (0) + κ − κ e−κt |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
ε(t) ≤ |
|
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
2 e−κtV (0) |
|
κ |
κ |
e−κt . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Из последнего неравенства следует экспоненциальная сходимость ошибки |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
управления к ограниченному множеству с границей |
|
= |
|
2 |
|
/ κ . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Алгоритм |
управления |
|
(1.6), основанный |
|
на |
|
алгоритме |
|
адаптации |
(2.7), также |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
обеспечивает устойчивость замкнутой системы и является робастным по отношению к внешнему возмущению. В то же время алгоритм (1.6), (2.7) позволяет парировать недостатки робастного алгоритма управления (1.6), (2.4). Так, при отсутствии внешнего
возмущения верхняя граница может быть снижена до нуля за счет обнуления
коэффициента σ (т.н. гибридная σ -модификация). Кроме того, для уменьшения |
нет |
необходимости в значительном увеличении γ , которое влечет за собой рост амплитуды |
|
управляющего воздействия. Снижение |
можно обеспечить путем уменьшения |
коэффициента σ . |
|
Порядок выполнения работы |
|
12
1. На основе данных, приведенных в таблице 1.1, провести моделирование адаптивной системы управления, полученной в лабораторной работе №1, в условиях действия на объект возмущения вида:
δ(t) = 0.25sin10t +0.5 .
По результатам моделирования построить три графика. На первом вывести x и xm , на втором u , на третьем θ~ =θ −θˆ .
2. Заменить алгоритм адаптации (1.9) на нелинейную обратную связь (2.4) и повторить эксперимент для трех разных значений коэффициента γ . Определить влияние
этого коэффициента на свойства системы. По результатам моделирования для каждого γ построить два графика. На первом вывести x и xm , на втором u .
3. В п.1 заменить алгоритм адаптации (1.9) на робастную σ -модификацию (2.7). повторить эксперимент для трех разных значений коэффициента σ . Определить влияние этого коэффициента на свойства системы. По результатам моделирования для каждого σ построить два графика. На первом вывести x и xm , на втором u .
13
Лабораторная работа № 3. Адаптивное управление линейным многомерным объектом по состоянию.
Цель работы: освоение принципов построения адаптивной системы управления многомерным объектом.
Методические рекомендации. До начала работы студенты должны ознакомится с принципом построения алгоритмов адаптации на основе стандартной модели ошибки с измеряемым состоянием [Ник Уш]. Лабораторная работа рассчитана на 2 часа.
Теоретические сведения. Рассмотрим пример типовой задачи адаптивного управления многомерным объектом с использованием эталонной модели. При этом воспользуемся принципами решения аналогичной задачи для объекта первого порядка (см. лабораторную работу №1).
Постановка задачи
Дан объект управления, модель которого описывается уравнениями вида
x = Ax +b u , x(0) |
(3.1) |
& |
|
y = C x , |
(3.2) |
где x Rn вектор состояния, y регулируемая переменная,
|
|
0 |
|
1 |
|
|
0 |
|
0 |
|
|
|
||
A = |
|
M |
|
M |
|
|
|||
|
|
0 |
|
0 |
|
− a |
0 |
− a |
|
|
|
|
1 |
|
0 |
|
L |
0 |
|
|
0 |
|
|
|
||
1 |
|
L |
0 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
, C = [1 0 L |
0 0], |
|||||
O |
O M |
|
, b = |
|
M |
|
|||||
|
|
|
|||||||||
0 |
|
L |
1 |
|
|
|
0 |
|
|
|
|
− a |
2 |
L − a |
n−1 |
|
|
b |
|
|
|
||
|
|
|
|
|
|
0 |
|
|
|
||
ai , i = 0, n −1 неизвестные параметры, b0 известный параметр.
Задача управления заключается в компенсации параметрической неопределенности объекта и обеспечении следующего целевого равенства:
lim(xM (t) − x(t))= lim e(t) = 0 , |
(3.3) |
|
t→∞ |
t→∞ |
|
где e = xM − x вектор ошибки управления, |
xM Rn вектор, генерируемый |
|
эталонной моделью: |
|
|
xM = AM xM +bM g , |
(3.4) |
|
& |
|
|
yM = CM xM , |
|
(3.5) |
где g(t) задающее воздействие,
14
|
|
0 |
1 |
0 |
L |
0 |
|
|
||||
|
|
0 |
0 |
1 |
L |
0 |
|
|
||||
|
|
|
|
|||||||||
A |
= |
M |
M |
|
M |
O |
|
M |
|
, b |
||
M |
|
0 |
0 |
0 |
L |
1 |
|
M |
||||
|
|
|
|
|||||||||
|
− a |
M 0 |
− a |
M 1 |
− a |
M 2 |
L − a |
M n−1 |
|
|
||
|
|
|
|
|
|
|
|
|
||||
|
0 |
|
|
|
0 |
|
|
|
|
, CM = [1 0 L 0 0], |
|
= |
M |
|
|
|
0 |
|
|
|
|
|
|
a |
|
|
|
|
M 0 |
|
|
Коэффициенты эталонной модели aM i ,i =1, n −1 строятся на основе метода
стандартных характеристических полиномов [Боб, Мир + Григ Лук].
Решение задачи
Предполагая параметры объекта известными, синтезируем регулятор, который обеспечит условие (3.3) с заданными динамическими показателями качества временем переходного процесса tп и перерегулированием σ .
Для синтеза регулятора сформируем ошибку слежения ε = xM − x , возьмем ее производную и проведем ряд алгебраических преобразований:
|
|
|
|
0 |
1 |
0 |
|||
|
|
|
|
0 |
0 |
1 |
|||
e |
= xM − x = |
|
|||||||
|
M |
M |
|
M |
|||||
& |
& |
& |
|
||||||
|
|
|
|
0 |
0 |
0 |
|||
|
|
|
|
||||||
|
|
|
− a |
M 0 |
− a |
M 1 |
− a |
M 2 |
|
|
|
|
|
|
|
|
|||
L |
0 |
|
x |
M 1 |
|
L |
0 |
|
|
|
|
|
xM 2 |
|
|||
O |
M |
|
M |
|
|
L |
1 |
xM n−1 |
|
||
|
|
|
xM n |
|
|
|
|
|
|
||
L − aM n−1 |
|
|
|
||
|
0 |
|
|
0 |
|
|
|
|
+ |
M |
g − |
|
0 |
|
|
|
|
a |
|
|
|
M 0 |
|
0 |
|
1 |
|
|
0 L 0 |
|
x |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
0 |
|
0 |
|
|
1 |
|
L 0 |
|
|
|
|
1 |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
− |
M |
|
M |
|
|
M |
|
O |
M |
|
|
|
M |
|
− |
M |
u ± |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
0 |
|
0 |
|
|
0 L 1 |
|
x |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
n−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
− a1 |
|
− a2 |
|
|
|
|
x |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
− a0 |
|
L − an−1 |
|
|
|
|
b0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
0 |
|
|
1 |
|
|
0 |
L 0 |
|
x |
|
0 |
|
|
|
1 |
|
|
|
0 |
|
|
L 0 |
|
e |
|
|
|||||||||||||||||
|
0 |
|
|
0 |
|
|
1 |
L 0 |
|
|
|
|
x |
1 |
|
|
0 |
|
|
|
0 |
|
|
|
1 |
|
|
L 0 |
|
|
1 |
|
|
||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
e |
|
|||||||||||||||||||||||
± |
M |
|
|
M |
|
O |
O |
|
|
M |
|
|
|
|
|
|
|
|
= |
M |
|
|
|
M |
|
|
|
M |
|
|
O |
|
|
M |
|
|
2 |
|
+ |
||||
|
|
|
|
|
|
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
M |
|
|||||||||||||||||||
|
0 |
|
|
0 |
|
|
0 |
L 1 |
|
|
|
x |
n−1 |
|
|
0 |
|
|
|
0 |
|
|
|
0 |
|
|
L 1 |
|
e |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n−1 |
|
||||
|
|
|
− aM 1 − aM 2 |
|
|
|
|
|
|
|
|
x |
n |
|
|
|
− aM 1 |
− aM 2 |
|
|
|
|
|
|
e |
|
|
||||||||||||||||
− aM 0 |
L − aM n−1 |
|
|
|
|
− aM 0 |
|
L − aM n−1 |
|
n |
|
|
|||||||||||||||||||||||||||||||
|
|
0 |
|
|
|
|
0 |
|
0 |
|
|
|
|
|
L |
|
|
|
0 |
|
|
|
|
x |
|
0 |
0 |
|
|
|
|
|
|||||||||||
|
|
0 |
|
|
|
|
0 |
|
0 |
|
|
|
|
|
L |
|
|
|
0 |
|
|
|
|
|
1 |
|
|
0 |
|
|
0 |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
M |
|
|
|
|
M |
|
|
M |
|
|
|
|
|
O |
|
|
|
M |
|
|
|
|
|
M |
|
− |
M |
u + |
M |
|
g. |
|
|
|
|
|
|||||
|
|
0 |
|
|
|
|
0 |
|
0 |
|
|
|
|
|
L |
|
|
|
0 |
|
|
|
x |
|
|
|
0 |
|
|
0 |
|
|
|
|
|
|
|
||||||
|
|
|
+ a − a + a − a + a |
|
|
L − a |
|
|
+ a |
|
|
|
|
n−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
− a |
|
|
|
|
|
|
|
xn |
|
|
0 |
|
a |
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
M 0 |
|
0 |
|
|
M 1 |
1 |
|
M 2 |
|
|
|
2 |
|
|
|
|
|
|
M n−1 |
|
n−1 |
|
|
|
|
|
|
|
|
M |
0 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
||||||||||||||
Полученный результат запишем в более компактной форме:
e& = AM e + k[(− aM 0 + a0 )x1 + (− aM 1 + a1 )x2 +... + (− aM n−1 + an−1 )xn + aM 0 g −b0u], (3.6)
15
где k = [0 0 |
|
L 0 1]T . Введем обозначения: |
|
||||||||
|
|
θ1 = (− aM 0 + a0 ), |
θ2 = (− aM 1 + a1 ), ..., θn = (− aM n−1 + an−1 ), |
||||||||
после чего перепишем (3.6) в следующем виде: |
|
||||||||||
|
|
|
|
|
& |
AM e + k[θ |
T |
x + aM 0 g −b0u], |
(3.7) |
||
|
|
|
|
|
e = |
|
|||||
где θT = [θ |
,θ |
2 |
,...,θ |
n |
] вектор, содержащий параметры матрицы |
A . |
|||||
1 |
|
|
|
|
|
|
|
|
M |
||
Нетрудно видеть, что, если |
|
|
|
||||||||
|
|
|
|
|
u = |
|
1 |
[θT x + aM 0 g], |
(3.8) |
||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
то ошибка слежения подчиняется закону вида |
|
||||||||||
|
|
|
|
|
e = AM e |
|
|
(3.9) |
|||
|
|
|
|
|
& |
|
|
|
|
|
|
и экспоненциально затухает со временем.
Однако в исходной постановке задачи параметры матрицы AM неизвестны. Следовательно закон (3.8) физически нереализуем. Заменим в 3.7 неизвестные параметры θ на оценки θˆ и получим настраиваемый закон управления:
u = |
1 |
[θˆT x + aM 0 g]. |
(3.10) |
|
b |
||||
|
0 |
|
|
|
Подставим последнее выражение в (3.7) и получим модель ошибок: |
|
|||
& |
|
~T |
x , |
(3.11) |
e = |
AM e + kθ |
|||
где θ~ =θ −θˆ вектор параметрических ошибок.
Аналогично решению, преведенному в лабораторной работе №1, выбирается функция Ляпунова:
V = eT Pe + 21γ θ~Tθ~ ,
где P положительно определенная симметричная матрица, удовлетворяющая уравнению Ляпунова
AT |
P + PA |
= −Q , |
(3.12) |
M |
M |
|
|
где Q произвольно выбранная симметричная положительно определенная матрица.
Далее берется производная от функции Ляпунова, проводятся простейшие алгебраические преобразования и выводится алгоритм адаптации аналогичный алгоритму
(1.9):
& |
θˆ(0) = 0 |
(3.13) |
θˆ = γxkT Pe , |
||
|
16 |
|
где γ > 0 алгоритм адаптации.
Таким образом, алгоритм адаптивного управления состоит из настраиваемого регулятора (3.10), алгоритма адаптации (3.13), в котором матрица P находится из (3.12).
Адаптивный регулятор (3.10), (3.13) обеспечивает:
1.Устойчивость замкнутой системы;
2.Асимптотическое стремление ошибки e к нулю;
3.Ограниченность сигнала θˆ . При этом вектор θˆ асимптотически стремится к θ ,
если вектор x удовлетворяет условию неисчезающего возбуждения:
t+∫Tx(τ)xT (τ)dτ >αI , |
(3.14) |
t |
|
где α > 0 , T > 0 постоянные величины.
Условие (3.14) эквивалентно условию наличия не менее k / 2 гармоник в векторе
x .
4. При увеличении коэффициента γ до определенного оптимального значения
увеличивается скорость сходимости параметрических ошибок θ~ к нулю. Дальнейшее увеличение γ приведет к снижению скорости сходимости ошибки и появлению
колебаний в системе.
Порядок выполнения работы
1. На основе заданных в таблице 3.1 значений времени переходного процесса tп и максимального перерегулирования σ сформировать эталонную модель в форме (3.4), (3.5). Построить график переходной функции модели, на котором показать tп и σ ;
2. На основе предположения, что параметры объекта известны построить и промоделировать систему управления с регулятором (3.8). Провести три эксперимента при следующих условиях:
использовать расчетные значения параметров объекта θ1 и θ2 ;
незначительно отклонить параметры объекта так, чтобы система не потеряла устойчивость;
отклонить параметры объекта так, чтобы система потеряла устойчивость.
По результатам каждого эксперимента построить траектории x(t) и xM (t) на одном графике и e(t) на другом.
3.Провести моделирование адаптивной системы управления с регулятором (3.10)
иалгоритмом адаптации (3.13). В ходе моделирования проиллюстрировать свойства 1-4 алгоритма управления. Для этого необходимо провести эксперименты при следующих условиях:
повторить три эксперимента п.п. 2 для фиксированного значения γ ;
используя расчетные значения параметров объекта, провести эксперимент с тремя различными значениями γ ;
провести один из пердыдущих экспериментов данного пункта при g(t) =1 .
17
По результатам каждого эксперимента построить траектории x(t) и xM (t) на
одном графике, e(t) на втором, θ~ =θ −θˆ на третьем. 4. Сделать выводы по каждому пункту работы.
18
Таблица 3.1. Варианты заданий
|
|
|
Коэфф. |
Время |
Максимальное |
|
|
|
Матрица |
перерегули- |
|
||||
Вар. |
передачи |
переходного |
Сигнал задания g(t) |
||||
A |
рование |
||||||
|
b0 |
процесса, tп |
|
||||
|
|
|
σ, % |
|
|||
|
|
|
|
|
|
||
1 |
0 |
1 |
1 |
0.16 |
0 |
2sign[cost] +3 |
|
|
|
||||||
|
1 |
1 |
|
|
|
|
|
2 |
0 |
1 |
2 |
0.3 |
0 |
cos t +3sin 2t +5 |
|
|
|
||||||
|
1 |
−1 |
|
|
|
|
|
|
0 |
1 |
|
|
|
t |
|
3 |
|
|
3 |
0.9 |
0 |
∫sign[cos 2τ]dτ |
|
|
4 |
0 |
|
|
|
0 |
|
4 |
0 |
1 |
4 |
0.2 |
0 |
sign[sin 0.5t] + 2 |
|
|
|
||||||
|
1 |
4 |
|
|
|
|
|
5 |
0 |
1 |
5 |
0.6 |
0 |
10 cos 0.5t + 2sin t +12 |
|
|
|
||||||
|
− 2 |
5 |
|
|
|
|
|
|
0 |
1 |
|
|
|
t |
|
6 |
6 |
0.3 |
0 |
∫sign[sin 0.5τ]dτ + 3 |
|||
|
|
||||||
|
9 |
−1 |
|
|
|
0 |
|
7 |
0 |
1 |
7 |
0.7 |
0 |
0.5sign[sin 0.7t] +1 |
|
|
|
||||||
|
− 2 |
3 |
|
|
|
|
|
8 |
0 |
1 |
8 |
0.1 |
0 |
5sin 0.5t + 4 cos 0.1t +8 |
|
|
|
||||||
|
2 |
2 |
|
|
|
|
|
9 |
0 |
1 |
9 |
0.9 |
0 |
sign[cost] +3 |
|
|
|
||||||
|
0 |
0 |
|
|
|
|
|
|
0 |
1 |
|
|
|
t |
|
10 |
|
|
1 |
3.5 |
15 |
∫sign[sin 2τ]dτ |
|
|
0 |
−1 |
|
|
|
0 |
|
11 |
0 |
1 |
2 |
0.6 |
15 |
sin 5t + 0.5cos 0.2t + 2 |
|
|
|
||||||
|
−1 6 |
|
|
|
|
||
12 |
0 |
1 |
3 |
0.9 |
15 |
2sign[sin 0.4t]+3 |
|
|
|
||||||
|
10 |
1 |
|
|
|
|
|
|
0 |
1 |
|
|
|
t |
|
13 |
|
|
4 |
0.4 |
15 |
3∫sign[sin 0.5τ]dτ +3 |
|
|
−1 |
5 |
|
|
|
0 |
|
14 |
0 |
1 |
5 |
0.2 |
15 |
2sin 0.2t +sin 0.1t +8 |
|
|
|
||||||
|
9 |
0 |
|
|
|
|
|
15 |
0 |
1 |
6 |
0.5 |
15 |
4sign[cos 2t]+5 |
|
|
|
||||||
|
−9 |
1 |
|
|
|
|
|
|
0 |
1 |
|
|
|
t |
|
16 |
|
|
7 |
0.5 |
15 |
3∫sign[cos 0.2τ]dτ +3 |
|
|
−7 |
6 |
|
|
|
0 |
|
19
17 |
0 |
1 |
8 |
0.45 |
15 |
sign[sin 0.5t] + 2 |
|
|
|
||||||
|
−1 7 |
|
|
|
|
||
18 |
0 |
1 |
9 |
0.15 |
15 |
0.8sin 2t + cos 0.8t + 2 |
|
|
|
||||||
|
−10 |
6 |
|
|
|
|
|
|
0 |
1 |
|
|
|
t |
|
19 |
1 |
0.7 |
15 |
∫sign[sin 0.3τ]dτ +1.5 |
|||
|
|
||||||
|
2 |
2 |
|
|
|
0 |
|
20 |
0 |
1 |
2 |
1.2 |
0 |
7sign[cos 0.9t] +8 |
|
|
|
||||||
|
−15 |
6 |
|
|
|
|
|
21 |
0 |
1 |
3 |
1.5 |
0 |
0.4sin 3t + cos 0.1t +1.5 |
|
|
|
||||||
|
11 |
1 |
|
|
|
|
|
|
0 |
1 |
|
|
|
t |
|
22 |
4 |
0.8 |
0 |
6∫sign[sin 0.1τ]dτ +9 |
|||
|
|
||||||
|
3 |
3 |
|
|
|
0 |
|
23 |
0 |
1 |
5 |
0.9 |
0 |
2sign[sin t] + 4 |
|
|
|
||||||
|
12 |
6 |
|
|
|
|
|
24 |
0 |
1 |
6 |
0.2 |
0 |
9 sin 0.2t +9 cos 0.1t +15 |
|
|
|
||||||
|
0 |
11 |
|
|
|
|
|
|
0 |
1 |
|
|
|
t |
|
25 |
7 |
0.5 |
0 |
4∫sign[sin 6τ]dτ +5 |
|||
|
|
||||||
|
13 |
2 |
|
|
|
0 |
|
26 |
0 |
1 |
8 |
0.9 |
0 |
4sign[cost] + 3 |
|
|
|
||||||
|
15 |
6 |
|
|
|
|
|
27 |
0 |
1 |
9 |
1.3 |
0 |
sin 0.1t + cos 5t + 2 |
|
|
|
||||||
|
13 |
−1 |
|
|
|
|
|
|
0 |
1 |
|
|
|
t |
|
28 |
1 |
1.6 |
0 |
9∫sign[sin 0.1τ]dτ +12 |
|||
|
|
||||||
|
4 |
4 |
|
|
|
0 |
|
29 |
0 |
1 |
2 |
0.75 |
15 |
3sign[sin 4t]+8 |
|
|
|
||||||
|
9 |
7 |
|
|
|
|
|
30 |
0 |
1 |
3 |
0.65 |
15 |
7 sin 0.3t +8cos t + 20 |
|
|
|
||||||
|
−1 |
9 |
|
|
|
|
|
20
