Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Авиационные радиолокационные системы. Часть 1.pdf
Скачиваний:
59
Добавлен:
05.02.2023
Размер:
2.02 Mб
Скачать

 

 

 

42

 

 

 

 

Отсюда при F ' = F

 

 

 

 

 

 

 

opt

 

 

 

 

 

D

³ D' .

 

 

 

 

(2.10)

opt

 

 

 

 

 

 

 

Соотношение (2.10) остаётся в силе и в

случае, если F ' £ F

 

. Это

 

 

 

 

 

opt

 

означает,

что

оптимальная

система

обеспечивает

 

наибольшую

вероятность

правильного

обнаружения(наименьшую

вероятность

пропуска) из всех систем, у которых вероятность ложной тревоги не больше, чем у оптимальной системы. Указанные условия составляют

содержание критерия

Неймана-Пирсона. Как и критерий идеального

наблюдателя, он является следствием критерия минимума среднего риска.

Достоинством его

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

вероятностей наличия или отсутствия цели(сигнала), что существенно при решении задач радиолокационного обнаружения.

2.2Оптимальное обнаружение детерминированных сигналов

2.2.1

Постановка задачи

 

 

 

 

 

 

 

 

На интервале [0,T ]

наблюдается аддитивная смесьr(t) сигнала

и

шума (1.15). Сигнал представляет детерминированную функцию времени и

известных

параметров (1.8). Статистические

свойства

помехиn(t) ,

представляющей

гауссовский

белый , шумсчитаются

известными

и

определяются

соотношениями (1.16).

Задача

правильного

обнаружения

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

к нахождению

оценки q* = q* (r(t)) ), обеспечивающего минимум

среднего

риска (или

максимум

весового

критерия). Алгоритм

 

оптимального

обнаружения должен определить математические операции, по которым для

каждой принятой реализацииr(t) можно найти

наиболее

рациональный

ответ о наличии или отсутствии сигнала. Далее необходимо определить

качественные

показатели обнаружителя, т.е. условные вероятности D и F ,

и выяснить путь реализации оптимального обнаружителя.

 

 

 

 

2.2.2

Методика решения задачи обнаружения

 

 

 

 

 

Для

решения

задачи

оптимального

 

обнаружения

воспользуемся

критерием

(2.9).

Чтобы

 

определить

оптимальное

решающее

правило

q* = q* (r(t)) ,

необходимо записать развёрнутое выражение вероятностей D

и F . Эти вероятности характеризуются условными плотностями вероятности

реализации

r(t)

при

 

наличии

и

отсутствии

сигнала

в.

Чтобысмеси

определить

 

плотности

вероятности, рассмотрим

простейший

метод

дискретного наблюдения, при котором берутся отсчёты функцииr(t) через равные отрезки времени Dt . При этом интервал наблюдения разбивается на

43

m равных элементарных интервалов, где m = T(Dt) . Введём осреднённые

за элементарный интервал Dt

значения колебания r(t) , сигнала s(t) и шума

n(t) :

 

 

 

 

 

 

 

 

 

 

 

 

1

 

ti

 

 

1

 

ti

ri

=

 

ò

 

r(t)dt; si

=

 

ò s(t)dt ;

Dt

Dt

 

 

 

t -Dt

 

 

t

-Dt

 

 

 

 

i

 

 

 

i

 

 

 

 

 

1

 

ti

 

 

 

 

ni

(t) =

 

ò n(t)dt .

 

(2.11)

Dt

 

 

 

 

 

t -Dt

 

 

 

 

Очевидно, что

 

 

i

 

 

 

 

 

 

 

 

 

 

 

ni

= ri - si .

 

 

 

 

(2.12)

Таким образом, рассматриваются системы случайных величин, которые полностью характеризуются многомерными( m -мерными) плотностями вероятности. В дальнейшем нас будет интересовать предельный случай, когда

Dt ® 0 . Определим совместную плотность вероятности для случайных

величин ni ,

i =1,..., m. Случайные

величины ni с

учётом условий(1.26) и

(2.11)

являются

гауссовскими

и

имеют

характеристики: < n >= 0;

 

 

 

 

 

 

 

 

 

