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

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

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

25

6

0

0

1

2

-7

26

0

0

13

0

5

6

27

13

5

0

0

5

-6

28

0

0

17

0

2

2

29

17

7

0

0

2

-2

30

0

0

8

0

2

10

31

8

0

0

0

2

-10

32

0

0

9

3

4

11

33

9

3

0

0

4

-11

34

4

7

0

0

1

-9

35

0

0

15

5

3

8

36

15

5

0

0

3

-8

37

0

0

18

0

10

5

38

0

0

10

1

2

10

39

10

1

0

0

2

-10

40

0

0

19

0

12

20

41

19

0

0

0

12

-20

42

0

0

16

1

4

10

43

16

1

0

0

4

-10

44

0

0

11

1

2

11

45

11

1

0

0

2

-11

Номер варианта выбирается по величине суммы последних трех цифр шифра зачетной книжки студента (например, три последние цифры 525, следовательно, их сумма равна 12, значит, в этом случае необходимо взять данные для варианта №12).

Литература

1.Ким Д. П. Сборник задач по теории автоматического управ- ления: Многомерные, нелинейные, оптимальные и адаптивные си- стемы: учебное пособие. М.: Физматлит, 2008. 328 с.

2.Аттетков А. В., Галкин С. В., Зарубин В. С. Методы оптими- зации: учебник для вузов. М.: Изд-во МГТУ, 2001. 440 с.

91

3.Афанасьев В. Н., Колмановский В. Б., Носов В. Р. Математи- ческая теория конструирования систем управления. М.: Высшая школа, 1998. 574 с.

4.Ванько В. И., Ермошина О. В., Кувыркин Г. Н. Вариационное исчисление и оптимальное управление: учебник для вузов / под ред. В. С. Зарубина, А. П. Крищенко. 2-е изд. М.: Изд-во МГТУ, 2001. 488 с.

5.Иванов В. А., Фалдин Н. В. Теория оптимальных систем ав- томатического управления: учеб. пособие для вузов. М.: Наука, 1981. 331 с.

6.Куропаткин П. В. Оптимальные и адаптивные системы: учеб. пособие для вузов. М.: Высш. школа, 1980. 287 с.

7.Пантелеев А. В., Бортаковский А. С., Летова Т. А. Оптималь- ное управление в примерах и задачах: учеб. пособие. М.: Изд-во МАИ, 1996. 212 с.

92

4. РЕШЕНИЕ ЛИНЕЙНОЙ ЗАДАЧИ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ МЕТОДОМ ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ

Цель изучение метода динамического программирования

Знания и умения, приобретаемые студентом

врезультате освоения темы

Врезультате освоения темы студент должен:

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

вания, получать решения в терминах выбранного метода.

Формируемые компетенции (или их части)

-ОК-1способность совершенствовать и развивать свой интел- лектуальный уровень;

-ОК-8способность самостоятельно приобретать с помощью информационных технологий и использовать в практической дея- тельности новые знания и умения;

-ПК-21способность выбирать оптимальные решения.

Актуальность темы

Метод динамического программирования был разработан почти одновременно с опубликованием принципа максимума. В основе метода динамического программирования лежит принцип опти- мальности, сформулированный Беллманом. Согласно этому прин- ципу, оптимальное управление определяется конечной целью управления и состоянием системы в рассматриваемый момент вре- мени, независимо от того, каким образом система пришла в это со- стояние. Таким образом, будущее оптимальное управление не зави- сит от его предыстории. Это значит, что любой конечный участок (от любой промежуточной точки до конца) оптимальной траектории является тоже оптимальной траекторией.

Теоретическая часть

Задачей оптимизации считается определение оптимальных управлений U0(t) и траекторий X0(t) из условия минимума (макси- мума) функционала

93

T

 

J = òf0 (t, X, u)dt; T = tк

(4.1)

tн

 

для заданных уравнений состояния объекта

 

&

(4.2)

X = AX + Bu = f (t, X,u) ,

а также начальных и конечных фиксированных значений X(tн) и X(T) и интервала tн t T при наличии ограничений вида X(t) Gx

и U(t) Gu . Здесь Gx и Gu заданные допустимые области для ко-

ординат состояния и управлений.

При этом вводится вспомогательная функция Беллмана

T

 

 

(t, X, u)dt] .

(4.3)

