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

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

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

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

Добавлен: 11.04.2024

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

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

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

Но

У —7

II PxyPxz ||>

Z 2

 

( I I I I J I K - ^ ' ( 3 ^ ' 0 r =

Рѵѵ

0

- 1

р - 1

0

 

^ УУ

 

0

р

 

0

р - і

 

 

 

 

' ZZ

где использовано равенство Pj/Z = 0, поскольку случайные векторы у и z независимы. Следовательно,

 

 

 

 

 

 

 

 

«/ — У

 

 

 

 

 

 

О

Р.,—1

г —

г

 

= * + л * * 3 " 1

(у - у ) + ( z - s ) =

 

 

 

 

= £ С ф ) + £ ( . * | г ) - * .

 

 

4. Если

г/ и z не обязательно

независимы, то

 

 

 

£(х|г/,

z ) = £ ( * | & z);

 

(3-50)

 

 

Е(х\у,

z)=E(x\y)

 

+E(x\z)—x,

 

(3-51)

где

z=zE(z\y).

 

 

 

 

 

 

Уравнение (3-51) сразу следует из свойств

2 и 3, до­

казанных

выше.

 

 

 

 

 

 

 

Чтобы

доказать (3-50), убедимся, что в силу уравне­

ния

(3-41)

 

 

 

 

 

 

 

 

 

Е(х\у.

г) =

х\+\\РіуРх

РУУ

РѴг -1

У — У

 

 

 

 

 

 

Ргу

Р гг

Z 2

 

Из уравнений

(3-36) — (3-38)

получим:

 

 

 

 

 

 

Руу

Руг

 

 

 

 

 

 

 

 

Ргу

Ргг

 

 

 

 

 

 

РУУ + РУУ

РугСРгУРш

-Р^\Р«гС

 

 

 

 

 

 

 

,—1

 

 

 

 

 

 

 

~СРгуР\Гуу

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

С

(PzZ

Pzy^uu Pyz)

* '

 

 

116


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

 

 

Е {х\у,

z)

=

X + {РхуР-1

 

+

Piyp->PyzCPzyP^

 

 

 

-

 

 

-

P*zCPZyP-1)

 

-

у)

! +

(Р«С -

P*yP-'PyzC)

 

(z

-

г) =

 

 

=

* + I V " 1

+

iP^Pyz

 

-

P*z)

CPzyP^\

 

X

 

 

 

 

 

 

Х(У-У)

 

+

(P*z -

P^Pyz)C(z—z).

 

 

 

(3-52)

 

С

другой

стороны,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E(x\y,

г) =

Е(х\у)