i

s 2 =< n2 >= N

o

(2Dt); < n n

j

>= 0 при i ¹ j . Тогда совместная плотность

ni

i

 

i

 

 

 

 

вероятности на основании независимости случайных величин ni , входящих в

совокупность, примет вид

m

w(n1 , n2 ,..., nm ) = w(n) = Õw(ni ) =

i=1

æ p N ö

-m 2

æ

 

1

m

ö

 

 

 

 

= ç

o

÷

 

exp ç

-

 

åni2 Dt ÷ .

(2.13)

Dt

 

No

è

ø

 

è

 

i=1

ø

 

Здесь и в дальнейшем для сокращения записей векторомобозначены

 

значений n = {n1 , n2 ,..., nm }. Знаком

m

 

 

совокупности

Õ

 

обозначено

 

 

 

 

 

i=1

 

 

произведение, в

данном

случае– одномерных

плотностей

вероятности

независимых случайных величин ni .

 

 

 

 

 

При отсутствии сигнала в смесиri

= ni , поэтому

условная

 

плотность

вероятности совокупности значений реализацииr

совпадает с

плотностью

вероятности помехи (2.13), куда вместо ni

следует поставить ri .

 

 

При условии действия суммы сигнала и шума, подставив значения ni из

выражения (2.12), можно

записать

условную

плотность

вероятности

совокупности значений реализации r в следующем виде:

 

 

 

wс.п (r) = wп (r - s) =

44

-m 2

= æç p No ö÷ è Dt ø

æ

 

1

m

ö

 

exp ç

-

å(ri

- si )2 Dt ÷ .

(2.14)

No

è

 

i=1

ø

 

При рассмотрении реализацииr в виде совокупности случайных величин решающее правило представляет функцию этой системы случайных

величин

q* (r) = q* (r , r ,..., r ) . По

условию, решающая функция может

 

 

1

2

m

 

D определяется при

принимать

только два

значения: 0 и 1. Вероятность

действии

в

смесиr(t)

и

сигнала, и

шума, т.е. плотностью вероятности

wс.п. (r).

В

случае

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

решающей

*

функцииq (r) условная

вероятность D принятия решения q* =1 определяется выражением

¥

D = ò q* (r)wс.п. (r)dr .

Аналогично при действии одной помехи и произвольной решающей функции q* (r) условная вероятностьF принятия решенияq* =1 определяется выражением

¥

F = ò q* (r)wп (r)dr .

Используя приведённые выражения, запишем взвешенную разность, входящую в формулу (2.9):

 

¥

 

 

 

 

 

D - lo F = ò q* (r)wп (r)(l(r) - lo )dx .

 

 

(2.15)

 

 

 

 

 

 

Здесь введено обозначение

 

 

 

 

 

l(r) = wс.п. (r) wп (r) ,

 

 

 

 

(2.16)

которое называется соотношением правдоподобия.

 

 

Согласно

критерию (2.9), оптимальной

считается

та система, которая

обеспечивает

наибольшее

значение

взвешенной

разности(2.15). Поэтому

необходимо

отыскать

такую

решающую

 

*

функциюq (r) , которая

максимизирует

подынтегральное выражение

при

всех

сочетаниях значений

r1 , r2 ,..., rm . Можно убедиться в том, что оптимальная решающая функция имеет вид

qopt*

ì1,

если l(r) > l

;

(r) = í

o

(2.17)

 

î0, если l(r) < lo ,

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

о наличии цели(сигнала), если отношение

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

превышает порог lo , и принимается решение

об отсутствии

цели (сигнала),

если отношение

правдоподобия оказывается

меньше порога

lo . Таким образом, оптимальный обнаружитель действует по следующему

45

алгоритму: по принятым значениямr1 ,..., rm вычисляется отношение

правдоподобия и производится его сравнение с порогом.

Отношение правдоподобия (2.16) с учётом выражений(2.13) и (2.14) можно представить в виде

æ

 

1

m

ö

2

m

 

l(r) = exp ç

-

åsi2Dt ÷exp

åri si Dt .

(2.18)

No

No

è

 

i=1

ø

i=1

 

