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

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

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

неблагоприятно распределение 1110 =w , а байесово решение ь• - мини­

0

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

зываются равномерными.

Критерий последовательной проверки гипотез Вальда. В рас­ смотренных критериях ограничивалось (фиксировалось) время приня­

тия решений Т.шбл или объем выборки у1, у2,.. ·J'k· Однако можно заранее объем выборки не фиксировать. При критерии Вальда область Г делится на три подобласти Г1, ГO и Г2 нижним Т.1 и верхним Тв порогами, как по­

казано на рис. 3.3 ,6:

б(у) = d0 , если Лкs Т.1 - в этом случае справедлива гипотеза Но и у

принадлежит области Г0 ;

·

 

 

б(у) = d1,если Лк~Тв -

в этом случае справедлива гипотеза Н1 и у

принадлежит области Г1 ;

 

 

 

б(у) = d2 , если

Т.1sЛкs Тв - в этом случае

принимается решение

продолжить наблюдение.

 

 

 

_

1v(y1...yk /0 ~ l)

 

 

здесь Лк-Л(у1,У2,, .. ,у") -~-"-----.

 

 

 

w(y1... yk /0 ~ 0)

 

 

Таким образом критерий Вальда двухпороговый:

 

lтн~Лк~тв.1

 

 

 

(3.8)

Пороги определяются

вероятностями D и

F: Тн~ l - D ,

Тв~ D .

Длительность наблюдений -

 

1-F

F

величина случайная. Критерий Вальда яв­

ляется оптимальным в смысле минимизации среднего времени наблю­

дения (обнаружения) по большему ансамблю экспериментов.

Сведение сложной гипотезы к простой. Если кроме параметров 0=1 и 0 = О, имеются другие:µ в пространстве .Q 1 для случая 0 = l и v в пространстве .Q2 для случая 0 = О с распределениями w(y/µ,l) и w(y/v,O) при известных w0(µ) и w0(v), то можно сформировать отношение прав­ доподобия, не зависящее от параметров ~L, v:

J1v(y / µ, l)w0 (p)dp

 

 

 

 

== 1-v(y / е=l)

n,

(3 .9)

Л(

 

/ 0 l

v)Jl,1'

 

у

,1-

'

w(y / 0 == О)

Jw(y / 11, O)w (v)dv .

 

 

 

 

 

 

0

 

 

 

 

 

 

П2

 

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

критериях оптимальности:

51

- однопороrовые критерии с фиксированным временем Т.,абл(п)

(рис. 3.4,а);

- двухпороrовый критерий с переменным временем Тнабл(п)

(рис. 3.4,6).

Устройство Л(у) преобразует

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

w(y/0) в распре­

деление w(Л). По­ роговые устрой­ ства (ПУ), назы­

 

 

 

 

 

 

 

ваемые

реле

или

 

 

 

 

 

 

 

компараторам и,

 

 

 

 

 

 

 

осуществляют

 

 

 

 

 

 

 

сравнение

 

 

Л с

 

6)

 

 

 

 

порогом

 

Т.

Гра­

Р11с. 3.4. Струк,урные схемы однопороговых (а) и двухпороговых фическая

интер­

(6) обнаруж11телей

 

 

 

 

 

 

претация

различ-

ных критериев приведена в таблице.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Критерий

 

 

Порог

Гщн\>Ическая интерпретация

Однопороrовый

 

q Со1 -Соо

 

wrл,tш,

 

Байеса

 

 

 

 

 

 

 

 

 

 

 

 

Р С10-С11

 

