 
        
        книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfСовместное распределение вероятностей последова тельностей дискретных переменных її, как для шума, так и для сигнала с шумом, при преобразовании нестацио нарного процесса определяется из формулы (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). Каждая точка на та кой дискретной плоскости представляет собой состояние цепи Маркова, и каждому состоянию ставится в соответ ствие вероятность попадания в это состояние, зависящая только от вероятностен в двух точках, находящихся на предыдущей строке — непосредственно под определяемой и слева от нее. Если цепь имеет поглощающие состояния,
то они выделяются на ри сунке значками — кру жочками н квадратиками. В дальнейшем для про стоты опускаем обозначе ния состояний и записы ваем только соответствую щие им вероятности. По глощающие состояния со ответствуют выходу за пороги • последовательной процедуры или достиже нию их.
Плоскость блужданий делится на три области: принятия решения о на личии сигнала, отклоне
О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
