Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации

.pdf
Скачиваний:
47
Добавлен:
25.10.2023
Размер:
14.41 Mб
Скачать

шумоподобных флуктуациях ЭОП. При совместных флуктуациях ЭОП расчет производится по формуле

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 дБ,

 

 

 

 

 

 

 

 

 

 

 

 

 

П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

Соседние файлы в папке книги из ГПНТБ