Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Медич Д. Статистически оптимальные линейные оценки и управление

.pdf
Скачиваний:
65
Добавлен:
25.10.2023
Размер:
13.81 Mб
Скачать

Уравнение (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)][х(\)-х(\))>}

= Р.

Далее, предполагается,

что {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

Соседние файлы в папке книги из ГПНТБ