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

книги / Основы теории оценивания с приложениями к задачам обработки навигационной информации. Ч. 1 Введение в теорию оценивания

.pdf
Скачиваний:
35
Добавлен:
19.11.2023
Размер:
26.49 Mб
Скачать

Используя

(1.3.45),

при

т =30°,

а? = 4 ,

02 = 1 >

К = \

найдем

значение дисперсии

Dp{z),

которое

будет

равно

2,62. С учетом вышесказанного можем сделать вывод о том, что

ошибка в направлении т =30°

представляет собой гауссовскую

ф.п.р.в. /р (р) - N (р;0,2,62)

(рис. 1.3.6).

Рис. 1.3.6. График ф.п.р.в. для р

при т =30°,D P (T) = 2,62

1.3.5. Ормогонализация случайных величин. Связь матрицы ковариаций и среднеквадратического эллипса

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

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

 

ТРХТТ = А = [Xj }

j = Т

л ,

(1.3.46)

где ТГТ =Е; 7,т =[/1

tj

/„], a

h j j j ,

у = 1 .л

- собст­

венные числа и собственные векторы матрицы Рх

 

 

P ' t j ^ j t j ,

 

 

(1.3.47)

причем tjîj=ôij.

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

Очевидно что, проведя ортогонализацию гауссовского вектора, имеющего ф.п.р.в. в виде

/х М = N ( г , х , Р х ) = ------—т -^----- — — ехр{- 0.5(х -

х)т(Рх )

1 (дг - .т)}>

 

(27r)"/2(det Рх)и2

 

 

 

для вектора х получаем плотность

 

 

 

 

/х (* ) =

-ехр< --о уZ _J

(*/

- X j Y

(1.3.48)

 

(2к)п/2л]Х1. л п

7=1

 

 

 

 

 

 

 

 

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

Рис. 1.3.7. Эллипс ошибок для двухмерного гауссовского вектора с зависимыми компонентами

Ясно, что в системе координат Â^Oxj, выбранной так, как это

показано на рис. 1.3.7, этому эллипсу будет соответствовать диа­ гональная матрица ковариаций вида (1.3.46). Переход от представ­ ления вектора в системе координат х2Ох\ к его представлению в системе координат х2Ох\, повернутой относительно оси Ох2 про­

тив часовой стрелки на угол т*, осуществляется с помощью мат­ рицы Г , определяемой как [38, с. 58]

COS Г

*

* "

cos(90- г)

sin(90 - г)

sin т

cos т

4= Sillг$

- sin(90- т) cos(90-r)

- cos г

sin г

-sin г

cos г

Таким

 

образом,

решение

задачи ортогонализации сводится

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

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

показать,

что

собственные числа матрицы ковариаций (1.3.38)

-)

А. 2

О

А-1 - а" и

=b~, Aj >Х2, представляющие собой квадраты от

малой и большой полуосей эллипса, и угол, определяющий ориен­

тацию этих осей,

задаются следующими соотношения [23, с. 93]:

 

CTf + а \

± i j ( < J f - e 5 ) 2 + 4 K 2

(1.3.50)

 

tg2t =

2К

 

 

 

 

2

2 5

 

 

 

 

а 2 ~

с т 1

 

 

 

t,

 

02 ~ и\ >

К >0;

 

т = — arctg

 

<у2 -с!2 >0,

К <0;

(1.3.51)

т + я,

сг? —а 1

я

а2 - а2 <0.

 

 

,+ 2-

 

 

Эти соотношения позволяют рассчитать параметры средне­ квадратического эллипса для заданной матрицы ковариаций.

Полагая Рх

b2

О

и используя (1.3.21), (1.3.23), (1.3.49),

О

а2

 

 

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

 

a

9

cos

9

x +b

9

sin

9

т

9

 

9

)cosTsinx

 

р х = р х J1

 

 

 

 