d(I

 

 

f

 

,'\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т di

 

 

 

 

Однопороrовый

 

 

 

!.!...

 

1'1Л!~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

максимума

 

 

 

р

 

[J

 

1'

 

л

 

апостериорной

 

 

 

 

 

 

 

 

 

 

 

 

 

da

!!.

 

J,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

BCf)OЯTIIOCTИ

 

 

 

 

 

 

 

р

 

 

 

 

 

 

 

Од11011ороrовый

 

 

 

1

 

W(ЛI~

 

 

 

 

максимума

 

 

 

 

 

 

 

 

../,..

- л

 

правдоподобия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d«1

1 d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

НейманаПирсона

f,v(y / O)dy =F

W(Л!t'>{;h~

 

Одно1юроrовый

U)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

dQ

 

 

r d,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Однопоро1·овый

(_<L)

с., -с••

W(Л)ro,\

 

 

МИIIИМi\КСНЫЙ

 

 

 

 

 

 

 

 

 

 

Р

мм С10-С11

 

 

 

 

 

 

 

 

 

 

 

 

<lo

T.w.,.

d1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Двухпоро1·овый

 

 

D

1-D

W(Л)оо

 

 

 

 

 

 

 

 

 

 

 

 

 

Вальда

Т=- Т=-

 

 

 

 

 

 

 

 

 

 

 

в

F

" 1-F

 

D

 

f

 

 

Л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

do

Tw

dJ ·тw

tl1

'

 

 

52

Задача обнаружения решается в каждом элементе разрешения. Способ просмотра элементов разрешения определяется выбранным ме­ тодом обзора пространства. Количество элементов разрешения зависит от величины области обзора или пространства обнаружения. Вероят­ ность ложной тревоги F и вероятность пропуска цели 1-D обыLJно за­ даются на все пространство обнаружения. Вероятность правильного не­

обнаружения во всем пространстве 1- F = F равна произведению веро­

ятностей правильного необнаружения во всех т элементах:

тt/1

F =ПF, =П() -F, ).

1=1

i=I

 

111

Если F,=~= ... =f'i=const, то F'=П(I-F,)=(1-Fi)m.

 

i=I

При F 1<<1

можно считать F=(1-Fi)m~l-mFi, или 1-F=F=

=mF1Таким образом, если задано F, то

IFi=Flm.j

Пример. Задано F = 10-3, число элементов разрешения по дально­

сти mu = 150км/150м = 1ООО, число элементов разрешения по азимуту

та=360°/0,36°= 1ООО, число элементов разрешения по углу места

тр=90°/О,9°=100. Общее число элементов разрешения т = mumam/J = 108

и Fi =1о-з 1108 = 10- 11

Что касается D, то вероятность правильного обнаружения равна

1- D, а эта последняя равна произведению вероятности пропуска цели в

одном элементе на вероятность правильного необнаружения во всех ос­ тальных т -1 элемен:rах:

D=] -[D Fi(m-l)] =] -(1- D )[1-(m- l)F ] ~ 1-(] - D ) =D ,

1

1

1

1

1

т.е.

ID=D1-I

Таким образом, D = D 1, а F= mF1.

3.3.Модели радиолокационных сигналов

Детер.11инирован11ый сиг11ш1, или сигна1 с полностью известны.wи

параметрами имеет вид

y(t) = 0и(/) + n(t), u(t) = Иш(l) cos[ro 0/ + \/f(/) - <р],

где U111(t) - амплитуда сигнала; ro 0 - несущая частота сигнала; \/f(/) - функция угловой модуляции; <р - начальная фаза колебаний.

53

Все это точно известные неслучайные величины. Считается, что мы знаем время запаздывания сигнала, эффективную площадь рассея­ ния цели, форму сигнала и все параметры его модуляции. Эта модель наиболее идеализирована.

Квазидетерминированные сигналы, или сигналы со случайными

параметрами имеют вид

y(t) = eu/1/(t,µ) + n(t)'

где µ - вектор случайных параметров сигнала.

Возможны два случая:

а) сигнет со случайной начальной фазой <р:

u(t,<p) = u,,,(t) cos[юоt + ч1(t) - ],

где <р - неизвестная начальная фаза, распределенная равномерно от О до

2n, т.е. w0(<p) = 1/(2л);

б) сигнал со случайной фазой <р и флуктуирующей аиmитудой аиш(t), где <р - неизвестная начальная фаза с распределением }Vо(<р)=1/(2л);

а - коэффициент флуктуации амплитуды с распределением

}V(a}=~exp{- а\}.

