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

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

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

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

Добавлен: 13.07.2024

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

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

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

Итак, пусть движение объекта описывается системой обыкновенных дифференциальных уравнений с постоянными коэффициентами

Рх = та (t) 4- ф,

(3-47)

где X — lxt (f), ..., л-,, (/)]' — п-мерный вектор координат объекта, и (t) — координата регулятора (управляющее воз­ действие), ф = [фі (0, ...,ф„ (/) 1' — n-мерный вектор внеш­ них возмущений, компоненты которого ф, (/) — стационар­ ные случайные процессы с нулевым математическим ожи­ данием и дробно-рациональной матрицей спектральных плотностей s^(w), Р и пг — матрица п X п и п-мерный вектор-столбец (матрица п X 1) соответственно, элементы которых Рі] (р) и т1 (р) — операторные полиномы от р.

Требуется определить уравнение регулятора в цепи об­ ратной связитак, чтобы функционал

е = (x'Rx) - f с (и2) -ф- (и2)

(3.48)

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

Считаем, что координаты объекта, необходимые для фор­

мирования закона управления, измеряются с

аддитивными

помехами, т. е. уравнение регулятора ищем в виде

 

w0(р) u ( i) = w (X - f

ф),

(3.49)

где ф = [фх (/),

..., ф„ (01' — вектор

погрешностей изме­

рения координат

объекта, компоненты которого (ошибки

Фг (0) — стационарные случайные процессы с

нулевым ма­

тематическим ожиданием и матрицей спектральных плот­

ностей 5 Ф (со).

 

 

Используя преобразование Лапласа

к уравнениям (3.47)

и (3.49), получаем

 

 

Р (s) x(s) = m (s) и (s) 4- ф (s),

и (s) =

w (s) [x (s) 4- ф (s)],

где

 

 

° ' (s)== w

S (s )-

(3-50)

Таким образом, задача сводится к определению вектора w такого, чтобы замкнутая система была устойчива и функ­ ционал (3.49) достигал минимума.


Согласно (3.28), решение задачи определяется формулой

w =

 

(с — s3) Д

(s)

~j- (k_— [Q— /-)-) D * tn

—1

 

 

X

 

 

к* (s) q* (s)

 

X

0b - - l 0 - l+ ) D ~ ' P

h* (s) q*(s)

 

 

(3.51)

где

 

 

 

 

 

 

 

 

 

 

 

 

Af, (s) =

det P,

 

 

 

(3.52)

 

 

 

 

 

 

 

 

 

 

 

n =

 

Ap (s) P~lm,

 

 

 

(3.53)

 

 

 

q (s) =

Ap (s) ß (s) +

an,

 

 

(3.54)

 

 

 

nh* (s)S )ti(s)h --

q4s)qis),

 

 

(3.55)

 

g* (s) g (s) =

ntRn +

Ap (s) (c — s2) Ap (s),

(3.56)

 

 

 

DD* — Sip -}- PSyP

 

 

(3.57)

k0

k-Sf-

 

k - =

 

 

1

n:l.R — h(s) а

P-'D,

(3.58)

 

h*(s) q* (s)

 

/0 +

/+ + /——

 

Л* (s) q* (s)

n*RSipP ,P , .

(3.59)

Здесь элементы вектора-строки а и скаляр

ß (s) — поли­

номы от

s,

 

удовлетворяющие

требованию аналитичности

матрицы Z =

Р ~ т

 

 

 

 

 

правой

полу­

ß (s) j вместе с обратной в

плоскости,

т.

е. полином

det Z = Ар (s) ß (s)

 

an =

q (s)

должен

быть

гурвицевым

полиномом.

 

 

 

Пусть полином g (s) имеет нули только в левой полупло­ скости. Тогда функция h (s), имеющая все нули и полюсы в левой полуплоскости, согласно (3.55), определится фор­ мулой

 

 

h(s)

g(s)

 

 

 

 

 

 

<?(«)

 

 

 

 

 

 

 

 

 

 

 

и соотношения (3.51), (3.58) и (3.59) можно записать

так:

