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

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

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

Применяя описанное выше правило нахождения ф.п.р.в. для функции от случайных величин и принимая во внимание тот факт, что

det dh{y) = det

coscp

-/isin tp

= A,

dyT

simp

Лсовф

 

можно получить следующее выражение для совместной плотности рас­ пределения амплитуды и фазы [44]:

 

(^costp-X])2

 

- \ 2

 

/ а,<|Мф) =

 

(i4sin<p-x2)

.,(1.3.11)

Â

+

Di

2K(D{D2)1/2 eXP j - i

 

A > 0 | ф |< %.

 

 

 

Интегрируя (1.3.11) по ф или

A , нетрудно найти плотности распре­

деления отдельно для амплитуды и фазы. Различным значениям парамет­ ров X], х2,D ],D2 будут соответствовать различные виды распределе­

ний [44, с. 122]. Так, в простейшем случае, когда 3cj=3t2 =0,.

D \ -D 2 =D, амплитуда и фаза между собой независимы, для амплитуды

получаем плотность распределения Рэлея, а для фазы - равномерную плотность в пределах от —71 до + к, т. е.

/ А(А) = Аехр{-(А2) / 2D)}/ Д Л>0;

(1.3.12)

1

ф е [ - л ,

+ 7 i],

 

/ Ф(ф) = 2п

1.3.13)

 

 

О

ф ё [ - 71,

+ п ] .

 

Обратное преобразование от полярных координат к декартовым осу­ ществляется с помощью (1.3.10). Таким образом, в условиях, когда со­ вместное распределение случайных величин А и (р имеет вид (1.3.11), случайный вектор, компоненты которого получены с использованием (1.3.10), будет гауссовским с плотностью распределения, задаваемой со­ отношением (1.3.7). В частности, когда амплитуда А, имеющая плотность распределения Рэлея, и равномерно распределенная фаза независимы

между собой, случайные величины

= A coscp и х2 = Asincp будут

центрированными независимыми случайными величинами с одинаковы­ ми дисперсиями. Данный результат представляется примечательным, поскольку гауссовские величины Х|,х2 формируются на основе нели­

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

Обратим внимание на то, что размерности векторов х и у в соотношении (1.3.6) предполагались одинаковыми. Вместе с тем часто возникает необходимость отыскания ф.п.р.в. случайного вектора у , связанного функциональной зависимостью с другими векторами, размерности которых не совпадают с размерностью у Например, нередко известна совместная ф.п.р.в. / Vx(v,x) для двух случайных векторов v,x, размерности п и ш и требуется найти ф.п.р.в. для m -мерного вектора у = g ](x,v). В этом случае целесообразно ввести составные векторы одинаковой размерности

’V

"у"

g|(x>v)'

H Z У"

V

У

X

X

-Z2_

X

и воспользоваться приведенным выше правилом для нахождения

ф.п.р.в. /zj,z2 ( У

^2 ) Для вектора

z = g(z).

Искомая плотность

fzj (^î ) - /у 0 ;)

далее

получается путем

интегрирования

/? b z2(? i,z 2) по

z2 - x .

Поясним

эту процедуру на следующем

весьма важном при решении задач оценивания примере.

♦ П р и м е р 1.3.4. Пусть задана ф.п.р.в. / v x (v, х) двух случайных векторов х , v размерности п и т и вектор у , связанный с этими век­ торами

y = s(x) + v ,

(1.3.14)

где s(x) - т -мерная вектор-функция известного вида.

Требуется найти ф.п.р.в. для вектора у

В данном случае вектор z целесообразно ввести в виде

У"

У

j(x) + V

У Z 2 ) + Z j '

У

X

X

. z2 .

При этом обратная функция z = A(z) будет определяться как

1

y

_

V

_

 

(N N

1

X

 

 

 

 

<>» l

X

1

_

i

Î

N.

Ю

 

N'

 

Z 2

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

 

= det

E - d s ( z 2 )

 

= 1,

 

l dzT

dz-y

 

 

 

0

-ly

 

 

 

 

 

 

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

 

 

 

 

 

 

 

f ü (z ) = fy,x (У> f \ , x

(^1 s(.z 2 )> ^2 ) — f \ , x (У ~

^0 (1.3.15)

Теперь для нахождения

искомой

ф.п.р.в. y= s(x) + v

достаточно

проинтегрировать стоящее справа выражение по аргументу х , т.е.

/z,(?1) = /у (T) = j Л,хO'-■*(*)»

 

(1.3.16)

Если векторы х и v независимы, то это выражение конкретизируется к виду

/у O0=J fv(y-s(x))fx(.x)dx.

В частном случае, когда у = х + v , используя (1.3.16), получаем

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

/ , М

= ] Л д О-- * ,* ) * .

(13.17)

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

 

/ у0 ') = |

Ш Ь Ь - х Ц х .

(1.3.18)

При нахождении ф.п.р.в. для суммы двух векторов в качестве

X + V

вектора z может быть также использован вектор z =

V

Тогда вместо (1.3.17), (1.3.18) можно получить (см. задачу 1.3.3) следующие аналогичные выражения:

f y ( y ) = \

ЛдО'.з'-1’)*’;

