Файл: Кондратьев, С. Л. Применение метода функционального моделирования для оценки помехоустойчивости систем связи.pdf

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

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

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

Добавлен: 21.10.2024

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

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

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

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

Можно показать, что если параметр у является случайной величиной, с известной совместной нормальной плотностью ве­ роятности W(u, у), то оценка максимального правдоподобия также совпадает со средним значением (теперь уже условным).

Какую же непрерывную обработоку u(t, у) в общем случае следует производить, чтобы получить максимально правдоподоб­ ную оценку различных параметров? Прежде всего заметим, что все предыдущие результаты были основаны на предположении нормальности закона распределения ошибок, а при нормальном шуме на входе линейной системы оценивания это было очевид­ ным. К сожалению, многие параметры, которые нас интересуют в теории связи, получаются путем нелинейной обработки сиг­ нала. Напомним, что для получения амплитуды сигнала необ­ ходимо использовать идеальный амплитудный детектор, выпол­ няющий операцию

 

A(t)=]/ а»(0+«*(*).'

 

(2.3.2)

Для получения фазы

используется фазовый

детектор

 

 

 

« W = a r c t g - i ^ L .

 

 

(2.3.3)

Производная

от ср (t)

дает

мгновенную

частоту,

также,

оче­

видно, являющуюся

нелинейной функцией,

т. е.

 

 

 

mlWj'{t)uW-*{t)~u(t)

 

 

 

( 2 з 4

)

 

2 т с

d t

и» (*)+

«•(<)

 

 

где и (£)—сопряженное по Гильберту

значение сигнала.

 

Поэтому

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

ошибок

оценки

не

обязано

6biTb

нормальным, если даже на вход вместе с регулярным сигна­

лом

поступает белый

гауссовский шум. Но

в

таком

случае

максимально правдоподобная оценка

может

оказаться

совсем

не оптимальной. Общего строгого

решения

задач

оценки при

нелинейных

процедурах

на

ограниченных

интервалах

наблю­

дения до сих пор не имеется.

 

 

 

 

 

 

 

 

При этом используются

процедуры

оценки,

которые

удов­

летворяют

какому-либо

критерию. Все они в

той или

иной

степени

исходят

из

того,

что

информация

об

интересую­

щем

нас

параметре

содержится

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

плотности

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

(7,) =

W (т/и).

 

 

 

 

 

 

 

 

59



К ним относятся:

 

 

 

 

 

 

 

1. Оценка

по минимуму

среднеквадратичной

погрешности,

по которой

оценка т* выбирается

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

чтобы

 

 

 

D (Д,) = j

(т-Т*)2 Wa

(T j )

rfT=min.

(2.3.5)

При

этом,

как и ранее,

а качестве оценки

используется сред­

нее

Значение

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

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

 

(2.3.6)

 

 

 

 

1* = !г - Л ( т М | . .

 

 

Выше было показано, что такая оценка является

оптимальной

при

нормальном

распределении. Убедимся,

что эта

обработка

удовлетворяет

 

минимуму

среднеквадратичной

Погрешности,

Возьмем первую

производную

от

дисперсии D (Д,)

и прирав-

няем нулю, а затем найдем вторую производную. Если вторая производная положительная, то действительно, такая оценка обеспечивает минимум среднеквадратичной ошибки. Читателю рекомендуется проделать это самому.

2.Вторая процедура основана на максимизации обратной (апостериорной) вероятности W(i,ju), т . е . на выборе такой оценки т*, которая доставляет максимум апостериорной ве­ роятности, или, что то же, In WYj/и).

3.Третья процедура основана на максимизации функции правдоподобия И^(и/т)=/(ч), т. е. на выборе такой оценки -f*, которая обеспечивает выполнение равенства:

max IT («/?) = max Щ и / г ) ,

(2.3.7)

либо, так как логарифмирование не смещает максимум, вместо этого равенства можно записать:

max In W(и/г) = max In W(lift*).

(2.3.8)

Обе последние процедуры не единственные,

более того, они

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

60


тем обстоятельством, что в асимптотике (т. е. при очень боль­ шом отношении сигнал/шум) распределенияамплитуды, фазы и частоты близки к .нормальным. В частности, обобщенный закон Рэлея становится практически нормальным.

Покажем связь между оценками по максимуму апостериор­ ной плотности и максимуму функции правдоподобия. Восполь­ зуемся для этого формулой Байеса. Функция плотности обратной (апостериорной) вероятности равна

 

 

 

 

 

 

 

 

 

 

 

