Файл: Ларин В.Б. Синтез оптимальных линейных систем с обратной связью.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 13.07.2024

Просмотров: 132

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

и задача сводится к определению W таким образом, чтобы первая вариация функционала (5.19) обращалась в нуль.

Первый этап решения сформулированной игровой зада­ чи, включая вычисление первой вариации функционала (5.19), целиком совпадает с соответствующим этапом реше­

ния

задачи

синтеза оптимальной

системы стабилизации

(§ 3,

гл.

1) и, согласно (1.76),

 

 

& =

 

/°°

 

 

- f

J

Sp {6Ф*

+ С271 (Л^*/г —

'- / с о

-GJjQ-'A) Р - 1] + 5 Ф [CD.Qr'G.GQ-1 +

+ р 7' ( W

-

^ Q r ’G.G) Q -1] 6Ф} ds,

где

 

 

 

G*G =

N^RN +

Ap (s) ЛДр (s),

 

Ap (s) =

det P,

N =

Ap (s) P _1M,

Q =

Ap (s) ß + ЛУѴ,

а полиномиальные матрицы Л и В должны удовлетворять требованию аналитичности в правой полуплоскости матрицы

'Р — М -1

Z =

А В

Следующий шаг при решении сформулированной задачи уже несколько отличается от решения, приведенного в гл. 1. Ввиду наличия отрицательных множителей Лагранжа мат­

рица Q7l G^GQ~l уже не будет положительно определенной и ее нельзя представить в виде произведения двух сопря­

женных матриц. Решение задачи

можно получить, предста­

вив матрицу QP’G^GQ-1 в виде

 

 

Q7lG*G<3-1 =

H J H , 1

(5-20)

где Н — вещественная матрица, аналитическая вместе с об­ ратной в правой полуплоскости, Т — симметрическая мат­ рица с постоянными элементами.

1 Факторизация (5.20), в отличие от приведенной в предыдущих главах, где Т = Е, имеет место для любой параэрмитовой матрицы, не имеющей вместе с обратной полюсов нечетной кратности на мнимой оси.

127


Матрица Ф, обращающая в нуль первую вариацию функ­ ционала (5.19) и имеющая полюсы только в левой полупло­

скости, запишется в виде Ф = —Я - 1 (Ко + К+) Г-1 , где аналитическая вместе с обратной в правой полуплоскости матрица Г получается в результате факторизации матрицы 5ф, т. е. = ГГ..., а матрицы Ко и К + определяются из оазложения

