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

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

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

Совместное распределение вероятностей последова­ тельностей дискретных переменных її, как для шума, так и для сигнала с шумом, при преобразовании нестацио­ нарного процесса определяется из формулы (1.7). Коэф­ фициент правдоподобия при проверке двух простых гипо­ тез с учетом соотношения (1.7) принимает вид

t

П її Р№

(*>

Л ( ц ) = p ( " l 6 , ) =

* = 0 a = 1

(2.188)

п п 'bw<*>

Логарифмируя (2.188), получаем

шл й=gs (А) щ

=s s < *> ^ <*

А = 0 а = !

А = 0 а = 1

 

(2.189)

Поскольку значения л ; а , а = 1 , г ,

для каждого /г =

= 0,/г—1 образуют полную группу несовместимых собы­

тий, вторая

сумма в (2.189) содержит для данного

/г =

0, п—

1

всегда

одно

слагаемое,

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

попаданию

случайной

величины щ в интервал

квантова­

ния

Д« .

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

Структурная схема

алгоритма (2.189)

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

на

рис. 2.36.

Дискретные

переменные

йи

i = 0,

1, 2, ... ,

с выхода преобразователя записываются в /п-канальный

регистр сдвига (Рг.

с)

с числом разрядов

п.

Каждое

выборочное

значение

нг- параллельным кодом

вводится

в регистр адреса постоянного запоминающего

устройства

( П З У ) , хранящего

весовые

коэффициенты

Wa

(k).

В каждом из п ПЗУ

хранится

по 2т

различных

ве­

сов Wi(k),

W^ik),...,

W.r{k). С помощью

параллельного

сумматора образуется сумма весовых коэффициентов, которая сравнивается с порогом в цифровом компарато­ ре ( Ц К ) . Такая схема используется при обнаружении сигнала с неизвестным временем прихода.

Если последовательности дискретных переменных ІЇі для обеих гипотез соответствуют стационарным случайISO

u(t)

П

и. (к)

 

Рг.с

 

 

 

 

 

 

п-1

п-2

0

 

 

1

 

 

Сдвиг

Щп-2) .

 

 

ЦК

Рис. 2.36. Структурная схема алгоритма (2.189).

НЬІМІ процессам

«(•/), то в соответствии с

(1.9)

 

 

іплй =sА .1 п -^=Е X J «•

( 2 Л 9 0 )

Структурная

схема,

 

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

 

алгоритму

(2.190) , содержит г счетчиков, определяющих

 

Яа ; г ПЗУ,

хранящих веса

Wа; г

умножителей

(XaWa)

и

 

параллель­

ный сумматор с г входами.

 

 

 

 

 

 

Более удобна при реализации другая,

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

интерпретация

алгоритма

(2.190)

 

 

 

 

 

 

 

Ш Л ( « ) = 5

£

da(k)Wa.

 

 

 

(2.191)

На рис. 2.37 приведена структурная схема

 

алгоритма

(2.191) . Дискретные переменные гТг- с выхода

преобра­

зователя (П) поступают

в регистр

числа

ПЗУ, храняще­

го весовые коэффициенты

Wa,

ia=l,

г. С помощью на­

капливающего

сумматора

образуется

сумма

(2.191), ко­

торая

сравнивается

с

порогом

с

либо в

 

цифровом

компараторе (ЦК) , либо

в сумматоре, если

в него был

введен

обратный код числа с.

При этом

отрицательные

весовые коэффициенты удобно хранить в ПЗУ в допол­ нительном коде.

При малых интервалах квантования Д'«<^; (0,5ч-0,7)а, где а—среднеквадратическое значение помехи или сиг-

181

alt)

ЦК

ПЗУ

Рис. 2.37. Структурная схема алгоритма (2.191).

нала с помехой, из (2.46) получим

 

р а

=

j" ю,

(ut) dui

^

wl