(а~-Ъ

 

.(1.3.52)

(a

9

 

 

9

 

 

 

 

9

9

x +b

9

cos

9

т

 

 

- 6 “)cosTsinT

а~ sin

 

 

 

 

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

не зависит от их ориентации. В частности,

+ а 2 =

~ . Этот

факт есть следствие того, что при ортогональных преобразованиях след матрицы не меняется. Таким образом, введение радиальной ошибки (1.2.24) в качестве количественной меры неопределенно­ сти местоположения представляется вполне оправданным.

Задачи к разделу

Задача 1.3.1. Пусть для гауссовского вектора

z = (xT,v T)T,

включающего два подвектора х и v

размерности

п и т , задана

ф.п.р.в.

 

 

 

 

(VX1

X

р х в

 

/« ( * ) = *

э

>

 

 

V

V

_ в Т

pV. J

 

VL J

 

 

 

Получите выражение для совместной ф.п.р.в. для составного вектора z = (х т ,у т )т , в котором у = s(x) + v , где s(x) - известная

вобщем случае нелинейная функция.

Ре ш е н и е . В данном случае в качестве z = g(z) имеем

"V ---1N1to1

X

X

Z1

.У.

s(x) +V

_J(Z!> +Z2_

таким образом, z = h(z) определяется как

Z = ’ z l"

X

 

 

X

 

 

 

X i___

_z 2 _ V

I

1

z i

_ z2

Принимая во внимание тот факт, что

О

 

 

 

Е

 

det

= det

ds(?i )

Р = 1

 

сГХ

,

 

dz.

JJ

 

и используя (1.3.6), можем записать

 

 

 

/

 

-Y

~px в

 

 

.V

 

TV

h (S') = Л,У (*, у) = Л ,г (й(х, V)) = N

V

_y-5(x)_’ V

B T

 

 

 

 

 

Задача 1.3.2. Используя результаты задачи 1.3.1, конкретизи­

руйте вид совместной плотности для вектора z = (xT,y T)T, когда

функция s(x) линейная, т.е. когда 5(х) = Н х . Уточните вид плот­ ности, когда х и V центрированные и независимые между собой.

Р е ш е п и е. При s(x) =Hx ф.п.р.в. плотность / х у(х, у) при­

нимает вид

 

 

 

 

X “1 X

_BT

в 11

Л , у (*. у) =fx.y (К*, v)) = N J_v-tfxj

1

V

H ,

 

_

5

 

Поскольку х и v центрированные и независимые между собой, можно записать

f г X

n "o'

/ х,у (х,у) = У

_oj

! j?

----1

px

о

о__1 ----1

cexpi - - Vv• (PxT lx + (y - tfv)V T ’ (>

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

Преобразуем выражение, определяющее показатель экспоненты

хт (Рх)~] х +(у - Hx)T(Pv)~l(у - Нх) =

= х г(Рх)~' X + y T(Pvy ]У - УГ(Р 'Т 1Нх - X1Я 1 (У'Т 1у + хтЯ т(/>'■ Г 1Нх

X

= [хту т]

~(pvy ]H

! OP1’)-1

Используя формулы обращения блочных матриц, получаем

*)-1 +ят(РТ'# j tfT(Pv)~r-1

\jP L l PXH T

(P V)~XH

I

(P v)_1

_H P X i Pv +HPxH r _

Отсюда следует, что

1

 

 

 

f r

 

_ P _ x_ _

j t____Р _ Х Н ____ Т

Л,У(x,ÿ) =N

3"°1

 

X

 

 

 

I b J

.o j

Н ТРХ ! P v +н р хн т1 '

Тот факт, что при s(x) =Hx плотность для составного вектора будет гауссовской, является также следствием сохранения гаус­ совского вида плотности при линейных преобразованиях гауссов­ ских векторов. Таким образом, для конкретизации гауссовской плотности достаточно найти выражение для математического ожидания и матрицы ковариаций составного вектора.

Задача

1.3.3.

Убедитесь

в том,

что для ф.п.р.в. вектора

у = х + V справедливо соотношение

 

 

 

f y ( y ) = \

fs',x(v , y

- v)d v >

которое при независимых х и v принимает вид

 

 

/ у О ) = j

f y ( v ) f x ( y - v)dv

П о я с н е н и е . При получении этих соотношений используйте

векторы z =

y

X + V

V

 

=

И Z =

 

 

V

X

X

 

Задача 1.3.4. Пусть с.в. y = x + v сформирована в результате суммирования двух независимых между собой с.в. х и v , распре­ деленных по равномерному закону в интервале [0, b].

Покажите, что ф.п.р.в. для с.в. у является треугольной и опре­ делятся в виде

Ату, при у G[о,б],

/ у ОО = <

(2b - у), при у е [b,2b] ,

(1)

О,

при у < 0,у > 2b .

 

Решение. Принимая во внимание тот факт, что

0,xg[0 ,b],

 

и используя (1.3.17), можем записать

 

/у (У)=J (*)/v~ x)dx = т \ / Л у~x)dx ,