(Уа 2а-а

Таким образом, u(t,<p,a) = aU111(t) cos[юof + ЧJ(f}- <р].

Модели сигналов охватывают случаи одиночных импу,пьсов и па­ чек импульсов. Пачки импульсов разделяют на пачки когерентных и не­ когерентных импульсов. Вид пачек показан на рис. 3.5.

,,,,,- -......

а)

6)

Рис. 3.5. Форма пачек импульсов при плавном (а) и ступенчатом (б) обзоре пространства

В пачке когерентных радиоимпульсов начальные фазы импульсов

коррелированы: (<pi<pk) -:t:- О . Пачка некогерентных радиоимпульсов со-

стоит из импульсов с независимыми начальными фазами (<p,<pk) = О .

54

Формирование

когерентных радиоим­

пульсов и особенности

их спектров показаны на

рис. 3.6. Стабильный за­

дающий генератор (син­

тезатор частот) форми­ рует колебание на часто­

те /j. После умножителя

частоты (УмЧ) колеба­

ния несущей частоты

fo=n.fj попадают на усиР11с. 3.6. Формирование коrере•m{ЫХ радиоимпульсов

литель мощности (УМ), где усиливаются и модулируются импульсами. Формирование некогерентных радиоимпульсов показано на рис. 3.7.

На генератор колебаний

ра-

м

f.

 

 

 

V

Or~

диочастоты, работающии в

~

 

режиме самовозбуждения, по-

м6

П ГJ

даются модулирующие

им-

_

пульсы и формируются мощ-

f 0

1

 

 

ные радиоимпульсы частотой

/о со случайнои начальной

фазой, что объясняется слуР11с. 3.7. Форм11рование неко,·ерентных р,щионм­ чайным характером началь­ пульсов

ных условий самовозбуждения.

Что касается помех: то в дальнейшем рассматриваются модели не­

коррелированной помехи - "белый" шум, коррелированной помехи -

"пассивная" помеха и "неrауссовой" активной помехи.

3.4.Синтез оптимальных обнаружителей

Обнаружение одиночного радиоимпульса с полностью извест­ ными параметрами на фоне «белого» шума. «Белый» шум складыва­

ется с сигналом аддитивно: y(t) = 0u(t) + n(t), 0 = О, 1, при этом

M{n(t)} =О,

R(т)= M{n(t)n(t- т)}= No б(т),

2

G(ш)= } R(т)ехр{- jшт}dт= : 0 ,

o-l = R(O) = cFJfG(ш)dш =оо

-мощность шума на входе приемника.

55

В пределах полосы пропускания приемника ci=N0Лf= =КшкТ0ЛJ;

где Кш - коэффициент шума приемника (Кш=ТJТ0+ +Тпрм!Тu-1, То=290К, kT0=4· 1о-21 Вт/Гц). Для прикидочных расчетов можно ориентироваться

на следующие значения коэффициентов шума.

Тип приемника:

Кш дБ

с параметрическим усилителем

..................................................... 2 -

4

сЛБВ ................................................................................................

4-8

с туннельным диодом ......................................................................

5 -

6

с балансным смесителем .................................................................

6 -

9

Пусть наблюдение входной реализации ведется дискретно через

интервалы времени Лt = Тна6/п в моменты времени tk:{t1,/2, •.. ,t1.,...,tn}. По­ скольку п = 2.fвТ, отсчеты в моменты tk независимы: y(tk) = Yk = + nk, k = 1,2,3,4, ... ,п. В этом случае «белый» шум имеет нормальный закон распределения вероятностей (0 = О) и плотность распределения вероят­ ностей (ПРВ) имеет вид

w(yk /0 = 0) =iv(yk /0) = d-:-

ехр{- Yi,,},

 

 

 

 

 

 

-v21UY

20--

 

 

 

 

 

 

8 силу независимости отсчетов совместная ПРВ

 

 

 

 

 

 

И\У1,У2,Уз,···,У11!0=O)=w(у/O)=п,v(yk/0)=

г;:-

ехр

{

 

11

.

- ~ LYk

 

,,

1

 

l

2 }

 

 

k=I

"ШО-

 

 

20-

