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

ОиАС-7

.pdf
Скачиваний:
11
Добавлен:
26.03.2016
Размер:
570.49 Кб
Скачать

Решая полученную систему из шести дифференциальных уравнений с помощью метода Рунге – Кутта, получим:

p0

53,4 109 ; p0

147 106 ;

p0

12,6 108 ;

 

11

 

 

12

 

13

 

(4.1.34)

p0

92 104 ;

p0

44 105 ; p0

 

116 106

 

22

 

23

 

33

 

 

 

Искомые параметры регулятора вычисляются на основе чисел (4.1.34):

c = – 0,126 107; с

2

= 0,44 104; c

3

= –116 103.

(4.1.35)

1

 

 

 

11

Аналитическое конструирование по

критерию обобщенной работы

12

В 1967 г. Александр Аркадьевич Красовский предложил упрощение процедуры АКОР с вычислительной стороны. Для этого в функционал (4.1.3)

J x Qx u u dt

0

вводится дополнительное слагаемое, с учетом которого функционал оптимизации принимает (в развернутой форме) вид

 

 

 

n

 

 

 

xi xj

J qij

0

ij 1

 

 

 

 

m

uk2

k 1

 

1

 

m n

v

 

 

4

xi

 

k 1 i 1

2

bik dt (4.1.36)

где квадратичная форма v=x'Px содержит положительно-определенную матрицу Р, являющуюся решением матричного алгебраического уравнения

PA+A'P+Q = 0.

(4.1.37)

Оптимальное управление определяется по-прежнему на основе формулы (4.1.13). Для того чтобы убедиться в этом, положим вначале n=m=1. Функционал (4.1.36)

примет вид

2

 

2

 

1

 

v

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

qx

u

 

 

 

 

 

x

b

 

dt

(4.1.38)

 

 

 

 

0

 

 

 

 

4

 

 

 

 

 

Подставляя в уравнение (4.1.6)

 

 

 

 

 

 

 

v

 

v

ax

1

 

v

2

qx

2

 

 

 

 

 

 

t

x

 

 

 

x

b

 

 

 

 

 

 

 

4

 

вместо qx2 выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

2

 

 

 

 

 

 

 

 

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

qx

 

 

 

 

x

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

получим вместо нелинейного алгебраического уравнения

 

2pa p2b2 + q = 0

(4.1.9)

 

линейное уравнение

 

 

 

2pa + q = 0

(4.1.39)

для определения коэффициента р квадратичной формы v = px2.

13

Таким образом, аналитическое конструирование по критерию обобщенной работы состоит в решении линейного алгебраического уравнения (4.1.37) и вычисления искомой матрицы С по формуле (4.1.13).

Уравнение (4.1.37)

PA+A'P+Q = 0.

называется уравнением Ляпунова. Оно имеет единственное решение Р>0, в частности, тогда, когда собственные числа матрицы А имеют отрицательные вещественные части. При этом условии нетрудно показать, что синтезированная система асимптотически устойчива. Действительно, в соответствии с прямым методом Ляпунова примем в качестве функции Ляпунова функцию v=x'Px > 0, вычисляя ее полную производную по времени, получим, что

dv/dt= –х' [Q + 2CC'] х < 0.

Функционал (4.1.36) называется критерием обобщенной работы. Это название связано с тем, что последнее слагаемое в (4.1.36) можно записать как

uопт uопт dt

0

который выражает собой «энергию» (обобщенную работу) оптимального

управления (uопт).

14

Пример 4.1.3. Рассмотрим следующую модель объекта управления:

x1 x2 ;

x2 x1 x2 u

Необходимо определить оптимальный закон управления по критерию

обобщенной работы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

v

 

 

2

 

 

 

 

 

 

 

 

 

 

J x12

x22 u2

 

 

 

 

 

bi

 

 

dt

4

 

x

 

0

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

Решение. В соответствии с процедурой синтеза v=x'Px, где матрица P определяется из решения матричного уравнения

PA+A'P+Q = 0.

p

p

0

1

 

0

1 p

p

 

1

0

 

11

12

 

 

 

 

11

12

 

 

 

0

p12

p22 1

1

 

1

1 p12

p22

0

1

 

 

 

p11

1,5;

p12

0,5, p22

1

 

 

 

15

Тогда функция v запишется как

v 1,5x12 x1 x2 x22

функционал обобщенной работы принимает вид:

 

2

2

 

2

 

1

2

 

J x1

x2

u

 

 

 

x1 2x2

dt

 

4

0

 

 

 

 

 

 

 

А закон управления, обеспечивающий минимум этого функционала записывается следующим образом:

u CT x PB T x 0,5x1 x2

16

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