
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfчанне относится и к рис. 2.29, |
2.30. |
На рис. 2.28 и 2.29 |
|||
приведены |
зависимости |
порогового |
сигнала |
однократ |
|
ных ~К=\ |
и двукратных |
К=2 |
систем |
ЧПВ без |
ограничи |
теля в функции от ширины спектра помехи для цифро вой и аналоговой системы. Сигнал и помеха представ ляли собой нормальные процессы с гауссовой функцией
корреляции. |
Рисунок |
2.30 |
иллюстрирует зависимость |
||||||
порогового |
сигнала |
|
•от |
радиальной |
скорости цели |
||||
q\d6 |
AfpT„=0,f |
AfsTn-0,05 |
|
||||||
|
|
||||||||
|
N |
|
п=Ю |
|
т^З |
|
|||
|
Л = 0,5 |
|
|
Г=Ю'г |
f |
||||
|
ч |
|
|
С ограничителем |
|||||
|
N |
|
|
|
|
\ |
X/1 У Без огра- |
||
|
ч |
\х^ ' |
/ |
/ |
|
||||
|
|
X |
v |
\ |
|
у /тчи |
теля__ |
> |
/ |
|
|
|
\ |
\'Анало зовая
О |
Я_ |
2Х_ Я |
ЇЯ |
5Я |
(р |
|
3 |
3 |
3 |
3 |
is |
Рис. 2.30. |
Пороговые |
сигналы |
ЦЧГЩ. |
|
160
(допплеровского набега фазы за период Та). Сигнал и помеха представляли собой нормальные процессы с га уссовой функцией корреляции.
2.4.6. Цифровой эквивалент оценочно-корреляционно го обнаружителя. Возьмем за основу алгоритм (2.13), произведя дискретизацию и квантование. Тогда
|
|
|
|
|
k |
- |
|
|
|
In Л (u„u8 |
%) = |
ггіЛ(/г) = |
£ |
V y V , |
|
(2.164) |
|||
|
|
|
|
|
ix=i |
|
|
|
|
где и — квантованные |
и |
дискретизированные |
|
выборки |
|||||
процесса |
u(t); 5 |
— оценка |
сигнала |
S[x(t),t] |
|
в |
момент |
||
времени |
t=v.At |
на основании |
выборок и, |
|
= |
1,2, |
|||
At полагаем равным периоду дискретизации |
цТя. |
|
|||||||
Техническая |
реализация (2.164) |
возможна |
с |
исполь |
зованием умножителя и накапливающего сумматора (•цифрового интегратора), т. е. достаточно проста. Одна
ко реализация блока оценки наталкивается на серьезные |
|||
трудности. Дело в том, что структура блока оценки за |
|||
висит от статистических свойств |
процесса S [ x ( r ) , |
і]. Так, |
|
если S[x(t), |
t] — стационарный |
гауссов процесс, |
то опе |
ратор оценки
t
3(f) = j A ( f , * ) a ( * — ( 2 . 1 6 5 )
о
представляет собой интеграл свертки, а его импульсная
переходная |
характеристика |
h(t, т) |
определяется из ре |
|||
шения |
интегрального |
уравнения |
Винера — Хопфа |
[7]. |
||
Если |
спектр |
процесса |
S(t) |
является или может |
быть |
аппроксимирован дробно-рациональным, то оценка для дискретного времени получается на выходе фильтра Калмэна[78] [см. (3.19)].
Цифровая интерпретация соотношений (2.165) и (3.19) требует многоразрядного аналого-цифрового пре образования процесса u(t) наряду с многоразрядным представлением функции <h(t, х) и коэффициентов филь тра 'Калмэ'на. В противном случае решение уравнения
(3.19) может стать неустойчивым |
[79], а |
решение |
(2.165)—сопровождаться большими |
ошибками. |
Эти во |
просы рассматриваются в § 4.6. Если оператор оценки 5
является нелинейным, то вопросы влияния |
квантования |
П-14Ю |
16] |
и дискретизации на устойчивость и точность решения со ответствующего нелинейного дифференциального (конеч но-разностного) уравнения остаются открытыми. Многоразрядность аналого-цифрового преобразования — это первое, что затрудняет цифровую реализацию алгоритма (2.13). 'Второе затруднение связано с обобщением алго ритма (2.13) па случай произвольных помех (2.16), по скольку теперь оптимальный приемник становится двухканальным. В чі. 2.6.2 рассматривается цифровой обнару житель, имеющий более простую техническую реализа цию.
2.4.7. Обнаружитель движущегося окна при зависи
мых отсчетах. После компенсации пассивных |
помех с по |
мощью систем ЧП'К дальнейшая обработка |
некомпен |
сированных остатков обычно сводится к |
возведению |
в квадрат, накоплению и сравнению с порогом. Такое накопление наиболее просто осуществляется при бинар ном квантовании остатков. Неполная компенсация помех (частичная декорреляция) приводит к статистической связи между нулями и единицами в бинарной последо вательности. Такая последовательность может быть опи сана цепью Маркова соответствующей связности. Опти мальная обработка таких последовательностей при обна ружении рассматривается в п. 2.6.2. На практике в каче стве накопителей могут использоваться более простые устройства: цифровые программные обнаружители (см. п. 2.7.2), цифровые накопители позиций (см. п. 2.7.1).
Если бы обработка остатков вычитания осуществля лась линейным аналоговым устройством, то алгоритм его функционирования можно было представить в виде
Уi = Y, |
flfcXi.k, |
(2.166) |
k=0 |
|
|
где xi, l=\, 2, ... , — последовательность |
напряжений |
с выхода системы ЧПК; Яд — весовая функция, завися щая от корреляционных свойств «сигнальной» и «помеховой» последовательности на выходе системы ЧПК-
Наиболее просто цифровой эквивалент (2.166) реали зуется в обнаружителе движущегося окна с бинарным квантованием последовательности {xi} и весовой функции Я (см. п. 2.3.3). При анализе эффективности обнаружи теля движущегося окна необходимо учитывать статисти ческую связь между отсчетами последовательности {xi}. 162
Пусть последовательность {хі} нулей и единиц 'ПОСЛИ бинарного квантования остатков вычитания образует односвязную однородную регулярную цепь Маркова. Для определения вероятности правильного обнаружения D и частоты ложных тревог >N±(c) необходимо найти финаль ные векторы для цепей Маркова, определенных на выхо де обнаружителя при наличии и отсутствии цели. Задача построения цепи Маркова для обнаружителя движущего ся окна рассмотрена в '[80, 81]. Переходная матрица та кой цепи является якобиевой, т. е. трехдиагональной. Определив финальные векторы по (2.236), с помощью (2.21), (2.66) и (2.67) найдем D, F и М+{с).
Ограничимся рассмотрением упрощенной методики анализа, когда число импульсов в пачке эхо-сигналов достаточно большое (/г^20), а коэффициент корреляции (1.14) цепи Сг не слишком велик. >В этом случае напря жение на входе 'порогового.устройства имеет одномерное нормальное распределение [82]
|
™ЛУ)=-77= |
е х Р |
- ' " ) * > |
(2.167) |
|
|
|
||
где |
р — вероятность превышения |
порога |
квантования; |
|
о2у = |
пр (1 — р) (1 + riii+1)l{ |
I — riti+1); |
п,і+1 |
— коэффициент |
корреляции (1.14). |
|
|
|
Тогда вероятность правильного обнаружения равна
Ь= 1 - Л ( ( с - я р 5 > „ 8 р ) .
ачастота ложных тревог [см. (2.66)]
где Fl (z) — функция Лапласа [1].
2.5. Ц И Ф Р О В Ы Е Н Е П А Р А М Е Т Р И Ч Е С К И Е О Б Н А Р У Ж И Т Е Л И
2.5.1. Проблема стабилизации уровня ложных тревог. В системах автоматической обработки радиолокацион ной информации весьма важной задачей является под держание среднего числа ложных целей на уровне, не превышающем допустимую границу. В противном слу-
П* |
163 |
чае возникает перегрузка памяти и увеличивается время на выделение полезной информации в ЭЦВУ вторичной обработки, возрастает число ложных траекторий. Причи нами, вызывающими увеличение вероятности ложной тре воги при первичной обработке, являются изменения ин
тенсивности внешних помех и нестабильности |
коэффици |
ентов усиления каскадов радиолокационного |
приемника |
и пороговых устройств. |
|
Стабилизация мощности шума на выходе амплитуд ного детектора радиолокационного приемника чаще всего осуществляется с помощью шумовой автоматической ре гулировки усиления (ШАРУ) . Рассмотрим одну из воз-
упч |
АД |
|
Начальный ввод |
|
|
KB |
Реверсивный |
|
|
счетчии |
|
|
|
|
|
|
|
,"0 |
|
|
|
а(6г) |
ФПКН |
|
|
|
Рис. 2.31. Цифровая система ШАРУ.
можных реализаций цифровой |
ШАРУ |
(рис. |
2.31). |
На |
пряжение с выхода амплитудного детектора |
(АД) |
по |
||
ступает на вход квантизатора |
( K B ) , порог квантования |
|||
«о установлен так, чтобы вероятности появления 1 |
и О |
|||
были равны, т. е. /?1=/?о = 0,5. |
Каждая |
1 увеличивает, |
а каждый 0 уменьшает число, записанное в реверсивном счетчике на единицу. В реверсивный счетчик вводится
начальный код x 0 |
= 2 n _ 1 , где п — число разрядов счетчи |
ка. Требуемая |
регулировочная характеристика УПЧ |
обеспечивается с помощью функционального преобразо
вателя код — напряжение |
(ФіШШ) . |
|||
Если |
рі = р0=0,Ь, |
то в |
счетчике |
будет находиться на |
чальный |
код x = x 0 |
= 2 n _ 1 . |
При pi>po |
содержимое счетчи |
ка увеличится, а при pi<po— уменьшится. Соответствую щее управляющее напряжение с выхода ФПКН и{о2) из менит коэффициент усиления УПЧ. Такая схема произ водит измерение мощности шума путем оценки вероят-
164
Мости pi = pN, которая для закона Релея (-2.33) равна
/?д, = ехр(— и\ /2а2 ).
Тогда мощность шума
°2 = — «о |
l2LNPN- |
Для оценки вероятности PN используется ее стати стическое определение
N
PN = |
T,dilN, |
i-l
где N — число повторения опыта, пропорциональное вре мени усреднения в реверсивном счетчике; di — значение і-го измерения (0 или 1 ) . Таким образом, оценка мощ ности шума равна
^ = - " о /2 1 п Яг
Управляющее напряжение и (а2 ), пропорциональное мощно сти шума а2 , изменяет коэффициент усиления УПЧ К так, чтобы компенсировать изменение мощности шума.
Такой метод стабилизации вероятности ложной тре воги является параметрическим, т. е. зависящим от рас пределения шума и связан с прямым или косвенным из мерением необходимого параметра, в данном случае — мощности шума. Требования к схеме ШАРУ достаточно жесткие. Так, например, при бинарном накоплении пачки из 50 эхо-сигналов увеличение мощности шума всего на 0,5 дБ по сравнению с расчетным увеличивает вероят ность ложной тревоги в 10 раз [83]. Отметим, что в дан ной схеме необходимо предусмотреть бланкирование сиг налов от местных предметов, расположенных в начале развертки дальности.
Широкие возможности для стабилизации уровня лож ных тревог открывают непараметрические обнаружители,
крассмотрению которых мы и перейдем.
2.5.2.Основные определения и примеры. Обычно не параметрическим (или независящим от распределения) называется обнаружитель [84], обеспечивающий постоян ную вероятность ложной тревоги при достаточно слабых ограничениях на статистические характеристики входных данных. Чаще всего это означает, что функциональная форма закона распределения шума (помехи) может быть
165
неизвестна. Непараметрические обнаружители, как пра вило, просты в технической реализации, но по сравнению с оптимальными параметрическими имеют более высокие пороговые сигналы.
Сравнение эффективности непараметрических обна ружителей и оптимальных параметрических можно про изводить по пороговым сигналам и с помощью критерия асимптотической относительной эффективности (АОЭ) Пптмена. Пороговые сигналы при малых выборках опре
деляются |
методом Монте-Карло. При больших |
выборках |
||
и слабых |
сигналах используется |
критерий |
АОЭ [см. |
|
(2.80)]. Относительная эффективность |
обнаружителя 1 |
|||
по отношению к обнаружителю 2 равна |
|
|
||
|
Єа,і=Лі(Л D, Wi)/n2(F, |
D |
о;,), |
(2.168) |
где «і и п2 — объемы выборок для обнаружителей / и 2 соответственно, обеспечивающие одинаковые вероятности ложной тревоги F и правильного обнаружения D; Wi — распределение выборочных данных.
Асимптотическая относительная эффективность есть
Е„л = lim е„л (F„ £),, /г,, л,), |
(2.169) |
« , - > 00 /12 ->00
где Fi, Di — фиксированные значения F и D. Естествен но, что АОЭ любого обнаружителя 2 по сравнению с оп
тимальным 1 при гипотезах Нй |
и Ну не больше |
единицы. |
|||
В свою очередь, при других Н0 |
и Hi, когда |
обнаружитель |
|||
/ не оптимален, |
АОЭ обнаружителя |
2 |
по отношению |
||
к обнаружителю |
/ может превышать |
1. Это |
означает, |
||
что квазиоптнмальнын обнаружитель 2 может |
оказаться |
менее чувствительным к отклонениям от исходных рас пределений, чем оптимальный.
В качестве примеров рассмотрим два непараметриче
ских |
обнаружителя. |
|
|
|
|
|
|
1. Знаковый непарам^трический обнаружитель. |
Имеем |
||||||
выборочный вектор |
и ==(«,, |
ипу, |
|
где |
и„,... |
— по |
|
следовательность |
независимых |
случайных |
величин. |
Опре- |
|||
д л им вероятность |
р = р(«г -^>0) = |
1 —Fl(Q), |
где |
F1(z)= |
|||
|
г |
|
|
|
|
|
|
— |
j"до,(х) dx — интегральная |
функция |
распределения. |
||||
—00 |
|
|
|
|
|
|
|
Рассмотрим две непараметрические |
гипотезы: Н0 при р = |
||||||
= 7г; Ні прир>72, |
Wi(x)—произвольная |
плотность. Ины- |
166
ми словами, при гипотезе Я 0 случайные величины ut есть независимые наблюдения шума с нулевой медианой, а .при гипотезе Я] из-за аддитивного сигнала медиана больше нуля. Решение в пользу 'гипотезы Н0 ИЛИ Hi вы носится путем испытания плотности вероятности выбо рочного вектора на симметрию, т. е.
td(Ui)>c, |
(2.170) |
где d(Ui) —функция единичного скачка |
(2.60а). Посто |
янство уровня i r = const обеспечивается |
при амплитуд |
ной нестабильности радиолокационного |
тракта и не |
стационарное™ помехи.
Приведем АОЭ знакового обнаружителя относительно линейного [84], образующего 'статистику
л
i=\
При обнаружении постоянного сигнала в произвольном
шуме с симметричной |
плотностью распределения Wi(x) |
||
|
E2ii |
= |
Ao2w2 (0). |
В частности, |
для гауссова |
шума Е2, І = 0,636. |
|
2. В двухканальных системах, когда сигнал присутст |
|||
вует в двух |
каналах, |
а их шумы статистически .незави |
симы, можно использовать непараметрический обнару
житель — коррелятор совпадения |
полярностей. |
Имеем |
||||
два |
выборочных вектора |
и= («і, и2, |
ип)т |
и |
У=(УІ, |
|
vz, |
• • v n ) T . При этом |
шумы в |
каналах |
независимы и |
имеют плотности вероятности, симметричные относитель но нуля. Алгоритм коррелятора совпадения полярностей таков:
f,d(iwi)2»c. |
( 2 . 1 7 1 ) |
-> -*• |
->• |
Если и, v, а также выборочный вектор сигнала 5 — нор мально распределенные последовательности, то опти мальный обнаружитель по критерию Неймана—Пирсона образует статистику [84]
("=i
167
Асимптотическая относительная эффективность корреля тора совпадения полярностей относительно оптимального
обнаружителя |
|
равна |
£ , 2 , і = 2/л;2 ^0,202. |
Однако при не |
||||||||
гауссовых |
шумах |
может оказаться £ 2 , і > 1 |
|
|
|
|||||||
|
Алгоритмы (2.170) и (2.171) легко реализуются сред |
|||||||||||
ствами цифровой |
вычислительной |
техники, |
однако их |
|||||||||
практическое |
применение при обработке |
|
радиолокацион |
|||||||||
ной |
информации |
весьма ограничено из-за |
требований |
|||||||||
к плотностям |
вероятностей и структуре системы. Так, ве |
|||||||||||
роятность ложной тревоги F не зависит от плотностей |
||||||||||||
вероятностей |
сигнала |
и шума, но требуется симметрия и |
||||||||||
равенство |
нулю медианы распределения шума. |
|
||||||||||
|
Значительный практический интерес представляют йе- |
|||||||||||
параметрические |
фазовые, ранговые |
и знаковые обнару |
||||||||||
|
|
|
|
|
|
|
жители, |
к |
рассмотрению |
|||
aft) |
-су UgO) |
|
|
|
|
которых мы и переходим. |
||||||
Кг Ф |
|
|
2.5.3. |
Цифровой |
фазо |
|||||||
|
|
|
|
|
|
|
вый |
коррелятор. Рассмо |
||||
|
|
|
|
|
|
|
трим |
фазовый |
обнаружи |
|||
|
|
|
|
|
|
|
тель |
когерентных |
сиг- |
|||
Рис. 2.32. Фазовый |
обнаружитель |
н а л о в |
/ р и |
с 2.2,2). Узко- |
||||||||
|
когерентных |
сигналов. |
|
Ч Г |
|
„ |
' „ |
,,/ 4 \ |
||||
|
1 |
|
|
|
|
|
полосный |
процесс |
и(г) |
|||
|
|
|
|
|
|
|
[см. (В . 1)] с выхода УПЧ |
|||||
с полосой |
<А</ подвергается |
жесткому |
|
двухстороннему |
ограничению, а затем осуществляется когерентная филь трация и сравнение выходного напряжения с порогом. Не рассматривая вопросы технической реализации фазового
обнаружителя, отметим некоторые |
его свойства [85—87]. |
|||
Если |
выполняются условия |
|
|
|
|
|
?2 <0,5, Л/»1/Г, |
(2.172) |
|
где q2— |
отношение |
сигнал/шум на |
входе |
ограничителя, |
а Т — длительность |
сигнала, то фазовый |
обнаружитель |
имеет свойства квазилинейной системы и его пороговые сигналы не намного (на 1—2 дБ) хуже, чем у оптималь ного когерентного обнаружителя. В то же время фазо вый обнаружитель является непараметрическим по отно шению к широкому классу плотностей шума Wi(u), кото рые должны обладать свойством симметрии [87]. Кроме того, при наличии перекрывающихся сигналов от не скольких целей все они разрешаются и не влияют друг на друга. Если лее условия (2.172) ие выполняются, то потери в пороговом сигнале достигают 6 дБ, а обнару житель теряет свойства квазилинейности, т. е. слабые
сигналы частично подавляются сильными и могут возни кать ложные цели {'88, 89].
При цифровой реализации фазового обнаружителя не обходимо перейти к квадратурным составляющим про цесса u(t) (В.'2). Для определения структуры дальней шей обработки (накопителя) запишем выражение для комплексной огибающей на выходе когерентного филь тра [43 ]
|
|
|
00 |
|
|
|
|
v(t) = -T |
§H*(t-z)ii(i)dx, |
|
(2.173) |
||||
|
|
|
—оо |
|
|
|
|
где H*[t)—комплексно-сопряженная |
огибающая |
им |
|||||
пульсной переходной |
характеристики |
фильтра; |
u(t) |
— |
|||
комплексная |
огибающая колебания и(t) |
(В . 5) . . |
|
|
|||
Перейдем от описания мгновенными значениями не |
|||||||
прерывных |
функций |
к описанию совокупности |
кванто |
||||
ванных значений в несколько последовательных |
момен |
||||||
тов времени. Тогда |
вместо u(t) |
и H*(t) |
будем иметь |
||||
|
|
и = (и1,и2, |
...,ип)т, |
|
(2.174) |
||
|
И* |
= |
{Н*1 ,Й*„...,Й*пу. |
|
(2.175) |
Здесь компоненты вектора — последовательные во вре мени значения, а знак « ~ » означает квантование.
Ограничившись интервалом, на котором подынтег ральное выражение отлично от нуля, и перейдя к сумми рованию с шагом Г д , получим (опуская множитель 1/2)
л-1 |
- * -* |
|
о(*) = Е Й\ІІ^Тл |
= Й**й'Тл, |
(2.176) |
где
I 0 при /г<ц;
ТА — период Дискретизации; Я * т — вектор, комплексно-со
пряженный и транспонированный к Н.
169