k=I

 

При 0 =1 отсчеты входной выборки Yk = щ+пk, поэтому щ =Yгlik и тогда

w(vtf0 = 1) = w(v/1) = г::l-- ехр{

(Yk -

Щ)

2

}

 

 

 

 

 

 

"2,ro-

 

2о-2

 

 

 

 

с учетом независимости отсчетов

 

 

 

 

 

 

 

1

 

 

 

п

 

1

ехр

{

1

 

п

 

2 }

w(i !0 =t)= Пw(yk/1)=

r---

-- L<Yk-щ)

 

 

 

 

 

 

'\J 21ro-

 

2

 

 

 

 

 

 

 

 

k=I

 

20-

k=I

 

 

Подставляя w(v! l) и w(v/0) в выражение для Л, получаем

Л=ехр

{

2

1

Lukyk --- 2LИi .

 

 

 

 

 

 

,

 

п

1

п

'}

 

 

 

 

 

 

 

 

О-

k=I

20-

k=I

 

 

 

 

 

 

 

Можно сравнивать Л с порогом решения Т, можно также сравни­ вать монотонную функцию от Л с такой же функцией порога. Напри­

мер, часто используют сравнение lnЛ~ lnT. Учитывая, что

1

11

1

п

111 л =-2 LЩYk --- 2 LИi

О-

k=I

20-

k=I

56

и обозначая а-2

1 п

получаем алгоритм оптимального

ln Т+

Lиf =ипор ,

 

2k=I

 

tl

обнаружения LllkYk ~ llnop•

k=I

Если перейти к непрерывному времени, то нужно устремить ин­ тервал дискретизации к нулю Лt О, при этом T,iaбл=const, n oo,

Nт•=~., I тнаб.,

-0 а-2 • Тогда lnЛ -2 Ju(t)y(t)dt-- f u2ct)dt.

2Лt

 

 

 

 

 

N0

о

 

N0

о

 

 

 

 

 

1;.;.б.1

 

 

 

 

 

 

 

Обозначим

f u(t)y(t)dt = =(t) и заметим, что это корреляцион­

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

1;136_,

 

 

 

 

 

 

ный интеграл, а

Ju2 (t)dt = Е -

энергия сигнала. Объединяя : 0 ln Т+

 

 

 

 

 

 

о

 

 

 

 

1

т

 

 

 

N

 

 

Е

 

 

 

 

11:16.,

 

 

 

 

 

 

 

+-

 

f и2(t)dt

=- 0 ln Т+- =ипор, получаем алгоритм обнаружения:

2

 

о

 

2

2

 

 

 

 

 

 

 

.,.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = набf.,

u(t)y(t)dt ;u110P •

 

 

 

(3.1 О)

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Структурная

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

схема корреляционно­

 

 

 

 

 

го обнаружителя (КО)

 

 

 

 

 

показана на рис. 3.8,а.

 

 

 

 

 

 

 

 

Этой

структуре

 

 

 

 

 

эквивалентна

структу­

 

 

 

 

 

ра

фильтрового обна­

 

 

 

 

а)

ружителя

 

(ФО)

 

 

 

 

z

(рис.3.8,6).

 

 

 

 

 

У(~, ОФ

 

 

 

 

 

 

 

--1

 

 

 

Импульсная

ха­

 

 

 

 

рактеристика

фильт­

 

 

 

 

 

ра,

максимизирующе­

 

 

 

 

 

го отношение мощно-

 

 

 

 

б

 

 

 

 

 

 

 

 

 

 

 

 

сти сигнала к мощно-

 

 

 

 

сти

 

шума q,

 

Р11с. 3.8. СтруКl)'ры карреляцнаннаrа (КО) (а) и фильтра­

 

является вага (ФО)

(6) абнаружrт:лей

зеркальным

 

отобра-

 

 

 

 

жением входного сигнала ГJ(t)

= и(Тс - t), поэтому выходной сигнал

имеет форму

57

т т

up1,,x (Те)= J17(~ - t)y(t)dt = fu(t - ~)y(t)dt =z .

оо

