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