
книги из ГПНТБ / Медич Д. Статистически оптимальные линейные оценки и управление
.pdfУравнение (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
Не исключается возможность, что процессы |
{w(k), |
|||||||||||
k^I) |
и {v(k+\), |
kŒl} |
взаимно |
коррелированы, |
т. е.: |
|||||||
E{[w(j+l)—w{j+l)][u(k+l)—v{k+l)]'} |
|
|
= |
|
S{k+l)öjk, |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
(4-41) |
для |
/, k = 0, |
1, . . . , |
где S(k+ |
1) — взаимная корреляцион |
||||||||
ная |
матрица |
размера |
рХт. |
В |
большей |
части |
дальней |
|||||
шей работы рассматривается только случай |
S ( é + 1 ) = 0 |
|||||||||||
для |
всех k. |
Однако |
иногда |
в |
рассуждениях, |
а |
также |
|||||
в предлагаемых |
задачах |
затрагивается |
и более |
общий |
||||||||
|
|
|
|
|
|
|
|
ѵ(к+1) |
|
|
|
|
|
х(к+1): |
|
Н(к+!) |
|
|
|
|
|
|
|
||
Рис. 4-6. Модель |
дискретной линейной |
измерительной |
системы. |
случай. В любом случае модель измерений имеет струк турную схему, приведенную на рис. 4-6.
Причины |
выбора линейной структуры |
уравнения |
|||
модели измерения (4-38) |
совпадают с причинами вы |
||||
бора |
линейного уравнения |
модели |
динэхмики системы |
||
(4-28). |
|
|
|
|
|
По |
тем же |
причинам, по каким |
процесс |
возмущения |
считается гауссовский и белым, ошибка измерения так же считается гауссовской белой последовательностью. В модели измерения предполагается, что ошибки изме рения, проделанного в какой-либо момент времени, не зависят от ошибок измерения, проделанного в любой другой момент времени. После решения в этом предпо ложении задачи фильтрации в гл. 5 будет показано, как изложенную теорию можно распространить на случай коррелированных ошибок измерения.
Предположение о независимости *(0) и |
{v(k+l), |
,/ге/} можно обосновать тем, что из физических |
сообра |
жений факторы, вызывающие ошибки измерения и колебания начального состояния, обычно считают не зависимыми. Вообще говоря, измерительная система яв ляется внешней для динамической системы, так что не точности оборудования, вызывающие ошибки измерения, функционально независимы от начального состояния.
Из общих соображений также естественно ожидать, что факторы, вызывающие возмущения системы, во
152
многих |
случаях |
|
не связаны |
с |
факторами, влияющими |
|||||
на ошибки измерения. |
|
|
|
|
|
|
|
|||
В качестве |
примера |
случая, когда предположение |
||||||||
о независимости |
процесса |
{v(k+\), |
k^I) |
от |
х(0) |
и |
||||
{w{k), |
kŒl} является |
разумным, |
рассмотрим |
еще |
раз |
|||||
сверхзвуковой |
пассажирский |
самолет. |
Предположим, |
|||||||
что он |
совершает |
рейс |
между |
Лос-Анджелесом |
и Нью- |
Йорком и над Сент-Луисом проводится сеанс допплеровской навигации. Поскольку действие электронной допплеровской системы не зависит от координат само лета при взлете, можно считать, что начальное состоя ние и ошибки допплеровских измерений, проделанных за время сеанса, являются независимыми. Кроме того, поскольку электронная допплеровская система является внешней для инерциальной навигационной системы, разумно считать, что ошибки измерения и случайные возмущения в навигационной системе независимы.
В то же время можно привести примеры, в которых возмущения в системе и ошибки измерения коррелированы. Например, в наземной РЛС слежения датчик положения может быть связан с осью привода антенны радиолокатора через соответствующую зубчатую пере дачу. Тогда случайные порывы ветра, воздействующие на антенну, не только являются возмущениями системы,, но и присутствуют в выходном сигнале датчика в виде шума измерения.
Как следует из уравнения (4-38), в начальный мо мент времени измерения не проводятся, т. е. z(0) отсут
ствует. Это не является |
существенным |
ограничением. |
||
Удобно |
предположить, что случайный |
вектор х(0) и |
||
случайные |
процессы |
{w(k), k^I] |
и |
{v(k+\), |
в уравнениях модели (4-28) и (4-38) имеют нулевые математические ожидания. При этом общность не те
ряется, поскольку |
математические ожидания |
x{k+\) |
||||
и z(k + l) |
можно |
легко |
вычислить, |
используя |
соотно |
|
шения |
|
|
|
|
|
|
х ( А + 1 ) = ф ( / г + 1 , |
k)x(k)+T(k+l, |
k)w(k); |
|
|||
|
I{k |
+ |
|
\)=H{k+\)x{k+\)+v{k+\) |
|
|
для k = 0, |
1, |
где *(0), {w(k), kŒl} |
и {o(k+l), |
А е / ) |
имеют ненулевые математические ожидания. Такое до пущение будет использовано при рассмотрении задачи оценки в гл. 5 и 6 и задачи управления в гл. 9.
Основная идея здесь заключается в том, что ненуле вые математические ожидания входят в описание систе мы как детерминированные составляющие, которые мож но легко вычислить в любое время. Следовательно, их можно считать второстепенными (но не обязательно пренебрегать ими) при анализе.
Гауссовская |
марковская |
последовательность |
второго порядка |
|
|
Рассмотрим |
случайный |
процесс {x(k), 6 = 0, 1 . . . } , |
описываемый линейным векторным разностным уравне нием второго порядка
х(£ + 2) = |
Ф( £ + 2, k+l)x(k.+ |
l) |
+ |
|
|
+ Ѳ(& + 2, k)x(k}+T{k |
+ 2, |
k)w{k) |
(4-42) |
||
для k = 0, 1, . , в |
предположении, |
что |
х(0) и |
х(\) |
имеют совместное гауссовское распределение с извест
ными параметрами: |
' |
|
|
£{х(0)] = х(0); |
£ { * ( 1 ) ] = * ( 1 ) ; |
||
Е {[х (0) - x |
(0)] [x (0) —x (0) П = Яоо; |
||
Е{[х(1)^х-(\)}[х(\)-х(\)]'} |
= Рій |
||
Е {[х(0)-х(0)][х(\)-х(\))>} |
= Р0і. |
||
Далее, предполагается, |
что {w(k), |
k = 0, Г . . . } — |
гауссовская белая последовательность с известными ма
тематическим "'ожиданием |
и |
корреляционной матрицей, |
|||||
независимая от'х(0) |
и |
х(\). |
|
w—77-вектор; |
|||
В уравнении |
' |
(4-42) |
х — п-вектор;' |
||||
Ф(& + 2, ß+T ) й |
Ѳ(& + 2, |
k) — матрицы |
размера |
п х « ; |
|||
T(k + 2, k) — матрица размера |
пхр. |
|
|
||||
В силу гауссовости х(0), х(1) и {w(k), |
k = 0, 1 . . . } и |
||||||
линейности уравнения |
(4-42) |
[x(k), k = 0, |
1 . . . ] , |
очевид |
но, является гауссовский процессом. Более того, по скольку уравнение (4-42) второго порядка, условная функция распределения вероятностей состояния систе мы x в любой момент времени при условии, что известны состояния системы в любые прошедшие моменты вре мени, зависит только от состояния системы в два пре дыдущих момента времени, т. е. процесс является мар ковским второго порядка, 154
Структурная схема системы |
показана на |
рис. 4-7, |
где указаны два блока задержки, |
поскольку |
уравнение |
(4-42) имеет второй порядок. |
|
|
Вместо того чтобы исследовать гауссовскую марков скую последовательность второго порядка, удобнее пре образовать ее в гауссовскую марковскую последова-
г(к+г,к)
Ѳ(к+2,к) |
|
ф.х<7г+2) |
|
Ф(к+г,к+1) |
|
вз |
х(к+1) |
ВЗ |
х(к) |
|
Рис. 4-7. Структурная схема модели системы с гауссов ский марковским процессом состояния второго порядка.
тельность, чтобы можно было использовать результаты, полученные ранее. Вначале введем новый п-вектор с по мощью соотношения
x{k+l)=y(k). (4-43)
Используя это обозначение, перепишем уравнение (4-42) в виде
y(k+l)=<b(k |
+ 2, k+l)y(k)+e\k |
|
|
+ 2, |
k)x(k) + |
|||
|
+ T(k + 2, |
k)w(k). |
|
|
(4-44) |
|||
Далее, вводя 2я-мерный- вектор |
|
|
|
|
||||
|
|
|
|
У |
|
|
|
|
|
|
|
i |
|
|
fi 1 .»'; |
|
|
запишем уравнения (4-43) и (4-44) |
в |
виде одного соот |
||||||
ношения |
|
|
|
|
|
|
|
|
j |
|
о |
|
|
|
/ |
|
|
Ѳ(6 + |
2, k) |
|
. Ф(к+ |
|
2, |
k+ |
1) |
|
+ |
О |
k) |
w{k), |
|
|
|||
Г ( 6 + |
2, |
|
|
|||||
где / — единичная |
матрица |
размера |
пу^п. |
|
155
Обозначая
Ф*(/г + 1, * ) = |
|
О |
|
/ |
|
|
Ѳ(6 + |
2, |
k) |
Ф(к + 2, |
k+ 1) |
||
|
Г * ( * + 1 , k) = |
О |
|
|
||
|
Г(/г + |
2, £) |
|
|||
можно записать объединенное уравнение в виде |
||||||
x*(k |
+ \) = Ф * ( & + 1 , |
A ! ) J C *(Ä)+r*(Ä+l, |
А)о>(£) |
|||
для /г = 0, |
1 .. . |
|
|
|
|
(4-45) |
|
|
|
|
|
||
По определению |
вектор |
|
|
|
||
|
|
|
х(0) 1 |
х(0) |
|
|
|
|
|
У(0) |
1 |
х(!) |
|
очевидно, является гауссовским случайным 2/г-мерным
вектором с математическим |
ожиданием |
|
|
|||||
|
|
X* |
(0) = |
*(0) |
|
|
|
|
|
|
г(і) |
|
|
|
|||
|
|
|
|
|
|
|
|
|
и корреляционной матрицей |
|
|
|
|
|
|||
|
|
Р(0) = |
Р |
Р |
|
|
|
|
|
|
' 0 0 |
' 0 1 |
|
|
|
||
|
|
|
Р' |
р |
|
|
|
|
|
|
|
^ 01 |
* п |
|
|
||
Поскольку |
{w(k), |
k = 0, |
1 |
. . . } — гауссовская белая |
||||
последовательность, независимая |
от х*(0), |
а |
уравнение |
|||||
(4-45) |
имеет |
тот же |
вид, |
что и |
уравнение |
(4-28), то |
||
[x*(k), |
k = 0, 1 . . . } — гауссовская |
марковская |
последова |
|||||
тельность. |
|
|
|
|
|
|
|
|
Аналогичная процедура была использована во второй |
||||||||
части примера 4-8. |
{x(k), |
k = 0, 1 . . . } здесь был пре |
||||||
Случайный |
процесс |
|||||||
образован в |
гауссовскую |
марковскую |
последователь |
ность ценой введения нового вектора состояния, имеюще го в 2 раза больше компонент, чем исходный вектор состояния. Второй возможный подход заключается в использовании функций распределения второго поряд ка. Как уже указывалось, первый подход был избран для того, чтобы воспользоваться полученными ранее результатами.
Теперь рассмотрим несколько иную ситуацию, в ко торой также возникает гауссовская марковская после-
156
Довательность второго порядка. Предположим, ЧТО си стема описывается двумя разностными уравнениями первого порядка
x(k+l)=a>(k+l, |
|
k)x(k)+T(k+l, |
|
k)w(k); |
(4-46) |
||||||
w(k+l)=e(k+l, |
|
k)w(k)+A(k+l, |
|
k)l(k) |
(4-47) |
||||||
для k = Q, 1, |
где |
X — п-вектор; |
w— р-вектор; |
| — |
|||||||
9-вектор. |
Матрицы |
Ф(/г+1, |
k), T(k+l, |
k), |
®{k+\, |
k) |
|||||
и A(k+\, |
k) имеют размеры |
nXn, |
nXp, |
pXp |
и |
pXq, |
|||||
соответственно. Пусть |
x — вектор |
состояния, |
a w— |
воз |
|||||||
мущение |
системы. |
|
|
|
w(0) |
|
|
|
|
|
|
Предположим, |
что |
х(0) |
и |
имеют |
совместное |
||||||
гауссовское распределение с параметрами |
|
|
|
|
|||||||
|
Е[х(0)] |
= х(0); |
E[w(0)] |
= w(0); |
|
|
|
||||
|
Е {[x(0)-x |
(0)][х(0)-x |
(0)]'} = Рхх(0) |
; |
|
|
|||||
E{[w(0)-w(0)]{w(0)-w(0)]'} |
|
|
= |
Pww(0); |
|
|
|||||
Е {[x(0)-x(0)]{w(0)-~w(0)]'} |
|
|
= Pxw(0). |
|
|
||||||
Кроме |
того, предположим, |
что |
(|(&), /г = 0, |
1 |
. . . } |
является гауссовской белой последовательностью, неза
висимой от |
х(0) |
и |
w(0), |
с известными |
математическим |
|||||
ожиданием и корреляционной |
матрицей. |
|
|
|
|
|||||
Очевидно, {w(k), |
& = 0, |
1 . . . } является |
гауссовской |
|||||||
марковской |
последовательностью. Следовательно, |
урав |
||||||||
нение (4-46), имеющее |
ту |
же форму, |
что |
и |
уравнение |
|||||
(4-28), описывает |
систему |
с |
возмущением |
в |
виде |
гаус |
совской марковской, а не гауссовской белой последова тельности. Иными словами, здесь рассматривается обоб щение модели (4-28) на случай, когда возмущение си
стемы |
|
является |
частным |
случаем |
коррелированного |
||||||||
случайного |
процесса. |
|
|
|
|
|
|
{x(k), |
|||||
Проще |
всего |
убедиться |
в |
том, |
что |
процесс |
|||||||
k = 0, |
1 . . . } представляет собой |
гауссовскую |
марковскую |
||||||||||
последовательность |
второго |
порядка, |
|
рассматривая |
|||||||||
структурную |
схему |
системы |
(4-46), (4-47), изображен |
||||||||||
ную |
на рис. 4-8. |
|
Система |
со входом | ( £ ) и выходом |
|||||||||
x(k+\), |
|
очевидно, |
второго |
порядка. |
Поэтому, |
так |
как |
||||||
х(0) |
и |
w(0) |
|
имеют совместное |
гауссовское |
распределе |
|||||||
ние, а |
{£(&), |
k = 0, |
1 . . . } — гауссовская |
белая |
последо |
вательность, независимая от начальных |
условий, случай |
ный процесс [x(k), k = 0, 1 . . . } имеет те |
же свойства, что |
157
и |
Процесс, описываемый |
уравнением |
(4-42), а именно, |
|||
он |
является |
гауссовской |
марковской |
последователь |
||
ностью второго порядка. |
|
|
|
|||
|
û(k+1,k) |
|
іи(к+і) |
|
||
|
|
|
|
|
||
|
|
|
+4> |
|
|
|
|
|
|
|
ш(к) |
|
|
|
|
|
Ѳ(к+1,к) |
53 |
|
|
|
\т(м,к)\ |
+ |
|
|
xfk+1) |
|
|
|
|
|
|||
|
|
|
|
|
||
|
|
|
|
х(к) |
БЗ |
|
|
|
|
|
Ф(к+ик) |
|
|
|
Рис. 4-8. Структурная схема системы (4-46), (4-47). |
|
||||
|
Наконец, |
легко показать, что процесс {x*(k), |
& = 0, |
|||
1 . . . } , где X* |
— (и+р)-мерный вектор |
вида |
|
представляет собой гауссовскую марковскую последо вательность.
Метод сведения марковских процессов высшего по рядка к простым марковским процессам, изложенный здесь, называется методом расширения вектора состоя ния ,[Л. 4-7].
4-3. МОДЕЛЬ С Г А У С С О В С К И М МАРКОВСКИМ П Р О Ц Е С С О М СОСТОЯНИЯ
В этом параграфе описана модель непрерывной линейной системы. Метод описания модели аналогичен методу предыдущего параграфа. Вначале формулирует ся понятие гауссовского белого шума, которое затем используется для построения модели системы. Как и в § 4-2, в заключение параграфа рассматривается дина мическая система с моделью состояния в виде марков ского процесса второго порядка.
158
Гауссовский белый шум
По аналогии с определением гауссовской белой последовательности назовем /7-мерный случайный про цесс {w(t), t^to} гауссовским белым шумом, если он является чисто случайным гауссовским процессом с из вестным математическим ожиданием
E[w{t)] = w{t) |
|
|
и матричной корреляционной функцией вида |
|
|
Е {[w (Z.)-w У) ] [w (т) -W |
(%)}'} = Q (t) ô |
(t-x), |
где to — известное начальное |
время; t, x'^to, |
Q{t)—не |
прерывная неотрицательно определенная матрица раз мера рХр; ô(t,—-г) —дельта-функция Дирака. Послед няя функция вводится по аналогии с символом Кронекера в выражении для матричной корреляционной функ ции гауссовской белой последовательности.
Приведем обоснование этой формулировки, принад лежащее Калману [Л. 4-7] и основанное на рассмотре нии предельного поведения кусочно-постоянной гауссов ской белой последовательности, в которой частота пере
ключения |
становится произвольно большой. |
Пусть |
{w(k), k = 0, 1 . . . } — гауссовская белая после |
довательность с нулевым средним и матричной корре
ляционной |
функцией |
E[w(j)w'(k)]=Q(k)-6jk; |
j , & = 0, |
||||
1, |
где |
расстояние |
между последовательными |
момен |
|||
тами отсчетов во множестве индексов |
времени |
At>0. |
|||||
Пусть |
/ обозначает непрерывное время, |
to |
соответствует |
||||
k — 0, |
а А соответствует k = N. Тогда |
t{ = to+'NAt или, что |
|||||
то же самое, NAt = t\—to- |
значения. N{w(N~>(t), |
||||||
Пусть |
для некоторого заданного |
||||||
to^.t^ti} |
обозначает |
кусочно-постоянную |
гауссовскую |
белую последовательность, одна из компонент которой изображена на рис. ,4-9. Заметим, что ее значение на каждом отрезке определяется на левой границе отрезка.
Отметим также, что математическое ожидание про
цесса |
{w(k), |
k = 0, |
1 . . . } считается |
нулевым просто из |
|
соображений |
удобства. |
|
|
||
Полагая /4 постоянным с ростом N, так что |
|
||||
|
|
|
NAt = ti—/o = const, |
|
|
исследуем поведение процесса {wiN) |
(t), t0<t<t1}. |
В пре |
|||
деле |
при N-—ѵоо |
интервал At—>-0, |
т. е. частота |
пере |
ключения становится произвольно большой-
159