Файл: Красовский, А. А. Фазовое пространство и статистическая теория динамических систем.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