Добавил:
Закончил бакалавриат по специальности 11.03.01 Радиотехника в МИЭТе. Могу помочь с выполнением курсовых и БДЗ по проектированию приемо-передающих устройств и проектированию печатных плат. Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Литература / бакулев радиолокация распозн

.pdf
Скачиваний:
30
Добавлен:
10.09.2023
Размер:
15.09 Mб
Скачать

Перейдем к непрерывному времени: Лt О, n oo, Лt dl, тогда плот­

ности распределения вероятностей переходят в гауссовы функционалы:

W(Y/0) F(Y/0).

Для 0=0 и 0=1 эти функционалы выглядят следующим образом:

F(У/0 = О)= К,,ехр{ ½'Т'JQ(/1,t2 )y(l1)y(t2 )dt1dt2 } ,

F(Y16 =1) =k,, ехр{-½тj'17'Q(t1,12)(у(11) - u(l1) ](y(t2)- u(12 )] dlIdl2}.

Уравнение обращения корреляционной матрицы становится инте­

гральным:

1~-~

JR(t1,t)Q(t,t2 )dt = д(t1 -t2 ),

о

при этом отношение правдоподобия имеет вид

лF(Y /0 =1)

F(Y /0 =0)

=ехр{тj''1'Q(lI , 12)y(tI )u(t2 )dt1dt2

'j''1'Q(t1, 12 )u(11)u(t2 )dt1dt2 }.

о

о

о о

Введем

весовой коэффициент

фильтра обработки W(t1)=

1;.;.б.1

f Q(t1.t2 )u(t2 )dt2 и получим алгоритм обнаружения

о

7~-~

f W(t)y(t)dt ;unop .

о

Перейдем в частотную область, для чего применим преобразова­ ние Фурье к левой и правой частям уравнения фильтра обработки:

W(to -f1}'/'тQ(t1,t2)u(to -lz)dtz.

о

С использованием интеграла свертки получаем

ел

ел

1;,ali1

fW(t)exp{-jwt)dt=

f[ f Q(t1,t2 )u(t0 -t2 )dt2 ]exp{-jwt}dt=k(Jw).

 

-а,

о

71

Полагая t0 - t2 = t, по теореме о спектре свертки имеем

т1Т1'Q(t1,l2 )u(l0 -t,)dt2 ]~p{-jwt} dt =

=00JQ(t1.t2 )exp{ - jюt}dt200fu(t0 -t2 )exp{ - jюt2}dt2 =

= P(jm )ехр{ - jmt0 }S*(.im) ,

(3.18)

гдеs• (}ю)= r/)fu(t)exp(iюt)dt .

-00

Таким образом, k(jro) = P(jro)S*(jro)exp{-jrot0}.

Вычисляя интеграл Фурье от уравнения обращения корреляцион­

ной ~атрицы, находим

P(jro )G(jro)=const,

где P(jro) = 00JQ(t1,t2)exp{ - jOJt2 }dt2 ;

G(jю) = 00JR(t1,t2 )exp{ - jюt2}dt2(3.19)

Решая уравнения (3 .18) и (3 .19) совместно и исключая P(jro ), при­

ходим

к уравнению

коэффициента передачи оптимального фильтра

(устройства):

 

 

 

 

К

1

 

} cexp{jwt0 )--1.!!! [

 

. (3.20)

( ·ш)=cs•Uro)ехр{- jrot0

1

 

 

 

s*( · )

 

 

от J

GUw)

No

l + Gп(jro)

 

 

 

 

 

No

На рис. 3.22,а показана структурная схема такого обнаружителя. Он состоит из последовательно включенных «обеляющего» фильтра с коэф-

1

фициентом передачи koб(jro) G ( . ) и фильтра, оптимального для

1+_д__1!!!___

No

обнаружения сигнала на фоне «обеленной» помехи с коэффициентом пе-

редачи k0nтCiro)=cexp{-Jюt0}S*(jm). На рис. 3.22,б показано прохожде­

ние спектральных составляющих через ••обеляющий" фильтр.

72

Обнаружение

произвольного

 

~

сигнала на фоне произвольной по­

 

 

~

мехи с независимыми значениями.

 

Gвх

Gоб

Входная реализация y(ti.)=0u(tk)+~(tk) =

 

 

а)

 

 

 

=0щ+~k- Плотность распределения

 

 

 

вероятности помехи с независимы­

 

 

 

ми значениями ~k обозначим wk(~),

 

 

 

где k=l,2,3, ... ,N - номер периода по­

 

 

 

вторения.

 

 

 

 

 

 

 

 

Отношение правдоподобия име­

 

 

 

ет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.21)

 

 

б)

 

 

 

 

 

Р11с. 3.22. СтруК1}'ра обнаруж:ителя ра­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

диосигналов на

фоне коррелированной

 

 

 

 

 

 

 

по~rехи (а) и спектральная картина обеле­

а

 

 

 

 

 

 

ния помехи (б)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

