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

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

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

6. СИНТЕЗ ОПТИМАЛЬНОГО ПИ-РЕГУЛЯТОРА

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

процесса в оптимальной системе и рассчитать прямые показатели качества.

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

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

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

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

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

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

Линейные стационарные системы с оптимальными регуляторами, синтезированными по вектору состояния, позволяют с течением вре-

мени свести к нулю влияние на выход объекта ненулевых начальных условий или кратковременных возмущений. Однако, если входные воздействия остаются постоянными или меняются с малой скоро- стью, то такие системы не обеспечивают равенство нулю статической ошибки регулирования. Для ликвидации этой ошибки регулятор должен включать в себя две составляющие. Одна из них должна за- висеть от фазового вектора, а другая от его интеграла по времени.

Пусть уравнение линейной стационарной системы имеет вид

 

 

&

= AX + BU;X(tH ) = XH .

(1)

 

 

X

Зададим критерий управления в следующем виде

 

J =

[X

T

QX + U

T

& T

&

 

ò

(2)

 

 

RU + U

SU]dt,

tH

где Q и R неотрицательно определенные матрицы; S положи- тельно определенная симметричная матрица.

51

Требуется найти такой вектор управления, который бы миними- зировал данный критерий. При этом будет считаться заданным па- раметр U(tH) = UH.

Введем новые переменные:

&

é

Х ù

, A1

éA Bù

, B1

é0

ù

,

V = U , Z = ê ú

= ê

0

0

ú

= ê ú

 

ëU û

 

 

 

ë

û

 

ëI

û

(3)

R 1 = S, Q1

 

éQ

0

ù

 

 

 

 

 

 

 

=

ê

 

R

ú.

 

 

 

 

 

 

 

 

 

ë0

 

û

 

 

 

 

 

 

 

 

Запишем уравнение движения и критерий управления в новых координатах:

&

(4)

Z = A1 Z + B1V ,

 

J = ò[ZT Q1 Z + V T R 1V ]dt .

(5)

t H

Используя процедуру метода динамического программирования, можно определить оптимальное значение вектора V:

 

V = -R1−1 × B1T × K × Z,

 

 

(6)

где К решение матричного уравнения Риккати:

 

 

&

T

−1

T

K − Q1 .

(7)

K = −KA 1 − A1

K + KB 1R 1

B1

Тогда оптимальное управление при заданных начальных услови-

ях определится следующим образом:

 

 

 

 

 

 

 

 

 

 

& 0

 

 

 

 

−1

é0ùT

éK11 K 21

ù

é X ù

 

U

= -S

 

 

êI

ú ×

êK

21

 

K

22

ú ×

êU 0

ú =

 

 

 

 

 

 

 

 

ë

û

ë

 

 

 

 

û

ë

û

 

S

−1K

21

X S −1 K

22

U 0

= C

 

X + C U .

(8)

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

Из исходного уравнения движения получим:

 

 

 

U = B

−1

 

&

− AX) = (B

T

B)

−1

B

−1 &

− AX).

(9)

 

 

 

(X

 

 

 

(X

Продифференцируем это уравнение по времени и подставим со- отношение для расчета оптимального управления:

52

&

0

= C1X + C 2 (B

T

B)

−1

B

−1 &

 

 

 

 

=

 

 

 

U

 

 

 

(X − AX )

 

 

 

C

2 (B

T

B)

−1

B

−1

&

 

 

 

 

 

T

B)

−1

B

−1

A)X

=

(10)

 

 

 