Ко + К+ + К - = (T~lH -'Q -lN:i;R - НА) Р -'Г (5.21)

(элементы матриц Ко — полиномы от s, К + — правильные дроби с полюсами только в левой полуплоскости, К - — правильные дроби с полюсами только в правой полупло­ скости).

Решение задачи, т. е. искомая матрица W, определится формулой

W = [(К0 + К+) Т~ХМ - Н ВГ' 1(К0 + К+) Т~'Р + НА].

(5.22)

Как и в задаче оптимальной стабилизации, матрица не зависит от конкретного выбора матриц А и В, удовлет­ воряющих требованию аналитичности матрицы Z вместе с обратной в правой полуплоскости. Доказательство этого

утверждения

полностью совпадает с доказательством, при­

веденным в гл. 1 .

 

 

Если элементы матрицы

константы, то решение рас­

сматриваемой задачи принимает вид

 

W = [(К0

+ К+) М — Н В ]-1[(К0 + К+) Р + НА],

(5.23)

Ко +

К + + К - = (T -'H : ]Q-'N*R - НА) Р-'

(5.24)

и не зависит от элементов матрицы Если к тому же параметры объектов позволяют выбрать

матрицы А и В в виде А = Е, В = 0, то процедура опре­ деления матриц Ко и К+ значительно упрощается. В этом случае Q = N и, согласно (5.20),

я = т-'я.-'я:1[NJRN + д; (S) ЛДр (5)1я - 1=

=T -'H -'R + T -'H -'P JA -'A M -'P ,

аматрицы Ко и К+ определяются из разложения

Ко + К + + К - = — Т -'Н -'Р^М -'А М '1.

(5.25)

Отметим, что, согласно результатам гл. 2, решение за­ дачи в виде (5.23) может быть использовано при решении

128


игровых задач с ненулевыми начальными условиями и ну­ левыми внешними возмущениями.

Примеры. I. Рассмотрим игровую задачу о сближении двух объектов. Движение первого (догоняющего) объекта описывается уравнением хг = их -f- ф1( а движение второ­

го (убегающего) — х2 = и2+ ф2.

Возмущения і ^ и ^ являются стационарными случайны­ ми процессами с нулевым математическим ожиданием и ма-

ѲІ

трицей спектральных плотностей 5ф =

, где Ѳх и

О

0 2J

Ѳ2— константы.

 

При ограничениях на управляющие воздействия

(и2) = к 2

требуется найти закон управления [ult и2У = W [хи х2У (матрицу W передаточных функций регулятора), обеспечи­

вающий минимаксное значение (min max) функционалу "і «і

е0= <(Xi — x2)z).

Решение задачи сводится к исследованию стационарных

точек функционала

е = ((хг х2)2) +

(и2) -f- %2 (и\), где

множители Лагранжа 1, >

0 и К2 <С 0.

 

Применительно к рассматриваемой задаче величины,

входящие в (5.17) — (5.19),

имеют вид

 

х = [х1, х2}',

u =

[ul t u2]',

Ф = [фІ;

фа]', Р =

diag(s; s),

М = diag {1;

1},

R =

 

А — diag

Ä,2}.

Приступим к решению задачи. Требование аналитич­

ности матрицы Z~l в правой полуплоскости удовлетворим, положив А = Е2, В = 0.

Так как Др (s) — s2, N = sE2, Q — sE2, то

G*G = N,RN + А; (S ) AAp (s) = -

'1 — M 2

1

s2

1

1 - M

T

 

 

 

 

' 1 — V

s

1 '

 

 

 

Q7%GQ~1

 

1 - K s \

 

 

 

1

 

 

 

Для факторизации последней матрицы в

виде (5.20),

т. е. для

определения аналитической

вместе

с обратной

в правой

полуплоскости и вещественной матрицы

Н и

9 3-582

129



симметрической матрицы Т, воспользуемся алгоритмом Дэ­ виса.

Детерминант

матрицы

т-1 п

 

 

равен X ^ s

2

|s

2

Q*.

G^GQ

 

 

Xj ~i~ Я3

, и, следовательно,

факторизация в виде (5.20)

KK

 

 

 

 

Яі 4- Х^

 

Л

 

 

 

 

 

 

 

 

возможна,

если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что соответствует

неравенству

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I Яа I

1|.

 

 

 

 

 

(5.26)

Нули определителя

матрицы

Н:

hx — 0,

2

=

 

k =

___ ■|уГ

К 4- Я,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

XjXa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Умножив последовательно матрицу Q. G^GQ

 

 

справа

на матрицы Тг и Тг,

а слева на

Д *

и Гг», где

 

 

 

 

 

 

 

 

s

 

Т —

 

1

Я і

+ х а

 

 

 

 

 

 

 

 

 

 

 

 

 

s-\-k

 

 

 

 

 

 

 

 

_1_

 

1 г

 

 

 

 

 

 

 

 

 

 

 

 

U

S

 

 

 

 

О

 

s -f- к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим

 

 

 

 

 

 

 

~j~ %2

 

 

 

 

 

 

 

 

Тп 1

—і/

 

-I— 1 г

 

^ 2

 

 

 

 

 

 

 

 

 

ТгТ2 =

 

Я2

Я#2

 

 

 

 

 

Таким

образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я2

 

 

 

 

 

 

 

 

 

 

 

Т =

1

+ ^ 2

 

 

 

 

 

 

 

 

 

 

 

 

* 2

 

Я2

1

%п

 

 

 

 

 

 

 

 

 

s +

%2

 

ь

 

 

 

 

 

 

Н =

ТГ1 1 =

Яг-Т Xj

 

 

Яі + х2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

Ja

 

 

s +

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

Теперь,

согласно

(5.25),

 

 

 

 

 

 

 

 

 

 

 

К0 +

К + + К - =

 

XjX2

s k

'

Я2

Х.2

'

 

*

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

— Я

Я

 

 

 

 

 

 

Ях (s — k)

 

Я2

(s k)

 

 

 

 

 

 

 

 

kJ

 

 

я2(s

 

Х„

k

 

 

 

 

 

 

 

 

 

 

 

X j

-f-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X,

 

 

 

 

 

 

130