(1.3.19)

/,(У ) = |

/ . M A Ü '- v y v

(1.3.20)

Определение ф.п.р.в. для суммы двух независимых векторов по известным плотностям для слагаемых называется композицией ф.п.р.в. [87, с. 47]. Из (1.3.18), (1.3.20) следует, что композиция представляет собой интеграл свертки.

С использованием (1.3.16) или (1.3.18) можно, в частности, по­ казать, что сумма двух независимых с.в., имеющих одинаковое, равномерное распределение, будет представлять собой с.в. с ф.п.р.в. треугольного вида (см. задачу 1.3.4), а при увеличении числа слагаемых это распределение будет стремиться к гауссов­ скому. В этом случае говорят, что происходит нормализации ф.п.р.в. Более общий результат, связанный с нормализацией ф.п.р.в. формулируется в виде центральной предельной теоре­ мы. Суть этой теоремы заключается в том, что при формировании случайной величины в результате суммирования независимых ме­ жду собой центрированных с.в. с одинаковыми дисперсиями дока­ зывается, что при увеличении числа слагаемых ф.п.р.в. такой с.в. стремится к гауссовской ф.п.р.в. [44]. Это весьма важный с прак­ тической точки зрения результат, поскольку он служит некоторым обоснованием часто используемого предположения о гауссовском характере тех или иных величин. Действительно, если некоторая величина, например ошибка измерителя, формируется в результате суммирования ошибок, обусловленных различными независимы­ ми между собой факторами, то такую с.в. можно с определенной долей приближения считать гауссовской с.в.

С учетом приведенных выше соотношений можно таюке убе­ диться в том, что сумма совместно гауссовских векторов является гауссовским вектором. Более общий результат формулируется так: при линейных преобразованиях гауссовских векторов гауссов­ ский вид ф.п.р.в. сохраняется. Весьма существенно при этом, что гауссовской должна быть совместная ф.п.р.в. преобразуемых век­ торов. В частности, если совместная плотность, стоящая под зна­ ком интеграла (1.3.16) не является гауссовской, то в общем случае плотность для суммы у = х + V не будет гауссовской. Примеры,

иллюстрирующие этот факт, приводятся, в частности, в [85]. Если преобразуемые векторы независимы, то достаточно, чтобы гаус­ совскими были плотности каждого из слагаемых (см. задачу 1.3.2).

1.3.3. Линейные преобразования случайных векторов

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

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

х = Тх,

(1.3.21)

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

Решение этой задачи легко получить, если воспользоваться со­ отношениями (1.2.9), (1.2.10). В первое подставляется g(x) = Tx, а

во второе g(x) = Т(х-х))(х-х)Тт . В этом случае можем записать:

Мхх = М хТх = ТМхх - Тх ;

Мх {(х-х)(х-х)т \=МХ{г(х-х)(Г(х-х))т }=МХ{г(х-х)((х-х))гт }=

= Ш х|>-х)((х-х))т}7’т

Таким образом, получаем

х = 7х,

(1.3.22)

Р* =ТРхТт

(1.3.23)

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

♦ П р и м е р 1.3.5. Пусть задан вектор z , включающий два подвек­ тора х и V размерности л и т , т.е.

z = ( x \ v T)T,

(1.3.24)

идля него определены математическое ожидание и матрица ковариаций

ввиде:

х

M Tz =z =

(1.3.25)

V

 

M z{ z - z ) ( z - z ) r = Рх В

(1.3.26)

Вт Pv

им ную корреляцию двух векторов.

Требуется найти математическое ож идание и матрицу ковариаций

вектора

X

' Е

0" X

.У.

Н

(1 .3 .27)

Е_ V

Принимая во внимание (1 .3 .22), (1 .3 .23), легко получить

z

X

(1 .3 .28)

 

 

 

Hx + v

р х

|

 

РХН Г +В

Р

 

 

(1 .3 .2 9 )

_ВТ +НРХ !

НРХН Т +НВ +В ТН Т + Р \

Из рассмотренного примера, в частности, следует, что для век­ тора

у = Нх + V,

(1.3.30)

его математическое ожидание и матрица ковариаций определяют­ ся как:

ÿ =m + v;

(1.3.31)

Ру =HPxH r +НВ + В ТН Т + P V

(1.3.32)

Очевидно, что если векторы х и v некоррелированы, то

 

Ру =НРХН Т +Р"

(1.3.33)

Если предположить, что в соотношении (1.3.21) вектор

х гаус­

совский, т.е.

 

f x(x) =N(x;x,Px),

(1.3.34)

то, как отмечалось в подразделе 1.3.2, вектор хтакже будет гаус­ совским, а его плотность будет иметь вид

h (X) = Щх; Тх, ТРхТ т).

(1.3.35)

♦ l i p it Me p 1.3.6. Пусть размерности x и x