S(t , X) = min[ò f0

U t′

 

 

Минимум функционала (4.1) при условиях (4.2) зависит от мо-

мента времени t/+

t, значения X t

и конечного момента времени T:

S(t′ +

t, T, X t ) = min[

T

(4.4)

ò F(t, X, u)dt] ,

 

U

t′+ t

 

где t/ + t произвольный промежуточный момент времени в ин-

тервале t/< t/ +

t < T;

X t = X(t′) +

X вектор координат состояния в момент времени

t/ + t.

 

Пусть t/ фиксированный момент времени, t малое положи- тельное число (t/ < t /+ t < T). Тогда с учетом (4.4) вспомогательная

функция

 

 

t′+ t

T

 

S(t, X) = min [

ò f0 (t, X, u)dt +

ò f0 (t, X, u)dt ] =

 

 

u

t′

t′+ t

(4.5)

 

t′+ t

 

 

= min [

ò f0 (t, X, u)dt + S(t′ +

t, T, X t )].

 

u

t′

 

 

 

Первое слагаемое в правой части выражения (4.5) с точностью

до малых величин более высокого порядка, чем

t, можно заменить

приближенным значением

 

t′+ t

 

ò f0 (t, X, u)dt = f0 (t, X, u) t .

(4.6)

t′

 

Второе слагаемое в правой части выражения (4.5), определяю- щее значение функции S(t′ + t, X + X) для любого момента време-

94

ни t/ + t, разложим в ряд Тейлора. Ограничиваясь линейными чле- нами относительно приращений xi и t и переходя к пределу, за- меним это слагаемое приближенным значением, если функция S(t, X) имеет непрерывные частные производные по всем t и xi(t):

 

 

 

 

 

 

 

n

S

 

S

 

 

 

 

 

 

S(t + Dt, X + DX) = S(t, X) + å

 

Dx i +

 

Dt =

 

 

 

x i

t

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

(4.7)

é n

S

dx

 

 

S ù

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

= S(t, X) + êå

 

 

 

+

 

úDt.

 

 

 

 

 

 

 

 

x i

 

 

 

 

 

 

 

 

 

 

 

ëi=1

dt

 

 

t û

 

 

 

 

 

 

 

 

 

На основании (4.5), (4.6) и (4.7) запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

é n

S

 

dxi

 

S

ù

S(t, X) = min{f0 (t,X, u)Dt + S(t,X) + êå

 

+

úDt},

xi dt

t

u

 

 

 

 

 

 

 

 

ëi=1

 

û

откуда после деления на

t всех членов и перехода к пределу при t

0 получим нелинейное дифференциальное уравнение Беллмана в частных производных относительно неизвестной функции S(t, X):

 

n

S

 

S

 

 

min {f0 (t, X, u) + å

 

fi (t, X, u) +

 

} = 0

(4.8)

x i

t

u

i=1

 

 

 

при условиях

 

 

 

 

 

 

xi (t) = fi (t, X, u) ;

 

 

 

&

 

 

 

 

 

(4.9)

x(t) ÎG x ,

U(t) ÎG u , t н £ t £ T .

 

 

 

Оптимальное управление U0(t) найдем в результате решения уравнения (4.8). Если координата управления ограничена U(t)ÎG U ,

то для внутренней точки области, определяемой множеством GU ,

условие (4.8) можно заменить функциональными уравнениями в

частных производных

 

 

 

n

S

 

 

S

 

ü

 

f0 (t, X,u) + å

 

 

fi (t,X, u) +

 

 

= 0;

ï

 

xi

t

 

 

 

 

i=1

 

 

 

ï

(4.10)

æ

 

n

S

 

ö

 

ý

 

 

 

 

 

 

ç

(t, X,u)

+ å

 

÷

 

ï

 

 

 

 

 

 

 

 

çf0

xi

fi (t,X,u)÷

= 0.ï

 

 

u è

 

i=1

 

ø

 

þ

 

Если S не зависит явно от времени, то St = 0 . В этом случае при

решении уравнений (4.10) для квадратичных функций f(…) функци- онала (6.1) и tн=0; tк=T=∞ А. М. Летовым предложено искать вспо- могательную функцию в виде квадратичной формы дифференциру- емой по всем координатам xi:

95

n

S(X) = åAij xi x j . (4.11)

i, j=1

Для линейных объектов функция S(X) является функцией Ляпу- нова при t → ∞. В результате решения (4.10) с учетом (4.11) найдем

оптимальное управление в функции координат вектора состояния U0(X), которое обеспечивает устойчивые процессы.

Таким образом, синтез оптимального управления методом дина- мического программирования ведется в следующем порядке:

1)из условия равенства частных производных по компонентам вектора управления определяется допустимый вектор U*(S);

2)используя зависимость U*(S), решается уравнение Беллмана при заданных краевых условиях и определяется функция S(X). Ре- шение для S обычно ищется в виде квадратичной формы:

 

 

 

 

S = XT × A × X ,

где А симметричная матрица.

S = (x1x 2 )

 

a11

a12

 