(3)

где

 

-,хе[у-Ь,у],

 

М у - х) = О

( 4 )

о,хе[)> -ь,у ].

 

Ясно, что ф.п.р.в. /у (у) отлична от нуля лишь для

v е [0,2б],

т.е. при >’ < 0 и у > 2Ь, / у(_у) = 0. Используя (3), (4), получаем

ЛО')=~т\ d x >

(5)

Ь~ Ja

 

где П = [у - Ь, >’] n [О,b].

Здесь первый интервал следует из (4), а второй из (2). Нетрудно убедиться, что

г

J'

[ dx = у, при V[о, ъ],

jdx = •

П| dx = 2b - у, при у е \b,2b\ .

у-Ь

Подставляя (6) в (5), получаем (1).

Задача 1.3.5. Заданы два центрированных, независимых между собой случайных вектора х и v размерности п и т с матрицами

ковариаций Рх и R , а также п и т -мерные векторы

77

 

 

х = АГу

у = Нх + V

 

где К

м Н - известные матрицы соответствующей размерности.

Для

составного

вектора

X = (xT,v 1,y 1,x )T найдите

матрицу

преобразования Т,

такую

что

Х - Т г , где z = (xT,v T)T

а также

матрицу ковариаций для вектора X = Тъ .

 

Ре ш е н и е. Матрица преобразования будет иметь вид

Е0

 

0

Е

 

 

 

Т =

Е

 

 

 

Я

 

 

 

к н

К

 

 

и таким образом

 

 

 

 

Р х

0 Р Х Н

Р Х Н Т К Т

0

R

R

R K

1

Р Х =

 

р у

 

Х

Н Р Х

R

Р У К

 

К Н Р х

K R

К Р }

К Р У К Т _

где Ру =НРХН Т +R.

Задача 1.3.6. Пусть компоненты л-мерного случайного вектора £, определяются как

£, =V,.+d,

где V,- - центрированные, некоррелированные между собой слу­

чайные величины, дисперсии которых равны г} i = ].n; d - цен­

трированная случайная величина, некоррелированная с v ( , i =\.n,

с дисперсией <3 2d .

Запишите выражение для матрицы ковариаций Ре вектора £. Дополнительно полагая, что все величины гауссовские, запишите выражение для ф.п.р.в. Считая дисперсии для всех v (- одинако-

2

о

---

выми, т.е. /-•

= г *•, / = 1 л , и, используя лемму об обращении мат­

риц, запишите выражение для матрицы, обратной Рг

Р е ш е н и е . Матрица ковариаций и ф.п.р.в. определяются как

Р£ = diag {/}2}+

, Д е) = N(e;О, Р Е),

где /)1ХИквадратная лхи-матрица, составленная из единиц. Используя соотношение (П1.1.50), получаем

2

?

г 1

1

£

2

>

____I

п*п

r E n +(TdInxn\

= -=-

ип

2 2

 

 

 

 

 

tlCTd+r

 

где Еп - лх/;-единичная матрица.

Задача 1.3.7. Найдите плотность распределения случайной ве­ личины y = axi + Рх2 , полагая, что Xj, х2 являются компонента­ ми гауссовского вектора с математическим ожиданием

х = (xj,х2 )т и матрицей ковариаций Рх =(<52 КЛ

, а а ,р -

 

известные параметры.

 

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

(y-ÿf

2<Jy = е >

2

у

где у н Оу определяются согласно (1.3.40), (1.3.41), т.е.

ÿ = ctx1 + Р х 2 ;

а 2 = а 2Р|"\ + р2Р22 + 2ар/Г

Контрольные вопросы

1.Почему решение задачи нахождения ф.п.р.в для с.в. у = g(x)

в условиях, когда функция g(x) имеет обратную функцию,

упрощается по сравнению со случаем, когда эта функция не­ однозначна?

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

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

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

иматрицы ковариаций для этих векторов?

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

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

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

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

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

Соседние файлы в папке книги