2)n 1v/yk - uk) -

ln ,v, (yk )] .

 

 

 

 

 

k:I

 

 

 

Разложим lш11~(yk- щ) в ряд по степеням щ:

 

 

 

 

 

VJ

(-li

· d;

(3.22)

lnw~(J.1k-Щ) = lnw~(yk) + L~-и~-;ln w,;(Yk),

 

 

 

 

i=I

l.

 

dyk

 

 

N

а) (

ti

d;

 

 

 

 

 

тоrдаZ=LL~и~-,ln1v,:(Yk).

 

 

 

 

 

k=I

i=I

l.

dy"

 

 

 

 

 

Объединив

в

(3.22) часть

сомножителей

в функцию j,'(yk)=

(-Ii

d;

 

 

получим

 

 

 

 

 

=-.,---; ln 1v,(yk),

 

 

 

 

 

l.

dy"

 

 

 

 

 

 

 

 

 

<n

 

 

 

 

 

 

 

 

Z = LZ; ;un011 ,

 

 

 

 

 

(3.23)

 

1=1

 

 

 

 

 

 

 

 

N

гдеZ; = LJ;(yk)иt.

k=I

Алгоритм показывает, что обнаружитель - многоканальное уст­ ройство с бесконечным числом каналов (рис. 3.22,а) ,в каждом из кото­ рых стоит блок нелинейной обработки (БНО), осуществляющий нели­ нейное преобразование f,(y). Число каналов стараются уменьшить, что

можно сделать с некоторыми потерями, если отношение сигнала к по-

73

На выходе АЦП заданы

условные плотности распреде­

ления вероятностей w( Ui0)= =w(Y,,,,tl0). Вычислим вероят­

ность появления единицы на

k-й позиции при наличии

только шума:

 

 

 

00

 

 

 

Рш1= J1v(U;l0 = O)dU; .'

 

 

li,.a

 

 

 

Такая же

вероятность

при наличии сигнала:

 

 

r,:,

 

 

 

 

 

 

Реш;= Jw(И;l0=t)dU ,

Р11с. 3.24. Иллюстрация работы цифрового обна­

 

1

ружителя: а - схема; б - ква1ттование; в - вероят­

h"s

 

ности превышения порога шумом Рш и сигналом с

где 1-РшFqш;;

1-Рсш;= qcшl

шумомрсш

-

 

вероятности появления нуля на i-й позиции.

Условные вероятности принятия случайной величиной Бk любого

из двух возможных значений (О, t ), показанных на рис. 3.24, б,

При статистически независимых наблюдениях

N

 

Б 1-Б

 

P(o1,02, .. ,,0N

/0 =О)= пРш;qШI 1

'

 

i=I

 

 

N

 

P(o ,o ,...,oN /0 =1) =п ~~;q~:~,

,

1 2

 

 

 

1=1

 

поэтомуЛ

N

(

~ш~

)ii, (

qcш_i

Jl-o,

= П

 

 

,

·

i=I

 

рш1

 

qш1

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

L Б;W,;Unop '

 

 

 

(3.24)

 

1=1

 

 

 

 

 

 

 

где W; - весовой коэффициент: W, = ln рсшiqш; .

 

 

 

 

 

 

 

 

рш1qсшi

Этот алгоритм соответствует структуре весового накопителя (ин­ тегратора). Если шум стационарный: Рш; = Рш, а пачка импульсов имеет

75

прямоугольную огибающую Реш;= Реш, то W, = W = const. В этом случае

алгоритм упрощается:

N

 

 

 

 

 

 

Lб,;h'

 

 

 

 

 

(3.25)

1=1

 

 

 

 

 

 

например, пусть распределение реализаций после детектора

и

·!'"

и2

}

 

 

w(U;i0=O)= ,,-; expt- ~2

,

 

 

1v(UJ 0 = 1) = И~ ехр~-

И,2

+ ~:; } /

0

( И;~сi) ,

а-

L

2сг

 

а

тогда вероятности для одной реализации Ршi и Реш;

r/)

Рш;= Jw(U, /0 = O)dU;,

hкв

Рсщi(U/0 =1) = NJw(U, /0 = l)dU, .

hкв

Если задана вероятность Рш;, то порог квантования можно найти из

соотношения

h.,=,,-J21n 1 .

PUI

При Рш, = Рш и Реш;= Рш статистика пачки имеет биноминальное

распределение:

N

D = ~ Ck p_k (I - р )N-k

,

_L.J

N сш

сш

k=l\1

 

 

 

где ct = (N )=

N!

- биномиальный коэффициент; h0 - наиболь-

k

k!(N-k)!

 

N

шее целое число, удовлетворяющее неравенству F-5: L ctP~(I-Pш)'v-k.

k=J,0

Понятно, что одинаковые вероятности D и F можно получить при различных сочетаниях N и Рш, т.е. при разных сочетаниях порогов /1..,8 и h0 Рассмотренный метод обнаружения с накоплением бинарно­