( « * J Лив , а =

ГТгТ

 

где « * а £

Диа

— средняя

точка

интервала

Дма . Тогда ве­

совые

коэффициенты будут

равны

 

 

 

 

 

W (k) =

ln-

777- =

In

.- 5,

v - •

(2 192\

Алгоритм

('2.189)

с учетом

(2.192)

соответствует

(2.55).

В

задачах

обнаружения

некогерентных сигналов на

фоне шумов при независимых отсчетах процессов наи­ большее распространение получило бинарное квантова­

ние. Распределения

вероятностей и алгоритмы

обработки

при бинарном

квантовании

получаются

по

формулам

(2.188) — (2.192)

при подстановке

a='0,I,

г=2.

Алгоритм,

соответствующий (2.191), имеет вид

 

 

In Л (и) = 2

К (k) W0

(k) +

dx (k) Wx (k)}.

 

Следует учесть, что теперь dB(k)=l

—dt{k).

Тогда

In Л (и) = §

* {W0 (k) + rf, (k) [W, (k) -

Wa (k)\) =

- q + S ' - w

?™тЛ

 

° q + S " - w y w -

fe=o ft=o

 

 

 

 

 

(2.193)

 

 

 

 

 

 

где величина

Q = S V ( * ) = const

ті ее следует отнести к пороговому уровню.

Реализация алгоритма (2.193) аналогична реализа­ ции,182 представленной на рис. 2.36, но вместо /и-каналь-

його регистра сдвига используется одноканальний и каж­ дый из весовых коэффициентов W(k) принимает одно фиксированное значение. По сравнению с алгоритмом (2.68) весовая функция теперь определяется не произ­ вольно, а наилучшим образом, хотя, как показывает со­ ответствующий анализ, эффективность подобных обнару­ жителей слабо зависит от весовой функции. Можно по­ казать i[95], что при обнаружении слабых некогерентных сигналов W{k)^G\, k = 0, п—1, т. е. весовая функция равна квадрату диаграммы направленности по мощно­ сти. При изменении 'интенсивности сигналов весовая функция будет также изменяться, хотя с практической точки зрения это не существенно.

Таким образом, для независимых выборочных значе­

ний сигнала и шума

рассматриваемый метод синтеза да­

ет те же

алгоритмы

обнаружения, что и

метод синтеза

цифровых

эквивалентов по аналоговому

прототипу (см.

§ 2.3).

 

 

 

Анализ эффективности некогерентных обнаружителей бинарно-квантованных сигналов проведен в § '2.3. Сле­ дует указать, что алгоритм (2.193) соответствует опти­ мальному обнаружению целей при отсутствии флуктуа­ ции и при' шумоподобных флуктуациях. В случае друж­ ных флюктуации пачки эхо-сигналов коэффициент правдоподобия определяется в соответствии с (2.56). В результате получается алгоритм, реализующий макси­ мально правдоподобное обнаружение (2.8) с весовой функцией, зависящей от интенсивности отраженного сиг­ нала. Однако такое усложнение системы обработки, как отмечалось в п. 2.3.3, практически не оправдано.

Перейдем теперь к вопросам обнаружения некоге­ рентных импульсных сигналов методом последователь­ ного анализа. Преимущества последовательной процеду­ ры Вальда по сравнению с процедурой Неймана — Пир­ сона хорошо известны [61]. Однако реализация алгоритма последовательного анализа в, многоканальных многоце­ левых Р Л С стала возможной при использовании фазиро­ ванных антенных решеток, программного обзора про­ странства и применении ЭЦВМ для обработки информа­ ции.

Последовательные процедуры в многоканальных си­ стемах можно разделить на три группы:

I . С игнорированием разрешающей способности.

183

2.С однократными пересечениями (независимым при­ нятием решения в каждом дальномерном канале и с по­ глощающими порогами). При этом в каждом канале вы­ числяются парциальные коэффициенты правдоподобия, значения которых сравниваются с двумя порогами. После достижения одного из них выносится соответствующее решение 'и испытания в данном канале прекращаются до окончания процедур в остальных каналах.

3.С многократными пересечениями порогов (иепо-

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

Алгоритм процедур первой группы [96] связан с обра­ зованием коэффициента правдоподобия по формулам (2.3), (2.4), где параметр и соответствует времени за­ паздывания отраженного сигнала Т. Практическое ис­ пользование подобных алгоритмов сильно ограничено требованиями наличия только одной цели и независи­ мостью дальномерных каналов. Кроме того, такая систе­ ма не обладает разрешением по дальности.

Процедуры второй группы обладают разрешающей способностью и могут использоваться при любом числе целей на одном направлении антенного луча. Ограничим­ ся рассмотрением процедур второй группы только с би­ нарным квантованием некогерентных сигналов. Лога­ рифм коэффициента правдоподобия в каждом дально­

мерном канале

с учетом

(2.189)

можно представить

в виде

 

 

 

4^f"

<«>=S \

+ (1 -

d (k)) In

In Л , (и) ) ) \

dС (k) In ^

 

 

fc=l L

 

 

 

 

 

 

 

 

ln \ ~ P s »

n,

PN

Л—РЫ

1 —PN

 

 

 

 

(2.194) где Ai — число единиц в выборке объема п. На каждом

—>

шаге 1пЛп(й) сравнивается с двумя порогами (2.19). Путем преобразования (2.194) можно перейти к весо­ вому варианту последовательной процедуры [61] с поро­ гами, заданными величинами

В

' = 0 '

A - = S =

. . K l - X w - ^ , ! •

< 2 ' 1 9 В >

184

 

 

 

 

с которыми на каждом шаге сравнивается величина

. . , , ~

— i n В

.

^ [ p s N ( l - p N ) / p N ( \ - p S N ) ] - i

 

^ • ( ^ - 1 ^ + — i n m - P s N w - p N ) ] —

х

п

 

 

п

 

п

 

 

Х ' £

d(k)-

J ] ( l - d ( * ) ) = r 0 + 2][d(*)r -

 

ft=l

 

/6=1

 

ft=l

 

 

 

 

 

-

(1 -d.{k))s\.

(2.196)

Первое

слагаемое

(2.196)

интерпретируется

как на­

чальный вес процедуры го, а

множители перед

сумма­

ми — как веса

г и s,

добавляемые при каждом

появле­

нии единицы или нуля

соответственно.

 

 

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

информация,

поступающая с

момента окончания

испы­

тания в отдельном

канале до завершения испытаний во

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

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

одновременно

во всех

каналах. Среднее число выборок для систем вто­

рой и

третьей групп будет определяться

наибольшим

средним числом выборок, требуемых для элемента с мак­

симальным временем

анализа. Теоретически

процедуры

с многократными пересечениями имеют некоторый вы­

игрыш по сравнению с

процедурами

второй

группы, од­

нако, как показывают

расчеты, этот

выигрыш очень мал.

Для анализа эффективности последовательных про­

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

весьма

простым и наглядным

методом

[97] . Последовательная

процедура от шага к шагу представляет собой

случай­

ное -блуждание, описываемое цепью Маркова.

Каждому

состоянию

Ах=Аы

такой

цепи

ставится в

соответствие

комбинация

из k единиц и / нулей,

которая

получается

после

п-го

испытания

 

данного

элемента

разрешения

(n = k-\-t).

Комбинации,

содержащие одинаковое

число

единиц и пулей, но с различным

их чередованием,

отно­

сятся к одному состоянию. Вероятность попадания

в со­

стояние

Аы

 

равна

р(Ам)

=С(п,

k)ph{\—p)n-k.

 

 

 

Число

С {її, k)

соответствует

возможным

способам

образова­

ния комбинации из к единиц и пк нулей.

Случайные блуждания, интерпретирующие последо­ вательную процедуру, представим на плоскости, где по оси ординат откладываем шаги испытаний, а по оси

185

Число „нулей.'"
Рис. 2.38. Плоскость случайных блужданий последовательной процедуры.

абсцисс — число нулей (рис. 2.38). Каждая точка на та­ кой дискретной плоскости представляет собой состояние цепи Маркова, и каждому состоянию ставится в соответ­ ствие вероятность попадания в это состояние, зависящая только от вероятностен в двух точках, находящихся на предыдущей строке — непосредственно под определяемой и слева от нее. Если цепь имеет поглощающие состояния,

то они выделяются на ри­ сунке значками — кру­ жочками н квадратиками. В дальнейшем для про­ стоты опускаем обозначе­ ния состояний и записы­ ваем только соответствую­ щие им вероятности. По­ глощающие состояния со­ ответствуют выходу за пороги • последовательной процедуры или достиже­ нию их.

Плоскость блужданий делится на три области: принятия решения о на­ личии сигнала, отклоне­

О1 2 3 * ния решения и продолже­

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

сигнал/шум <72 среднее число выборок многоканальной процедуры при отсутствии сигнала HN было •минимальным. Считаем, что решения принимаются отдельно в каждом канале и что сигнал от цели присутствует в данном ка­ нале все время анализа пли не присутствует совсем. При поиске оптимальных областей конкретная структура обнаружителя не задается, 'Поскольку в результате ана­ лиза можно получить либо усеченную процедуру, либо процедуру с постоянными или переменными порогами, либо весовой вариант испытаний. Однако при использо­ вании неусеченного последовательного анализа стремле­ ние к.уменьшению EN ведет к увеличению среднего числа выбцрок при нал'ичи сигнала HSN, которое может стать

186

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

Рассчитаем вероятности F и D для примера, пред­ ставляющего усеченную на пятом шаге последовательную процедуру (рис. 2.38). Вероятность ложного обнаруже­ ния в элементе разрешения равна

Отсюда по заданной величине F можно найти вероят­ ность превышения порога квантования м0 шумовым вы­ бросом pN. По известной величине «о и вероятности PN из таблиц или графиков {98] находится вероятность psx [см. (2.47), (2.42), (2.43)]. Вероятности правильного об­ наружения и пропуска цели определяются как

D = РІК + 3P3SN

4SN +

<?SN '

 

M = \ - D = q*SN+2pSN

q3SN+5plN

g3Sf/

(2.198)

Вероятности пересечения или достижения порогов при отсутствии сигнала от цели за і шагов рг определяются следующим образом:

pl

=

0,p2 = irN,pl

= pa +

p3lt.

 

р, =

р3

+ Зр3, qN +

2pNq3N ,P i =

1.

(2.199)

Вероятности пересечения порогов за і шагов при отсут­ ствии сигнала в М элементах разрешения рм, І вычис­ ляются по формуле

Рм.і=Р"-

i = l - 2 5 -

(2-20°)

Вероятность пересечения порогов на і-м шаге

•*М.І = РМ,І-РМ.І-І-

С 2 - 2 0 1 )

