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

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

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

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

Добавлен: 11.04.2024

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

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

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

показано, что эти операции вновь приводят к

гауссов-

ским

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

 

 

 

 

 

 

 

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

что х — гауссовский

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

тор с математическим ожиданием х

и

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

матрицей

Рхх,

а А — матрица размера

тХп.

Тогда

слу­

чайный m-вектор

у=Ах

 

является гауссовским

вектором

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

ожиданием у=Ах

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

матрицей РѴУ=АРХХА'.

 

Доказательство

этого утвержде­

ния

легко получить

 

с

помощью

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

функций.

Согласно уравнению

(3-32)

характеристиче­

ская

функция

X имеет вид:

 

 

 

 

 

 

 

 

 

ІР, (s) =?Е (e'x's)

 

=

ехр ( jx's

-

s'Ps^.

(3-44)

Характеристическая

функция

у по определению

 

где г — m-вектор.

Ь(г)

 

=

Е(е'у'г),

 

 

 

(3-45)

 

 

 

 

 

 

 

 

 

 

 

Так как у —Ах,

то

 

 

 

 

 

 

 

 

 

 

 

(г) =

Е [ехр / {Ах)'

г) =

Е [ехр jx'

{А'г)} =

<?х

(А'г)

=

 

 

=

ехр

jx'A'r

 

г'АРА'г

 

 

 

 

 

 

=

ехр

j (Ах)'

г

\г-г' (АРА') г

 

 

 

и исходное утверждение доказано.

 

 

 

 

 

Теперь допустим, что х и у—гауссовские

случайные

n-вектор и m-вектор соответственно

с совместной харак­

теристической функцией (3-34). Пусть

А и В — действи­

тельные матрицы

размера

рХп

и рХт.

Покажем,

что

случайный

р-вектор

z=Ax+By

является

гауссовским

вектором с математическим ожиданием z=Ax

+ By и

корреляционной матрицей

 

 

 

 

 

 

 

 

Рzz=АРХХА

-\-АРХѴВ'-\-ВРѴХА

-\-ВРл1уВ'.

 

Доказательство в основном совпадает с доказатель­ ством, приведенным выше. Пусть % — случайный (п + т)- мерный вектор вида

% =

Пусть, кроме того, q означает

(п + т)-мерный

вектор,

а С

обозначает матрицу \\АВ\\

размера

рХ(п+т).

Тогда

z = C\.

 

 

I l l


Согласно

уравнению

(3-32)

 

 

 

 

 

 

 

: ехр

 

1

 

(3-46)

 

 

 

 

 

 

где

 

 

р

 

 

 

 

 

 

Р =

РхУ

 

 

 

 

 

1 XX

 

 

 

 

 

 

РУх

Руу

 

 

 

Характеристическая функция z имеет вид

 

 

 

<?z(s)=E(eiz's),

 

 

 

(3-47)

где s — р-вектор.

 

 

 

 

 

 

Так как

2 = С £ ,

из уравнений

(3-46)

и (3-47)

сле­

дует, что

 

 

 

 

 

 

 

Ь (s) = Е [ехр j (Ck)' s] =

E [ехр j?

(Сs)] =

<р. (Сs)

=

=

ехр ß'(C's)

lT{C's)'P(C's)]

=

 

 

— ехр

j(CÏ)'

s--^-s'

 

{CPC')s^.

 

Это означает, что z является гауссовским случайным вектором с математическим ожиданием

 

z=Œ^\\A

В\

=

Ах +

Ву,

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

матрицей

 

 

 

PZZ

= CPC^=\\А

В

Р*х

РхУ

А'

В'

 

 

 

 

 

Рух

РѴѴ

=

АРХХА'

+

АРхуВ'

+

ВРухА>

+

ВРууВ',

что и требовалось доказать.

 

 

і

Вывод

формулы

для

плотности

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

С помощью обратного преобразования (3-21) здесь будет показано, что плотность распределения (3-33) со­ ответствует гауссовской характеристической функции (3-32).

Пусть X — гауссовский случайный n-вектор с харак­ теристической функцией (3-32). Предположим, что ма­ трица Р положительно определена. Тогда матрице Р со­ ответствует единственная положительно определенная 112


матрица А = Рм, т. е. такая матрица, что А2

= Р (см.

[Л. 3-6]).

 

Теперь рассмотрим случайный я-вектор

 

у = Р-Ч2(х-х),

(3-48)

где матрица Р _ І / 2 является матрицей, обратной Р 1 / 2 . Слу­ чайный вектор у, очевидно, гауссовский, поскольку он получен из X с помощью линейного преобразования. Ясно, что

Е (У) =

0; Руу = Е (уу') = P-Wppr-W

=

; .

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

 

 

 

?ѵ (s) =

exp ( Y s's\ = exp I - - i

 

-

s2 J.

 

4

1=1

'

 

 

 

 

Подставляя этот результат в уравнение (3-21), полу­ чаем:

 

 

 

 

 

00

 

00

 

 

 

 

 

 

 

— 0 0 — 0 0

 

 

 

X

 

 

 

 

 

 

 

 

Так как

 

»=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j "

exp {—'msi

 

Ys2j

dsi

=

 

 

—00

 

 

 

 

 

 

 

 

 

00

 

 

exp ( d s *

 

 

=

J exp (r-iyiSi)

 

=

 

 

—00

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

-

-yiß

 

 

 

 

 

 

 

 

I

cos (t/t-s,-) exp ( —

Y

d s

=

V2ve

 

i = V%

 

то отсюда

следует, что

 

 

 

 

 

 

 

 

 

1

 

exp

 

 

 

 

 

 

î(y)

= Ѵ{2к)п

--ТУ'У)-

 

 

Из уравнения (3-48)

 

 

 

 

 

 

у'у = (je

-

je)' р-'12р-1'2

(Х-х)

=

(х-

х)' Р'1 (х— X)

8—85

113


й искомая плотность распределения имеет вид

f{x) =

V (2п)" ехр

• (х %)' Р~1

х) =

 

 

-^(х-х-уР-Чх-х)

Ѵ ( 2 я ) " | Р |

где J = \ P \ -1/2

якобиан преобразования (3-48).

Независимость и

корреляция

В § 3-4 показано,

что два независимых случайных

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

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

важное

и полезное свойство гауссовского распределения легко

доказать.

 

 

 

Пусть X и у соответственно п- и m-мерные

гауссов-

ские

некоррелированные

случайные

векторы.

Тогда

ху

и уравнение (3-34)

принимает

вид:

 

(s,

r) =

ехр//(*'s

+

y ' r ) - 4 -

 

О

Р,

 

 

 

 

 

 

 

 

 

 

 

ѵѵ

 

 

ехр / {x's

+

у'г)

-

4 - (s'Pxxs +

r'Pyyr)

 

=

ехр ( jx's

 

к- s'Pxxs)

ехр ( jy'r

2

r'Pyyr ) =

 

 

 

 

 

=

 

?x{s)fy(r),

 

 

 

что и требовалось доказать. Это свойство

можно

также

получить

для

совместной

плотности

 

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

f(x, у)

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

и Рѵѵ

поло­

жительно

определены.

 

 

 

 

 

 

Сформулируем два основных результата, полученных

выше,

в виде

теоремы.

 

 

 

 

 

 

Теорема 3-2.

 

 

 

 

 

 

 

 

1)

Независимые

 

гауссовские

случайные

векторы не-

коррелированы

и

наоборот.

 

 

 

 

2)

Линейное

преобразование

и линейная

комбинация

гауссовских

случайных

векторов

являются

гауссовскими

случайными

векторами.

 

 

 

 

 

 

1U

 

 

 

 

 

 

 

 

 

 

 


Свойства

гауссовского условного

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

ожидания

 

 

 

 

 

В заключение параграфа получим некоторые свой­

ства условного

математического ожидания

гауссовских

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

фундаменталь­

ное значение для задач оценки.

 

 

 

1. Е(х\у)—случайный

вектор, компоненты

которого

являются линейной комбинацией компонент у.

 

Это

свойство следует

непосредственно

из

уравнения

(3-41)

и п. 2 теоремы 3-2.

 

 

 

 

2. Разность

x—Е(х\у)

не зависит от случайного век­

тора, полученного с помощью любого линейного преоб­

разования

вектора

у.

 

 

 

 

что х — п-вектор,

Чтобы

показать

это, предположим,

у — г-вектор,

а

А — произвольная

матрица

размера

mXr.

Случайные

векторы £ = х—Е(х\у)

и t, = Ay

явля­

ются

гауссовскими. Очевидно,

 

 

 

 

 

 

 

 

1 = Е®

= х-Е[Ех(х\у)]

=

0;

 

 

 

 

 

 

 

Z=

E®

=

Ay.

 

 

 

 

Отсюда

ясно, что

 

 

 

 

 

 

 

 

 

 

Р К = Я [ ( 5 - Г ) ( С - Q ' ] =

 

 

 

— Е{[х — Е (х\у)} (у -

у)' А'} =

Е {[(x - х)

-

 

-

Р*уР^

ІУ -

П

(У-Т>'}А'

= (р-

Л , ) Л ' = О,

и искомый

 

результат

сразу

следует

из п. 1

теоремы

3-2.

 

у

и z

 

 

 

 

где z — m-мерный,

ay —

3.

Если

независимы,

r-мерный случайные

векторы, то

 

 

 

 

 

 

Е(х\у,

z)=E(x\y)+E(x\z)—x.

 

 

(3-49)

Для доказательства свойства определим | как (г + пг)- мерный вектор вида

Тогда из уравнения (3-41) имеем:

Е (х\у, г) = Е (х\І) = x - f РАР-' (* _ 6).

U5