X + (C1 − C 2 (B

 

 

 

=

 

&

+ C 4 X,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C3 X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

C 3 = C 2 (B T B) −1 B −1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C 4 = C1 − C 2 (B T B) −1 B −1A.

 

 

 

Проинтегрируем выражение (10) по времени и получим:

 

или

 

U0 (t) = U(tH ) + C3 (X(t) − X(t H )) + C4 ò Xdt

 

 

 

U 0 (t) = C 5 + C 3 X (t ) + C 4 ò X (t )dt .

 

 

 

 

 

(11)

Таким образом, оптимальное управление будет реализовано в виде ПИ-регулятора.

Структурная схема системы с таким регулятором будет иметь вид:

 

 

u

Х ′ = АХ + ВU

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C5

 

 

 

C3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C4

 

ò X ( t)dt .

 

 

 

 

 

 

 

 

Пример. Определить оптимальное управление объектом, задан- ным уравнением:

T y& (t) + y(t) = ku (t)

или

x = ax + bu ,

(12)

&

 

где Т и k постоянная времени и коэффициент усиления объекта; a = − T1 ; b = Tk ; х = y(t); u = u(t)

впроцессе перехода объекта из фиксированного начального х(0) = 1

вфиксированное конечное х(∞) = 0 состояние из условия минимума функционала:

53

 

 

 

 

2

 

2

 

& 2

 

 

 

 

J = ò

(x

+ u

 

)dt.

(13)

 

 

 

 

+ u

Решение

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем обозначения:

 

 

 

 

 

 

 

 

 

éxù

&

=

éa bù

 

 

é0ù

 

é1 0ù

é1 0ù

Z= ê ú, V = u, A1

ê

ú, B1 =

ê ú,Q1

= ê ú, R1

=S= ê ú. (14)

ëuû

 

 

ë0 0û

 

 

ë1û

 

ë0 1û

ë0 1û

На основании (6) записываем линейный закон оптимального управления:

V = -S

−1

ék11

k12

ù

éx

 

ù

= -k12 x - k 22 u

0

&

0

. (15)

 

×[0 1]× ê

 

ú

× ê

0

ú

 

= u

 

 

 

ëk 21

k 22 û

êu

ú

 

 

 

 

 

 

 

 

 

 

ë

 

û

 

 

 

 

 

Коэффициенты k12 и k22 находим из матричного уравнения Рик-

кати (7), учитывая, что K& (t) =0:

é0 0ù

ék11

k12

ù

éa bù

-

éa 0

ù

ék11

k12

ù

ék11

k12

ù é0ù

´

ê

ú

= -ê

 

 

 

ú

× ê

ú

ê

ú

× ê

 

ú

+ ê

 

ú × ê

ú

ë0 0û

ëk21

k 22 û

ë0 0

û

 

ëb 0û

ëk 21

k 22 û

ëk 21

k 22 û ë1

û

(16)

 

 

ék

 

k

 

ù

é1

0ù

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

´ S−1[0 1]× ê

11

 

12 ú

- ê

 

ú.

 

 

 

 

 

 

 

 

 

 

 

 

ëk 21

k 22 û

ë0 1û

 

 

 

 

 

 

 

 

 

 

Перемножив матрицы из уравнения (16), запишем три алгебраи- ческих уравнения:

ì

- 2ak 11 + k122

- 1 = 0;

 

 

ï

- bk 11 - ak 12

+ k 21 k 22 = 0;

 

í

(17)

ï

- 2bk 12 + k 222 - 1 = 0.

 

î

 

Решая эти уравнения, найдем k11, k12 и k22.

 

Из уравнения (12) выразим u:

 

 

 

 

 

 

 

 

 

 

u =

1

&

 

 

a

x

 

 

 

(18)

 

 

 

 

 

 

 

 

 

 

 

b

x -

 

 

b

 

 

 

и подставим (18) в (15):

 

 

 

 

 

 

 

 

 

 

 

 

1

 

a

 

 

&

0

= -k12 x

- k22 (

&

x).

(19)

 

 

 

 

u

 

b

x -

b

 

 

 

 

 

 

 

 

 

 

Проинтегрировав выражение (19) по времени с учетом началь- ных условий, получим:

54

 

 

u0

 

 

k

22

 

 

k

22

a − k

12

b

k

22

 

 

 

(t) = −

 

 

x(t) +

 

 

 

 

 

ò x(t)dt +

 

 

 

 

b

 

 

 

b

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

tH

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u 0 (t) = K 1 x (t) + K 2

 

 

 

 

 

 

 

 

 

 

 

 

ò x (t)dt + K 3 ,

 

(20)

 

 

 

 

 

 

 

 

 

 

 

t H

 

 

 

 

 

 

 

 

 

где K1

= −

k 22

, K2

=

k 22a − k12b

,

K3

=

k 22

.

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

b

 

 

 

 

 

 

 

 

b

 

 

 

Таким образом, получен оптимальный ПИ-закон управления.

Перечень используемого оборудования

Перечень используемого оборудования приведен в лаб.работе №1.

Указания по технике безопасности

Указания по технике безопасности приведены в лаб.работе №1.

Указания по порядку выполнения работы

Работа выполняется на ПЭВМ каждым студентом самостоятельно.

Задание в виде таблицы значений ординат переходной функции объекта (кривой разгона) каждый студент получает от преподавателя.

1.По заданным в таблице 1 лабораторной работы № 1 ордина- там кривой разгона объекта управления построить кривую разгона.

2.Провести касательную к точке перегиба и определить пара- метры передаточной функции объекта.

3.Определить закон оптимального управления объектом, задан-

ным уравнением Ty& (t) + y(t) = ku (t) , из условия минимума функ-

 

2

 

2

&

2

 

ционала J = ò

(x

+ u

)dt при заданных начальных условиях

 

 

+ u

 

0

х(tH) = 1 и u(tH) = 0.

4. Построить график кривой переходного процесса в оптималь- ной системе и определить прямые показатели качества.

Содержание отчета

Отчет о выполнении лабораторной работы должен содержать: 1. Исходные данные в виде таблицы ординат кривой разгона.

55

2.График кривой разгона, на котором приводится аппроксима- ция и определение параметров объекта k и Т.

3.Решение задачи аналитического конструирования регулятора.

4.Структурную схему оптимальной системы управления и схе- му, составленную с использованием пакета MATLAB.

5.График кривой переходного процесса.

6.Значения прямых показателей качества процесса.

7.Анализ полученных результатов.

Контрольные вопросы

1.Как решается задача синтеза оптимального ПИ-регулятора?

2.Какой вид имеет матричное уравнение Риккати?

3.Какой критерий оптимальности использован при выполнении лабораторной работы?

4.Какой закон оптимального управления объектом был получен при выполнении лабораторной работы?

Литература

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

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

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

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

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

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

7.Пантелеев А. В., Бортаковский А. С., Летова Т. А. Оптималь-

ное управление в примерах и задачах: учеб. пособие. М.: Изд-во МАИ, 1996. 212 с.

56

7.СИНТЕЗ ОПТИМАЛЬНОГО РЕГУЛЯТОРА ДЛЯ ОБЪЕКТА С ЗАПАЗДЫВАНИЕМ

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

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

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

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

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

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

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

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

Наиболее распространенной является задача определения опти- мальных уравнений U0(t), обеспечивающих оптимальные процессы в функции времени при заданных начальных Y(t0) и конечных Y(tк) зна- чениях выходных переменных. Оптимальные управления U0(t) объ- ектом и их реализацию различают для двух основных случаев: без применения обратных связей (разомкнутая система) и с применением обратных связей (замкнутая система). В разомкнутой оптимальной системе управляющее устройство УУ по заданным граничным усло- виям формирует оптимальный закон управления (рис. 1, а). При этом

координаты объекта О не используются в формировании закона управления. В замкнутой оптимальной системе на вход УУ подаются соответствующие сигналы с выхода объекта (рис. 1б).

57

а)

Y(t0)

 

 

U(t)

 

 

Y(t)

 

УУ

 

О

 

 

 

 

 

Y(tk)

 

 

 

 

 

 

 

 

 

 

 

б)

 

