
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfшумоподобных флуктуациях ЭОП. При совместных флуктуациях ЭОП расчет производится по формуле
Dc^D(a0)w;(a0)daa, |
(2.65) |
|
|
6 |
|
где D(ao)—условная |
вероятность |
правильного обнару |
жения нефлуктуирующего сигнала, принявшего данное значение амплитуды а0 ; Wi(a0)—плотность вероятности флуктуации амплитуд на выходе квадратичного детек тора, определяемая из (2.34) с подстановкой o" = as и
При неизвестном времени прихода пачки эхо-сигналов вероятность правильного обнаружения также опреде ляется по формулам (2.62) и і(2.65). Для определения вероятности ложной тревоги FT (2.24) на интервале на блюдения Т„ (ожидания прихода пачки эхо-импульсов) необходимо определить частость ложных тревог N+(c). Эту задачу можно решить, по крайней мере, двумя спо собами. Во-первых, в соответствии с формулой ('2.26)
N+ (с) = F (с)/т (с) « F {с)1Тл = F (с) |
(2.66) |
поскольку .при высоком лороге с, т. е. при малых значе ниях вероятности ложной тревоги, средняя длительность выброса обнаружителя движущегося окна над порогом х(с) с большой вероятностью не превышает шага ди скретности по оси времени, т. е. Т„.
Во-вторых, используя соотношение (1.60),
|
N+(c)=AfxBM, |
(2.67) |
где хнм=Тн+тм, |
а величины хм |
и тн определяются в со |
ответствии с (1.55) или (1.59). |
Соотношение'(2.67)—точ |
|
ное, а (2.66) —приближенное, |
но зато более простое. |
|
2.3.4. Весовой бинарный обнаружитель некогерентных |
||
сигналов. Структура весового |
бинарного обнаружителя |
соответствует алгоритму (2.58) при бинарном квантова
нии элементов последовательности |
v.j, / = 1, 2 , . . . и от |
|||
сутствии квантования |
(или при многоуровневом |
кванто |
||
вании) весовой функции |
G2 |
|
|
|
yvJ^TPfi^ |
с; |
6 = 0 , 1 , 2 |
r=l7M |
(2.68) |
110
Получение аналитических |
соотношений |
для |
/ г = Я ( с ) и |
||||||||||
D=D(c, |
q2) |
связано |
с интегрированием |
характеристиче |
|||||||||
ской функции, которая имеет вид [58] |
|
|
|
|
|||||||||
|
|
|
|
|
|
л/2 |
|
1VQ2 |
|
|
|
|
|
|
|
|
Ш = |
І [ ( |
Р * |
'+іГ> |
|
|
|
(2-69) |
|||
|
|
|
|
|
|
/=1 |
|
|
|
|
|
|
|
где <7=1—р; |
р — вероятность |
превышения порога кван |
|||||||||||
тования, |
определяемая |
|
для |
одного |
шума |
из |
(2.64), |
||||||
(2.49), а для |
сигнала |
с |
шумом — и з |
(2.47) и |
(2.43); |
v — |
|||||||
произвольная |
вещественная .переменная. |
|
|
|
|||||||||
Кумулянты распределения, соответствующего харак |
|||||||||||||
теристической функции (2.69), определяются по |
форму |
||||||||||||
ле [1] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N s ^ w l |
|
' |
|
|
( 2 - 7 0 ) |
||
и с учетом (2.69) |
равны |
|
|
|
|
|
|
|
|
||||
|
|
|
|
л/2 |
|
|
|
л/2 |
|
|
|
|
|
|
|
|
|
л/2 |
|
|
|
|
|
|
л/2 |
|
|
Ks = 2pq(q |
— p) |
£ |
Gj, |
Х 4 = |
2pq (<?2-j- p 2 — dpq) £ |
G®, |
|||||||
|
|
|
|
M i |
|
|
|
|
|
|
/=i |
|
|
|
x e = |
2 W ( 9 |
- p ) ( ^ |
+ p * - 1 0 W ) 2 |
О ) 0 . |
|
(2.71) |
||||||
|
|
|
|
|
|
|
|
|
/ і |
|
|
|
|
Определяя |
p |
для смеси |
сигнала с |
шумом, |
следует |
||||||||
учитывать, что pj=p(Gj), |
|
т. е. зависят от диаграммы |
на |
||||||||||
правленности |
антенны. |
Поэтому |
приближенно |
можно |
|||||||||
считать, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
л |
|
|
|
|
|
|
|
|
|
^ « ( 1 / « ) Е Ps»t- |
|
|
|
|
|||||
|
|
|
|
|
|
|
1=1 |
|
|
|
|
|
|
Для |
расчета зависимостей |
F = i f (с) |
и D=D(c, |
q2) |
при |
известных кумулянтах (2.71) полезно использовать ряд
Эджворта |
(см. приложение |
1, алгоритм 3). Зависимость |
|
D=D(c, |
qz) |
можно определить методом Монте-Карло на |
|
Э Ц В М |
(см. приложение 1) |
. В случае совместных флук |
|
туации ЭОП цели вероятность D определяется в соответ |
|||
ствии с |
(2.65). |
|
Ш
Техническая |
реализация |
алгоритма (2.68) |
возможна |
|
с использованием |
регистра |
сдвига, |
регистров, |
хранящих |
весовые коэффициенты G2, |
/ = 1 , а, |
сумматора |
и цифро |
|
вого компаратора. |
|
|
|
2.3.5. Оптимизация обнаружителей бинарно-квантован ных сигналов. Эффективность цифровых обнаружителей.
Структурная схема |
одноканального |
обнаружителя би |
|
нарно-квантованных |
сигналов (рис. 2.13) содержит одно |
||
уровневый аналого-цифровой преобразователь |
напряже |
||
н и е — код (квантизатор), накопитель |
(весовой |
или рав- |
|
u(t) |
Нанопитель |
ПУ |
|
|
|
1 7
Рис. 2.13. Одпоканальиый обнаружитель бинарно-квантованных сиг налов.
новесный), и пороговое (решающее) устройство (цифро вой компаратор). Порог квантования щ называется ана логовым, а порог решающего устройства с — цифровым. Вероятности D(u0, с, qz) и Е{ий, с) суть функции обоих порогов. Возникает задача выбора наилучшей в опреде
ленном |
смысле комбинаци |
порогов |
«о и с. Эта |
задача |
в силу |
своей актуальности |
ставилась |
и решалась |
рядом |
исследователей с различных позиций преимущественно для одноканального равновесного обнаружителя («/г из п») и нефлуктуирующего сигнала.
Различные критерии выбора порогов и0 и с можно разбить на 4 группы: 1) информационные критерии, ми нимизирующие потери информации при квантовании, 2) критерий, максимизирующий отношение сигнал/шум при квантовании, 3) критерий относительной асимптоти ческой эффективности, 4) критерий, минимизирующий средний риск либо пороговый сигнал при обнаружении.
Перейдем к характеристике названных критериев.
1. Информационные критерии связаны с максимиза цией количества информации или ее верхней грани в би нарной случайной величине 4=1, 0 относительно нали чия на входе квантизатора смеси сигнала с шумом или одного шума. При этом используются различные инфор мационные меры. Информационная мера Шеннона (1.124) в нашей задаче определяет количество информа-
112
ции в бинарной случайной величине d относительно слу чайной величины 6 — дискретного случайного парамет ра (1.127), при этом
d |
_ { |
1 с |
вероятностью |
J 0, |
и0 ), |
(2 |
72) |
|
\0 с |
вероятностью р(0|8, |
и0). |
|
|
||
Величина |
/(0, |
J w0 ) зависит от и0. |
Оптимальное |
зна |
чение и0 соответствует максимальному количеству инфор мации, содержащемуся в d относительно 9, а значит, и минимуму среднего риска при квантовании [см. п. 1.7.4. (1.129)]. С учетом соотношений (1.124), (2.72) и (1.133) необходимо решить уравнение
У(8, d ) = m a x I ^ |
^ |
|
р (6<) р {d5 | б,-, и0) X |
|||||
|
|
\ 1=0 1=0 |
|
|
|
|
||
X In , |
Р ^ ^ а Л |
] . |
|
(2.73) |
||||
|
Ё |
P$i)p(dj\h. |
и.) |
|
|
|||
|
;=о |
|
|
|
|
|
|
|
Решение (2.73) для различных распределений |
Wi(v^) и |
|||||||
ш'і'(овя) дает для р = <7 = 0,5 |
окончательный |
результат |
||||||
Wi{vN |
= u0Opt) |
=wi(vSN=Uoopt), |
|
(2.74) |
||||
что эквивалентно |
|
|
|
|
|
|
|
|
оо |
|
|
|
|
u„ |
Л |
|
|
! |
|
|
|
|
|
|
|
|
j |
о», (»Л ) Лі Л |
+ |
J « » i ( o 5 A , ) d o s w J . |
(2.75) |
||||
Аналогичный |
результат |
получается при |
использова |
нии критерия минимума среднего риска при квантовании [59]. К сожалению, такой подход не позволяет минимизи ровать средний риск решения.
Поскольку максимизация информационной меры (1.124) эквивалентна минимизации среднего риска реше
ния, в соответствии |
с (1.124) |
запишем уравнение |
|
/(8, " ) = m a x j £ £ p |
(Є*) р (щ \ 6,, и.) X |
|
|
|
JfaO /=1 |
|
|
Х 1 п |
JP ("J.I Bt. До) |
(2.76) |
|
|
|
||
£ |
Р (9») jo («і I 9i. "о) |
|
8—1410 |
113 |
где |
Р(Щ\ЬІ, |
иа) = р'фі. «о)9п "г '(9г. "oh |
(2-77) |
|
|
Р$)=рЪ(Ъь |
«o )<7n -'J (0b |
«о), |
|
(2.78) |
|
а |
р(0і, |
ы0 ) = 1—q(Qi, |
«о) |
определяются |
по |
формулам |
|
(2.48), |
(2.64) при 0о |
и в |
соответствии |
с |
(2.47), (2.43) |
||
при 0і. |
|
|
|
|
|
|
|
|
Решение уравнения (2.76) с учетом (2.77), (2.78), так |
||||||
же |
как |
и более простого |
уравнения |
(2.73), |
возможно |
численными методами. Отметим, попутно, что в задачах' радиолокационного обнаружения байесов подход, как правило, не используется из-за отсутствия априорных ве роятностей р и q. В работах [60, 61] при выборе опти мального порога и0 использовались информационные ме ры Вальда и Фишера.
2. Критерий максимума отношения сигнал/шум при
квантовании |
[62] |
|
|
|
<L = |
V*(PSN |
- PNWPN^-PN) |
(2-79) |
|
используется |
для |
биномиальных распределений |
сигнала |
|
с шумом и |
шума после |
квантования (2.62), (2.63) при |
3. Критерий асимптотической относительной эффек тивности Питмена (АОЭ) дает возможность установить оптимальные дороги и0 и с в бинарном обнаружителе, а также оценить предельную (асимптотическую) эффек тивность этого обнаружителя при п—>-со. Он основан на следующем определении эффективности выбранного кри терия обнаружения:
|
|
|
|
ЧГ"1* |
1*4.1 |
, |
(2.80). |
|
|
|
|
, |
|
||
|
|
П-*оо " |
|
|
|
||
где |
т1 |
[vSN \—среднее |
значение |
суммы |
сигнала и |
шума; |
|
|
|
п |
|
|
|
|
|
a 2 [ f w |
] — дисперсия |
шума; а — амплитуда сигнала. |
|
||||
|
71 |
|
|
|
|
|
|
Асимптотическая |
относительная эффективность |
одно |
|||||
го обнаружителя относительно другого |
определяется как |
||||||
|
|
Я « = |
е_Ю/е,(«>п). |
|
(2-81) |
||
При сравнении обнаружителей по критерию (2.81) |
пола |
||||||
гают, что в обоих случаях одинаковы величины |
Ft—Fz, |
||||||
Di=Dz, |
a! = fl2. В качестве эталона берется обнаружитель |
||||||
114 |
|
|
|
|
|
|
|
с числом уровней квантования г\—>-оо (2.61). Известно, что его эффективность как при наличии, так и при от сутствии флуктуации ЭОП равна ez(vn) = \.
Альтернативным методом в случае близких гипотез является следующий. В соответствии с выражением (2.206) для скалярного параметра X и дискретного рас пределения Р(и) имеем
где
О in Р(и | К)
Х=1о
— информация Фишера, т. е. информация, содержащая-
ся в выборке и относительно параметра К, определенная в точке ЧК=%о.
Характеристики обнаружения в данном случае пол ностью определяются отношением сигнал/шум на входе
порогового |
устройства |
|
|
|
|
|
|
||
|
|
_/», [Z1] — m1[Z0 |
|
|
|
|
|||
|
|
|
|
t |
|
|
|
|
|
Максимизация qa |
эквивалентна |
максимизации |
функ |
||||||
ции Tn{h) |
по порогу квантования «о: |
|
|
|
|||||
Л Л . |
«оорі) = п'ах/я (Ав 1 |
и0)=ш-Іп(Ій, |
"о)=° - |
||||||
|
|
|
«о |
|
0 |
U ° |
|
|
|
При оптимизации обнаружителя (2.60) необходимо |
|||||||||
максимизировать |
|
ei(vn) |
по и0, |
а зная «оopt и асимптоти |
|||||
ческое выражение |
для F, найти с = /г0. Решение |
этой за |
|||||||
дачи дает следующие окончательные результаты [63]: |
|||||||||
|
«о opt = 1,785 a, |
PN opt = 0,2035, |
|
|
|||||
|
lim |
|
= |
0,2035, £ S 1 |
= 0,648. |
(2.82) |
|||
|
n-»oo |
|
|
|
|
|
|
|
|
Отношение |
сигнал/шум |
в обнаружителе |
бинарно-кванто |
||||||
ванных сигналов |
необходимо |
увеличить |
в 1/ У0,648 |
раз |
|||||
при заданных Д |
F, п. Это эквивалентно потерям |
в поро |
|||||||
говом сигнале 0,944 дБ, |
|
|
|
|
|
|
|||
8» |
|
|
|
|
|
|
|
П5 |
Отметим теперь некоторые особенности рассмотрен ных критериев. Критерии, минимизирующие потери ин формации или риска при квантовании, не связаны с кри терием минимума риска при вынесении решений. Урав нение, минимизирующее риск решения (2.76), можно ре шить лишь численными методами. Критерий асимптоти ческой относительной эффективности, равно как и ин формационные критерии (1.139), (1.140), а также крите рий максимального отношения сигнал/шум справедливы лишь при больших выборках: п—»-оо.
Поэтому более целесообразным па практике для ко нечных выборок является подход, основанный на мини мизации потерь квантования при решении, т. е. при об
наружении, |
реализуемый численными методами [64, 65]. |
|
4. Численный |
метод оптимизации. Рассмотрим обна |
|
ружитель |
(2.60), |
эффективность которого определяется |
в соответствии с (2.02), (2.63). Выбор порогов квантова ния и0 и решения с возможен из решения двух уравне ний:
D = max{£>(u0 , |
с, |
q2)} |
при ?a = |
<£i n , |
|
|
||
|
"о |
|
|
|
|
|
|
(2.83) |
|
|
|
|
|
|
|
|
|
|
F(u0, |
с) — const. |
|
|
|
|
||
Ниже |
приводятся |
результаты |
численных |
решений |
||||
уравнений |
(2.83), выполненных |
при |
F= |
10~ 2 н - Ю - 8 |
п= |
|||
= 5-^-50, при аппроксимации диаграммы |
направленности |
|||||||
функцией |
(2.41). |
|
|
|
|
|
|
|
Вычисления производились по следующей методике. |
||||||||
Вначале определялись семейства зависимостей F'(iio, |
с) = |
|||||||
= const для различных |
« 0 |
и с и трех |
названных |
обнару |
жителей. При этом для обнаружителя (2.60) использопались таблицы кумулятивного биномиального распре
деления |
[66], |
для |
|
обнаружителя |
(2.68) —кумулянты |
||||
(2.71) и ряд |
Эджворта |
(см. приложение 1, алгоритм 3), |
|||||||
для |
обнаружителя |
«k/k—І» — метод |
Монте-Карло. Да |
||||||
лее, по порогам «о и с, |
соответствующим заданной веро |
||||||||
ятности F, методом Монте-Карло определялись вероят |
|||||||||
ности |
D—D(UOJ |
С, qz), |
при этом выбирались те значения |
||||||
ы0 и с, для которых |
выполнялись условия (2.'83). |
||||||||
Как показали результаты численного анализа (64, 65], |
|||||||||
оптимальный |
цифровой |
порог |
с определяется |
числом Л |
|||||
и типом |
флуктуации |
ЭОП. в |
свою |
очередь, |
оптималь |
ный порог квантования «о определяется заданной ве«
т
|
|
|
|
|
Т А Б Л И Ц А |
2.1 |
||
|
|
|
Тип бинарного |
обнаружителя |
|
|
||
Tim флуктузцні! ЭОП |
|
,fe из п" |
|
.klk-l' |
|
|||
|
|
|
весовой |
|
|
|
||
Отсутствие |
флук |
|
"о.5 |
|
|
|
||
туации |
|
|
|
|
2 |
|
|
|
Шумоподобные |
|
1,5 V n 0 l 5 |
|
2 |
|
|||
Совместные |
|
|
|
"о. 5 |
|
1.5 |
|
|
|
|
|
2 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
(=1 |
|
|
|
|
|
П р и м е ч а н и е . Оптимальное значение |
с = fe для |
обнаружителя |
||||||
„А/А — /" |
при |
шумоподобных флуктуациях |
равно с = 1 . |
Однако |
с |
|||
целью защиты от несинхронных импульсных |
помех |
целесообразно |
||||||
взять с = |
2. |
|
|
|
|
|
|
|
роятностью |
F |
и числом /г и не зависит |
от q2. |
В |
табл. 2.1 |
приведены оптимальные значения порога с для трех названных выше бинарных обнаружителей. Зная опти*
мальное значение |
с при |
заданной |
вероятности |
F, |
||
можно определить |
величину |
и0 из уравнения F(uo, |
с) = |
|||
=F. |
Такая зависимость для |
обнаружителя (2.60) |
задана |
|||
формулой (! 2.63), а |
для обнаружителя |
( 2 . 6 8 ) — в |
виде |
|||
ряда |
Эджворта (см. приложение 1, алгоритм 3) и куму-* |
лянтов распределения (2.71). С помощью соотношений
(2.48), (2.49), (2.60) определяется |
величина. pN 0 pt, обее* |
|
печивающая выполнение условий |
(2.83). |
|
На рис. 2.14 приведены зависимости Pivopt(-F), вычис |
||
ленные в соответствии с алгоритмом |
(2.60), для трех ти |
|
пов флуктуации ЭОП (я = 6 , ГО, |
30, |
50). Оптимальное |
значение с выбиралось в соответствии с табл. 2.1, а соот ветствующее PN — из формулы (2.63). Величина поро га квантования и0 определялась из (2.64), а вероятность
PSN — W3 формул (2.47) |
и (2.43). В точках, обозначенных |
||
на графике кружочками, приведены пороговые |
значения |
||
qz при £) = |
0,5. Как следует из графика, в общем случае |
||
Pivopt отличается от |
= 0,2036 как для слабых |
(—2 д Б ) , |
|
так и сильных сигналов ( + 2,5 дБ) и только в |
частных |
||
случаях pN |
opt:M),2Q35, |
|
|
II?
Для оценки эффективности цифровых обнаружителей необходимо построить характеристики обнаружения. Ме тодики расчета характеристик обнаружения для бинар
ных |
равновесного |
(2.60) и весового |
(2.68) обнаружите |
|
лей |
рассмотрены |
в пп. 2.3.2 и '2.3.3, |
а вопросы |
оптими |
зации порогов щ |
и с — в п. 2.3.4. Проанализируем эф |
|||
фективность весового аналогового |
накопителя |
(2.57). |
||
Для этого необходимо получить зависимости F(c) |
и D(c, |
|||
q2). |
Вначале найдем плотность вероятности для взвешен |
ной суммы шумовых |
напряжений |
(2.57). Характеристи |
|
ческая функция для |
амплитуды |
шумового |
напряжения |
на выходе квадратичного детектора имеет |
вид (1] |
Є І ( » ) = ( l—2fo(T2 )
Для напряжения шума, умноженного на весовой коэффи циент G2h, характеристическая функция равна [67]
Ьк (v) = 6 , (G[v) = (l - 2 W G V ) - I .
Рассмотрим случай симметричной весовой' функции, что всегда выполняется на практике, и будем считать, что число импульсов в пачке четно, а импульсы попарно симметричны относительно максимума диаграммы на правленности. Последнее ограничение не является стрі тим и не снижает общности получаемых результатов.
118
Обозначим git — (З2/, а2 , тогда характеристическая функция напряжения на выходе весового сумматора запишется в виде
6(о) = П ( 1 - 2 ш £ э - ) - 2 . |
(2.84) |
а искомая плотность вероятности в виде
со
Перейдем к интегрированию на комплексной плоско сти, сделав замену переменных — i v = S :
|
|
і со |
|
|
|
|
|
|
|
л/2 |
|
|
1 |
Г |
|
е' |
|
|
|
|
|
resf t . |
|
|
2пг |
|
п/2 |
|
|
|
|
|
|
|
|
|
- / с о |
Д ( l + 2 S f t ) » |
|
|
* = і |
|
|||||
Подынтегральная |
|
функция |
|
имеет |
полюс |
кратности |
2 в |
||||
точке |
Sj = — l / 2 g j . Применив |
формулу |
вычисления |
вы |
|||||||
чета |
относительно |
/г-кратного |
полюса |
|
|
||||||
|
|
1 |
|
d n _ 1 |
|
|
|
|
|
|
|
|
r e s f c = ( д _ |
i)[ |
|
|
К |
( 5 ) |
( S |
— S * ) l f a s f c « |
|
||
получим |
|
п , - |
|
|
л/2 |
|
|
|
|
||
|
|
|
|
|
і |
e - |
w ^ x |
|
|||
|
a). ((/)== |
— ^ r e s f |
e |
= ^ |
|
||||||
|
2 ^ |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
*=1 |
|
|
|
|
|
|
|
|
|
X 2 |
^ |
|
|
|
- |
- |
|
|
|
|
|
|
/ є * |
g h |
П |
C - ^ ) 3 |
|
|||||
|
|
|
|
ye |
|
|
* |
|
|
(2.85) |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
2 e* |
П |
|
( |
J - |
^ ) |
2 |
|
|
Здесь множества M, Mh, |
Mh} |
определяются следую |
|||||||||
щим образом: |
|
|
|
|
|
|
|
|
|
|
|
|
М = 1, 2, ... , л/2, Mk=M—{k}, |
|
A l f c i = M f c — { / } , |
|
119