(с- s

2) Др (s)

+ ( / е _ - / 0~ / + ) D-'m

 

 

w =

 

•X

 

S* (S)

 

 

 

 

 

 

(£_ — /„ — /+) О“ 1/3

1

'h-R

 

(3.60)

8* (s)

 

 

 

 

 

 

h +

At- + /_

 

1

 

 

 

(3.61)

 

g* (S) 'h-RSqP. , ; Ö ,

,

 

 

 

 

 

 

 

(3.62)

93


г д е

 

 

1

giß)

(3.63)

g*(s) n*R

9(s)

V.

 

 

В предыдущем параграфе была доказана независимость

решения задачи (матрицы W и ет іП) от произвола в

выборе

матриц А и В. Однако формул для W и ет |П, явно не зави­ сящих от А и В, в общем случае получить не удалось. В за­ даче стабилизации объекта одним управляющим воздей­ ствием удается выразить решение только через исходные данные, т. е., как и в § 1 гл. 2, получить формулу для w, явно не содержащую полиномы щ (s) и ß (s).

Действительно, как видно из формул (3.60) — (3.63), полиномы а £ is) и ß (s) неявно входят лишь в выражение для вектора и_, т. е. если возможно получить соотношения для вектора и_, явно не зависящее от а £ (s) и ß (s), то вектор передаточных функций регулятора w также определится формулой, не содержащей a t (s) и ß (s). Поскольку вектор и_ (3.63) совпадает с соответствующим вектором (2.15) в § 1 гл. 2, для элементов которого правомочны формулы (2.17)

 

 

 

м

bj іщ) + Cj (gt)

 

N

d j jgt)

 

 

« /- (s) =

2

+

E

 

 

 

S — gt

 

s — gt

 

 

/=1

 

 

 

r = M - H

 

 

 

 

(/ =

1, 2...........n),

 

 

где постоянные коэффициенты bj {gt), Cj (gt)

и dj (gt) имеют

вид (2.18), (2.23) и (2.19) и не зависят

от полиномов щ (s)

и ß(s), то элементы

вектора

и_ (3.63) определятся эти­

ми же формулами и тоже не

зависят от полинонов а,1(s)

и

ß (s).

 

 

 

 

 

 

 

 

 

Примеры.

 

 

 

 

 

 

 

 

I. Пусть Р — а — s, т = М = const,

5ф (со) = Яф, 5 Ф(со)=

 

 

 

 

е =

с2 (и2) +

(и2).

 

 

 

Согласно (3.52),

(3.53),

 

 

 

 

 

Ар (s) = det Р =

Р (s) — а — s,

п =

Ар (s) Р~'т = М.

Полагая

ß (s) = 0,

a

(s) =

получаем q (s) = Ap(s)ß(s) +

+

а /г =

1. Из (3.56)

и

(3.57)

имеем

 

 

 

 

g* (s) g (s) =

+ s) (с2 — s2) (а — s),

94


т. е.

 

 

 

8

(s) = + s) (а +

s),

 

 

 

 

 

 

 

 

 

 

 

 

DD* = 4

+

(а2 -

s2) ( - b t ) 2 =

( - ^ ) 2 (d2 - s2)

 

(d = | /а 2 + Г > 0 ) , т. e. D = - ^ - ( d + s).

 

Согласно (3.61)—(3.63),

 

 

 

 

 

 

 

 

 

 

 

 

(с - f

s) (а 4 - s)

__

 

2 а

(с +

а)

 

 

 

 

 

 

М

s)

 

М

(а —

s) ’

 

 

 

k—=

 

2 а

(с +

о)

X (d +

s)

 

 

 

 

 

 

М

(а —

s)

 

 

 

 

 

 

 

(а +

с) (д + d)

/0 =

/+

= /_ =

0.

 

 

M X .

 

 

(а — s)

 

 

 

 

 

 

 

 

 

 

И, наконец, подставив необходимые величины в (3.60),

(3.30),

(3.29),

получим

 

 

 

 

 

 

 

W ~

 

 

 

 

