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

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

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

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

Добавлен: 11.04.2024

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

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

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

Показать, что оптимальная оценка х =

А°у

вектора х

имеет

вид:

 

 

 

 

 

 

 

 

 

 

х =

Е(ху')[Е(уу')]-'

 

 

у.

 

 

Иными словами, показать,

что

 

 

 

 

 

 

А° = Е{ху')[Е{уу')]-і

 

=

 

РхуР~\у

 

в предположении, что матрица Руу

несингулярна.

 

 

Показать также,

что

корреляционная

матрица ошибки

оценки

X = X — X имеет вид

 

 

 

 

 

 

 

 

 

р

р

__ р

р— 1

р

 

 

*

— ~ — ' X X

'

х У г

ци

1

УХ'

 

 

 

XX

 

 

 

а

а

 

 

 

3-14. Какой физический смысл имеют результаты решения задач 3-11 и 3-13, если функция плотности распределения f(x, у) —гауссов- ская?

Г л а в а ч е т в е р т а я

ЭЛЕМЕНТЫ ТЕОРИИ СЛУЧАЙНЫХ ПРОЦЕССОВ

И ПОСТРОЕНИЕ МОДЕЛЕЙ СИСТЕМ

Целью настоящей главы является построение мо­ делей систем, для которых впоследствии будут сформу­ лированы и решены задачи оценки и управления. Основ­ ной метод главы заключается в объединении результа­ тов гл. 2 и 3. Здесь используется вероятностное описа­ ние векторов начальных условий, возмущений системы и ошибок измерения в дискретной и непрерывной моде­ лях линейных систем из гл. 2. Такое описание можно было ввести непосредственно в начале главы. Однако предварительное изложение некоторых идей теории слу­ чайных процессов способствует лучшему пониманию ко­ нечных результатов.

Теория случайных или стохастических процессов изу­

чает явления, управляемые вероятностными

законами,

в которых время или какая-либо другая

величина

является параметром. Строгое и законченное

изложение

теории выходит за рамки данной книги, поэтому здесь формально рассмотрены только аспекты теории, сущест­ венные для наших целей. Читатель, интересующийся более подробным изложением теории случайных процес­ сов примерно на том же уровне, на каком она рассма­ тривается в этой главе, найдет более полезными учебни­ ки Парзена [Л. 4-1] и Папулиса [Л. 4-2]. Более тради-

120


ционное с инженерной точки зрения изложение

теории

можно найти, например, в

книгах Лэннинга и

Бэттина

[Л. 4-3], Давенпорта и Рута

[Л. 4-4], Ли [Л. 4-5],

Джейм­

са [Л. 4-6] и др.

 

 

После изложения в § 4-1 теории случайных процес­ сов в § 4-2 и 4-3 вводятся две модели систем. Эти моде­ ли, так называемые гауссовская марковская последова­ тельность для дискретных линейных систем и гауссовский марковский процесс для непрерывных линейных систем, с успехом используются в инженерной практике моделирования систем, подверженных влиянию случай­ ных начальных условий, возмущений и ошибок измере­ ния. Причина заключается в том, что эти модели оказа­ лись достаточно простыми для аналитических и вычис­ лительных целей и в то же время настолько содержа­ тельными, что с их помощью удалось получить сущест­ венные результаты.

4-1. ЭЛЕМЕНТЫ ТЕОРИИ СЛУЧАЙНЫХ ПРОЦЕССОВ

Случайный процесс

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

Случайным

процессом

называется

семейство

случай­

ных векторов

{x(t),

.1),

зависящих

от параметра

t, все

значения

которого

принадлежат некоторому множеству

индексов

I.

 

 

 

 

 

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

121


Множество индексов / может быть абстрактным. Однако выше уже говорилось о времени как о параме­ тре случайного процесса. В дальнейшем будут исполь­ зоваться только два множества индексов времени. Пер­ вое является множеством дискретных моментов времени l={tk: k = 0, 1, . . . } , где th<tk+u и моменты времени не обязательно равно удалены друг от друга. Для упро­

щения обозначений заменим tu на k. Тогда

множество

индексов

примет

вид / = {6:& = 0, 1, . . . } .

В

этом

случае

процесс

называется случайным

процессом

с

дискретным

временем.

 

 

 

 

 

 

 

Другое используемое здесь множество индексов вре­

мени

включает

интервалы

на

оси времени,

такие как

/ = {/ : 0 ^ / ^ : Т) или I = {t : t^to).

Тогда

процесс

назы­

вается

случайным

процессом

с непрерывным

временем.

Понятие «семейство случайных векторов» в определе­ нии случайного процесса особенно важно. Оно означает, что случайный процесс состоит из множества или из ансамбля случайных векторов, определенных на множе­ стве индексов. Ансамбль может включать в себя как счетное, так и несчетное множество элементов. Для на­ ших целей представляет интерес только последняя си­ туация.

Комбинация понятий времени и ансамбля подразу­ мевает, что случайный процесс фактически является

функцией двух переменных. Это можно записать с по­

мощью

обозначения

[x(t,

ю), IŒI,

C O Œ Q } ,

где Q

выбо­

рочное

пространство

рассматриваемого

эксперимента.

При постоянном значении t x(t,

•)

является

вектор-

функцией, определенной на выборочном пространстве Q,

т. е. случайным вектором. Если

же

ш постоянно,

то

x(-, со)—вектор-функция времени,