Соотношение для коэффициента передачи согласованного фильтра получаем по формуле

ш

k(jш) = f!]{t)exp{ - jcot}dt =ехр{-jшt}S*(jш),

при этом формируется максимальное отношение сигнала к шуму:

q111ax = 2EIN0 (при реШ!ьных шумах в диапазоне частот О оо q = EIN0 ).

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

онным интегралом, зависит от времени запаздывания tu и расстройки по

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

только по скорости.

В КО на выходе инерционного фильтра (ИнФ) будет нарастающее

напряжение в момент fu+t'и, равное z(т11), в то время как в ФО на выходе

возникает радиоимпульс, по форме совпадающий с корреляционной функ­

цией входного сигнала. Для устранения колебаний внутри огибающей ра­ диоимпульса на выходе ОФ ставят детектор огибающей (ДО) (рис. 3.9, а-г).

Для определения величины вероятностей правильного обнаруже­

ния D и ложной тревоги F надо знать плотности распределения вероят­ ностей величины z на входе порогового устройства при 0 = О и 0 =). Ес­

ли 0 = О, то на входе - толы<о шум n(t), поэтому y(t)=n(t). Операция ин­

тегрирования является линейной, и(t) - детерминированная величина.

Следовательно z(t) будет иметь то же распределение вероятностей, что и n(t), т.е. нормальное распределение с параметрами

M{z/0 = О}=О,

o/=M{z2}-(M{z} )2=M{z2 }= 7~абf., 1;r.абf., u(t)u(t')M {n(t)n(t')}dtdt' =

 

 

 

 

 

 

о

 

о

 

 

 

Т

,

7' ~

N

 

 

 

 

N

Т

 

N Е

нао.,

нао.,

 

 

 

 

наб.,

=

f

f u(t)u(t')-fб(t-t')dtdt'=-f

 

f

и2ct')dt=T·

 

о

о

 

 

 

 

 

 

 

о

 

Таким образом,

 

 

 

 

 

 

 

 

w(x/0 =0) =

]

ехр

{

.,.2

}

.

 

 

 

г;::-

 

- ~

 

 

 

 

 

 

v2tr6z

 

 

26=

 

 

 

 

 

58

При 0 =1 меняется лишь среднее значение из-за того, что u(t) - де­

терминированная величина:

М{z 10 =1} =М{rJu(t)[u(t) +n(t))dt}} =7'и'(l)dt=Е.

 

 

 

1

Следовательно,

 

 

 

J'(t)/ tl)'(t)U(tч ИвФ 1

W(z/0 = 1) =

 

 

=-1-ехр{ (z--E}

2

'Гша-=

~

 

U(t)

 

 

 

Отсюда можно

вы­

а)

числить искомые

вероят­

 

ности Fи D:

 

 

 

в

Рис. 3.9. Форма С11ГШ1ЛОВ в коррСЛЯЦИОIIНОМ (а).

фильтровом (6) и фильтровом с детектором оп1баю­ щей (в) обнаружителях

где h = u110/z; Ф(h)- интеграл вероятности.

Аналогично вычисляется вероятность правильного обнаружения:

D =

 

 

 

{ (

=-

Е2)2

} d==

а,

{

,,

r-y-

 

f Б

ехр

 

 

 

~ fехр

 

_У~

 

 

 

21ra-?

 

 

2а-.,.

 

"27r

 

2

 

 

11пор

-

 

 

 

-

 

_,,

 

 

 

-;,_,,

h-,Jц;;;;

ехр

{

21,1

 

 

 

 

 

 

1

 

_J;

JY=1-Ф(h-.jq.,~),

 

 

 

где 'lmax=E/No.

Обнаружение квазидетерминированного радиоимпульса. «Бе­

лый» шум аддитивно складывается с квазидетерминированным сигна­

лом, у которого начальная фаза неизвестна: u(t) = U,u(t)cos[ro 0!+\/f(l)-<p].

Здесь <р - неизвестная начальная фаза, распределенная равномерно в пре­ делах от О до 21r, с плотностью распределения вероятностей w(<p) = 1/21r.

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

59