
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfГ Л А В А В Т О Р А Я
О Б Н А Р У Ж Е Н ИЕ РАДИОЛОКАЦИОННЫХ СИГНАЛОВ ЦИФРОВЫМИ УСТРОЙСТВАМИ
2.1.В В Е Д Е Н И Е
Вданной главе рассматриваются вопросы синтеза структуры и анализа эффективности цифровых устройств обнаружения радиолокационных сигналов на фоне по мех. В настоящее время известны оптимальные и квази оптимальные алгоритмы обнаружения как для непрерыв ных сигналов, так и для последовательностей. С матема тической точки зрения нет никакой проблемы при пере ходе от непрерывных алгоритмов к их цифровым экви валентам: чем меньше шаг квантования и интервал вре менной дискретизации, тем ближе свойства цифрового эквивалента к свойствам его аналогового прототипа. Однако на практике при построении Э Ц В У обработки информации такой прямой подход не всегда целесообра зен.
Во-первых, уменьшение шага квантования и интерва ла временной дискретизации приводят к существенному усложнению аппаратуры ЭЦВУ: увеличивается разряд ность чисел, представимых в ЭЦВУ, растет емкость ЗУ, повышается требуемое быстродействие узлов ЭЦВУ, что ведет к снижению надежности ЭЦВУ и увеличению га баритов, потребляемой мощности, стоимости. Во-вторых, хорошо известно, что ряд обнаружителей обладает весь ма высокой эффективностью уже при двух-четырех уров нях квантования, т. е. при наиболее простой технической реализации. К их числу относятся: некогерентные бинар ные обнаружители (см. § 2.3), когерентные фазовые обнаружители (см. п. 2.4.2), устройства спектроанализа [39] и ряд других.' При этом потери в .пороговом сигнале или времени накопления за счет квантования не превы шают 2—3 дБ . С другой стороны, увеличение шага вре менной дискретизации до величины интервала корреля ции обнаруживаемого случайного процесса при обнару жении на фоне узкополосного шума ухудшает пороговый
70
сигнал всего на 3 дБ по сравнению с непрерывными ана логовыми накопителями. Последние технически нереализуемы, если время накопления составляет десятки или сотни секунд [40].
Оптимальной структурой Э Ц В У обнаружения счи тается та, которая при заданной сложности технической реализации обеспечивает минимум потерь порогового сигнала по сравнению с аналоговым алгоритмом, синте зированным методами теории статистических решений. Нахождение оптимального, в этом смысле, решения за дач дискретизации и квантования в общем виде затруд нительно. Д а ж е первая задача не решается однозначно: существуют различные методы синтеза линейных им пульсных фильтров (см. гл. 4). Поэтому в рамках метода синтеза ЭЦВУ, который мы называем синтезом структу ры по аналоговому прототипу, в общем случае не удает ся однозначно определить оптимальную структуру и па раметры ЭЦВУ .
Выбрав структуру и параметры цифрового обнаружи теля, необходимо определить его эффективность, т. е. вычислить пороговые сигналы или построить характери стики обнаружения. Сравнение эффективности цифрово го эквивалента (технически реализуемого) и его анало гового прототипа (обычно технически нереализуемого) позволяет определить проигрыш в -пороговом сигнале как плату за упрощение. Если оказывается, что требуемое число уровней квантования велико (128—51'2), а шаг временной дискретизации мал, то статистические свой ства сигналов « помех, преобразованных в цифровую форму, мало отличаются от статистических свойств не прерывных сигналов и помех, а цифровой эквивалент аналогового алгоритма оптимален в той же мере, что и его прототип. Если же обнаружитель работает при малом числе уровней квантования (2—16), то статистические свойства преобразованных в цифровую форму сигналов и помех будут существенно отличаться от статистических свойств непрерывных сигналов и помех. Поэтому циф ровой эквивалент аналогового прототипа не будет опти мальным для сигналов и помех с новыми статистиче скими свойствами.
Отсюда возникает другой подход к синтезу структуры цифрового обнаружителя: определив статистические свойства сигналов и помех, преобразованных в цифровую форму, синтезировать оптимальный алгоритм, используя
71
методы теории статистических решений. Как будет пока зано в п. 2.6.2—2.6.4, при большом числе уровней кван тования, а также для слабого сигнала получаемая структура совпадает с цифровым эквивалентом анало гового прототипа. При малоразрядном преобразовании и произвольном сигнале синтезируемая структура отли чается от известных ранее п обладает рядом интересных свойств. В частности, в рамках выбранной марковской аппроксимации дискретных распределений подобные устройства являются структурно-инвариантными по от ношению к широкому классу сигналов и помех, стати стические свойства которых влияют лишь на параметры устройства. Такой результат является следствием того, что синтезированные устройства (см. п. 2.6.3) соответст вуют табличным алгоритмам, а с помощью таблиц мож но задавать любые функции.
Наряду с отмеченными выше двумя методами синтеза существует и третин, который уместно назвать эвристи ческим. При этом используются характерные особенности сигналов и помех, преобразованных в цифровую форму, известные из опыта. Эвристические структуры характер ны тем, что они просты в технической реализации, но получаются не в результате статистического синтеза и поэтому выбор того или иного варианта в значительной мере случаен.
2.2.О П Т И М А Л Ь Н Ы Е А Л Г О Р И Т М Ы
О Б Н А Р У Ж Е Н И Я Р А Д И О Л О К А Ц И О Н Н Ы Х
СИ Г Н А Л О В
2.2.1.Синтез оптимальных и достаточных приемников. Из теории статистических решений известно [69], что структура оптимального устройства обнаружения детер минированных сигналов на фоне помех в случае двух простых гипотез (Но — присутствует только помеха, НІ — присутствует сигнал с помехой) определяется коэффи циентом правдоподобия:
|
|
Л(и) = |
юя (ы|в1 )/гии (и|в0 ), |
(2.1) |
|
где wn |
{и | б,), wn (и | б0) — функция правдоподобия смеси |
||||
сигнала |
с |
помехой и помехи |
соответственно |
для выборки |
|
|
|
-у |
|
|
|
объема |
п\ |
и = (и1,и„ |
...,ип)т |
— выборочный |
вектор. |
72
Структура оптимального приемника в случае М про стых альтернативных гипотез определяется совокупно стью М—1 коэффициентов правдоподобия:
|
Лг -(u) = wn{u \ bi)fwn (иI 80 j, і = 1,/W — 1. |
(2.2) |
|||
Обычно |
техническая |
интерпретация |
соотношений |
||
(2.1) и |
(2.2) |
довольно сложна. Поэтому на практике |
пе |
||
реходят |
от коэффициента |
правдоподобия |
к его информа- |
ционному эквиваленту J (и) —достаточному приемнику, имеющему ту же эффективность, что и оптимальный, и, как правило, более простую структуру. Такой переход
осуществляется |
чаще |
всего |
логарифмированием А(и) и |
|
|
|
-> |
моделированием |
тех |
членов |
1 п Л ( « ) , которые учитывают |
взаимодействие сигнала с помехой [57]. |
|||
В .реальных |
условиях сигналы, отраженные от целей, |
являются либо квазидетермииированными, т. е. при известной форме колебания содержат ряд неизвестных параметров (амплитуда, частота, фаза и т. д . ) , либо стохастическими. Эхо-сигналы могут быть стохастиче скими по двум причинам: во-первых, вследствие излу чения шумоподобных ('чаще всего узкополосных) колеба ний, во-вторых, из-за шумоподобных статистических свойств отражающей поверхности цели и среды распро странения, модулирующих излученные сигналы. В свою очередь, помеха также содержит ряд неизвестных пара метров. .В таких ситуациях мы имеем дело со сложными параметрическими гипотезами, когда необходимо выне сти решение о принадлежности выборки к распределению данного класса, определяемого одним или несколькими параметрами. Для простоты будем считать, что парамет ры помехи известны. В общем случае полезный сигнал
можно записать |
в виде |
' |
|
||
|
|
|
S{t) = |
S(t,aJ), |
|
|
-> |
|
|
|
|
где |
а = (а,, а2 , . . . ) т — в е к т о р неизмеряемых |
параметров; |
|||
-у |
|
|
|
|
|
Я = |
(Я 1 Д 2 |
) т —вектор измеряемых параметров. |
|||
|
Приведем |
классификацию |
функций, |
выполняемых |
|
устройством |
обнаружения [42]. |
|
|
||
|
1. Простое |
обнаружение — проверка |
двух альтерна |
тивных гипотез, обнаружение полностью известного сиг нала S(>t) на фоне помехи с известными статистическими свойствами.
73
'2. Сложное обнаружение — проверка сложных аль тернативных гипотез, обнаружение сигнала с неизвестны-
ми параметрами S(t, а) без их измерения.
|
|
-> |
3. Простое |
измерение —обнаружение ,гсигнала |
5 ( ^ Д ) |
|
-» |
|
и измерение его |
неизвестных параметров Я. Эта |
задача |
может быть сведена к задаче проверки многих простых
альтернативных |
гипотез. |
|
S(t,a,X) |
||
[ |
4. Сложное |
измерение —обнаружение сигнала |
|||
с |
измерением |
его неизвестных |
параметров Я. Эта |
задача |
|
эквивалентна |
задаче проверки |
многих сложных альтерна |
|||
тивных гипотез. |
|
|
|
Впрактике радиолокации всегда осуществляется сложное измерение.
Вслучае сложных альтернативных гипотез оптималь ного обнаружителя, как правило, не существует. Исклю чение составляют, например, обнаружитель узкололосных сигналов с неизвестной фазой на фоне белого гаус сова шума |[43], а также байесовский оценочно-корреля ционный обнаружитель [27], рассматриваемый далее.
При проверке сложных гипотез важным понятием является равномерно наиболее мощное решающее пра вило [44]. Равномерно наиболее мощным (т. е. опти мальным) называется такой обнаружитель, который при заданной вероятности ложной тревоги F = const обеспечи
вает |
максимум |
вероятности |
правильного обнаружения |
|||||
£>=max{£>t (Si)}, |
где |
X — множество |
индивидуальных |
|||||
индексов |
обнаружителей. |
Несмещенные |
обнаружители |
|||||
образуют |
более |
узкий |
класс |
и обеспечивают |
D(S)^F |
|||
для |
всех |
|
и D(S)^\F. |
для всех 5 е Я 0 . |
Соответст |
венно оптимальный несмещенный обнаружитель обеспе чивает D—<max{Di(S)} для класса несмещенных обнару-
жителей.
Чаще всего используют обнаружители, обладающие свойством локальной оптимальности, т. е. обеспечиваю щие выполнение условий D = raax {£>i(S,)} при F = const
н,->я0
(оптимальность для слабого сигнала). В дальнейшем основное внимание будет уделено локально оптимальным обнаружителям. При этом будет использована теория проверки простых гипотез. Иногда бывает полезно осу-
74
ществить предварительное преобразование (редукцию)
выборочного |
U |
и |
параметрического |
01Э(Я,а) |
про |
|||||
странств и искать |
решение только |
в классе |
обнаружите |
|||||||
лей, использующих такие |
преобразования |
(см. л. 1.7.4). |
||||||||
Пусть а —вектор |
неизмеряемых |
параметров, |
заданный |
|||||||
в области |
Л, |
a wm(a) |
— плотность |
вероятности |
распре- |
|||||
деления а. |
Тогда |
структура |
оптимального |
приемника оп |
||||||
ределяется |
в соответствии |
с |
выражением |
|
|
|||||
|
|
Л(и) = < Л ( и | а ) > „ т Д |
= |
|
|
|||||
= |
f A ( w , a)da= |
f Л(u |а)oi m (а)do . |
(2.3) |
АА
Если интеграл (2.3) не вычисляется аналитически, то можно заменить его конечной суммой, т. е. осуществить интегрирование методом прямоугольников. Для скаляр ного параметра а имеем:
М |
(/ + 1)До |
м _^ |
Л (и) « £ Л (и |а,) |
f |
tw, (a) da = 2 Л (« | аг-) рг-, (2.4) |
1=1 |
|'Да |
1=1 |
где Ла=А/М. В этом случае получается многоканальная структура с числом каналов М. Используя более точные интерполяционные методы (метод трапеций, метод Симпсона), для получения той же точности можно ограничить ся меньшим числом каналов.
|
|
|
|
-> —> |
Для вычисления~(2.3) |
можно |
также |
рчзложить |
Л(ц|а) |
в ряд и, ограничившись |
некоторым отрезком ряда, про- |
|||
ц?вести почленное усреднение в соответствии с |
(2.3). |
|||
Е^ли параметры, характеризующие гипотезу, измеряе |
||||
мые, то |
|
|
|
|
А(и) = < А ( и | Я ) > - = |
f A(u|l)a»i(X)dt, |
(2.5) |
||
где Я — вектор измеряемых параметров, |
заданный в |
облас |
ти Ф. В этом случае возможны следующие приближения,
приводящие к построению .приемников совместного |
обна |
ружения — измерения. |
|
1. Замена интеграла (2.5) произведением |
|
А ( " ) = Л г а а > , Я 0 ) Ф 0 , |
(2.6) |
75
где |
Ф 0 — эквивалентная |
величина |
области |
Ф, |
Соотноше |
|||||||
ние |
(2.6) |
основано |
на |
том, что |
в |
окрестности |
некоторого |
|||||
значения |
параметра |
Я. величина |
Л (и, Я0 ) |
с |
большой веро- |
|||||||
ятностью |
значительно |
превосходит |
значения |
Л(ы, Я) |
||||||||
вне |
этой |
области. |
Этот |
способ |
формирования |
статистики |
||||||
Л (и) = max Л (и, Я) |
по |
методу |
максимума . коэффициента |
|||||||||
правдоподобия |
сводится |
либо |
к |
многоканальной |
системе |
|||||||
(т. е. системе |
параллельного поиска), либо |
к |
одноканаль- |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
-> |
ной (т. е. последовательной), когда параметр Я имеет смысл времени запаздывания. В обоих случаях имеется
возможность получения оценки Я параметров Я.
Алгоритм формирования Л (и) путем отбора максималь ного выходного эффекта одного из М каналов (2.6) про ще в реализации, чем алгоритм усреднения (2.4). Мож но показать, что
1пЛ(ц) = 1пЛт а х й?) — I n A f . |
(2-7) |
Иными словами, платой за упрощение алгоритма является уменьшение выходного эффекта достаточного приемника, т. е. увеличение порогового сигнала. Однако на практике это увеличение невелико: при УИ = 104 поро говый сигнал увеличивается всего на 2 дБ [45].
2. Применение к статистике (2.5) обобщенной теоре мы о среднем [13] приводит к выражению
|
Л (и) = |
f Л (и | Я) wt |
(Я) S, |
= Л (и\ Я* (и)), |
(2.8) |
||||
|
|
Ф |
|
|
|
|
|
|
|
где |
Л * ( ы ) с Ф - |
|
|
|
|
|
|
|
|
Выражение |
(2.8) |
соответствует задаче |
проверки |
слож |
|||||
ных |
гипотез. |
В этом |
случае |
можно |
использовать |
крите |
|||
рий, |
максимального |
правдоподобия |
[6], |
в |
соответствии |
||||
с которым необходимо образовать |
статистику |
|
|||||||
|
|
т а х ш п ( « | Я * ( « ) , |
01)/вим |
(и [ б0 ), |
|
|
|||
|
|
|
|
-> ->• |
|
~¥ |
-+ |
-> -> |
|
которая дает |
[оценку |
Я (и). 'Положив Я*(ы) = Я(ы), |
полу |
||||||
чим |
значение''выражения (2.8). |
|
|
|
|
76
Оценка |
Оценка |
|
я/35 |
Ключ |
|
Ми) |
||
|
||
ПУ |
|
|
|
»о |
|
Рис. 2.1. Самонастраивающаяся система совместного обнаружения - |
||
измерения. |
|
|
Следует отметить, что максимально |
правдоподобные |
обнаружители могут иметь структуру и параметры, со впадающие с оптимальными байесовскими обнаружите
лями [45]. |
|
|
|
|
(2.8), |
|
|
|
|
|
Структура, |
соответствующая |
изображена |
на |
|||||||
рис. 2.1. Она состоит из двух частей: |
|
|
|
|
||||||
— системы |
максимально правдоподобной |
оценки |
пара- |
|||||||
метров |
Я по выборке и; |
|
|
|
|
|
|
|
||
— системы |
со |
структурой, |
соответствующей |
точно |
||||||
известному сигналу |
S(t,Xn), |
изменяющей |
свои параметры |
|||||||
|
|
|
|
^» -> |
|
|
|
|
|
|
в соответствии с оценкой |
Я (и) |
и |
воздействующей |
на по |
||||||
роговое |
устройство (ПУ). При |
этом обучающая выборка |
||||||||
должна |
быть |
классифицирована -и соответствовать гипо |
||||||||
тезе Н\К Таким образом, алгоритм |
(2.8) |
есть |
самонастра |
|||||||
ивающаяся система совместного |
обнаружения—измерения. |
|||||||||
Поскольку оценки максимального правдоподобия при |
||||||||||
весьма |
слабых |
ограничениях на вид плотности wn(u\X, |
б,) |
|||||||
являются асимптотически |
эффективными, |
асимптотически |
||||||||
несмещенными |
и |
асимптотически |
нормальными |
и |
схо- |
дятся к истинному значению параметров Я (и) —+ЯИ [46], можно считать, что обнаружитель (2.8) при наличии классифицированной обучающей выборки является асим птотически эффективным и при увеличении объема вы-
Если выборка и соответствует гипотезе Я 0 , получаем псевдо-
оценку Я и параметров X. В этом случае потребуется дополнительное рассмотрение поведения системы.
77
боркн (п—>-оо) он эквивалентен обнаружителю полно стью известного сигнала. Однако при малом объеме вы борки п эффективность такого обнаружителя может ока заться невысокой.
""При наличии нескольких целей, разрешимых по пара-
метру Я, целесообразно построить многоканальную систе-г му, осуществляющую многоальтернативные решения: обна ружение с оценкой неизвестного параметра. Выходные
эффекты |
всех М каналов, т. е. Л (и, Яг-), і = |
1, М обраба |
тываются |
независимо. В тех случаях, когда |
отраженный |
сигнал содержит несколько неизвестных параметров; воз можны различные комбинации рассмотренных выше ме тодов формирования статистики. Конкретные алгоритмы
обнаружения и их цифровые эквиваленты будут |
рассмо |
||||
трены ниже. |
|
|
|
||
Необходимость в различных интерпретациях [форму |
|||||
лы (2.4), |
(2.6)—і(2.8)j |
связана с невозможностью анали |
|||
тического |
вычисления |
интеграла (2.3). В тех же случаях, |
|||
когда |
такое вычисление возможно, структуры |
получае- |
|||
|
|
|
|
-+ -+ |
|
мых |
обнаружителей |
явно зависят как от Л(м|а), так и |
|||
от wm(a). |
В то же |
время известны и структурно-инва |
риантные обнаружители: оценочно-корреляционный об наружитель [26, 27, 47] и цифровой обнаружитель, рас сматриваемые в пп. 2.6.2, 2.6.3.
Оценочно-корреляционный обнаружитель, предложен ный впервые Р. Л . Стратбновичем и Ю. Г. Сосулиным [26, 47], обладает рядом замечательных свойств и допус кает широкие обобщения (27]. Рассмотрим задачу про верки сложной гипотезы НІ против простой альтернати вы Я 0 :
Hl:u(t)=S(x*(t), t]+n(t),
Но : u(t) =яі(0>
где п (t) — белый гауссов шум с нулевым средним: п (t) — О
ифункцией автокорреляции n{t)n(t-\--z) = b{t), 8(т) —
дельта-функция, 5 — заданная функция |
своих аргументов, |
|||
в частности S[x* |
(t), t] — х* |
(f) — произвольный, |
немарков |
|
ский и негауссов, |
процесс. |
S[x*(t),t] |
имеет |
известную |
78
плотность вероятности и ограниченную энергию
І_____
JS 2 [x*(t),t] |
Л < о о . |
о |
|
В частном случае jc*(f) = a = (a1 , а„,... aO T )T —вектор
случайных параметров, тогда сигнал S(t„ а) является квазадетеїрминироваиньїм. В общем случае .процесс x*(t) может быть аппроксимирован с любой степенью точности компонентой некоторого многомерного марковского про-
цесса [48]: x(t)={x(t), |
Xi{t), ... , x ( _ i ( i ) } T . Такая'замена |
чаще всего является |
приближенной и требуемая точность |
аппроксимации достигается увеличением числа компо нент / ч. Сделав замену x*(t) =x(t), можно показать [47], что оптимальный байесовский приемник образует стати стику
|
|
|
|
|
t |
|
|
|
t |
|
|
|
|
|
In Л [и (1)) = |
j S (х) и (х) di |
— - і - J З 2 (т) dx, |
(2.9) |
|||||||||
где |
|
|
|
|
о |
|
|
|
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S(t) = |
§ S[x(t), |
t]wt(x, |
Л:,, |
|
Xi_lt |
t)dxdx1... |
|
dxt_x, |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.10) |
|
S2 (£) = |
j S2 [x (t), t] Wi (x, x1,..., |
Xi.,, |
t) dxdxt... |
dxt |
_г |
|||||||
— |
оценки |
сигналов |
S[x(t), |
>t] и 53[.v(f), |
t] |
по |
критерию |
||||||
минимума |
среднеквадратической |
ошибки |
(при гипотезе |
||||||||||
Hi); |
Wi(x, |
ХІ, ..., |
Xi-i, |
t) |
—апостериорная |
плотность рас- |
|||||||
|
|
|
|
|
|
-+ |
|
|
|
|
|
|
|
пределения вероятностей x(t), |
определяемая |
уравнением |
|||||||||||
нелинейной фильтрации [47]. |
|
|
|
|
|
|
|
||||||
|
Оценки S(t) |
и S°(t) |
зависят от наблюдаемой |
реализа |
|||||||||
ции u(t), |
которая |
содержит |
сигнал |
только |
при гипотезе |
||||||||
#_. |
При |
гипотезе |
Я 0 |
S(t) |
и S3 (f) |
суть |
псевдооценки. |
||||||
Оценка сигнала |
S{x\(t), |
t] в случае марковской |
аппрокси |
||||||||||
мации по |
минимуму |
ереднеквадратичеекого |
отклонения |
||||||||||
носит причинный характер и может |
вычисляться рекур- |
||||||||||||
|
Ч Число компонент / |
многомерного |
марковского |
процесса |
зави |
сит от порядка стохастического дифференциального уравнения систе мы, формирующей процесс x(t) из белого гауссова шума (см. п. 1.7.1.).
79