U(Y,Xвх)

 

Y(t)

 

 

 

 

Xвх(t)

УУ

 

 

О

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 1. Структурная схема оптимальной системы управления

Время τ между началом скачкообразного изменения входной ве-

личины и началом изменения выходной называют чистым или транспортным запаздыванием.

Рассмотрим АСР, структурная схема которой представлена на рис. 2.

U

 

X

объект

 

 

 

 

 

ε g(t)

регулятор

Рис. 2. Структурная схема АСР

В момент времени t = 0 на регулятор поступает ступенчатое воз-

действие

ì 0, при t < 0, g(t) = í

îG,при t ³ 0.

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

a × х(t) + b× x(t) + x(t) = K × U ×(t - t),

(1)

&&

&

 

где α и β – постоянные коэффициенты; К коэффициент усиления объекта; τ – чистое запаздывание.

За критерий оценки качества переходного процесса примем функционал

58

¥

2

 

2

&

2

 

 

J = ò[e

(t) + A

(t )]dt ,

(2)

 

 

× e

 

t

где e(t) – отклонение регулируемой величины от заданного значе- ния; А коэффициент веса, величина которого выбирается в соот- ветствии с требованиями к переходному процессу.

e(t) = g(t) - x(t).

(3)

Нижний предел интегрирования в (2) выбран равным t, так как на