квантованных импульсов по любым k реализациям из N соответствует

так называемому обнаружителю tипа «k из N».

76

Обычно выбирают порог h0 из

 

соотношения h~ 1,5 fii. Проигрыш

аналоговому обнаружителю не пре­

вышает 1,5 - 3 дБ. Пример схемы

бинарного обнаружителя (обнару­

житель в «скользящем» окне) пока­ зан на рис. 3.25.

Графики сигналов в различных

точках схемы приведены на

рис.3 .26,б. В режиме обзора про- Рис. 3.25. Бинарный обнаружwrель, рабо-

странства на вход поступают вид~о- тающий 8 «скользящем окне»

импульсы пачки, промодулирован-

ные по амплитуде вследствие движения диаграммы направленности ан­

тенны (т./ на рис.3.26,а). Селектор дальности (Сел R) пропускает на

квантователь (Кв) импульсы це-

 

1tAJJJ.lli-~nop

 

 

 

лп только с определенного эле­

 

 

мента разрешения по дальности.

 

,, л1ТГiъ 1:

 

Квантователь работает по алго­

 

 

ритму: если иk~hкв, то на выхо­

 

3

t-,-,-гт,-г--

 

де появляется стандартный им­

 

 

 

 

 

 

 

пульс о,= 1, а если И,<hкв, то о;=О

 

 

6)

 

(т.2 на рис. 3.26,а).

а

 

 

 

 

 

 

Таким образом, на регистр Р11с. 3.26. Схема квантоватеня (а) и графики

сдвига (РС) и реверсивный счет­ сипщ1юв (б) в точках /, 2 и 3 устройства

чик (РСч) подается последова- рис. 3.26

тельность единиц и нулей в п периодах повторения, где п - число импуль­

сов в пачке (т.3 на рис. 3.26,а). Регистр сдвига имеет число ячеек, равное п,

и управляется тактовыми импульсами (ТИ). С выхода Кв стандартные им­ пульсы попадают на вход РС и на суммирующий вход РСч. На вычитаю­

щий вход РСч импульсы подаются с выхода РС. Код числа накопленной

пачки стандартных импульсов сравнивается в пороговом устройстве (ПУ)

с порогом h0, после чего выносится решение об обнаружении цели.

3.5. Эффективность систем обнаружения сигналов

За критерий эффективности систем обнаружения при их сравнении

можно выбрать функцию f, связанную с потерями C j{С). Чаще других

используются:

1) условный r или средний r риск (лучше та система, у которой

риски 1·, r меньше);

2) кривые (.-rарактеристики) 0611аружения D = I (F,CJ) (лучше та

система, у которой кривая для F = const лежит левее);

77

Пороговой лющ11остью называют минимальную мощность сигнала

на входе приемника, при которой он обнаруживается с заданными веро­ ятностями правильного обнаружения и ложной тревоги. Обычно, если

решение принимает оператор, то берут D = 0,5, а если решение прини­

мает автомат, то D=0,9; 0,99; 0,999; .... В то же время F=J0-3 .•• 10-20.Зная

D и F, по характеристикам обнаружения определяют qnop так, как это

показано на рис. 3.28.

Таким образом, Рrюр относят ко входу приемника, а qnop - ко входу

порогового устройства или к выходу приемника. В общем случае связь

Pnop=/ (qпор) неизвестна. Однако при оптимальной обработке qвых=С/шах=Ев/N0 с учетом того, что спектральная плотность реального

шума в области положительных частот равна N0 Аналогично для опти­

мального приемника можем считать, что l/rюp=En0/No, поэтому Enop=q.ю11No.

При обнаружении одиночного импульса Enop=Pnoptи, следовательно,

pnop = qnop kшkT .

(3.26)

Если обработка не оптимальна, то это обстоятельство можно

учесть с помощью коэффициента потерь v = qmaxfq :

 

kшkTv

(3.27)

 

pnop = Чпор --- ·

 

 

 

 

Учитывая, что tи =1 /!:if, а Ршв."( = kшkТЛ/, можно представить соот­

ношение (3.27) в виде Pnop =k8Ршвх•

где k0= q 110Pv - так называемый ко­

эффициент видимости сигнала на фоне шума.

При обработке пачки когерентных радиоимпульсов Eriop=Pnoptиn

(Pnop - пороговая мощность одного импульса; п -

число импульсов в

пачке), поэтому

 

 

 

 

kшkTv

 

(3.28)

 

_

 

 

рпор - qnop --- ·

 

 

 

Ми

 

 

 

 

 

 

Выражение (3.28) можно представить с учетом сужения полосы

пропускания приемника при накоплении п импульсов:

pnop = kвРшвх,

k kT

где Рш =_ш_ =kшkT4[.

Ми

При обнаружении пачки некоrерентных радиоимпульсов накапли­ ваются огибающие радиоимпульсов после детектора, в котором возника­

ют потери. Эти потери зависят от отношения мощностей сигнала и шума

79