книги из ГПНТБ / Красовский, А. А. Фазовое пространство и статистическая теория динамических систем
.pdf150 |
РЕ Ш Е Н И Е Ф П К -У РА В Н Е Н И Я МЕТОДОМ РЯ Д О В [ГЛ . XII |
других по величине, то подбором а можно, по крайней
мере для начальных распределений, близких к (3.111), добиться высокой вероятности нахождения главного эк стремума и малой вероятности прихода к неглавному экстремуму. В этом случае данный алгоритм поиска может применяться без дополнительных процедур. Если
же экстремумы близки по величине, то согласно получен ному решению вероятности прихода в окрестности этих экстремумов будут сопоставимы по величине. В этом слу чае путем статистической обработки можно выполнить анализ найденного экстремума и в случае его отличия от искомого, осуществить переброс начальной точки поиска в соседний «квадрат» (заранее выбранный район поиска). Подобные операции «глобального» поиска широко извест ны [3.20], [3.13], [3.21].
Заметим, что, помимо задачи поиска экстремума, урав нениям (3.110) и изложенной теории можно дать совер шенно другую интерпретацию. Действительно, уравнения
§ 3.4] |
СТАТИСТИЧЕСКОЕ И ССЛЕДОВАН ИЕ СИСТЕМ |
151 |
(3.110) при п — 3 представляют собой уравнения движе
ния броуновской заряженной частицы в электростати ческом поле с потенциалом F. В этом случае xv x2, ха — геометрические координаты частицы, кх = ка = ка —ве
личина, пропорциональная заряду частицы и обратно пропорциональная коэффициенту вязкого трения (силой инерции в сравнении с силой Стокса пренебрегаем),
Sx = iS2 = |
S a —■ спект |
дх. |
|
|
|
ральные плотности ,флу- |
|
|
|||
ктуационного теплового |
|
|
|
||
воздействия среды. |
|
|
|
||
Решение (3.111) оз |
г- ■^)— |
— *{х]— |
Ц |
||
начает, что |
плотность |
||||
распределения g частиц |
|
|
|
||
соответствует экспонен |
|
|
|
||
циальной |
функции по |
|
|
|
|
тенциала. |
Это может ис |
|
Рис. 3.6. |
|
|
пользоваться, |
напри |
|
|
|
|
мер, при рассмотрении электростатических улавливающих фильтров (электроулавливания), электроокрашивания, электросепарации (о данных технологических процессах см., например, [3.22]).
Система экстремального регулирования с синхронным детектированием при случайных сигналах поиска. Из вестная структура системы экстремального регулирова ния с синхронным детектированием изображена на рис. 3.6. Если объект безынерционный и описывается квадра тичной формой, то уравнения системы запишутся в виде
П |
|
А±1 = — кх 2 ajk(Axj + 8х,)(кхк’ 4- 8хк) Sx— *iB, (3.124) |
|
i = 1 , 2 , . . ., n, |
|
где Дач = x\ — Xi3 — отклонения от |
экстремума, ач — |
выходные величины интегрирующих |
звеньев, xia — ско |
рости дрейфа эстремума, 624 — поисковые составляющие, кi — коэффициенты усиления каналов, a.jh — постоянные
величины. Будем считать, что поисковые колебания пред ставимы с помощью формирующих фильтров типа аперио дических звеньев:
Tlabti + bxt = £*п, i — 4, 2, . . ., re, (3.125)
152 |
РЕ Ш Е Н И Е Ф П К -У РА В Н Е Н И Я МЕТОДОМ РЯ Д О В (ГЛ . III |
гДе £гп — некоррелированные белые шумы, Тщ — задан
ные постоянные величины. Компоненты скорости дрейфа точки экстремума также полагаем случайными и подчи няющимися уравнениям вида
|
Тig&io “Ь ±ia = iija. |
|
(3.126) |
|
где |
— белые шумы, не коррелированные |
с |/ п, Т 1Я— |
||
постоянные. Используем обозначения: |
|
|
||
Xi = |
Axi -f- 6Xi, ^j+n = |
®i+2n = |
i |
1 , 2, . . .,ra. |
|
|
|
|
(3.127) |
Тогда уравнения (3.124), (3.126) можно представить в виде
пЗп
+ |
S |
a>rc? |
|
S |
a i)klx ix kx l = |
|
(3.128) |
||||
где |
3=1 |
|
|
j,k,l=1 |
|
|
|
|
|
||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
при |
|
i |
га, |
/ |
= |
i + |
га, |
|
|
|
т,„ |
|
|
||||||||
|
при |
|
i |
га, |
/ |
= |
i + |
2га, |
|
||
|
1 |
|
|
(3.129) |
|||||||
|
1 |
|
при |
w< i = |
/ < |
2га, |
|
||||
|
|
|
|
||||||||
|
?’in |
|
|
|
|
|
|
|
|
|
|
|
7’. |
|
при |
2га<^ i = |
7 ^ |
|
|
|
|||
|
1э |
|
|
|
|
|
|
|
|
|
|
и все остальные ац = О, |
|
|
|
|
|
|
|
||||
a»3ki = |
М м |
|
при |
i, ft, Z< |
«. |
/ = |
1 + ге» |
(3.130) |
|||
a»k3l = |
|
|
|
|
|
|
|
|
|
i3ltI |
|
a iklj — aHlk — |
a ilkj = |
a il)k |
|
|
|||||||
и все остальные |
а^ь; = 0, |
|
|
|
|
|
|
||||
|
|
|
ф— при |
i < |
|
|
|
||||
|
|
|
■'in |
|
|
|
|
|
|
(3.131) |
|
li = |
|
^in |
при |
га ^ |
^ |
^га, |
|||||
|
Т- |
|
|
|
|
|
|
|
|||
|
|
|
**п |
|
|
|
|
|
|
|
|
|
|
|
|
- |
при 2га < |
i < |
3n- |
|
|||
|
|
|
‘ ■* |
гэ |
|
|
|
|
|
|
|
Уравнения (2.14), записанные для iV = |
т- е- Для пРеД" |
5 3.4] СТА ТИСТИЧЕСКОЕ ИССЛЕД О ВА Н И Е СИСТЕМ 153
ставления искомой логарифмической плотности полино
мом четвертой степени, |
в данном случае имеют вид |
|||||
|
|
|
3п |
|
|
|
А } |
2 ^ p iA p |
2 |
(4-pgi "4“-ApiAy) — 0 , |
|
||
|
P =1 |
|
P, 4=1 |
|
|
|
|
3n |
|
|
|
|
|
■A-ik — 2 (a P iA p k |
+ apk-Api) |
|
|
|||
3d |
P=1 |
|
|
|
n |
|
|
|
|
|
|
||
2 ‘-*P4(3-4p4i/c_f‘ % A pik^q~\- |
A p iA q k ) = 6 2 ®PpilO |
|||||
P. 4=1 |
. |
|
|
p=1 |
|
|
|
3/i |
|
|
|
|
|
A ik i |
2 |
iflpiApiti "Ь О'ркАрЦ -(- ®p;^4pi/i) |
|
|||
|
P=1 |
|
n |
|
|
|
|
|
n |
|
|
|
|
|
— 3 ^ (IpilclAp 2 |
S p q |
[ЗЛрЛдЦ.; |
(3.132) |
||
|
|
P=1 |
P .4 = l |
|
|
|
|
3n |
-j- {A p jA giii "l- ApfcAgii 4" ApiAqin)] = |
0 , |
|||
|
|
|
|
|
|
|
Aikla—^ i(apiApkii-\"apkApiis-\-<lpiApii!i-\- apsApiki) |
|
|||||
|
p= 1 |
|
|
|
|
|
ЗП |
|
|
|
|
|
|
2 |
(apiklAps + |
dpiksApi "f" |
"b O’pklsApi) |
|
||
P=1 |
|
|
|
|
|
|
|
2 |
*^434!A pi Aqftis “i“^4p k A g iis~\~A p iA q jfo A ' |
|
|||
P.4 = 1
-f A PsA q Ucl-\-
I
-j-(i4pi(t^ g(a4-Api/c4 Qil+ ^ p sk/leij+
+ А р ц А д /и + А Р1 ,А а1к |
A plsA g ik)] = 0. |
Если начальное распределение является симметричным, так что
A j (0) = A i = О, А 1ы (0) = А ш = О,
эти уравнения имеют решение, при котором
Af (0 — 0, A{hl (f) — О,
154 |
РЕ Ш Е Н И Е Ф П К -У РА В Н Е Н И Я МЕТОДОМ РЯ Д О В [ГЛ . Ш |
a A ik, A ikls определяются уравнениями
A ilc —2
v—i
( а Р ^ Р к +a P k A p i) —
зп |
Зп |
|
2 S p q |
( 3 A p qi]( -|- A p i A q i i ) — |
62 a PPiki |
|
||
|
P, 9=1 |
|
|
|
P=1 |
|
|
3n |
|
|
|
|
|
A i k l s |
— 2 ( a P i A p h l s |
~\~a p k A p i i s + d p i A p i b s + |
(3.133) |
|||
|
P=1 |
3n |
|
|
|
|
|
|
|
|
|
|
|
+ a PsApikl) |
S |
(a Pikl-Aps + a piksApi + d p iit A p b ^ |
|
|||
|
|
P= 1 |
|
|
|
|
|
|
|
3n |
|
|
|
Ч- |
®PklsApi) |
|
2 |
^Pq (ApiAqftte -{- |
A p kA q j;s -f- |
|
|
|
P.<I=1 |
|
|
||
+ ApiAqifo -f- ApsAqikl) = 0.
Дальнейшие конкретизация и упрощение уравнений воз можны на основе учета выражений (3.129)—(3.131).
Как и в задачах, рассмотренных в двух предшествую щих параграфах, коэффициенты A ih, A ikis разделим на
три группы. К первой группе относятся коэффициенты, соответствующие только основным фазовым координатам
xt — Ах* -f бxt, i = 1 , 2, . . ., п.
Для этой группы коэффициентов согласно (3.129), (3.130), (3.132)
ЗП
к |
2 |
'Spqi’^Apqiic -f- Aj^Aq^) = |
0, |
|
|
p,q=i |
|
|
|
|
Зп |
^ p q (A p iA q itf s -j- A p kA q U a -(- |
|
|
|
2 |
(3.134) |
||
|
p,q= l |
A PsAqikl) = |
0, |
|
|
|
+ ApiAqixs + |
||
i, k, l,s^.n.
Ко второй группе относятся коэффициенты, имеющие индексы п, т. е, соответствующие поисковым состав
§ 3.4l |
СТАТИСТИЧЕСКОЕ ИССЛЕД О ВА Н И Е СИСТЕМ |
155 |
ляющим bxi и компонентам скорости дрейфа т1ь. Для
этой группы коэффициентов согласно тем же выражениям
A-ik (&ц "Ь &кк) Ajk
3n
— |
2 ^Pq i^Apqiis + ApfAqx) = |
0, |
P.4=1 |
|
|
A\kls — (aii 4" akk + |
all 4“ ass) ^ikls |
(3.135) |
3n |
|
|
|
|
|
2 *^Pe (A p iA q k is -f" AgfrAqiis -j- |
|
|
P .«=l |
“1" ApiAqifcg -f" ApsA.qift[) = |
0, |
|
||
3n '^ i,k ,l,s'^> n .
Третью группу коэффициентов составляют величины A ih, A ihls, у которых часть индексов ^ п, а другие ин дексы |> п. Уравнения для этой группы коэффициентов
оставим в прежнем виде (3.133).
Анализ полученных уравнений в общем виде затруд нен. Основной способ исследования — численное инте грирование этих уравнений. Самое простое, но наиболее грубое приближение получается при учете лишь квадра тичных членов, т. е. когда решение ищется в виде нор мального распределения. В этом случае
П |
зп |
A-ik 2 f a p i A p k ~Ь Q p k A p i) |
2 ^ p q A p i A q k = |
Р=1 |
P,q—1 |
|
3n |
|
— в 2 ®ppik> (3.136) |
|
p=i |
i, к = 1 , 2, . . 3п.
Эти уравнения в матричной форме имеют вид
А — Аа — атА — A SA = р, |
(3.137) |
где |
|
А — I A th ||, а — || aih ||, S = | S ih ||, |
р = || 62aPJ><h |
p = i
156 |
РЕ Ш Е Н И Е Ф П К -У РА В Й Е Й И Я М ЕТОДОМ РЯ Д О В [ГД . I l l |
— матрицы размерности 3n X Зп. Если перейти к обрат ной матрице — матрице вторых моментов
М = - А - 1,
то вместо (3.137) получаем
М + аМ + МоА - М рМ = S, |
(3.138) |
т. е. матричное уравнение Риккати для ковариационной матрицы. Система в рассматриваемом приближении со храняет работоспособность, если решение, соответствую щее установившемуся
состоянию
аМ у + .
A-MyCi1 — МурМу ~ S,
(3.139)
устойчиво. Линеаризуя уравнение (3.138) в ок рестности установивше гося состояния, полу чаем
АМ + ( а - М у р ) А М +
+ A M ( a * — p M v ) = 0 .
(3.140)
Для того чтобы решение АМ = 0 этого линейно
го уравнения было устойчиво, необходимо и достаточно, чтобы характеристические числа матрицы а — М ур имели
отрицательные действительные части.
Для одномерной системы при отсутствии дрейфа эк стремума, как видно из выражений (3.129)—(3.131),
|
|
0 |
&ic<ii |
<7l2 |
>7*2 |
а — • т |
г |
л и |
п |
||
, р = fciaii |
0 1, |
Л = |
т>2 |
||
0 |
тв |
|
|
пп2 |
|
|
|
|
|
1 п |
лл |
Уравнение (3.138) в скалярной форме в данном случае
§ 3.4] СТАТИСТИЧЕСКОЕ И СС ЛЕД О ВА Н И Е СИСТЕМ 157
имеет вид
М ц — тг,— М i2 |
|
|
ц = — , |
|
|
|||
|
|
|
|
|
|
Т1 |
|
|
М у ь |
-1------ тт, |
A f 22 |
"1-----77; |
Л ^12 |
|
|
|
|
|
1 п |
|
п |
|
|
|
|
(3.141) |
|
|
|
— ^1аи (Af?2 -|- М\\Муъ) = |
?,2 |
||||
|
|
|
|
|||||
■^22 Н--- ™--- -^22 — 2&1<Хц.М12М 22 — |
j |
|
||||||
Для |
установившегося |
состояния |
|
|
|
|||
2 |
- - |
-- |
|
|
«Ус |
’ |
|
|
|
■-Л^1/12 — |
2/с1а 11М щ12М 1/11: |
ф2 |
|
|
|||
|
|
|
|
|
П |
|
|
|
----М у Ъ2 -] |
77, |
М у 12 |
^ 1 а 11 |
{ M v l i |
~1~ М |
у ц М y w ) |
j |
|
J n |
|
Jn |
|
|
|
|
|
J п |
— M y i 2 — 2 к г а п М v i 2M У22 = |
j • |
|
|
|
||||
” |
|
|
|
|
■*П |
|
|
|
На рис. 3.7 представлены кривые изменения диспер
сии М п = (A#!*-}- ба^)2 во времени при различных зна чениях постоянных х = /с1а117’п, т] = S'JTn, получен
ные путем численного интегрирования уравнений (3.141). Видно, что увеличение спектральной плотности случай ного поискового сигнала, так же как и увеличение коэффи циента усиления («форсирование» процесса), неизбежно приводит к потере устойчивости.
ГЛА&А IV
РАВНОВЕСНЫЕ РАСПРЕДЕЛЕНИЯ В СИСТЕМАХ С КУСОЧНО-ЛИНЕЙНЫМИ ХАРАКТЕРИСТИКАМИ.
ФПК-УРАВНЕНИЕ И СИНТЕЗ ДИНАМИЧЕСКИХ СИСТЕМ
В данной главе рассматриваются относительно разно родн ы е во п р о сы статистической теории динамических
систем. Сначала выводятся условия существования ку сочно-нормального равновесного распределения в систе мах с кусочно-линейными характеристиками при наличии шумов.
Затем рассматриваются постановки задач синтеза си стем управления и оценивания и пути решения этих задач. Часть затрагиваемых вопросов носит проблемный характер и не имеет исчерпывающего решения.
§ 4 . 1 . Равновесные распределения в системах с кусочно-линейными характеристиками
Предположим, что re-мерное фазовое пространство ди намической системы может быть разбито на конечное чис ло областей C?v, в каждой из которых характеристики системы линейны:
П |
|
/|ч)(хъ . . . ,xn,t) = a(i ] + 2 aikXk, |
(4.1) |
k=l |
|
где a?\ a$ — заданные постоянные или функции време
ни. Области Gv в совокупности занимают все фазовое пространство и разделены заданными непрерывными ги перповерхностями
/'V (жц . . . |
, хп) = Cjx = const. |
(4.2) |
В дальнейшем в качестве таких граничных гиперповерх ностей чаще всего будут рассматриваться гиперплоскости
П |
|
(4.3) |
2 |
= с\>- |
|
г=1 |
|
|
§ 4.1] РА В Н О В Е С Н Ы Е Р А С П Р Е Д Е Л Е Н И Я В СИСТЕМ АХ |
159 |
Поверхности (4.2), (4.3) будут именоваться поверхно стями переключений.
Будем полагать, что шумы, действующие на систему, одинаковы для всех областей фазового пространства.
Тогда уравнения движения в v-й области будут иметь вид
П
|
4" ai ^ + 2 aik x k = £i, |
(4.4) |
|
/с=1 |
|
где |
— белые шумы. Системы данного вида обычно на |
|
зываются системами с кусочно-линейными характери стиками. Подклассом этих систем являются релейные
системы вида
П
x i 4" + 2 aikx k — £ii (4-5)
<с= 1
в которых при переходе из одной области в другую меня-
(v)
ется скачком «постоянная составляющая» щ .
Другим подклассом кусочно-линейных систем являются системы с переменной структурой, исследованные С. В. Емельяновым и его учениками [4.1]—[4.5]. Для этого
подкласса ^ = 0, а граничными поверхностями служат
гиперплоскости, проходящие через начало координат. Часто используемая аппроксимация характеристик эле ментов систем управления в виде ломаных также при водит к уравнениям вида (4.4). Достаточно очевидно, что использование концепции фазового газа для систем с ку сочно-линейными характеристиками правомерно.
Ввиду локального характера ФПК-уравнения оно справедливо внутри каждой из областей Gv:
д In j |
|
^ |
Ду) д In р(ч) |
|
|
|
|
|
|
||
Ш |
|
А |
и |
дх. |
|
|
|
|
|
|
|
|
|
i=i |
|
1 |
|
|
|
|
|
|
|
|
|
___ i_ |
Vi |
о / &1In |
|
d In |
d In |
_ |
|
||
|
|
2 |
.“ |
^ \ |
дх.дх. |
"T" |
dxt |
dx, |
J |
|
|
|
|
|
t,j= i |
' |
г |
j |
|
5 |
i |
/ |
|
_ |
v |
|
a In pW |
|
|
+ 2 |
aikxk) |
a in P<v) |
|||
~ |
^ |
dx, |
~ |
dt |
|
|
dx. |
||||
|
1=1 |
* |
|
|
|
|
|
k=l |
|
|
|
|
|
|
(P In />(v) |
|
d In p (v) |
d In pW ) |
-^i |
(v) |
|||
|
|
г,1=1 |
|
dx,dx, |
|
~ |
dx, |
|
dx. j ~ |
2 л |
a ii > (4.6) |
|
|
|
y |
|
|
|
|
|
l |
|
|