одинаковы, а матрица

T в (1.3.22) не вырождена. Убедимся, что если х

- гауссовский, то х

также будет гауссовским.

 

Принимая во внимание равенства [8] и используя (1.2.11). можно за­ писать:

/х(*) =

deter-1)

 

(2jt)"/2(detJP'Y)1/2

expj- 0.5((.v - T x ) f (Г т)“ 1 (Рх Г 1Т~1(х - 7х)}.

Т \ - 1 .

- 1ПТ*—1

(ТРХТ Г)~] = ( Т т) ~' (Рху

det(rrvr T) = (detr)2 detPJC detT-1 = I / del Г

Таким образом, легко убедиться в справедливости (1.3.35), поскольку

/х(*) =

г ехр{-0.5(х - х ) г (ГРхТ тГ 1

(2jt)"/2(det77>AT T)1/2

Если предположить, что в примере 1.3.5 совместное распреде­ ление двух векторов х и V гауссовское

fr

п

./x.vO',v) = 7V

X

1

VLV-*

X

рх

— “

____

I

I

в

pri (1.3.36)

то вектор z , определяемый согласно (1.3.27), также будет гауссов­ ским, т.е.

M Z) = N (Z ; Ï , P - \

(1.3.37)

а параметры этой ф.п.р.в. будут задаваться соотношениями (1.3.28), (1.3.29).

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

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

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

Задан двухмерный случайный вектор х = (Х[,Х2 )т с математи­ ческим ожиданием ,т = (л-!,,v2) и матрицей ковариаций

Р

о f

К

 

(1.3.38)

 

К

о \

Требуется найти математическое ожидание и дисперсию ска­

лярной случайной величины у , связанной с х как

 

y = ^/jX[ + ci2^ 2

(1.3.39)

где d\, <72 ~ известные числа.

 

Задача такого типа уже рассматривалась в примере 1.2.1. Заме­ тим, что ее решение также может быть легко получено с использо­

ванием (1.3.22), (1.3.23), если в качестве Т принять

T =(dA,d2) В

результате имеем

 

 

y —d\X\ +с/2%2,

(1.3.40)

oÿ = d [ erf

—d\d2K

(1.3.41)

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

проекции вектора х на некоторое произвольное направление, за­ даваемое единичным вектором (рис. 1.3.4)

е т = (sin т, cos т)т

(1.3.42)

Рис. 1.3.4. Определение проекции вектора x = (x i,x ? )T

на единичное направление

Поскольку величина, определяющая с учетом знака длину про­

екции, может быть представлена в виде

 

р = х, s in t + x 2 c o st,

(1.3.43)

для ее математического ожидания и дисперсии можно записать:

р =х: sinx + x2 COST;

(1.3.44)

Dp(т) = ст2 cos2 т + а 2 sin2 т + sin 2т

(1.3.45)

В частности, если х = (х15х2)т - центрированный вектор, то

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

Обращаем внимание на то, что дисперсия (1.3.45) может быть представлена в виде квадратичной формы

D0(t) = (sin т cos х)Р'

yCOS Xj

При изменении направления х единичного вектора (1.3.42) очевидно, что меняется и значение дисперсии.

Интересной представляется задача нахождения таких значений углов, при которых дисперсия принимает минимальное и макси­ мальное значения. Фактически это есть известная задача нахожде­ ния минимального и максимального значений квадратичной фор­ мы на единичной окружности [8]. Ее решение определяют собст­ венные числа матрицы ковариаций (1.3.38), характеризующие максимальное X] и минимальное Х2, ^ 1 - Х2 значения диспер­ сии, и соответствующие им собственные векторы. Если изначаль­ но матрица ковариаций диагональная, то направления, соответст­ вующие наибольшему и наименьшему значениям дисперсий, сов­

падают с направлениями координатных осей, при этом

Xj - а 2,

Х 2 = Ь 2

 

 

 

Из сказанного выше также

следует, что

если

вектор

x = (x j,x 2)T, задающий ошибку

местоположения

на плоскости,

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

скую ф.п.р.в. с математическим ожиданием (1.3.44) и дисперсией (1.3.45). Знание ф.п.р.в. для этой величины позволяет в полном объеме описать свойства ошибки вдоль заданного направления, что весьма важно при решении ряда задач обработки навигацион­ ной информации.

П р и м е р

1.3.7. Пусть компоненты двухмерного центрированно­

го гауссовского вектора х = (x i,X 2) T описывают ошибки местоположе­

ния

на плоскости

относительно некоторой заданной точки, и матрица

ковариаций этого вектора имеет вид: Р,V = ^

м~ Требуется найти

ч 1

V

ф.п.р.в. н ее параметры для с.в., определяющей длину проекции в на­ правлении прямой, расположенной под углом 30° относительно оси л~2

(рис. 1.3.5).

01 Y

!

•to .10

Рнс. 1.3.5. График ф.п.р.в. двухмерного гауссовского вектора и соответствующего ей эллипса

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