интервале времени от 0 до t мы не можем воздействовать на регули- руемую величину, а следовательно и на ошибку регулирования (3).

ε(t) ≡ Gппр0 < t ≤ τ.

(4)

Граничные условия задачи

 

ε(τ) = G;ε(∞) = 0.

(5)

Необходимо с учетом (4) и (5) получить такой закон управления U = L(ε), чтобы соответствующая ему функция e(t) доставляла ми- нимум функционалу (2).

L – искомый оператор регулятора.

Для нахождения оптимальной e(t) решим вариационную задачу. Составим в соответствии с (2) функцию f:

f= e2 + A2 × e&2

изапишем уравнение Эйлера

fe - dtd × fe& = 0,

которое после взятия производных сводится к виду

А2 ×&e& - e = 0.

Решение уравнения Эйлера с учетом (4) и (5) будет иметь вид

ì

 

-

t-t

 

ï

 

A

,при t ³ t,

e(t) = íG ×e

 

 

ï

G,при 0 < t £ t.

î

Применив к (7) преобразование Лапласа, получим

e(p) =

G

-

G

 

×e-t×p.

p

p×(A×p +1)

Согласно (3) и (8):

 

 

 

 

 

G

 

 

Х(p) =

 

 

× e-t×p .

p × (A × p + 1)

 

 

 

 

(6)

(7)

(8)

(9)

59

Имея в виду, что

 

 

 

 

 

&

= 0,

 

(10)

Х (+0) = Х (+0)

 

и учитывая (9), из (1) получим

 

 

 

U(p) =

G × (a × p2

+ b × p + 1)

.

(11)

K × p × (A × p + 1)

 

 

 

С учетом (8) и (11) определим передаточную функцию опти-

мального регулятора

W0

(p) = U(p)

=

a × p2 + b ×p +1

.

(12)

 

p

e(p)

 

K ×(A × p +1- e-t×p )

 

 

 

Если передаточную функцию регулятора

Wp (p) = p2 ×+b××p +1 K A p

охватить обратной связью с передаточной функцией

W

(p) = W (p) ×[1- e-t×p ]=

 

К

×[1- e-t×p ],

(13)

 

 

 

о.с.

1

p2

+b× p +1

 

 

 

 

 

то можно реализовать структуру оптимального регулятора (12). Схема, реализующая оптимальный регулятор, имеет вид

 

 

U

 

 

 

 

К × е-t×р

 

 

 

 

 

 

 

 

X

 

 

 

 

 

a × р2 + b × р +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е-t×р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a × р2 + b × р +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wp(p)

 

 

 

 

 

 

 

 

 

U

 

 

 

 

 

 

 

 

ε

 

 

 

G(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a × р2 + b × р +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К × А × р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3. Схема, реализующая оптимальный регулятор

60

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