Файл: Ларин В.Б. Синтез оптимальных линейных систем с обратной связью.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 13.07.2024
Просмотров: 149
Скачиваний: 0
|
|
|
|
|
|
|
|
«4-1 |
|
|
|
|
|
|
|
|
|
'X |
|
M S) |
£ |
(s) + |
S PA«A (S) 2 ( - |
1)" |
8! |
|
|||||
|
|
|
|
|
|
|
fc=i |
|
|
|
|
|
|
||
|
|
|
ll-l |
* |
|
|
» /«, |
/z -|- 1 \ |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
,'n, |
n + |
1\ |
|
|
|
ft=i |
|
|
|
M s ) |
Ä * W f f ( * ) « l Af |
n + 1J(*) + |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
(1 - 1 |
|
|
|
|
|
|
|
|
|
|
+ |
(— |
l)fe |
' nk (s) 2 |
|
гсПі (s) (— |
1)' m ^ |
* |
j J |
(s) |
|
|||||
|
|
( - |
i),l+1 |
d-i |
|
Wt (s) |
|
. |
|
In, n -]- |
1\ |
||||
- lim |
|
|
|
|
|
|
|||||||||
_ |
. |
, |
, |
|
|
|
g*(s) |
+ |
Pi |
( - » |
“ |
l i . e + l |
(S)' |
||
S-HI, |
/"> |
,l + |
1 |
|
|
|
что совпадает с (2.50) с точностью до обозначений, так как
(__ ])“+Ч-і |
т |
11, п + |
1 |
||
/п, |
п + |
|
|||
1 |
|
і, и + |
( Р 7 % , |
||
U . |
п + |
I |
(S) |
1 |
|
|
|
||||
|
|
|
|||
где (Р 71 Ь)і — і-я |
компонента вектора Р 7 1Ь. |
Таким образом, есть основания утверждать, что задачи синтеза оптимальных систем стабилизации объектов при стационарных случайных внешних возмущениях могут быть сведены к соответствующей задаче аналитического кон струирования регуляторов большей размерности, если по грешностями измерений фазовых координат объекта можно пренебречь.
Наконец, докажем для случая произвольного числа управляющих воздействий еще одно утверждение о связи задачи аналитического конструирования регуляторов при ненулевых начальных условиях с задачей оптимальной ста
билизации. |
|
|
|
|
|
|
||
|
Пусть для объекта, движение которого описывается ка |
|||||||
нонической системой дифференциальных уравнений |
х = |
|||||||
— Ах -f- Ми -f- ф, закон |
управления, обращающий в нуль |
|||||||
первую |
вариацию функционала е = (x'Rx) + |
(и'Си) |
на |
|||||
классе |
устойчивых |
замкнутых |
систем, имеет вид и — Wx, |
|||||
где |
X = |
[х1%хг, |
|
хпУ, |
и = [иѵ иг...........ит]', ф = [фі. |
|||
ф2, |
. . . , |
фл]', |
=- |
[5,7 (со)], |
Sij (со) = const, |
А, М, R, С |
78
и W — матрицы п X п, п X т, а X п, т X т к т X п
соответственно, элементы которых константы. Утверждается, что этот закон управления обеспечивает
обращение в нуль первой вариации функционала в соот ветствующей задаче аналитического конструирования регу ляторов при ненулевых начальных условиях, а именно, если движение объекта описывается системой дифферен
циальных уравнений |
у = |
Ау + |
Мѵ, у (0) = у0, то |
закон |
|||
управления |
|
|
|
|
|
|
|
|
|
|
|
и — Wy |
|
(2.60) |
|
обращает в нуль первую вариацию функционала |
|
||||||
|
|
I |
j {y'Ry + |
v'Cv) dt. |
|
|
|
|
Стационарное значение функционала е можно записать |
||||||
в виде |
|
|
/со |
|
|
|
|
|
|
|
|
|
|
|
|
|
« = |
(x'Qx) |
= — $ |
Sp |
ds, |
|
|
|
|
|
|
—/оо |
|
|
|
а функционал |
/ при |
законе управления |
(2.60) — в |
виде |
|||
|
|
оо |
|
|
/оо |
|
(2.61) |
|
I = |
Jy'Qydt = |
- ± j - |
j y0FX.,QFxy0ds, |
|||
|
|
0 |
|
J |
—joo |
|
|
где |
Q = R + W'CW, |
Fx = |
(sE — А — MW)~X. |
|
Используя возможность циклических перестановок мат риц-сомножителей под знаком Sp, перепишем функционал
(2.61) так: |
|
/оо |
/оо |
1 = |
2л/ |
J |
y'oFx*QFxy<>ds = |
j |
Sp (y0Fx*QFxy0) ds = |
|
—/со |
|
y' oo |
|
|
|
со |
|
|
||
|
Sp г |
/ |
-J |
|
|
|
/J |
|
|||
|
ds = - |
j J Sp (F№QFXS) ds, |
|||
|
— /оо |
|
|
|
— /со |
где элементы матрицы S = |
y0yQ— константы. |
Но так как матрица Ц7, получаемая в результате ре шения вариационной задачи (обращающая в нуль бе), не зависит от матрицы 5ф, когда элементы последней кон станты х, то и первая вариация б/ также обратится в нуль, что и требовалось доказать.1
1 См. § з гл. 1.
79
ГЛАВА 3 СИНТЕЗ СИСТЕМ СТАБИЛИЗАЦИИ ПРИ ИЗМЕРЕНИИ КООРДИНАТ С ПО МЕХАМИ
В рассмотренных выше задачах предполагалось идеальное измерение координат объекта, необходимых для формиро вания закона управления. На практике, однако, довольно часто уровень погрешностей при измерениях тех или иных координат настолько значителен, что ими пренебрегать уже нельзя. В связи с этим представляется целесообразным учесть погрешности в измерениях при постановке и реше нии задач синтеза оптимальных систем стабилизации, что и будет сделано в настоящей главе для случая, когда коор динаты объекта измеряются с помехами, предполагаемыми стационарными случайными процессами. Сначала будет приведено общее решение задачи при т управляющих воз действиях (§ 1), а затем рассмотрен частный случай, когда
т = 1. Метод |
решения аналогичен методам, |
изложенным |
|
в § 3 гл. 1 и § |
1 гл. 2. |
|
|
|
§ 1. |
Р Е Ш Е Н И Е З А Д А Ч И П Р И т У П Р А В Л Я Ю Щ И Х |
|
|
|
В О З Д Е Й С ТВ И Я Х |
|
Пусть движение объекта, |
как и в § 3 гл. 1, |
описывается |
системой обыкновенных дифференциальных уравнений с постоянными коэффициентами
Рх — Ми-\- ф |
(3.1) |
(см. (1.51)).
Требуется определить уравнение регулятора в цепи об ратной связи так, чтобы замкнутая система объект -)- регу
лятор была устойчива и функционал |
|
е — {x'Rx) + (и!Си) |
(3.2) |
достигал минимума. |
|
80
Считаем, что необходимые для формирования закона управления координаты объекта измеряются с аддитивными помехами, которые предполагаются стационарными слу чайными процессами, т. е. уравнение регулятора ищем в виде
W0u = W (х + ф), |
(3.3) |
где Ф = Іфі ( 0 . •••. Фп ( 0 1' — вектор |
ошибок измерения |
координат объекта, компоненты которого (ошибки ф, (/)) —
стационарные случайные процессы с нулевым |
математиче |
|||||||||
ским |
ожиданием |
и |
матрицей |
спектральных |
плотностей |
|||||
S ,, |
(©)• |
|
|
|
|
|
|
|
|
|
|
Используя преобразование Лапласа к уравнениям (3.1) |
|||||||||
и |
(3.3), получаем |
|
|
|
|
|
|
|
||
|
|
|
Р (s) X(s) = М (s) и (s) + ф (s), |
|
(3.4) |
|||||
|
|
|
и (s) = W (s) [X (s) + |
ф (s)], |
|
(3.5) |
||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W (s) = Го“ 1(s) W (s). |
|
(3.6) |
||||
|
Введем |
матрицы |
передаточных |
функций |
F f, F f, |
F f |
||||
и F f |
между координатами системы х (s) и и (s) и внешними |
|||||||||
возмущениями ф (s) |
и помехами измерений |
ф (s): |
|
|||||||
|
|
|
X (s) = Fx (s) ф (s) + F f (s) ф (s), |
I |
|
|
||||
|
|
|
U(s) = |
F f (s) ф (s) - f F f (s) ф (s), |
J |
|
|
|||
которые, согласно |
(3.4) и (3.5), определятся формулами |
1: |
||||||||
|
|
|
F * = ( P - M W ) - \ |
) |
|
|
|
|||
|
|
|
|
F f |
= W ( P ~ MW)-\ I |
|
|
|
||
|
|
F f = |
(P — MW)~XMW = |
(P — MUT)“ 1P - |
F, |
|
||||
|
|
F f = |
W (P — Л4Й7)-1 P. |
|
|
|
|
|
||
|
Поскольку замкнутая система должна быть устойчивой |
|||||||||
и, |
следовательно, |
элементы матриц |
F f, F,f, |
F f и F f |
не |
должны иметь полюсов в правой полуплоскости, то, исполь зуя преобразование Фурье (полагая s = /со), функционал
1 В дальнейшем аргумент s опускается.
6 3 - 5 8 2 |
81 |
|
(3.2) можно записать в следующем виде:
|
/ ©о |
|
е = j . |
j Sp l(F lR F x -I- Ff.C Ff) S,„ + |
|
|
— /оо |
|
+ |
(F?,/?Ff + F f.C F f) S,p]ds. |
(ЗЛО) |
Таким образом, задача сводится к определению матрицы № такой, чтобы замкнутая система была устойчивой, а функ ционал (ЗЛО) достигал минимума.
Вариации бW при минимизации функционала (ЗЛО) сле дует ограничить так, чтобы соответствующие вариации
бF f, 6F,f, |
бF f и бF f функций F f, F f, F f и F f имели бы |
полюсы в |
левой полуплоскости. |
Так как, согласно (3.8) и (3.9),
F f = F fF — Е,
(3.11)
F f = F fF ,
то из физической реализуемости функций F f и F f и их вариаций следует физическая реализуемость функций Ff и F f и их вариаций.
Матрицы F f и Fu удовлетворяют уравнению связи
F F f — M Ff = Еп |
(3.12) |
(это уравнение получается после подстановки (3.7) в (3.4)), следовательно, п (т + п) их элементов могут быть выра жены через пт независимо варьируемых функций.
Как и в § 3 гл. 1, введем матрицу Ф размера т X п следующим образом:
ф == /[F f + BFu, |
(3.13) |
где Л и В — полиномиальные матрицы т X п и т х т,
соответственно. Тогда все элементы матриц F f и F f могут быть выражены, согласно (3.12) и (3.13), через тп свободно варьируемых элементов матрицы Ф:
F f = |
F -1 + MQ-1(Ф — AP~l), |
(3.14) |
F f = |
Ap (s) Q~l (Ф — AP~l), |
(3.15) |
где |
|
|
|
Дp (s) = det P, |
(3.16) |
82