Чтобы осуществить переход к непрерывному наблюдениюr(t) ,

необходимо

в выражении(2.18)

перейти

к

приделуDt ® 0 . При этом

отношение

правдоподобия (функция

 

 

переменных r1 ,..., rm )

перейдёт в

функционал отношения правдоподобия:

 

 

 

 

 

æ

1

T

ö

æ

 

2

T

ö

 

l(r(t)) = exp ç-

 

òs2 (t)dt ÷exp ç

 

 

òr(t)s(t)dt ÷.

(2.19)

No

 

No

 

è

0

ø

è

 

0

ø

 

Функционалом называют переменную величину, значение которой зависит от вида функции. Оптимальное решающее правило(2.17) для непрерывного наблюдения записывается следующим образом:

qopt*

ì1,

åñëè

l(r(t)) > l ;

 

(r(t)) = í

åñëè

o

(2.20)

 

î0,

l(r(t)) < lo .

 

Правило сохраняется при переходе к монотонным функциям от l(r(t)) . В частности, применив логарифмическую функцию, получим

1nl(r(t)) = -

E

+

2

z(r(t)) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

No

 

 

No

 

 

 

 

 

 

 

 

где E = òs2 (t)dt - энергия сигнала;

 

 

 

 

 

 

0

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z(r(t)) = z(T ) = òr(t)s(t)dt

 

 

 

 

 

 

(2.21)

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

есть

корреляционный

интеграл. С

 

учётом

этих

соотношений

 

оптимальное решающее правило принимает следующий вид:

 

 

qopt* (r(t)) =

ì1, åñëè

z(T ) > z

ï

;

 

 

 

(2.22)

 

í

 

åñëè

 

 

 

 

 

 

 

 

 

 

î0,

z(T ) < zï ,

 

 

 

 

 

где порог zп

= (1nlo + E No ) No 2 .

 

 

 

 

 

2.2.3 Корреляционный оптимальный обнаружитель

 

 

В соответствии с алгоритмом обнаружения(2.21), (2.22) структурная

 

схема

оптимального

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

может

быть

представлена

ви

корреляционного

 

приёмника с

пороговым

устройством(рисунок 2.1). В

 

z(T )

 

 

 

 

46

 

 

 

 

 

пороговом

 

устройстве (ПУ)

производится

сравнение

значения

корреляционного

интеграла в

момент

ожидаемого

окончания действия

сигнала T

с порогом zп

и принимается решение о наличии или отсутствии

цели (сигнала).

Начало

интегрирования и

его

окончание

совпадают по

времени

с

началом

и

окончанием

ожидаемого

sсигнала(t) , что

обеспечивается

устройством

синхронизации(УС).

Это

же

 

устройство

синхронизирует работу генератора опорного сигнала (ГОС) для коррелятора.

Работа

корреляционного

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

поясняется

временными

диаграммами в точках 1-4 оптимального обнаружителя представленными на рисунке 2.2. Диаграммы рисунка 2.2, а соответствуют наличию в смеси r(t) и сигнала, и шума, диаграммы на рисунке 2.2, б – наличию в смеси только

шума. Сравнение напряжения, вырабатываемого на выходе интегратора с порогом производится в момент времениt = T , при этом величине является случайной и зависит от реализации шума на интервале[0,T ].

При наличии в смеси r(t) сигнала и шума величину z(t) можно представить в виде сумме:

T

T

z(T ) = òs2 (t)dt + òn(t)s(t)dt = E + zсл .

0

0

Второе слагаемое zсл представляет случайную величину. Её среднее значение равно нулю. Действительно, с учётом условия (1.26) имеем

T T

< zсл >=< òn(t)s(t)dt >= ò< n(t) > s(t)dt = 0 .

(2.23)

0

0

 

Дисперсия случайной величины zсл определяется выражением

T T

< zñë2 >= òò< n(t1 )n(t2 ) >s(t1 )s(t2 )dt1dt2 =

0 0

TT

=(N0 2)òòd (t1 - t2 )s(t1 )s(t2 )dt1dt2 .

0 0

Здесь учтено, что автокорреляционная функция белого шума определяется d - функции (1.25), получаем