Среднее число испытаний для шумовых элементов раз­ решения .

.(2.202)

І = І

187

Алгоритм работы устройства состоит• в накоплении первых трех единиц на пяти позициях, при этом выдается сигнал обнаружения; сброс осуществляется по появле­ нию на первых двух позициях двух нулей подряд пли по накоплению трех нулей па пяти позициях.

2.6.3. Синтез цифровых обнаружителей и анализ их эффективности при зависимых отсчетах [164—166]. За­ висимость между отсчетами квантуемых процессов соот­ ветствует ситуации, когда период временной дискретиза­ ции меньше времени корреляции процессов: Г д < т л . При­ чины существования соотношения Tn<Xh таковы.

1.Период дискретизации Гд выбирается не из соот­ ношения Гд^І/2/max, а из других соображений, например из условия однозначной дальнометрии в импульсной ра­ диолокации: Гд = Г п = '2 /?щах/С.

2.Часто с целью уменьшения ошибок восстановления днскретизированной информации (см. п. 4.2.2) период Г д

приходится

брать в несколько раз

меньше, чем l/2/mas-

3. Ширина спектра помехи AfP

может

быть неизвест­

на и Ді/р max</max, а период дискретизации

выбран из со­

отношения

Г д = 1 / 2 / ш а х .

 

 