×

 

x1

 

= a11x12 + 2a12 x1x 2 + a 22 x 22 ;

 

 

 

 

 

 

a12

a 22

 

 

 

x 2

 

 

3) используя функцию S(X), получаем зависимость U(X).

Вопросы и задания для самостоятельного решения Пример. Определить оптимальное управление объектом, заданным

уравнением

&

= a 0 x1

+ b0 u1 ,

Ty& (t) + y(t) = kU (t) или x1

где T и k постоянная времени и коэффициент усиления объекта

a0

= -

1

; b0

=

k

из

условия

минимума

функционала

T

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J = ò[q1x12 + r1u12 ]dt , где q1 > 0, r1 > 0

весовые коэффициенты.

 

0

 

 

 

 

 

 

 

 

Вспомогательную функцию зададим в виде

S(X) = A11x12 .

Для данного случая на основании (4.10) запишем функциональ-

ные уравнения Беллмана

2

 

 

2

+ (a 0 x1 + b0 u1 )S / x1

ü

q1x1

+ r1u1

= 0;ï

2r u

 

+ b

 

S / x

 

= 0,

ý

1

0

1

ï

1

 

 

 

 

þ

96

Откуда u

0

(t ) = −[b

0

/( 2r

 

)]∂S / ∂x

1

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или с учетом S(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u 0 (t) = −[b

0

/( 2r

)]A

11

x

1

(t) = k

o.c

y(t) ,

 

 

 

где ko.c = −A11b0 / r1 .

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

найдем производную ∂S / ∂x1

 

 

Для определения ko.c

и подставим

в первое уравнение системы Беллмана:

 

 

 

 

 

 

 

 

 

q

1

x 2 + (b

2

/ r )A

2

 

x 2 + 2a

0

A

11

x 2

− 2(b

2

/ r )A

2

x 2

= 0 ,

 

 

 

1

 

 

0

 

1

11

 

1

 

 

 

 

 

 

1

 

0

1

11

1

 

или b 2 A 2

 

− 2a

0

r A

11

 

− q

r

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

11

 

 

1

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда A

11

= [r (a

0

±

 

 

a 2

+ b 2 q

1

/ r

 

)] / b 2 .

 

 

 

 

 

 

 

1

 

 

 

0

 

 

0

 

 

1

 

 

 

 

0

 

 

 

 

 

 

Подставив значение А11 в выражение для ko.c , получим выражение k o.c = (−a 0 a 02 + b02 q1 / r1 ) / b0 .

Литература

1.Ким Д. П. Сборник задач по теории автоматического управ- ления: Многомерные, нелинейные, оптимальные и адаптивные си- стемы: учебное пособие. М.: Физматлит, 2008. 328 с.

2.Аттетков А. В., Галкин С. В., Зарубин В. С. Методы оптими- зации: учебник для вузов. М.: Изд-во МГТУ, 2001. 440 с.

3.Афанасьев В. Н., Колмановский В. Б., Носов В. Р. Математи- ческая теория конструирования систем управления. М.: Высшая школа, 1998. 574 с.

4.Ванько В. И., Ермошина О. В., Кувыркин Г. Н. Вариационное исчисление и оптимальное управление: учебник для вузов / под ред. В. С. Зарубина, А. П. Крищенко. 2-е изд. М.: Изд-во МГТУ, 2001. 488 с.

5.Иванов В. А., Фалдин Н. В. Теория оптимальных систем ав- томатического управления: учеб. пособие для вузов. М.: Наука, 1981. 331 с.

6.Куропаткин П. В. Оптимальные и адаптивные системы: учеб. пособие для вузов. М.: Высш. школа, 1980. 287 с.

7.Пантелеев А. В., Бортаковский А. С., Летова Т. А. Оптималь- ное управление в примерах и задачах: учеб. пособие. М.: Изд-во МАИ, 1996. 212 с.

97

5. РЕШЕНИЕ ЗАДАЧИ АНАЛИТИЧЕСКОГО КОНСТРУИРОВАНИЯ РЕГУЛЯТОРОВ МЕТОДОМ ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ

Цель изучение методики решения задач аналитического констру- ирования регуляторов методом динамического программирования.

Знания и умения, приобретаемые студентом

врезультате освоения темы

Врезультате освоения темы студент должен: знать метод динамического программирования;

уметь выбирать необходимый математический метод исследо-

вания, получать решения в терминах выбранного метода.

Формируемые компетенции (или их части)

-ОК-1способность совершенствовать и развивать свой интел- лектуальный уровень;

-ОК-8способность самостоятельно приобретать с помощью информационных технологий и использовать в практической дея- тельности новые знания и умения;

-ПК-21способность выбирать оптимальные решения.

Актуальность темы

Существует традиционная схема синтеза систем регулирования, согласно которой, имея уравнение объекта управления (ОУ) и пока- затели качества регулирования (время регулирования, перерегули- рование), выбирается регулятор. Если при выбранном регуляторе не удовлетворяется заданное качество, то синтезируются корректиру- ющие цепи. В выборе корректирующих цепей собственно и состоял синтез системы регулирования. Данная схема синтеза получила большое распространение, а методы синтеза непрерывно совершен- ствуются. На основе результатов этого метода сконструированы и

выпускаются промышленные регуляторы с различными законами (П, ПИ, ПИД и т. д.). Методика синтеза распространяется на линей- ные, нелинейные системы, системы с запаздыванием, системы, под- вергавшиеся действию случайных возмущающих воздействий.

98

Однако появились новые идеи, которые позволяют по-другому подойти к задаче синтеза систем. Эти идеи оформились в теории

аналитического конструирования оптимальных регуляторов (АКОР). Сущность ее состоит в том, что закон изменения управля- емой координаты и управляющего воздействия задается аналитиче- ски в виде некоторого функционала J(x, u) , который нужно мини-

мизировать подбором управляющего воздействия. Зная алгоритм управления U(t) или U(x) , конструируют регулятор.

Теоретическая часть

Пусть объект управления представлен нормальной системой уравнений

dx

1

 

= x 2 ;

ü

 

 

ï

dt

 

 

ï

dx 2

 

 

ý

 

 

ï

 

 

 

= a1x1 + a 2 x2

+ u.ï

dt

 

 

 

þ

Примем критерий оптимальности в виде

J(u) = ò (q1x12 + q2 x22 + u 2 )dt .

0

Требуется определить закон изменения U = f(x1, x2).

Составим функциональные уравнения Беллмана

q1x12

 

2 x 22 + u2

 

S

 

 

S

ü

+ q

+ x2

 

+ (a1x1 + a 2 x2

+ u)

 

= 0;ï

x1

x2

 

 

 

 

 

 

 

ï

 

S

 

 

 

 

 

 

 

ý

2u +

 

= 0.

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

þ

Определим u из второго уравнения системы (5.3)

u = -

1

 

S

.

2

 

 

 

x2

(5.1)

(5.2)

(5.3)

Чтобы определить закон изменения u как функции координат xi, нужно найти функцию S. Подставим в первое уравнение системы (5.3) выражение для u из (5.4):

q1x12 + q2 x22 + x2

S

+ (a1 x1 + a2 x2 )

S

=

1

æ

S

ö2

 

ç

÷ .

(5.5)

x1

x2

4

x2

 

 

 

è

ø

 

99

Получим нелинейное дифференциальное уравнение в частных производных. Для линейных объектов управления и квадратичных критериев оптимальности функция S представляет собой положи-

тельно определенную квадратичную форму вида

S = A11 x12 + 2A12 x1x 2 + A 22 x 22 .

Для определения коэффициентов этой формы составим следую- щую систему:

S

 

= 2A11x1

+ 2A12 x 2

 

ü

 

 

 

; ï

 

x1

 

 

 

 

 

ï

(5.6)

S

 

 

 

 

ý

 

= 2A12 x1

+ 2A 22 x 2

ï

 

x 2

 

.ï

 

 

 

 

 

þ

 

Подставим S в (5.5) и, приравняв нулю коэффициенты при со-

xi

ответствующих степенях x1, x2, составим алгебраические уравнения

для А11, А12 и А22:

2

ü

 

q1 + 2a1A12 - A12 = 0;

ï

 

q2 + 2A12 + 2a 2A22 - A222 = 0;

ï

(5.7)

ý

2A11 + 2a1A22 + 2a 2A12 - 2A12 A22

ï

 

= 0.ï

 

þ

Из системы (5.7) определим

A12 = a1 ± a12 + q1 ;

A22 = a 2 ± a 22 + q2 + 2A12 ;

A11 = 12 (-2a1A22 - 2a2A12 + 2A12A22 ) .

Из возможных решений А11, А12 и А22 выбираем только положи- тельные значения, так как для получения устойчивого решения функция S должна быть положительно определенной. Зная А11, А12 и А22, определяем закон управления:

 

1

 

S

 

1

 

2

 

u = -

 

= -

(2A12 x1

+ 2A22 x2 ) = åpi xi .

(5.8)

2

 

x2

2

 

 

 

 

i=1

 

В отличие от классического вариационного исчисления и прин- ципа максимума метод динамического программирования позволя- ет получить решение сразу в замкнутой форме, то есть в виде U = f(xi), что является его достоинством.

100

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