Файл: Кондратьев, С. Л. Применение метода функционального моделирования для оценки помехоустойчивости систем связи.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