Аналитическое конструирование оптимальных регуляторов. Практикум для лабораторных и практических занятий
.pdf
Производную по u от фигурной скобки приравниваем к нулю и находим оптимальное управление u* (t) в функции от ∂s / ∂x 2
|
|
|
|
¶ |
|
|
|
{...}= 2ru* + |
|
¶s |
= 0, |
|
|
|
|
|
|||||||||||||||||||
откуда |
|
|
|
¶u |
¶x2 |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
¶s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
u * = - |
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
(12) |
||||||||||||||||
|
|
2r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
¶x 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Подставив это значение в уравнение Беллмана, получим |
|||||||||||||||||||||||||||||||||||
|
2 |
2 |
|
1 |
|
|
|
|
|
∂s |
|
2 |
|
|
|
∂s |
|
|
|
|
|
1 |
|
|
∂s |
|
2 |
|
|||||||
x |
|
+qx - |
|
|
|
( |
|
|
|
|
) |
|
|
+ |
|
|
|
|
|
x |
|
- |
|
|
( |
|
|
) |
|
= 0 |
|||||
|
4r |
∂x |
|
|
|
∂x |
|
|
2r |
∂x |
|
|
|||||||||||||||||||||||
1 |
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|||||
или, приведя подобные члены, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
2 |
2 |
|
1 |
|
|
|
|
|
¶s |
|
|
|
2 |
|
|
|
|
¶s |
|
|
|
|
|
|
|
|
|
|
||||||
x1 |
+ qx2 |
- |
|
|
|
|
|
( |
|
|
|
|
) |
|
+ |
|
|
|
|
|
|
x2 = 0. |
|
|
|
(13) |
|||||||||
|
4r |
¶x2 |
|
|
¶x1 |
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Решение уравнения Беллмана (13) будем искать в виде квадра-
тичной формы |
|
S(X) = k11 x12 + k12 x1x 2 + k 22 x 22 . |
(14) |
Вычислив производные и подставив их в уравнение Беллмана, путем приравнивая к нулю всех коэффициентов при фазовых коор- динатах, получим систему уравнений для определения коэффици- ентов k11, k12 и k22:
ì |
k2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ï1- |
12 |
|
= 0; |
|
|
|
ì |
|
= 2 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ï |
4r |
|
|
|
|
|
|
|
ïk |
12 |
|
|
r ; |
||||||||||
ï |
|
|
|
|
1 |
|
|
2 |
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
- |
|
|
= 0; |
22 = |
|
r(q + 2 r ); |
|||||||||||||||
íq + k12 |
|
|
k22 |
ík |
|
||||||||||||||||||
r |
|||||||||||||||||||||||
ï |
|
|
1 |
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
ï |
|
|
k12 |
× k22 = |
ïk |
|
= |
q + 2 |
r |
. |
|||||||||||||
ï2k11 - |
r |
0; î |
11 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
î |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Следовательно, функция Беллмана принимает вид:
S(X) = 
q + 2
r × x12 + 2
r × x1 × x 2 + 
r(q + 2
r ) × x 22 ,
а оптимальное управление
u* (X) = -(k1x1 + k2 x2 ),
(15)
(16)
(17)
21
где |
k1 |
= |
1 |
|
; k2 = |
q + |
2 r |
. |
||
|
|
|
|
|
|
|||||
|
|
r |
||||||||
|
|
|
|
r |
|
|||||
Структурная схема оптимальной системы управления представ- лена на рисунке 1а, схема моделирования, составленная с использо- ванием пакета MATLAB, – на рис. 1б.
u |
ОУ |
х |
|
|
|
|
|
|
УУ
а)
∞
б I = ò(х12 + qх22 + ru2 )dt → min
0
Рис. 1. Структурная схема оптимальнойсистемы управления
22
Перечень используемого оборудования
Перечень используемого оборудования приведен в лаб.работе №1.
Указания по технике безопасности
Указания по технике безопасности приведены в лаб.работе №1.
Указания по порядку выполнения работы
Исследовать оптимальную систему управления, структурная схема которой представлена на рисунке 1. Объект управления опи-
сывается системой в нормальной форме
x1/ = x 2 ; x 2/ = u,
где u = -(k 1x1 + k 2 x 2 ), |
k1 |
= |
1 |
, |
k2 = |
1 |
(q+ 2 |
|
). |
|||
r |
||||||||||||
|
|
r |
|
|||||||||
|
|
|
||||||||||
|
|
|
|
r |
|
|
|
|
||||
Постоянные r и q являются коэффициентами квадратичной фор-
мы в интегральном критерии
∞
I = ò(x12 + qx 22 + ru 2 )dt.
0
Их численные значения приведены в таблице 1. Номер варианта задается устно преподавателем.
Исследование проводится в следующей последовательности:
1.Рассчитать для заданного варианта коэффициенты k1 и k2 (этот пункт выполняется при подготовке).
2.Собрать структурную схему системы управления (рис. 1) с использованием пакета MATLAB.
3.Установить начальные условия х(0) = x1(0) = 1 и х/(0) = x1/(0) =
x2(0) = 0.
4.Определить перерегулирование и время регулирования при каждой паре постоянных q и r.
5.Найти среди исследованных значений коэффициентов k1 и k2
такие значения k1* и k *2 , при которых время регулирования полу- чается минимальным, а перерегулирование – максимальным.
Содержание отчета
Отчет о выполнении лабораторной работы должен содержать:
1. Исходные данные, расчетные формулы, структурную схему системы управления и схему моделирования.
23
2.Решение задачи синтеза оптимальной системы управления методом динамического программирования.
3.Таблицы со значениями k1 и k2, соответствующими заданному варианту постоянных r и q.
4.Графики переходных процессов, полученные при каждой па- ре постоянных r и q.
5.Выводы о зависимости прямых показателей качества tp, σ и
характера переходных процессов от коэффициентов квадратичного
критерия.
Контрольные вопросы
1.Что такое оптимальная автоматическая система управления?
2.Что такое ограничение, краевые условия и критерий опти- мальности?
3.Как ставится задача синтеза оптимальной автоматической си- стемы управления?
4.Как формулируется принцип оптимальности?
5.Как записывается уравнение Беллмана и что оно выражает?
6.Как определяется оптимальный алгоритм управления линейным объектом с интегральным квадратичным критерием оптимальности?
Литература
1.Ким Д. П. Сборник задач по теории автоматического управ- ления: Многомерные, нелинейные, оптимальные и адаптивные си- стемы: учебное пособие. М.: Физматлит, 2008. 328 с.
2.Аттетков А. В., Галкин С. В., Зарубин В. С. Методы оптими- зации: учебник для вузов. М.: Изд-во МГТУ, 2001. 440 с.
3.Афанасьев В. Н., Колмановский В. Б., Носов В. Р. Математи- ческая теория конструирования систем управления. М.: Высшая школа, 1998. 574 с.
4.Ванько В. И., Ермошина О. В., Кувыркин Г. Н. Вариационное исчисление и оптимальное управление: учебник для вузов / под ред. В. С. Зарубина, А. П. Крищенко. 2-е изд. М.: Изд-во МГТУ, 2001. 488 с.
5.Иванов В. А., Фалдин Н. В. Теория оптимальных систем автома- тического управления: учеб. пособие для вузов. М.: Наука, 1981. 331 с.
6.Куропаткин П. В. Оптимальные и адаптивные системы: учеб. пособие для вузов. М.: Высш. школа, 1980. 287 с.
7.Пантелеев А. В., Бортаковский А. С., Летова Т. А. Оптималь- ное управление в примерах и задачах: учеб. пособие. М.: Изд-во МАИ, 1996. 212 с.
24
ПРИЛОЖЕНИЕ
Таблица 1
Значение постоянных r и q интегрального квадратичного критерия
|
|
Значение постоянных r и q интегрального |
|
|||||
№ |
|
|
квадратичного критерия |
|
|
|||
варианта |
Эксперим. 1 |
Эксперим. 2 |
Эксперим. 3 |
Эксперим. 4 |
||||
|
r |
q |
r |
q |
r |
q |
r |
q |
0 |
0.04 |
0.05 |
1 |
1 |
10 |
10 |
50 |
50 |
1 |
0.04 |
0.09 |
2 |
5 |
8 |
8 |
40 |
40 |
2 |
0.04 |
0.08 |
0.2 |
0.2 |
4 |
4 |
20 |
20 |
3 |
0.04 |
0.07 |
0.1 |
0.1 |
2 |
2 |
10 |
10 |
4 |
0.02 |
0.008 |
0.3 |
0.3 |
3 |
3 |
30 |
30 |
5 |
0.03 |
0.03 |
0.6 |
0.6 |
6 |
6 |
60 |
60 |
6 |
0.05 |
0.04 |
0.5 |
0.5 |
5 |
5 |
25 |
25 |
7 |
0.06 |
0.06 |
0.4 |
0.4 |
14 |
14 |
56 |
56 |
8 |
0.07 |
0.07 |
0.7 |
0.7 |
7 |
7 |
70 |
70 |
9 |
0.08 |
0.4 |
0.8 |
0.8 |
12 |
12 |
80 |
80 |
10 |
0.09 |
0.6 |
0.9 |
0.9 |
9 |
9 |
90 |
90 |
11 |
0.05 |
0.001 |
4.0 |
4.0 |
16 |
16 |
64 |
64 |
12 |
0.01 |
0 |
0.6 |
0.9 |
8 |
10 |
32 |
32 |
13 |
0.02 |
0.01 |
0.8 |
0.6 |
9 |
8 |
45 |
45 |
14 |
0.03 |
0.07 |
0.9 |
0.1 |
5 |
4 |
55 |
55 |
15 |
0.04 |
0.06 |
0.4 |
0.2 |
10 |
2 |
75 |
75 |
16 |
0.05 |
0.05 |
0.2 |
0.4 |
2 |
1 |
26 |
26 |
17 |
0.06 |
0.04 |
0.1 |
0.3 |
1 |
9 |
15 |
15 |
18 |
0.07 |
0.1 |
0.3 |
0.5 |
6 |
3 |
36 |
36 |
19 |
0.08 |
0.3 |
0.5 |
0.7 |
3 |
4 |
27 |
27 |
20 |
0.09 |
0.5 |
0.7 |
0.8 |
4 |
5 |
62 |
62 |
21 |
0.02 |
0 |
4.0 |
2.0 |
18 |
18 |
72 |
72 |
22 |
0.03 |
0 |
6.0 |
3.0 |
12 |
16 |
86 |
86 |
23 |
0.04 |
0 |
9.0 |
4.0 |
20 |
22 |
94 |
94 |
24 |
0.05 |
0 |
8.0 |
5.0 |
15 |
15 |
95 |
95 |
25
3. СИНТЕЗ ОПТИМАЛЬНОГО ПД-РЕГУЛЯТОРА
Цели: получить оптимальный закон управления методом дина- мического программирования; исследовать качество процесса регу- лирования в системе с оптимальным ПД-регулятором.
Формируемые компетенции (или их части)
-ОК-1 – способность совершенствовать и развивать свой интел- лектуальный уровень;
-ОК-8 – способность самостоятельно приобретать с помощью информационных технологий и использовать в практической дея- тельности новые знания и умения;
-ПК-21 – способность выбирать оптимальные решения;
-ПК-39 – способность проводить математическое моделирова- ние систем автоматизации с использованием современных техноло- гий проведения научных исследований.
Теоретическая часть
Сущность теории аналитического конструирования регуляторов состоит в том, что закон изменения управляемой координаты и управляющего воздействия задается аналитически в виде некоторо- го функционала J(X,U), который нужно минимизировать подбором управляющего воздействия. Зная алгоритм управления U(t) или U(X), конструируют регулятор.
Особое значение при аналитическом конструировании регулято- ров принимает выбор критериев оптимальности. Общих критериев оптимальности для любых объектов управления быть не может. Каждый объект управления имеет свои критерии качества. Крите-
рий качества должен полностью соответствовать наилучшей работе объекта и в то же время быть достаточно простым для того, чтобы
вариационная задача не оказалась слишком сложной или вообще аналитически неразрешимой.
В качестве критериев оптимальности можно выбрать минимум квадратичной ошибки
26
∞ |
|
J1 = ò x2dt |
(1) |
0 |
|
или минимум абсолютной ошибки
J 2 = |
∞ |
|
||
ò |
x |
dt . |
(2) |
|
|
0 |
|
|
|
Однако оценки (1) и (2) не учитывают того, что в системах имеет место ограничение мощности управляющего сигнала, а также его насыщение при определенных значениях ошибки. Эти ограничения можно задать в интегральном виде:
J3 |
∞ |
U 2dt £ |
ü |
|
||
= ò |
L;ï |
|
||||
|
0 |
|
|
|
ï |
(3) |
|
∞ |
|
|
|
ý |
|
|
|
|
|
ï |
|
|
J4 |
= ò |
|
U |
|
|
|
|
|
|
||||
|
dt £ L. ï |
|
||||
|
0 |
|
|
|
þ |
|
Таким образом, ставится задача обеспечения минимума (1) или
(2) при заданных интегральных ограничениях (3).
Как известно, минимизация критериев J1 и J2 при условии (3) от-
носится к решению задачи на условный экстремум вариационного исчисления. Поэтому функционалы (1) и (2) при наличии ограниче- ний на управление и на координаты приводится к одному функцио- налу. Например,
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
J(U) = ò( |
x |
+ c |
U |
)dt, |
(4) |
|||||||
или |
0 |
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||
J(U) = ò(x 2 + c |
|
U |
|
)dt, |
(5) |
|||||||
|
|
|||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
J(U) = ò(x2 + cU |
2 )dt, |
(6) |
||||||||||
или |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
. |
2 |
|
|
|
|||||
J(U) = ò(x2 |
+ cU2 |
+ c1 U |
|
+ ... + cr (U(r) )2 )dt, |
(7) |
|||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
где сi – множители Лагранжа.
27
Ограничение управляющего сигнала приводит к ограничению всех координат объекта. Поэтому системы, синтезированные по этим оценкам, обеспечивают плавные процессы с небольшим пере- регулированием.
Применение квадратичных оценок дает в результате расчетов линейный закон управления
n |
|
U = å pi xi , |
(8) |
i=1 |
|
поэтому системы, удовлетворяющие минимуму |
|
∞ |
|
J(U) = ò (x 2 + cU 2 )dt, |
(9) |
0 |
|
называются линейными оптимальными системами, а определение закона управления, доставляющего минимум J(U), – аналитическим конструированием регулятора.
Аналитическое решение задачи синтеза оптимальных систем выполняют методами теории оптимального управления. В данной работе произведем аналитическое конструирование регулятора ме- тодом динамического программирования.
Пусть объект управления представлен нормальной системой уравнений:
dx |
1 |
|
= x2 , |
ü |
|
|
|
ï |
|||
dt |
|||||
|
|
ï |
|||
dx 2 |
|
|
ý |
||
|
|
ï |
|||
|
|
|
= a1x1 + a 2 x 2 |
+ u.ï |
|
dt |
|
||||
|
|
þ |
|||
Примем критерий оптимальности в виде:
tk
J = ò(x2 + u2 )dt.
tн
Требуется определить закон управления U = f (x1, x2). Составим функциональные уравнения Беллмана:
(10)
(11)
28
x12 + u2 + x2 |
¶s |
+ (a1x1 + a2x2 + u) |
¶s |
||
|
¶x2 |
||||
|
|
|
¶x1 |
||
2u + |
¶s |
= 0. |
|
|
|
|
|
|
|
||
|
¶x2 |
|
|
|
|
Из второго уравнения системы (12) определим u:
u = - |
1 |
× |
¶ s |
. |
2 |
|
|||
|
|
¶ x 2 |
||
= 0;üï
ï
ý
ï
ï
þ
(12)
(13)
Чтобы определить закон управления u как функцию координат хi, нужно найти функцию S. Подставим в первое уравнение системы (12) выражение для u из (13):
2 |
|
¶s |
|
|
¶s |
|
1 |
æ |
¶s |
ö2 |
|
+ x2 |
+ (a1x1 |
+ a2x2 ) |
= |
ç |
÷ |
|
|||||
x1 |
¶x1 |
¶x2 |
4 |
ç |
¶x2 |
÷ . |
(14) |
||||
|
|
|
|
|
è |
ø |
|
Получим нелинейное дифференциальное уравнение в частных производных.
Для линейных объектов управления и квадратичных критериев оптимальности функция S представляет собой положительно опре-
деленную квадратичную форму вида
S = A x 2 |
+ A |
x x |
2 |
+ A |
22 |
x 2. |
(15) |
|
11 |
1 |
12 |
1 |
|
2 |
|
||
Для определения коэффициентов этой формы составим следую-
щую систему: |
|
|
|
|
|
|
|
|
¶s |
|
|
|
|
ü |
|
|
|
|
= 2A11x1 + A12 x 2 |
; ï |
|
||
|
¶x1 |
|
|||||
|
|
|
|
|
ï |
(16) |
|
|
¶s |
|
|
|
|
ý |
|
|
|
|
|
|
ï |
|
|
|
¶x2 |
|
= 2A22x 2 |
+ A12x1.ï |
|
||
|
|
|
|
|
þ |
|
|
Подставим ¶s / ¶xi в (14) и, |
приравняв к нулю коэффициенты |
||||||
при соответствующих степенях x1 и x2, составим систему алгебраи- ческих уравнений для определения А11, А12 и А22:
4 + 4a |
2 |
= 0; |
ü |
|
|
1A12 - A12 |
ï |
|
|||
A12 + |
|
2 |
= 0; |
ï |
(17) |
2a 2A22 - A 22 |
ý |
||||
2A11 + 2a1A22 + a 2A12 - A12 A22 |
ï |
|
|||
= 0.ï |
|
||||
|
|
|
|
þ |
|
29
Из системы (17) определим: |
|
|
|
|
|
|
|||||
|
|
A12 |
= 2a1 ± 2 |
|
a1 |
2 + 1; |
|
(18) |
|||
|
|
|
|
|
|
|
|
|
|||
|
A 22 |
= a 2 ± a 2 |
2 |
+ A12 ; |
(19) |
||||||
A11 = |
1 |
(−2a1A22 − a 2 A12 + A12A22 ). |
(20) |
||||||||
2 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
||
Из возможных решений А11, А12 и А22 выбираем только положи- тельные значения, так как для получения устойчивого решения функция S должна быть положительно определенной. Зная А11, А12 и А22, определяем закон управления:
|
1 |
|
¶s |
|
1 |
2 |
|
u = - |
× |
= - |
(A12x1 + 2A22x2 ) = åpi xi . (21) |
||||
2 |
¶x2 |
2 |
|||||
|
|
|
i=1 |
В отличие от классического вариационного исчисления и прин- ципа максимума метод динамического программирования позволя- ет получить решение сразу в замкнутой форме, т. е. в виде u = f (xi), что является его достоинством.
Однако применение этого метода также требует большого объе- ма вычислений, связанных с решением системы алгебраических уравнений для определения коэффициентов Аij квадратичной фор- мы S. В аналитическом виде это не всегда возможно сделать и по- этому приходится применять метод подбора значений Аij. Как видно из (21), оптимальная по минимуму функционала (11) система долж- на реализовывать ПД-закон регулирования.
Структурная схема оптимальной системы управления представ- лена на рис. 1а, схема моделирования, составленная с использова- нием пакета MATLAB, – на рис. 1б.
u |
ОУ |
х |
|
|
|
|
|
|
УУ
а)
30