W(u)

 

 

 

(2.3.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

W(u/i)—функция

 

правдоподобия

(если.ее

рассматривать

как

функцию

параметра

т), W (и)

и

^(ч)—априорные

плот­

ности

указанных

величин.

 

 

 

 

 

 

 

 

 

Левая часть показывает распределение параметра при усло­

вии,

что на

входе

имеется

некоторое

значение

напряжения

и.

Следовательно, можем записать:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.3.10)

 

Рассмотрим

два

случая:

 

7 = 7 , —постоянная величина,

 

 

а)

Пусть

f

есть

некоторая

а

поэтому

плотность

вероятности

W(f )

есть 8-функция. Тогда

W ( Y / w ) = £ I

^ ( И / Т ) . Т . е. если

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

яв­

ляется' единичной функцией, то определение

оценки по

мак­

симуму

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

вероятности

 

 

 

 

 

и

функции

правдоподобия

отли­

 

 

 

 

 

чается

на

постоянный

множитель,

 

 

 

 

 

поэтому их использование дает одит

 

 

 

 

 

наковый результат.

 

 

 

 

 

 

 

 

 

 

 

б)

Пусть

априорная

плотность

 

 

 

 

 

вероятности

 

W(f)

равномерна

в

 

 

 

 

 

области

 

Г (f е Г).

Прономируем

 

 

 

 

 

область

Г к

единице:

^-

— 1;

Рис. 22.

 

 

 

 

 

 

 

- i - ->

+

1 .

 

Будем рассматривать случай, когда отношение

сигнала

к

 

шуму

велико,

но

тогда

функция

 

W(«/r)

очень

«узка»

(см.рис. 2.2), т.е. имеет

малую

дисперсию.

 

 

•Отсюда видим, что априорная плотность, если она равно­

мерна,

не

влияет

на качество

оценки,

т. е. все

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

функцией правдоподобия. Более того, если даже истинное рас­ пределение W(y) неравномерное, но изменяется медленно, то с хорошим приближением можно не считаться с неравномер­ ностью при условии, что дисперсия параметра и при данном у невелика (т. е. функция правдоподобия очень узка). Однако уже она будет тогда, когда отношение сигнала к шуму оченьвелико. В пределе, когда ^ ° с > ^ э „ , можем делать предположение,

61


что априорное распределение равномерное, так как функция правдоподобия «вырезает» лишь ее узкий участок.

Таким образом, если хотим получить хорошую оценку, а сле­ довательно малую ошибку, то при неизвестном априорном рас­ пределении вполне справедливо предполагать его равномерным (за это расплачиваемся большим с^У^п)-

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

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

Wn (Д„ Д2, .. . , d ^ - J ^ J Y ^ 2 4 ' 2 , (2.3.11)

если за /V считать число отсчетов по Котельникову. N=2FTn, будет функция правдоподобия, рассматриваемая как некоторая функция параметра т (при условии, что сигнал является по­ стоянным, а шум «белым» гауссовским):

 

 

 

 

 

 

 

О

 

 

 

 

Логарифм

этого выражения

 

 

 

 

 

 

(2.3.12)

 

 

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I n L

frMn(

i - a ) * )

-

Д - j * i u

W-c&

f ) l 2

^ -

( 2 - З Л 3 )

 

 

 

 

 

 

 

о

 

 

 

 

 

Взяв производную

по

исключим

I n ^ ( у ^ -

j/vjи

н а и

Д е м

значение оценки. Покажем это на примерах.

 

 

 

 

1.

Найдем

оценку

амплитуды

сигнала

при

условии, что

^ У ^ п

велико:

,

 

 

 

 

 

 

 

 

 

 

 

 

 

u(t)=c(t,

 

т ) + » ( 0 .

 

 

 

 

Искомый параметр—амплитуда

а0. Запишем сигнал в виде

 

 

 

 

c(t, т )=а 0 с а (0>

 

 

 

(2.3.14)

где

( 0 = 1 -cos (о)^+<р)—сигнал

с единичной

амплитудой,

т. е.

пронормированный к единице по aQ=^AQ

и известной

фазе ».

Нужно

оценить а0,

наблюдая

напряжение и на интервале Т.

Воспользуемся

выражением

для функции правдоподобия:

 

62