Файл: Олянюк, П. В. Оптимальный прием сигналов и оценка потенциальной точности космических измерительных комплексов.pdf

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

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

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

Добавлен: 23.10.2024

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

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

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

лов и для одиночного сигнала со случайной начальной фазой равна

 

 

w[q!y) = Kw(q) exp ( -

Э / Л д

/„ ( 2 Z ' / V 0 ) .

( I I I . 1.1)

Из

этой

формулы

видно,

что зависимость

плотности

веро­

ятности

от

величины оцениваемых параметров носит слож ­

ный нелинейный характер . В общем случае

закон

распреде­

ления отличается от нормального.

 

 

 

 

Однако

в случае приема сильных сигналов, который

имеет

наибольшее

практическое

значение, аргумент

бесселевой

функции

в

формуле

( I I I . 1.1)

значительно

превышает

еди­

ницу,

и

возможна

аппроксимация

этой

функции зависи­

мостью

 

 

 

 

 

 

 

 

 

 

 

 

 

I0(x) =

e-v : |

2кх

,

 

 

 

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

МЧ У) = K{w (q) exp (2Z.'/V0 ) exp ( -- 3,W 0 )

 

(III.1.2)

Д л я небольших

отклонений параметров

движения от

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

р а з л о ж е ­

нием показателей

экспонент формулы (III.1.2)

в

ряд Тей­

лора в окрестности априорных значений. Производя это раз ­

ложение,

получаем

 

 

 

 

 

 

 

 

 

2 (q)

- у

Э (q)

=

Z (qa ) -

i - 3 (q a )

 

Z\ ( q e ) -

 

 

 

 

 

(ft — Я ад

+

S

 

 

 

 

 

 

 

 

 

 

 

2

i.j

 

 

 

 

 

 

 

 

 

{Я l — ЯаМЯ]

— Яа)> + •

 

(III.1.3)

Здесь

q a

априорное

значение

 

вектора

оцениваемых

пара ­

метров;

qai

i-я с о с т а в л я ю щ а я

этого

вектора;

Z'L

( q j —

значение

первой производной

АК Ф по

составляющей

вектора

 

оцениваемых

параметров

точке

q = q a ;

Z'ij (Ча ) з н а ч е н и е

второй производной

АК Ф по г-\\

и /-й

составляющим

в этой ж е точке.

 

 

 

 

 

 

 

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

57


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

условная плотность

вероятности

приема

данной

реализации

и з о б р а ж а е т с я

гауссовой

кривой и распределение

является

нормальным .

Этот

случай

з а с л у ж и в а е т более обстоятельного

рассмотрении,

так

как удается

получить

простые

аналитиче ­

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

априорным значениям параметров движения .

 

 

 

 

 

Примем, что ошибки

априорных

данных подчиняются нор­

мальному

 

закону

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

и

 

плотность

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

и з о б р а ж а е т с я формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

(q) =

С е х р

- • j

( q - q Q ) T

B - ' ( q - q n )

 

 

 

(III.

1.4)

где q a

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

 

вектора

априорных

з н а ­

чений

параметров

движения;

В 0

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

матрица

ошибок

априорных

данных.

 

 

 

 

 

 

 

 

 

 

 

 

 

П о д с т а в л я я

(III.1.3)

и

(III.1.4)

в

 

(III.1.2), получаем

сле ­

дующее векторное

соотношение:

 

 

 

 

 

 

 

 

 

 

 

 

 

w(q

у) =

Кг

exp j -~

(q — q n ) T

B j ' f a -

q0 1 —

 

 

-L[3(qe )-2Z(qf l )J+J-(q-qe )

 

Z ' ( q J -

Э'

(qa )

+

N 0

 

 

 

 

 

 

 

Л'о

 

 

 

 

 

 

 

 

 

 

 

 

 

+ -L(q_qeyr

Z " ( q „ ) -

y 3 " ( q j

( q - q „ )

,

 

( Ш . 1 . 5 )

которое в развернутом виде может

быть

переписано

так:

 

(?) .

ft

. •

• • . 1тМ

=

Къ

exp

j - ~Y.

 

ill

Qal)

Baii

 

~

Я a

 

 

 

 

 

 

 

 

 

 

(

 

 

* /,/

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

а1

,

q„2, • • • ,

am

)—2Z{q„

u

q

a 2

, . . ,

О ] ' ! "

Wo

Л'о

\Э(д

O

 

 

 

 

iv

+

Qai)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

58


1

71

 

 

Zu(</a\

ЯаО,

qam)

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

(9J —

Яа})

При этом формулу дл я

апостериорной

плотности вероят­

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

 

 

 

МФ) = К->

Г

1

Л

Ч ) т В - 1 (q — q)

( I I I . 1.6)

ех р ]~

(q -

где q — значение

вектора

параметров, при котором дости­

гается максимум

апостериорной

плотности

вероятности.

Так как апостериорное распределение симметрично, то q представляет собой оптимальную оценку, соответствующую обеим упомянутым ранее функциям потерь. Следовательно,

Л

вектор q можно называть вектором уточненных значений па­

раметров

движения .

Буквой

В

обозначена

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

матрица результирующих ошибок измерений.

Это

квадрат ­

ная

матрица,

размерность

которой

 

определяется

размер ­

ностью вектора

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

параметров .

 

 

 

 

Произведя

соответствующие

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

получаем

 

 

 

 

 

 

 

 

 

Л

 

 

 

д л я

матрицы

В

и

вектора-столбца

q

следующие

соотно­

шения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В - 1 =

в - 1

Z " ( q a ) - - | 3 " ( q 3

)

 

(III.1.7)

 

 

 

 

 

 

О L

 

 

 

 

 

 

 

 

 

А

 

 

2

В Z ' ( q e ) - y 9 ' ( q n )

 

 

(Ш.1.8)

 

 

 

 

 

 

 

 

которые

могут

оыть т а к ж е записаны

в виде

 

 

 

 

 

11 =

 

I

I

- -

 

 

 

1

 

Э], ( q e )

 

*Г?

 

 

 

 

 

 

 

 

 

 

 

 

о I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Zj(<\a)

II

 

^\\Э}(Ца)

