Файл: Медич Д. Статистически оптимальные линейные оценки и управление.pdf

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

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

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

Добавлен: 11.04.2024

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

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

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

является гауссовским и для его описания требуется мат­ ричная корреляционная функция P(k, / ) , а не P{k). Чтобы обосновать приведенное здесь описание, покажем,

что

знания матрицы P{k)

достаточно для

определения

 

По аналогии с уравнением (4-30) имеем:

 

 

 

 

k

 

 

x{k) =

<è{k, / ) * ( / ) +

Ц

Ф(*. /)Г(г,

i—\)w{i~\)

для

k >• /.

Отсюда ясно,

что

 

 

+

S Ф(*. і)Г(і,

t - l ) £ t œ ( i - l ) 3 c ' ( j ) l .

 

 

/ = / + і

 

 

 

Заменяя в уравнении (4-36) k на /, транспонируя обе части полученного уравнения, умножая их слева на wii1) и усредняя, получаем:

Е [w іі -1)х'

(/)] = Е [w іі

-1)х'

(0)] Ф' (/, 0)

+

+

£

E[wii-

\)w'il~

\)}Y'il,

/ _ 1 ) Ф ' ( / ,

/).

Как и ранее, E[wii1)Ж'(0)]

= 0

для

всех і = 1 , 2, ...

Кроме

того,

индекс і здесь

 

принимает

значения

1 + 2, ...,

k.

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

 

 

 

 

 

 

 

 

 

E\0ii—\)w'il—1)]

 

= 0

 

 

для рассматриваемых значений г и /. Поэтому £ [ й ? ( » - 1 ) Г ( / ) ] = 0.

Это значит, что выражение для Pik, j) сводится к

Pik, j)=Oik, })РЦ),

откуда ясно видно, что для описания случайной после' довательности достаточно знать ее корреляционную матрицу.

Короче

говоря, гауссовость случайной

последователь-

ности

(x(k),

kŒl}

требует для ее описания знания пары

[xik),

Pik,

f); k,

J Œ I } . Марковость позволяет свести это

требование к определению {x(k)t Pik),

ksi}.

148


Уравнение (4-37) можно интерпретировать как опи­ сание изменения во времени степени неопределенности системы, поскольку корреляционная матрица характе­ ризует меру отклонения распределения случайного со­

стояния системы

x(k) от его математического ожидания.

В силу этой

интерпретации уравнения (4-37), второе

из приведенных описаний гауссовской марковской после­

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

будет

играть

преобладающую

роль во

всей дальнейшей работе.

 

 

 

 

Модель легко обобщить, включив в нее управляющее

воздействие:

 

 

 

 

 

 

x(k

+ l)=<D(k+l,

k)x{k)+V{k+\,

k)w(k)

+

 

 

 

+ W(k+l,

k)u(k),

 

 

где слагаемое

^(k+l,

k)u(k)

описано ранее.

В

пред­

положении, что

управляющая

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

{u(k);

k = 0, 1 . . . } известна

или может быть сформирована по

желанию,

случайный

процесс

{x(k)\

& = 0, 1,

. . . } вновь

является

гауссовской

марковской

последовательностью.

Это следует из того, что присутствие известного управ­ ляющего воздействия просто добавляет детерминирован­ ную составляющую ко входному сигналу системы и,

следовательно, воздействует

только на

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

ожидание процесса

{x(k),

/> = 0, 1,

. . . } .

Тогда

 

x(k+l)

= 0(k

+ \,

k)x{k)+T(k

+ \,

k)w(à)

+

 

 

+ W(k+l,

k)u(k)

 

 

для k = 0, 1

... Поэтому

 

 

 

 

 

x(k+\)—x(k+l)=0(k+l,

 

 

k)[x(k)—x(k)]

+

 

+ Г ( Л + 1 ,

 

k){w{k)—w{k)l

 

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

При определении условного математического ожида­ ния, используя знание управляющего воздействия, по­ лучаем:

E[x(k+\)\x{k),

и( £ )]= . ф(А+1,

k)x(k) +

+ V(k + \,

k)ïo(k)+4f{k+\t

k)u(k),

где k — 0, 1 . . . Условная корреляционная матрица имеет тот же вид, что и раньше.

149



Уравнение (4-28) представляет собой модель линей­ ной дискретной системы с возмущающей функцией в ви­ де гауссовской белой последовательности, независимой от начального состояния системы. Эти допущения заслу­ живают некоторого комментария.

Во-первых, предполагается, что модель линейна. Это оправдывается тем, что использование линейных моде­ лей в инженерных исследованиях оказалось достаточно плодотворным и методы анализа линейных систем хо­ рошо изучены, в то время как для нелинейных систем таких методов в общем случае нет.

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

рующие

векторы

возмущения и начального состояния,

не связаны

и не

взаимодействуют.

Например,

можно

утверждать,

что случайные

внутренние

возмущения инер-

циальной

навигационной

системы

сверхзвукового

пассажирского самолета

(вызванные

неточностью гиро­

скопов,

акселерометров

и связанной сними электроники)

не зависят

от начальных

координат

(широты,

долготы

и высоты), подаваемых на вход навигационной системы непосредственно перед взлетом. Однако на практике встречаются ситуации, в которых это предположение не­ допустимо.

В качестве примера рассмотрим связь между случайными порывами ветра, воздействующими на кос­ мический корабль во время запуска и подъема в атмо­ сфере, и начальными условиями, т. е. координатами пусковой установки. Очевидно, роза ветров зависит от расположения пусковой установки. Тем не менее здесь

будут исследованы только случаи первого

рода.

В-третьих, предположение о том, что компоненты

вектора

возмущения имеют

в каждый рассматриваемый

момент

времени совместное

гауссовское

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

основанона центральной предельной теореме (теоре­ ма 3-1). Здесь неявно допускается, что макроскопический случайный процесс возмущения системы является сум­ мой большого числа, независимо действующих микроско­ пических случайных процессов. Это же справедливо и для модели начального состояния х(0). Можно также обосновать такое допущение тем, что, оно оказалось

150


полезной идеализацией при Моделировании многих фи­ зических явлений.

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

Измерение

Обратимся к построению модели измерения. Пред­ полагается, что она имеет вид

z{k+\)

= H{k+\)x(k+\)

+v(k+l)

(4-38)

для kŒl, где

z — m-вектор

измерения;

H — матрица

размера тхп,

связывающая вектор состояния и вектор

измерения; ѵ — m-вектор ошибки измерения. Предполагается, что случайный процесс ошибки из­

мерения является m-мерной гауссовской белой после­

довательностью {v(k + \), k<=I} с математическим

ожи­

данием >•

 

 

 

 

и матричной корреляционной

функцией

 

 

 

' ( Ер

(/ + 1)' —Щ + 1 ) I M * + I )—v\k + 1 ) ] ' } = / ? (А + 1 ) ôj k ,

известными для всех / и k из /, где ôjh — символ

Кроне-

кера;

R(k+l)—неотрицательно

определенная для

всех

k^I

матрица размера mXtn.

 

процесс {v(k +

Также предполагается, что случайный

+ 1), k.1} не зависит от х(0)

для всех

& Œ / . Следова­

тельно,

 

 

 

 

 

£{[х(0)—x{0)][ü(k+\)—

»(ft + l)]'} = 0.

(4-40)

151