Файл: Медич Д. Статистически оптимальные линейные оценки и управление.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 на /, транспонируя обе части полученного уравнения, умножая их слева на wii—1) и усредняя, получаем:
Е [w іі -1)х' |
(/)] = Е [w іі |
-1)х' |
(0)] Ф' (/, 0) |
+ |
|||||
+ |
£ |
E[wii- |
\)w'il~ |
\)}Y'il, |
/ _ 1 ) Ф ' ( / , |
/). |
|||
Как и ранее, E[wii—1)Ж'(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