4. При

проверке двух гипотез,

соответствующих рас­

пределениям с неперекрывающимися спектрами, выберем Г д = l/2/max, где /maxмаксимальная частота в спектре процесса, расположенного в области более высоких ча­ стот. Тогда для процесса, спектр которого находится в области более низких частот, отсчеты будут статисти­ чески зависимыми.

В соответствии с п. 1.7.1 последовательность дискрет­

ных

переменных Гц,

1 = 1 , 2, ... , принимающих значения

ха,

<х = 1, г, можно

аппроксимировать регулярными це­

пями Маркова

с г состояниями и связностью, определяе­

мой

из (1.81)

или (1.83). В общем случае связность це­

пи Маркова будет различной для обеих гипотез Я 0

и # i .

Ясно, что при

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

сле­

дует учитывать наибольшую связность vm ax, полагая, что обе гипотезы соответствуют распределениям со связ­ ностью Vmax-

Для выявления общих закономерностей и удобства сопоставления с обнаружителем при независимых отсче­ тах (2.189) рассмотрим синтез обнаружителя, когда по­

следовательности

дискретных переменных, соответствую­

щие

гипотезам Н0

и Ни

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

собой

односвяз-

ные

(v = l ) регулярные

неоднородные

цепи

Маркова

188

 

 

 

 

 

с г состояниями. Коэффициент правдоподобия для рас­ пределении типа (1.10) имеет вид

II—1 г

п

Л ( и ) = I '

/г=1

В=1

(2.203)

И — 1

Г

 

/»(2|8,)

п

 

 

/г=1 р=1

 

о = 1

В результате логарифмирования (2.203) получаем алго­ ритм оптимальной обработки

In Л (и) = J] ^

(0) + 5 '

X с/ар (/г)

(/г), (2.204)

а = 1

ft=l

а,В

 

где

 

 

 

WJ0)=^(P0aSNIP0«N)-

К?

{k)F=b(pa%SN(k)/Pe9N{k)).

Ha рис. 2.39 представлена структурная схема алго­ ритма (2.204). Ее отличие от схемы рис. 2.36, соответст­ вующей алгоритму (2.189), состоит в том, что т:перь весовые коэффициенты Wa? (k) зависят от значений ха и

в два смежных момента времени. Исключение состав­

ляет вес

Wa[0), зависящий только

от значения

ха. Соот-

u(t)

й(к)

 

 

 

 

n-1

n-2 P Z ' C

1

О

 

Сдвиг

; . .

. і ' \

ШО)

 

 

 

_J

 

 

 

тт

т т

 

 

 

z

 

ЦК

Рис. 2.39. Структурная схема алгоритма (2.204).

ветственно ПЗУ, хранящие весовые коэффициенты, долж­ ны иметь 2т-разрядный регистр адреса и хранить в об­ щем случае 2 2 т различных весовых коэффициентов. В остальном обе схемы аналогичны.

189

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