+ Е(у\г')-х

 

=

 

 

 

 

 

 

=

* +

P*»P~l

(У-У)

 

+

Р „PZl

-

 

E{z)\.

 

Ho

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£(z) =

£ [ z - £ ( z | t , ) ] =

0;

 

 

 

 

 

 

 

 

 

P

„ =

E{(x~x)[z-E{z)}>}

 

 

=

 

 

 

 

 

 

 

 

 

 

xz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

E{(x-x)[z-z-PzyP~l

 

 

 

 

(y -

y)]'}

=

 

 

 

 

 

 

 

 

 

 

P*z

 

PxyPyy

Pyz-

 

 

 

 

 

 

 

Кроме

того,

 

 

 

 

PzyP~x

 

 

 

 

 

 

 

 

 

P^

z z

=

E{[(z-

 

z) -

{y -

y)} [(z

-

z

)

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

PryP'l

{У -

У)}'}

 

=

Pzz

-

PzyP^Pyz-

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

(x\y.

z) =

X +

 

P^P-1

(y -

 

jr) +

 

 

 

+

(/>«

-

P*yP~lpyz)

(Pzz

-

V "

'

Py r ) " 1

[Z

Я

( ф ) ]

=

 

 

 

=

* +

 

 

(У -

 

 

+

(P«

- P*yP~lPyz)

 

X

 

 

 

 

 

 

X C [ ( z - 2 ) - P n , P - , 0 / - ^ ) ] =

 

 

 

 

=

* +

 

l * V £

+ ( Л ^ Ч *

-

P*z) CPzyP-1]

(У-У)

+

v

 

 

 

 

+ ( ^ z - - P x y P - 1 P y 2 ) C ( z - 2 ) .

 

 

(3-53)

Из уравнений (3-52) и (3-53) сразу следует уравнение (3-50).

117


Заметим, что если бы все случайные векторы имели нулевые математические ожидания, то многие преобра­ зования существенно упростились. Поскольку математи­ ческое ожидание случайного вектора изображается сред­ ней точкой евклидова пространства, около которой рас­ пределены точки, изображающие значения случайного вектора, часто бывает удобно перенести начало коорди­ нат с помощью формулы % — х—x и проводить все опера­ ции со случайным вектором g, имеющим нулевое мате­ матическое ожидание. Ненулевое математическое ожида­ ние можно ввести вновь в любое удобное время. В дальнейшем обычно будет использоваться именно та­ кой подход.

З А Д А Ч И

К ГЛ. 3

 

 

 

 

 

 

 

 

 

 

 

3-1. Вывести уравнение (3-3).

 

 

 

 

 

 

 

 

3-2.

Как

выглядит функция

f(x\y)

в

примере

3-2

для

случаев

ï = + l и у

=—П

 

 

 

 

 

 

 

 

 

 

 

 

3-3.

Изобразить функцию f(xi,

xz)

из примера 3-3

для Хі = х2 и

отметить на рисунке х.

 

 

 

 

 

 

 

 

 

 

 

3-4. Доказать четыре основных свойства условного математиче­

ского ожидания,

приведенные на стр. 98.

 

 

 

 

 

 

3-5.

Вывести

уравнения

(3-36) — (3-38).

 

 

 

 

 

 

3-6.

Пусть

x — гауссовский

случайный п-вектор с

математиче­

ским ожиданием х и корреляционной матрицей

Рхх.

Показать, что

вектор у=Ах+Ь,

 

где А—матрица

размера тХп,

а Ь — постоянный

гсг-вектор, является гауссовский

случайным

вектором

и определить

его математическое ожидание и корреляционную

матрицу.

 

3-7..

С помощью гауссовской плотности распределения

(3-33),

используя

определение Е(х)

и Е[(х—х)(х—г)1],

показать, что Е(х) =

= х и Е[(х—х)

(х—х)'] =

Р.

 

 

 

 

 

 

 

 

 

 

3-8.

Пусть

x — гауссовская

случайная

величина

с нулевым мате­

матическим

ожиданием и

дисперсией

Р =

о"2 >0.

Показать,

что А-й

момент

x,

определяемый

соотношением

р-ь = £ ( * * ) ,

где k—\, 2,

имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 X 3 X

•• • X

{k — 1) о* для четных k;

 

 

 

 

 

 

 

 

 

 

 

 

для

нечетных

k.

 

Указание.

 

Использовать

характеристическую

функцию случайной

величины и вывести соотношение

 

 

 

 

 

 

 

 

~ds* «=0

где }=Ѵ—1. Этот результат в n-мерном случае называется теоремой моментов [Л. 3-4].

3-9. Если x — гауссовский случайный п-вектор, то следует ли отсюда, что ^-вектор, включающий любые k компонент вектора х,

m


при k=\,

2, ..., h—i также является гауссовским случайным

векто­

ром?

 

 

 

 

 

 

 

 

 

 

 

3-10. Показать с

помощью

функции

плотности

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

(3-35), что два некоррелированных

гауссоЕСких

случайных

вектора

независимы в предположении, что матрицы РХх

и Руѵ

положительно

определены.

 

 

 

 

 

 

 

 

 

Остальные задачи имеют своей целью

иллюстриро­

вать

возможность

применения

некоторых

результатов

главы для получения оптимальных оценок

случайного

вектора.

 

 

 

 

 

 

 

 

 

3-11.

Пусть X и у

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

скими ожиданиями п- и m-мерные соответственно.

 

 

Предположим,

что

вектор у

известен,

и

обозначим через jT=

= h (у)

 

оценку х,

где

 

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

л-мерная

вектор-функция

вектора у; т. е. значение

Л(-) известно при любом заданном

значе­

нии у.

Поскольку ошибка

оценки,

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

здесь

как х

—х,

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

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

ошибку L, определяется выражением х = Е(х]у). Указание: Использовать тот факт, что

и

положить

 

£* [g (*)] =

£« {£*[g

(*)ЫЬ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g(x)

~?'&

=

(х—х)'

{х — х).

 

 

 

 

 

Этот результат показывает, что для квадратичного критерия ка­

чества наилучшей или оптимальной

оценкой

является

условное мате­

матическое

ожидание.

 

 

 

 

 

 

 

 

 

 

 

 

3-12. Как выглядит выражение для оптимальной оценки в задаче

3-11, если X и у

независимы?

 

 

 

 

 

 

 

 

 

 

 

3-13. Пусть

X — «-мерный,

а

у — /«-мерный

случайные

векторы

с

нулевыми

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

ожиданиями

и известной

совместной

плотностью

распределения f(x,

у).

Предположим,

что

значение у

известно и нужно определить

оценку х

вида

х=А°у,

где

А0—матрица

размера пхпг. Предположим

также, что А0

требуется

отыскать среди

множества всевозможных матриц А размера пХпг

таким

образом,

чтобы минимизировать

критерий качества

 

 

 

 

 

 

 

 

 

 

L =

 

E[(x-Ay)'(x-Ay)).

 

 

 

 

119