Файл: Красовский, А. А. Фазовое пространство и статистическая теория динамических систем.pdf

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

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

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

Добавлен: 23.10.2024

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

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

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

§ 4.2] Ф П К -У Р А В Н Е Н И Е И С И Н ТЕ З СИСТЕМ У П Р А В Л Е Н И Я

179

где Ь изменяется от <; t2 до t2- Сами оптимальные в

смысле минимума функционала (4.47) управления имеют вид

П

 

щ —— (А\ -(- 2 AiicXft

2 Aftix^Xi - ) - ...) (4.57)

Jt=1

К1=1

Этот ряд при t = t2 по условию сходится.

Для объектов без прямых нелинейных связей и обоб­ щенно консервативных объектов выражения для коэффи­ циентов оптимальных управлений принимают вид

A i

( t ) =О,

 

 

 

 

 

 

 

 

П

 

 

 

 

 

АН(0 =

— 2

Р**»* {h, t)

(<г, t) ,

 

 

 

v,H=l

 

 

 

 

 

 

 

n

(|

n

 

) A p v ( t r) +

 

А ц к

(0 = — 2

j*{ 2

 

 

 

 

v,9,t=l

t

p=1

 

 

 

+

flpvt (f) Ap\x (t') +

flpvji (t') Ape (<')]} X

 

 

 

 

X Ulvi (t\ t)

(t\ t) Wek (t\ t) dt',

Aijkl (t) —

 

 

 

 

 

 

=

2

 

 

(t2, t) W p . j ( t 2 , t) Wtk(t2, t) X

@.58)

 

v, 9.,*, 11= 1

n

 

f,

n

 

 

x Wni ( t2, t )

 

oo A p t ( t )-(-...

2

 

J { 2

 

 

«.И1.* .4 = 1

t

P=1

 

 

 

 

...- ( -

apvp.t (t') App (<')] +

 

 

n

 

 

 

 

 

 

 

+

2

[apen(0 -4pvn(0 +•••+ ЯрмЦ (^0 ^4рец(<')]) X

 

p=i

(t\ t) ww (t’, t) W t k

(t’, t) wnl (f, t) dt’,

 

 

X W t i

Если невозмущенное состояние исходного

объекта при

ut =

0

асимптотически

устойчиво по Ляпунову, т. е.

wik (t2, t) -*■ 0

при t ->

00, то для ограниченных по моду­

лю

коэффициентов

 

 

 

 

&ijkо • • •



180

РА В Н О В ЕС Н Ы Е РА С П Р Е Д Е Л Е Н И Я В СИСТЕМАХ [ГЛ . IV

все величины A i} (t),

A ijk (t), . . . стремятся к нулю при

t ->• —

оо. Это является предпосылкой быстрой сходимости

ряда (4.57) при tx

t2 для данного случая.

Синтез систем оценивания. Пусть процесс, координаты которого х 1п должны быть оценены с максимальной воз­

можной точностью, описывается уравнениями

(4.59)

гДе fin— известные функции, | г — белые шумы. Коорди­ наты х 1а измеряются с помощью датчиков, имеющих ошибки или шумы xi ш:

Ъя = Ъп + Хт, i = l , 2 , . . . , n .

(4.60)

В задаче фильтрации Калмана непосредственно измери­ мыми считаются некоторые линейные комбинации коорди­ нат Xia в совокупности с шумами. Однако предполагается

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

%im4* fim (#l.im • • • >#nin> t) =

i = 1, 2, . . . , П. (4.61)

Требуется построить систему оценивания (фильтр) вида

Х{ф - |- f \Ф (#1 ф , • . .

, # п ф , ^1д>

• • ■ » *^ПД> 0 ж

6 ,

i =

1 , 2 , . . . ,

и,

(4.62)

выходные величины которой x t^ в определенном смысле

наименьшим образом отличались бы]от истинных значений координат хщ контролируемого процесса по истечении

определенного времени наблюдения. При линейных функ­ циях /in./iffl) /гф квадратичном критерии приближения и известном нормальном начальном распределении коорди­ нат эта задача имеет решение в виде фильтра Калмана

(5.11], [5.12].


§ 4.2) Ф П К -У РА В Н Е Н И Е И С И Н ТЕ З СИСТЕМ У П РА В Л Е Н И Я

181

Представим систему уравнений (4.59)—(4.62) в виде

■^in 4" /in (^-Xnj • • • ) З-шп О *=

ijfl,

 

 

X-im 4* fim (®1ш> • • • »%пш, О =

?4пц

 

 

•^Хф4“ /гф (*1ф, ••* ) Хпф, ^ln 4" *^1ш» •••1 £пп +

,

(4.63)

 

4" Хпт, 0 — О,

 

 

г= 1, 2, . .. , п ,

изапишем для этой системы ФПК-уравнение в виде (1.13):

ПП

3 In р

 

ZJ

d in р

 

 

 

3 In р

 

 

 

dt

 

in -^ Г ------- Z) Jim

toim

 

 

 

 

i=l

. tn

 

t=i

 

 

 

 

 

S

4

d In p

1

n

сП /

3* In p

,

3 In p

d In p \

 

 

 

 

Г

*

 

ij \dx.ndx,n

+

~ Э 4 ~

3*,п )

 

i » l

 

i,i—1

 

'

гП 7П

 

гП

'

 

 

 

1 v

cm / 9* In р

, 3 In р d In р \ _

 

 

~ ~ Г £

1

 

 

 

+

 

^ ш Г

 

 

 

 

-

I

 

$

Г

+ ^

+ $<*«“ ») '

где Si} — спектральные плотности шумов | гп, 61™— спек­

тральные плотности шумов £г ш, причем шумы | jn счита­ ются некоррелированными по отношению к шумам £<ш. Обозначим векторы

Ха — (Хщ, . . . , хпв),

j

 

хш = (^1Ш) ■• • I я-пш)»

I

(4.65)

Хф == ( % ф ) ■ ■ • 1 Хпф)'

}

 

Если бы при заданном начальном распределении

р ( х в , Хф, Хт , t-у) = P i (Хв , Хф, х ш)

удалось найти такие функции / [ф, при которых распреде­

ление хп, Хф в момент времени t2 >

ty

00

 

Рг (*m Хф, <2) = \ р (хп, Хф,

хш, t2) dXja

(здесь интегрирование ведется по пространству ошибок датчиков) приближалось бы в определенном смысле


182 РА В Н О В ЕС Н Ы Е РА С П Р Е Д Е Л Е Н И Я В СИСТЕМАХ [ГЛ . IV

наилучшим образом к 6-распределению б (жп —£ф), то зада­ ча синтеза оптимальной системы оценивания была бы ре­ шена. Однако в такой общей постановке решение задачи встречает существенные трудности и в настоящее время не известно.

Представляя заданные функции / гп, /*ш и неизвестные функции / гф указанных аргументов в форме полиномов или рядов, можно решать уравнение (4.46) описанным выше методом рядов. При этом получается бесконечная система уравнений типа (2.11), в которой коэффициенты athi • ••», относящиеся к функциям / гп, / гш, заданы, а от­

носящиеся к / гф — свободны, с учетом, однако, связей, накладываемых составом аргументов этой функции. Если ограничиться конечным числом членов разложения текущей логарифмической плотности, то указанная система обыкновенных дифференциальных уравнений становится замкнутой. Далее необходимо в пространстве коэффици­ ентов разложения In р назначить критерий или функцио­

нал точности оценивания. Если начальное распределение задано, т. е. известны начальные значения коэффициентов разложения In р, то после всего этого задача сводится

к выбору варьируемых параметров ацц..., из условия экстремизации заданного функционала. Эта задача в принципе разрешимая, но трудная. Одна из трудностей заключается в том, что желаемое распределение является обычно концентрированным относительно ягф — а:<п (ма­ лые ошибки оценивания). Это может вызывать плохую сходимость разложения In р в ряд, необходимость учета

большого числа членов ряда и чрезмерную громоздкость системы уравнений типа (2.14). В связи с этим возникает мысль о решении той же задачи методом моментов, т. е. использовании линейных уравнений моментов типа (2.77). При концентрированном распределении разложение ха­ рактеристической функции (2.76) быстро сходится и можно ограничиться учетом моментов невысокого порядка. Одна­ ко начальное распределение может быть не концентри­ рованным по всем фазовым координатам. В этом случае возникают те же трудности, связанные с громоздкостью уравнений моментов.

Можно поставить задачу синтеза оптимальной систе­ мы оценивания как задачу аналитического конструирова­ ния. В этом случае синтезируемые функции /*ф без огра­

§ 4.2] Ф П К -У РА В Н Е Н И Е И С И Н Т Е З СИСТЕМ У П РА В Л Е Н И Я

Ig b

ничения общности представляем в векторном изображении в форме

/ф = / п (Хф, t) — Нф (Хф, Ха + Хт, t)

( 4 . 6 6 )

и рассматриваем Иф как синтезируемые управления. Такая задача решается в работе [4.8]. Трудность здесь заключается прежде всего в том, что имеет место неполная степень непосредственной наблюдаемости: в управлениях могут использоваться лишь выходные сигналы синтези­ руемого фильтра и датчиков. Для неполной степени наблю­ даемости в общем случае разработан лишь приближенный метод аналитического конструирования [4.8].

Кроме этих сложных путей синтеза существуют более простые, но, конечно, менее совершенные пути. Так, можно задаться равновесным распределением, найти условия его существования и связь параметров равновес­ ного распределения с параметрами фильтра. Далее пара­ метры фильтра выбираются так, чтобы обеспечить желае­ мый характер и параметры равновесного распределения. Рассмотрим этот путь, задавшись нормальным централь­ ным относительно Да;, = а^ф — хщ, Хщ, х ш равновесным

распределением

П

In р = Aq-|—2~ | [^-i/cXjaXfca "f" ^-ihXimXkm “Ь i, к=1

- f - А ы (х^ф Xia) ( х кф X ka)

2 A i k XiaXjan - f -

+ 2Afi^Xia (хкф— a:(in) + 2Аи^хцп (хщ ^(cn)] (4.67)

и представлением функций f ia, fm

в форме полиномов или

рядов

 

пп

/in = 2

“Ь

2

^ijkXinXkn “Ь

3—1

 

3, к=1

п

 

 

 

 

 

 

 

 

+ 2

®ijklXjnXknXin +

п

 

п

з, к, 1=1

 

 

( 4 . 6 8 )

/»ш — 2^I

Q'ijXjai -f-

21

&ijkXjinXкш ~f"

j= l

;

3,4—1

 

 

 

+

2

Q-ijkt XjmXкшХ1ш H * • • ,

з, к, 1=1