З д е с ь Z" (q a )=

|| Z"u(qa)

\\ -—

матрица

вторых

производных

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

интеграла

по

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

п а р а м е т р а м

59.



д в и ж е н ия

(производные рассчитываются в точках, соответ­

ствующих

априорным д а н н ы м ) ; Z ' ( q 0 ) = || Zj(q„) || вектор-

столбец первых производных корреляционного интеграла по

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

п а р а м е т р

а м в тех ж е

точках.

Формулы

( П Ы . 7 ) и

( I I I . 1.8)

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

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

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

возможность устранения к а ж у щ е й с я

неоднозначности

кро­

ется в избыточности измерительной

информации

с

одной

стороны, и в использовании априорных данных — с другой.

Остановимся

теперь на

свойствах корреляционной матри­

цы результирующих

ошибок измерений

( I I I . 1.7).

Если

на­

звать

элементы

матриц,

обратных

корреляционным

матри ­

цам

ошибок, мерами

точности, то

смысл

формулы

( I I I .

1.7)

можно выразит-ь следующим образом . Мера точности резуль ­

татов измерений равна сумме мер

точности

априорных дан ­

ных и произведенных измерений.

Очевидно, что если точ­

ность априорных данных слишком

мала, то

первые слагае ­

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

ветствовать точности априорных данных.

 

В формулах (III.1.7) и (III.1.8)

фигурируют

производные

от корреляционного

интеграла

в

точке, соответствующей

априорным данным . Корреляционный интеграл

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

собой функцию, у б ы в а ю щ у ю по

мере увеличения различия

между априорными

и истинными

значениями

п а р а м е т р о в

60