представляющая

со­

бой одну из

возможных

реализаций

или

выборочных

функций

случайного процесса.

 

 

 

 

 

Если выборочное пространство Й дискретно, то часто

вместо

слова

процесс

используется

термин

цепь.

Поэто­

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

Теперь приведем несколько примеров случайных про­ цессов.

122


Пример 4-1. Рассмотрим скалярный случайный процесс [x(t), t^3?0}, где x(t) = Vo sin t и V0—непрерывная случайная величина, рав­ номерно распределенная на интервале [—1, lj, т. е.

 

 

 

 

 

 

 

 

 

I

0

в остальных точках.

 

 

 

 

Ансамбль

здесь

состоит

из

несчетного

множества

элементов,

a

{x(t),

t^tts],

очевидно,

является случайным

процессом

с

непре­

рывным

временем.

Ансамбль

 

,2

 

 

 

 

 

процесса

можно

представить

+ 7

 

 

 

 

 

в виде рис. 4-1.

 

 

постоянно­

 

 

 

 

 

t

 

Для

некоторого

О

 

I

р Ъ ^ - '

i

 

го значения

? =

/ i ^ 0 ,

x(tt)

бу­

 

3sz

 

-1

 

 

дет,

очевидно,

случайной

вели­

•л jI

 

 

 

чиной. Любая

синусоида явля­

+ J

 

 

t

ется

выборочной

функцией

про­

- ч

I

 

 

i

цесса.

 

 

 

 

 

 

 

про­

О

\

I

/

 

 

 

 

Если дискретизовать

-1

-

i

 

 

 

 

цесс,

рассматривая

только

мно­

 

 

 

 

жество

моментов

 

времени

 

 

 

 

 

 

 

{tk

:/„ = Ая/4,

*=0,

1, . . . } , то

 

 

 

 

 

 

 

в

результате получится

случай­

 

 

 

 

 

 

 

ный процесс с дискретным

вре­

 

 

 

 

 

 

 

менем.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4-2. В качестве вто­

 

 

 

 

 

 

 

рого

примера

рассмотрим

ска­

 

 

 

 

 

 

 

лярный

случайный

процесс

с

 

 

 

 

 

 

 

дискретным

временем

(к),

 

 

 

 

 

 

 

k=0,

1,

. . . } , где для

каждого

 

 

 

 

 

 

 

k

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

что

у(к) —

Рис. 4-1. Представление случайно­

гауссовская

случайная

величи­

на с нулевым средним. Возмож­

го

процесса

в

виде ансамбля.

 

ную

зависимость

дисперсии

у (k)

от

 

 

 

 

 

 

k учтем, записывая плотность рас­

 

 

 

 

 

 

пределения

y{k)

в

виде

 

 

 

 

 

 

 

 

 

1

Î(У,

к):

(k)

Ѵ2л

 

 

 

где

a2(k)

=Е[у*(к)].

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

что

значение y(k)

для заданного k

статистически не зависит от его зна­ чения в любой другой точке множе­ ства индексов.

Две типичные выборочные функ­ ции процесса показаны на рис. 4-2.

Этот процесс называется чисто случайной гауссовской последова­ тельностью. Прилагательное гауссовский используется с целью подчерк­ нуть тот факт, что при известном k распределение значений y{k) по ан­ самблю гауссовское, а определение

?

3 "

о

Рис. 4-2. Выборочные функ­ ции дискретного случайного процесса.

123


«чисто случайный» указывает на статистическую независимость во времени элементов во всех выборочных функциях. Термин «последо­ вательность» возникает в силу того, что любая выборочная функция представляет собой просто последовательность случайных перемен­ ных, определенных в дискретные моменты времени.

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

Описание случайного процесса

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

Рассмотрим случайный процесс с дискретным време­ нем {x(k), kŒ\}, где / = {&:&=1, ..., N}, a N— некото­ рое постоянное положительное целое число. В этом слу­ чае естественным описанием процесса будет совместная

функция

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

вероятностей

N

случайных

/г-векторов

х(\),

 

 

x(N),

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

для любых

п-векторов

хх, ...,

 

xN

 

 

 

 

 

 

 

 

F(x\

 

х»)=Р{х(\)<^х\

 

x(N)^x»l

(4-1)

 

Эквивалентное описание можно ввести с помощью

совместной

плотности

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

вероятностей

 

 

f(xl

 

xN)

=

• - j

г-^-я

 

ж '

(4 -2)

 

,

ѵ

 

;

 

дх\...дх1п...дх^

...dxNn

ѵ

'

где

х'{і=\,

 

п\

/ = 1

 

N) является

і-й компонен­

той

n-вектора xj

в

предположении,

что

указанные

nN

частных производных существуют, или через совместную характеристическую функцию

 

 

 

 

 

N

x'

(t)

si

(4-3)

 

9(s\ ....

s")--

 

ехр / £

для любых и-векторов s'(t =l , ...,

N),

где

/=]/"—1-

 

Если множество индексов является дискретным, но

бесконечным, т. е. I = {k:k=l,

2 . . . } ,

или же непрерыв­

ным, ограниченным

или

неограниченным,

т. е. / = {^:0^

ssZt^zT}

ИЛИ I={t

: t^sO}, то

подобное

 

описание

не бу­

дет столь естественным. Однако с практической

точки

зрения

представляется

разумным,

что

случайный

про-

124