Герасимов Адаптивное управление
.pdfСанкт-Петербургский государственный научно-исследовательский университет информационных технологий, механики и оптики
Герасимов Д.Н., Никифоров В.О.
Синтез адаптивного управления линейными системами
Санкт-Петербург
2014
В пособии приводятся методические рекомендации к выполнению лабораторных работ по курсу “Адаптивное и робастное управление”. Лабораторные работы направлены на изучение методов построения адаптивных наблюдателей состояния линейных параметрически неопределенных объектов и методов синтеза следящих систем управления по состоянию в условиях параметрической неопределенности модели объекта управления, а также генераторов задающего и возмущающего воздействий.
Пособие служит для выполнения лабораторного практикума при подготовке инженеров по специальности 220201.65 «Управление и информатика в технических системах», бакалавров и магистров по направлению: 220200.62 «Автоматизация и управление».
2
Содержание |
|
Введение....................................................................................................................................... |
4 |
Лабораторная работа №1: Принципы построения систем адаптивного управления |
|
невозмущенными объектами...................................................................................................... |
5 |
Лабораторная работа №2: Принцип построения систем адаптивного и робастного |
|
управления возмущенными объектами..................................................................................... |
9 |
Лабораторная работа № 3. Адаптивное управление линейным многомерным объектом |
|
по состоянию. ............................................................................................................................. |
14 |
Лабораторная работа № 3. Адаптивное управление линейным многомерным объектом |
|
по состоянию. ............................................................................................................................. |
14 |
Лабораторная работа № 4. Робастное управление линейным многомерным объектом по |
|
состоянию. .................................................................................................................................. |
21 |
Лабораторная работа №5: Параметризация модели объекта управления (способ №1) .. |
23 |
Лабораторная работа №6: Синтез адаптивного наблюдателя состояния линейного |
|
объекта........................................................................................................................................ |
28 |
Лабораторная работа №7: Параметризация модели объекта управления (способ 2). |
|
Адаптивное управление объектом по выходу. ....................................................................... |
30 |
Лабораторная работа №8: Адаптивноеуправление линейным объектом по выходу – |
|
метод расширенной ошибки. .................................................................................................... |
34 |
Лабораторная работа №9:Адаптивный наблюдатель внешнего возмущения |
|
(Никифоров). .............................................................................................................................. |
37 |
Лабораторная работа № 10. Адаптивная компенсация внешнего возмущения. .............. |
38 |
Лабораторная работа № 11: Адаптивное слежение за неизвестным задающим |
|
воздействием. ............................................................................................................................. |
39 |
Лабораторная работа № 12: Адаптивное слежение за неизвестным задающим |
|
воздействием с компенсацией неизвестного возмущения..................................................... |
40 |
3
Введение
Целью работ является освоение методов адаптивного и робастного управления, основанных на аппарате функций Ляпунова. В ходе выполнения работ предусматривается изучение в качестве примеров простейших алгоритмов адаптивного и робастного управления объектами первого порядка, на основе которых последовательно синтезируются блоки регулятора, обеспечивающие компенсацию параметрической неопределенности модели многомерного объекта и моделей генераторов задающего и возмущающего воздействий. В отдельную работу вынесен метод построения адаптивного наблюдателя вектора состояния параметрически неопределенного объекта.
Перед выполненим лабораторных работ студентам необходимо освоить базовые курсы теорий линейных и нелинейных систем.
Для выполнения лабораторных работ необходимо предварительно изучить теоретический материал, относящийся к работе и освоить пакет программ MatLab версии не ниже 7.01 со средой моделирования Simulink версии не ниже R14.
Выполнение каждой работы включает проведение компьютерного моделирования, создание в печатной форме отчета (пояснительной записки) и защиту работы. Защита предполагает индивидуальную дискуссию преподавателя со студентом, в ходе которой выявляется качество приобретенных студентом знаний и навыков.
Отчет необходимо оформлять согласно правилам ГОСТ 2.105-95 ЕСКД . Отчет должен содержать:
1.Титульный лист, выполненный согласно образцу, представленному в приложении 1;
2.Цель работы;
3.Основные теоретические выкладки: формулы, диаграммы, схемы и т.п. с пояснениями;
4.Расчетную часть, в которой согласно порядку выполнения работы приводятся расчеты на основе формул либо листингов программ с пояснениями, результаты моделирования в виде схем и графиков. В обязательном порядке подписываются схемы, таблицы, листинги программ, скриншоты и т.п, а также оси графиков и диаграмм;
5.Выводы по выполненной работе.
4
Лабораторная работа №1: Принципы построения систем адаптивного управления невозмущенными объектами
Цель работы: освоение принципов построения систем адаптивного управления на основе решения задачи слежения выхода объекта первого порядка за эталонным сигналом.
Методические рекомендации. До начала работы студенты должны ознакомится с анализом устойчивости нелинейных систем методом функций Ляпунова. Лабораторная работа рассчитана на 2 часа.
Теоретические сведения. Рассмотрим решение задачи адаптивного слежения, начиная с постановки, и заканчивая структурной схемой закмнутой системы.
Постановка задачи
Дан объект, представленный моделью вида:
x =θ x +u , |
(1.1) |
& |
|
где x переменная состояния объекта, u |
управляющая переменная, θ |
неизвестный постоянный параметр.
Цель управления заключается в компенсации неопределенности θ и обеспечении следующего целевого равенства:
lim(xm (t) − x(t))= limε(t) = 0 , |
(1.2) |
|
t→∞ |
t→∞ |
|
где ε = xm − x ошибка управления, xm |
переменная состояния эталонной модели |
|
вида: |
|
|
xm = −λ xm + λ g |
|
(1.3) |
& |
|
|
где g эталонный сигнал, λ > 0 параметр, задающий время переходного процесса
эталонной модели.
Решение задачи
Положим, что параметр θ известен. Тогда, взяв производную ошибки ε по времени, с учетом (1.1) и (1.3), получим:
ε = xm − x = −λ xm + λ g −θ x −u . |
(1.4) |
||
& |
& |
& |
|
Для экспоненциального стремления ошибки к нулю должно выполняться равенство:
ε& = −λε ,
которое с учетом (1.4) позволяет вывести закон управления
ε& = −λ xm + λ g −θ x −u = −λε ,
или
−λ xm + λ g −θ x −u = −λε = −λxm + λx .
5
Откуда следует |
|
u = −θx −λx + λ g . |
(1.5) |
Пусть теперь, как в исходной постановке задачи, параметр θ |
неизвестен. Тогда |
для реализуемости закона управления (1.5) заменим величину θ на ее оценку θˆ : |
|
u = −θˆx −λx + λ g . |
(1.6) |
Таким образом, задача синтеза регулятора сводится к задаче нахождения функции θˆ(t) , которая обеспечит устойчивость замкнутой системы и цель управления (1.2). Для
нахождения θˆ(t) проведем ряд алгебраических преобразований и воспользуемся методом функций Ляпунова. Подставим последнее выражение в модель (1.1):
x& =θ x −θˆx −λx + λ g =θ~ x −λx + λ g ,
где θ~ =θ −θˆ параметрическая ошибка. Сформируем модель ошибки управления:
ε& = x&m − x& = −λ xm + λ g −θ~ x + λx −λ g = −λε −θ~ x
или |
|
|
|
|
|
~ |
|
|
|
|
|
|
|
||
|
|
|
& |
|
x. |
|
|
|
|
(1.7) |
|||||
|
|
|
|
ε = −λε −θ |
|
|
|
|
|
||||||
Выберем функцию Ляпунова в форме |
|
|
|
|
|
|
|
|
|||||||
|
|
|
1 |
2 |
|
|
1 |
~2 |
|
|
|
|
|||
|
|
|
|
V = 2 ε |
|
+ |
|
θ |
|
, |
|
|
(1.8) |
||
|
|
|
|
2γ |
|
|
|
||||||||
где γ > 0 параметр учетом (1.7) возьмем ее производную по времени: |
|
||||||||||||||
|
1 |
|
1 ~~& |
|
~ |
|
|
|
1 ~ & |
2 ~ |
~ 1 |
& |
|||
& |
|
|
|
2θθ = ε(−λε −θ |
x)− |
|
ˆ |
|
|
ˆ |
|||||
V |
= 2 |
2εε + |
2γ |
γ θθ = −λε −θ xε −θ γ |
θ . |
||||||||||
|
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
Исходя из условия асимптотической устойчивости замкнутой системы или положения равновесия ε* = 0 , что эквивалентно цели (1.2), алгоритм, генерирующий функцию θˆ(t) имеет вид:
& |
(1.9) |
θˆ = −γxε . |
Нетрудно видеть, что производная функции Ляпунова будет отрицательной:
V& = −λε2 < 0 ε ≠ 0.
Алгоритм, формирующий оценку θˆ называется алгоритмом адаптации, а параметр γ коэффициентом адаптации.
Таким образом, алгоритм адаптивного управления, обеспечивающий цель (1.2), представляется настраиваемым регулятором (1.6) и алгоритмом адаптации (1.9). Алгоритм наделяет замкнутую систему следующими свойствами:
1.Устойчивость замкнутой системы;
2.Асимптотическое стремление ошибки ε к нулю;
6
3. Асимптотическое стремление θˆ к θ , что следует из свойства 2 и выражения (1.7). Следует отметить, что, если величины θˆ и θ векторные, то в общем случае для стремления θˆ к θ необходимо выполнение условия неисчезающего возбуждения [].
4. При повышении до определенного оптимального значения коэффициента γ
повышается скорость сходимости параметрической ошибки θ~ к нулю. Дальнейшее увеличение γ приведет к снижению скорости сходимости ошибки и появлению
колебаний в системе.
Порядок выполнения работы
1. На основе заданных в таблице 1.1 значений параметров объекта (1.1) и эталонной модели (1.3) построить неадаптивную систему управления на базе регулятора (1.5). Провести моделирование системы в условиях скачкообразного изменения параметра θ на величину равную 10-и и построить два графика. На одном вывести x и xm , на другом u .
2. Повторить эксперимент для адаптивной системы управления, замкнутой регулятором (1.6) с алгоритмом адаптации (1.9). Коэффициент адаптации γ выбрать
экспериментально. Дополнительно построить график параметрической ошибки θ~ =θ −θˆ .
3.Провести моделирование адаптивной системы управления при разных коэффициентах γ .
4.Сделать выводы на основе сравнения результатов моделирования.
Таблица 1.1. Варианты заданий
|
Параметр |
Параметр |
Сигнал задания |
|
Параметр |
Параметр |
|
|
эталонной |
|
эталонной |
||||
Вар. |
объекта |
модели |
|
g(t) |
Вар. |
объекта |
модели |
|
θ |
λ |
|
|
|
θ |
λ |
|
|
|
|
|
|
||
1 |
6 |
1 |
|
sin 2t |
16 |
7 |
1 |
2 |
2 |
2 |
|
cos 4t |
17 |
8 |
2 |
3 |
4 |
3 |
sign[sin 3t] + 2 |
18 |
9 |
3 |
|
4 |
1 |
4 |
|
2 |
19 |
1 |
4 |
5 |
5 |
5 |
|
2 sin 5t |
20 |
2 |
5 |
6 |
8 |
6 |
|
1 |
21 |
3 |
6 |
7 |
0 |
7 |
|
7 sin t |
22 |
4 |
7 |
8 |
9 |
8 |
|
5sin 0.5t + 6 |
23 |
5 |
8 |
9 |
10 |
9 |
|
sign[cost] +3 |
24 |
6 |
9 |
|
|
|
t |
|
|
|
|
10 |
15 |
1 |
∫sign[sin 2τ]dτ |
25 |
7 |
10 |
|
|
|
|
0 |
|
|
|
|
11 |
0 |
2 |
sin 5t + 0.5cos 0.2t + |
26 |
8 |
1 |
|
|
|
|
|
|
|
|
|
12 |
8 |
3 |
2sign[sin 0.4t]+3 |
27 |
0 |
2 |
|
13 |
6 |
4 |
3∫t |
sign[sin 0.5τ]dτ |
28 |
1 |
3 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
14 |
7 |
5 |
2 sin 0.2t +sin 0.1t + |
29 |
2 |
4 |
|
|
|
|
|
7 |
|
|
|
Сигнал задания g(t)
cos t +3
4
cos 3t + 2
3sin t
0.5cos8t +1
sin 2t +6
9
cos t + 4
9 sin 0.5t +12
sin 8t + 2
10
cos 2t + 2
sign[cos t] +8
3sign[sin t]+ 4
|
|
|
|
|
|
|
|
15 |
11 |
6 |
4sign[cos 2t]+5 |
30 |
3 |
5 |
3sin t +10 |
8
Лабораторная работа №2: Принцип построения систем адаптивного и робастного управления возмущенными объектами
Цель работы: освоение принципов построения систем адаптивного и робастного управления на основе решения задачи слежения выхода возмущенного объекта первого порядка за эталонным сигналом.
Методические рекомендации. До начала работы студенты должны ознакомится с анализом устойчивости нелинейных систем методом функций Ляпунова. Лабораторная работа рассчитана на 2 часа.
Теоретические сведения. Рассмотрим пример задачи адаптивного и робастного слежения выхода возмущенного объекта за эталонным сигналом. Приведем несколько решений поставленной задачи. При этом воспользуемся результатами, приведенными в лабораторной работе №2.
Постановка задачи
Дан объект, представленный моделью вида
x =θ x +u +δ , |
(2.1) |
&
где δ ограниченное возмущение, удовлетворяющее неравенству δ(t) ≤δ .
Цель управления заключается в обеспечении следующего целевого неравенства:
xm (t) − x(t) |
|
≤ , |
t ≥T , |
(2.2) |
|
где , T точность работы системы управления и время ее настройки соответственно, xm (t) эталонный сигнал, генерируемый моделью (1.2). Предполагается, что параметры
и T можно изменять в соответствии с требованиями, предъявляемым к системе. Приведем три решения поставленной задачи.
Решение №1
В качестве решения расмотрим просинтезированный в работе №1 регулятор (1.6) и (1.9) и аналогичным образом построим модель ошибок:
& |
~ |
(2.3) |
ε = −λε −θ x −δ. |
Далее проведем анализ устойчивости замкнутой системы с помощью функции Ляпунова (1.8) с учетом последнего выражения и алгоритма адаптации (1.9).
|
1 |
|
|
1 ~~& |
|
|
|
|
|
~ |
|
|
|
|
1 |
~ & |
|
|
|
2 |
|
~ |
|
|
~ 1 |
|
|
|
|
|||||
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ˆ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
& |
+ 2γ 2θθ |
= ε(−λε −θ x −δ )− γ θθ = −λε −θ xε −δε +θ γ γxε = |
|
||||||||||||||||||||||||||||
V = 2 2εε |
|
|||||||||||||||||||||||||||||||||
= −λε |
2 |
− |
δε = − |
λ |
ε |
2 |
− |
λ |
ε |
2 |
−δε |
± |
1 |
δ |
2 |
= − |
λ |
ε |
2 |
|
|
λ |
ε + |
1 |
|
2 |
1 |
δ |
2 |
≤ |
||||
|
|
|
|
|
|
|
|
|
− |
|
|
|
δ |
+ |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
2 |
|
|
|
2 |
|
|
|
|
2λ |
|
|
|
|
2 |
|
|
|
|
2 |
|
2λ |
|
|
2λ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
≤ − λ2 ε2 + 21λ δ 2 ≤ − λ2 ε2 + 21λ δ 2
9
Из полученного неравенства следует асимптотическое стремление ошибки ε к некоторому ограниченному множеству, определяемому верхней границей сигнала
возмущения δ и параметром λ . При этом точность системы управления может быть увеличена путем увеличения λ . Однако из приведенного анализа не следует
ограниченности сигнала θˆ . Если продолжить анализ и предположить, что в случае, когда переменная x отделена от нуля (например, при решении задачи стабилизации), а ошибка ε уходит в окрестность и остается в ней на постоянном ненулевом уровне, то
θ&ˆ = −γxустεуст = C = const ,
откуда следует
θˆ = Ct
и неограниченный рост оценки θˆ с течением времени.
Таким образом, представленный регулятор (1.6) и (1.9) в общем случае не является робастным по отношению к внешнему возмущенению, что мотивирует необходимость модификации алгоритма (1.9).
Решение №2
Представим модификацию алгоритма (1.9) в следующей статической форме:
θˆ = −γxε , |
(2.4) |
Статическая зависимость гарантирует ограниченность сигнала θˆ в случае, если ошибка ε ограничена. Проанализируем устойчивость замкнутой системы управления со статической связью (2.4) и регулятором (1.6). Выберем функцию Ляпунова
V = |
1 |
ε2 |
(2.5) |
|
2 |
|
|
и возьмем ее производную. Учитывая (2.4) и модель ошибки (2.3), проведем алгебраические преобразования:
& |
1 |
2εε |
= −λε |
2 |
|
~ |
|
|
|
|
λ |
ε |
2 |
− |
λ |
|
ε |
2 |
− |
ˆ |
|
|
|
|
|
||||||||||
V = |
2 |
|
−θ |
xε −δε |
= − 2 |
|
4 |
|
(θ −θ)xε −δε = |
||||||||||||||||||||||||||
|
|
|
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= − |
λ |
ε2 − |
λ |
ε2 − |
δε ± |
1 |
δ 2 −(θ +γxε)xε = |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
2 |
2 |
λ |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
= − |
λ |
ε |
2 |
− |
|
|
λ |
ε |
+ |
1 |
|
2 |
1 |
|
δ |
2 |
± |
θ2 |
−θxε |
2 |
ε |
2 |
= |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
δ |
+ |
|
|
|
|
|
−γx |
|
|
|
|
||||||||||||||||
|
2 |
|
|
|
|
|
|
2 |
|
|
|
2λ |
|
|
2λ |
|
|
|
|
4γ |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
= − |
λ |
ε |
2 |
− |
|
|
λ |
ε + |
1 |
|
2 |
1 |
|
δ |
2 |
|
|
|
|
θ |
|
+ |
γ xε |
2 |
+ |
θ |
2 |
≤ |
|||||||
|
|
|
|
|
|
|
|
δ |
+ |
|
|
|
|
− |
|
|
|
|
|
|
|
||||||||||||||
|
2 |
|
|
|
|
|
|
2 |
|
|
|
2λ |
|
|
2λ |
|
|
|
|
|
2 |
γ |
|
|
|
|
|
4γ |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
≤ − |
λ |
ε2 + |
1 |
|
δ 2 |
+ |
θ2 |
= −λV |
+ |
1 |
δ 2 |
+ θ2 |
|
= −λV + , |
|
|
|
|
|||||||||||||||||
|
2 |
|
|
|
|
2λ |
|
|
|
|
4γ |
|
|
|
2λ |
|
|
|
|
4γ |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
|