Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Аналитическое конструирование оптимальных регуляторов. Практикум для лабораторных и практических занятий

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
974 Кб
Скачать

Производную по 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 + 2r × x12 + 2r × x1 × x 2 + r(q + 2r ) × 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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]