2 а (а -f- с) (а +

d)

 

 

 

 

М

(s2 +

(2a +

c +

d)s +

[2a(a +

c +

d) +

cd]}

 

 

У°°

 

 

 

 

 

 

 

 

 

 

 

ßmin =

у

j* k - (S)

(— S) ds =

 

(а - f c f (а - f d f Sv.

 

 

— /o o

 

 

 

 

 

 

 

 

 

 

 

При отсутствии помех в измерении координатых (X -*оо, S(p-> 0), значения w и ет іп совпадают с соответствующими значениями в примере III § 2 гл. 1.

II.Движение объекта описывается дифференциальным

уравнением первого порядка х -J- ах =

и

ф.

 

Требуется найти закон управления в цепи обратной

связи (передаточную функцию регулятора w) и =

w (х -1- ф)

так, чтобы при устойчивой замкнутой

системе

объект +

+ регулятор дисперсия х была минимальной, т. е. критерий качества имел виде = (х2).

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

ческим ожиданием и спектральными

плотностями 5ф (со) =

=

Л 2/(ц2 +

со2)

и 5 ф(со) = Д2/(ѵ2 +

со2).

 

 

 

В

принятых

обозначениях Р

s) = а - •s, m =

1.

 

Аналитичности матрицы Z —

P

— tri

вместе с

обрат-

 

а

ß

ной

в

правой

полуплоскости добьемся,

положив

ß = 0,

а

=

1,

т. е.

q — det Z = 1.

 

 

 

 

95


 

Согласно

(3.52) —

(3.57),

п =

1,

g — 1,

 

 

 

 

 

£>Д> =

ф _______(а2— s3) В2

 

 

 

 

откуда

 

 

*

(I2 -

 

 

V2 —

S3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

=

^52 ~l~ ^is

 

 

 

 

 

где

 

 

 

 

 

(H +

s ) ( v + s )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я0 =

У 5 2а2(.і2 +

ЛV s > О,

 

 

 

 

 

Вг =

У ß 2 (а2 +

У ) +

Л* +

2ВВ0>

0.

 

 

Из

(3.61) — (3.63)

имеем

 

 

 

 

 

 

 

 

 

lo + l+ +

1- =

 

 

В2(а — s) (|х — s)

В0)

'

k—— V—=

0.

 

 

(ѵ +

s) (Äs3 — BlS+

 

Следовательно,

/0 =

0,

 

 

 

 

 

 

 

 

 

 

 

/+

= - ■

 

5 2 +

ѵ) (р + у)

 

 

 

 

 

 

(ѵ +

s) (ßv3 +

ßxv -f- ß0)

 

 

 

И окончательно,

согласно

(3.60),

 

 

 

 

 

 

 

 

w

(s) =

 

4V2 +

yts +

y„

 

 

 

 

 

 

 

s) (fl +

 

 

 

 

 

 

 

где

 

 

 

 

( ( i +

 

V) ( p

+ V)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vo = ац (a +

v) (p - f

v ) ------ § -

(v2 +

 

+ -З в -j

,

 

Vi = ifl +

v) (a + P) (P +

v) —

 

(v2 +

 

V + A - J ,

 

 

 

Va =

V (а + p) + ац —

v ---- •

 

 

 

Рассмотрим случай, когда внешнее возмущение исче­

зает (Л -*■ 0, ß

 

0). Здесь возможны два

варианта.

Пусть

а >

0, т. е. объект устойчив,

тогда у0, ух, у2 -> 0 (так

как

ß 0 —V Вар,

Вх ->

В {а +

р))

и,

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

w -> 0,

(х2) ->■ 0.

 

 

 

 

 

 

 

 

0),

то

В0 -> Вщі,

Ву -■>-

 

Если объект неустойчив (а <

- > ß (H — а), a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а-)- V

ф о

(х2) ф 0 .

 

 

 

 

 

 

 

 

 

 

 

 

'

 

 

 

 

 

Казалось бы, что и в этом случае

w должно быть

рав­

но

нулю, но, положив

w =

0,

мы

не удовлетворили

бы

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

96