
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfКомплексный |
|
вектор принимает значения |
|
|||
и; = |
(0,0,...,0), и] = |
{0,0 |
«,), |
|
||
и2 — (0, 0, |
...,//.,, ц . , ) , / . ^ ^ ( H j , «2 , |
...,ип), |
||||
|
" ^ + |
| 1 = ("ft. "ft + i . |
ик+п)- |
|
|
|
Представим |
и и Я * в алгебраической |
форме |
||||
Я* |
= Н' — Ш", и = и ' + и?' |
(2.177) |
||||
и подставим ('2.177) |
в (2.176). В |
скалярном |
произведе |
нии комплексных векторов |
(2.176) выделим действитель |
|||
ную и мнимую части, |
выполнив |
подстановку |
(2.177): |
|
v (k)=5' |
( / з |
> + |
# ' > ' w |
^ + |
Выражение (2.178) представляет собой алгоритм сколь
зящего |
весового суммирования. В частных случаях |
|
(2.178) |
упрощается: |
|
а) |
u(t)—комплексный, |
H(t)—вещественна, |
v ( k ) = 1 S |
' |
#V"'W |
(2 -1 7 9 ) |
|J.=:0 |
Ц=0 |
|
|
б) ы (£) — вещественный, |
Я * (£) — комплексна, |
|
в) ы (tf) — вещественный, H (t) — вещественна,
Y >) = 2 # > W |
' |
(2-181) |
|
|j.=0 |
|
|
|
Здесь положено 7 ^ = 1 . |
|
|
|
Соотношения (2.178) — (2.181) |
и определяют |
структу |
|
ру обработки. Если излучаемый |
сигнал — монохромати- |
170
ческий, а допплеровское смещение частоты сигнала от сутствует, тогда Н(/) вещественна, a u(t) — комплексный из-за незнания начальной фазы. Структурная схема, соответствующая (2.179), приведена на рис. 2.33. В каж дом квадратурном канале осуществляется бинарное
cosy
Рг. с і.
|
|
|
Г |
Т |
Т " |
г |
f,lt) |
г д |
В |
|
ФП |
|
|
|
|||
|
|
|
Е2 |
|
|
|
|
|
|
|
ЦК |
|
|
|
|
рг.с |
2 |
Рис. 2.33. Цифровой фазовый коррелятор.
квантование и временная дискретизация с. периодом Тя выходных напряжений фазовых детекторов.
|
Дальнейшая обработка'производится в регистрах сдви |
|||||||||||
га |
(Рг.с). Каждый |
из |
весовых |
коэффициентов |
|
Н„Нй,... |
||||||
. . . , # „ |
представлен |
двоичным |
кодом. |
Единица |
с |
выхода |
||||||
/е-го разряда рггистра сдвига являзтся |
командой тіа |
вве |
||||||||||
дение |
веса |
Нк |
в |
параллельный |
/г-входовый |
сумматор. |
||||||
.Функциональный |
пр .образователь |
(ФП) |
выделяет |
огибаю |
||||||||
щую Ё (£)==]/"Яа (£) + |
Са (Л) |
или ее |
квадрат |
|
Ez(k) |
= |
||||||
= |
А2 (/г) -|- С 3 |
(/г). |
Для |
упрощения |
функционального преоб |
разователя иногда производят отбор максимального зна
чения Ж(к) |
или C(k). С этой целью в цифровом компара |
||
торе |
(ЦК) |
поочередно сравниваются с порогом A{k) |
и |
C(k). |
Результатом этого сравнения является наличие О |
||
или I на выходе ЦКСледствием такого упрощения яв |
|||
ляются дополнительные потери в пороговом сигнале, |
не |
171
превышающие 0,7 дБ. Соответствующий анализ показы вает [91], что цифровой фазовый коррелятор (ЦФК) при бинарном квантовании в случае выполнения условии (2.172) обладает квазилинейностью, но по сравнению с фазовым обнаружителем дает потери в пороговом сиг нале — 2 дБ .
Известно применение ЦФК для сжатия фазоманипулироваиного (ФМ) сигнала [92]. В этом случае весовые коэффициенты Hh, k—\, п принимают значения 1 или О и схема обработки упрощается. В соответствии с кодом ФМ сигнала в каждом разряде регистра сдвига напря жение снимается с выхода либо Qu, либо (?/,. Суммирова ние теперь можно производить с помощью параллельного счетчика [93]. Отметим, что при неизвестной допплеровской частоте эхо-сигнала структура накопителя опреде ляется алгоритмом (2.178).
2.5.4. Цифровые ранговые обнаружители. Ранговая ин-
формация, |
содержащаяся в выборочном векторе |
и=(щ, |
|||
« 2 , |
..., ип)т, |
определяется вектором рангов R = |
(Ri, |
R2, ... |
|
|
|
|
|
|
—> |
..., |
Rn)*, где Ri—ранг |
выборки Ui. Для получения R не |
|||
обходимо |
расположить |
выборки ии и.% ..., ип |
в |
порядке |
возрастания их значений слева направо. Порядковый но
мер каждой |
выборки И; в ранжированной |
(упорядочен |
||||
ной) последовательности и является ее рангом Ri. |
На |
|||||
пример, если |
и = ( 4 , 9, 8, 2, |
1) т , то |
# = ( 3 , |
5, |
4, 2, |
1) т . |
При этом предполагается, что вероятность |
равенства |
|||||
двух выборок равна нулю, т. е. все выборки |
имеют |
не |
||||
одинаковые ранги. |
|
|
|
-> |
|
|
|
|
|
|
|
|
|
Определим |
относительный |
ранговый |
вектор |
/ выборки |
||
и 94]. Положение наибольшей |
выборки • обозначим /,, |
счи |
тая слева направо. Позицию второй наибольшей выбор
ки среди оставшихся п—1 обозначим h |
и т. д. |
Наимень |
шая выборка с необходимостью будет |
иметь |
положение |
|
-> |
|
1п=1. Для рассмотренного выше примера /—(2,2, 1, 1,1) т . Можно также ввести вектор абсолютных рангов [84], расположив выборки в порядке увеличения их абсолют
ных значений:
k , l < K l < - < K n r
Перейдем теперь к рассмотрению структур ранговых об наружителей.
17
Обнаружитель Вилкоксона [84] является линейным ранговым и основан на сравнении суммы абсолютных рангов положительных выборок с порогом
|
t Вь>С |
|
(2Л82) |
где |
|
|
|
|
О |
при |
« t - f t < 0 , |
|
k |
при |
" i f t > 0 |
и \ut |
I << І «г21 <C-"<C| uin\- Такой |
обнаружитель используют, |
|
когда |
шум имеет симметричную |
плотность Wi(u) с нуле |
вым средним, а сигнал с шумом имеет среднее значение, большее нуля. При наличии сигнала положительные вы борки станут встречаться чаще и будут иметь большие абсолютные ранги. Асимптотическая относительная эф
фективность |
обнаружителя |
Вилкоксона |
по |
отношению |
|||||
к линейному при гауссовом шуме равна |
Ezi=0,95, |
а при |
|||||||
ЛЮбОМ Шуме С СИММетрИЧНОЙ ПЛОТНОСТЬЮ £7 21 = 0,864. |
|||||||||
Для технической |
реализации более |
удобна |
другая, |
||||||
эквивалентная |
форма |
алгоритма |
(2.182), |
связанная |
|||||
с парными суммами [94]: |
|
|
|
|
|
|
|||
|
£ |
gk= |
S |
£ |
<*("< + |
щ)>с, |
|
|
(2.183) |
где d(Ui-\-Uj) |
=d(Uh) |
определяется |
из (2.60a). |
|
|
||||
Обнаружитель Вилкоксона может быть выполнен на |
|||||||||
элементах вычислительной |
техники и |
его |
реализация |
сходна с рассматриваемым ниже обнаружителем, реали зующим тау тест Кендалла [94].
Тау тест Кендалла, называемый еще критерием трен да, является нелинейным ранговым и основан на обра
зовании статистики |
|
|
|
|
|
||
Г |
=2 |
JJ d ( щ - |
и < ) = |
Щ ± - |
J] |
к. |
(2.184) |
|
/=1 |
!=i+i |
|
» = i |
|
|
|
Иными |
словами, каждое |
значение ы,-, i—\,n |
|
сравнивает |
|||
ся со всеми |
остальными |
щ, j = |
i - j - 1. п |
и |
производится |
||
подсчет |
тех |
ситуаций, |
когда |
U j > u t - . |
|
Статистика Т% |
173
при гипотезе Н0 асимптотически нормальна с параметрами
m, 1^1 = - £ - ( / * - 1 ) ,
о 3 [ Г т ] = - 4 - ( я - 1 ) ( 2 / г + 5):
Алгоритм (2.184) может использоваться в «обнаружите лях движущегося окна», поскольку учет каждого нового значения u„+i связан с простым рекурсивным вычисле нием
7\ = 7\ + Е < / ( и п + 1 - И і ) .
равно как и учет старого значения щ, .покидающего «дви жущееся окно» —
(=2
Таким образом, в стационарном режиме обзора в мо мент времени tu при k~>n необходимо образовать стати стику
|
п |
п— I |
T(k) = T(k— 1) + |
Е d(uk — uh_i) — Ц d(u,t_n — ик.г). |
|
|
1=1 |
(=1 |
|
|
(2.185) |
Использование |
алгоритма |
(2.185) при обработке ра |
диолокационной информации связано с тем обстоятель ством, что пачка эхо-сигналов в Р Л С обзора промодулирована диаграммой направленности. Благодаря этой мо дуляции при наличии точечной цели образуется харак терный ранговый порядок (тренд). Сравнение критерия
тренда |
для когерентного |
(гауссовы распределения шума |
|
и сигнала с шумом) и некогерентного |
приема с плотно |
||
стями |
(2.34) и (2.43) |
и оптимальных |
обнаружителей |
приведено на рис. 2.34. Эти кривые для малых значений п получены методом Монте-Карло при аппроксимации диаграммы направленности антенны гауссовой кривой. При п—мх> определена асимптотическая относительная эффективность.
Приведенные результаты справедливы для нефлуктуирующей пачки эхо-сигналов. При шумоподобных
174
|
|
* |
|
|
|
/6 |
|
|
St |
|
|
256 |
|
|
n |
|
|
|
|
||
Рис. 2.34. Потери |
в пороговом |
сигнале |
тау-обнаружителя |
Кендалла |
|||||||||||||||||
|
|
по сравнению с оптимальным |
обнаружителем. |
|
|
|
|
||||||||||||||
флуктуациях |
требуется |
большее число |
импульсов |
|
в пач |
||||||||||||||||
ке, |
|
чтобы |
потери |
|
были |
допустимыми. |
Как |
|
следует из |
||||||||||||
рис. |
|
|
|
|
|||||||||||||||||
|
|
2.34, |
применение такого |
обнаружителя |
|
целесообраз |
|||||||||||||||
|
|
когда |
|
||||||||||||||||||
но, |
|
|
число импульсов |
в пачке |
велико |
(/г |
= 5 0 - М 0 0 ) . |
||||||||||||||
В |
противном |
|
случае потери |
могут |
быть весьма |
|
|
значи |
|||||||||||||
тельными. |
Следует |
отметить |
сложность |
технической |
реа |
||||||||||||||||
лизации критерия |
|
тренда, |
так |
как |
она |
связана |
|
с |
|
необ |
|||||||||||
|
|
|
|
|
|
ходимостью многоразрядного аналого-цифрового преоб
разования |
выходного |
напряжения |
приемника |
и |
многока- |
||||||||||||||||||||
нальностыо |
|
|
(число |
|
каналов |
равно |
|
|
числу элементов |
раз |
|||||||||||||||
решения |
по |
дальности). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
2.5.5. Модифицированный |
|
знаковый |
обнаружитель. |
||||||||||||||||||||||
импульсной |
|
Р Л С |
обзора |
|
с разрешением |
по |
дальности |
||||||||||||||||||
В угловым |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
и |
|
координатам |
|
при накоплении пачки |
бинарно- |
||||||||||||||||||||
квантованных |
-сигналов |
|
можно |
|
использовать |
непарамет |
|||||||||||||||||||
рический |
|
метод |
обнаружения. |
|
Вместо |
квантизатора |
|||||||||||||||||||
в этом случае применяется модифицированный |
знаковый |
||||||||||||||||||||||||
обнаружитель |
{83], |
|
изображенный |
|
на |
рис. 2.35. |
Видеона |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
на |
|||||||||||||||
пряжение с выхода |
|
амплитудного детектора поступает |
|
||||||||||||||||||||||
вход линии |
|
задержки |
с |
отводами. |
Число |
отводов |
т== |
||||||||||||||||||
= 5 - И 0 , |
задержка |
|
х между |
|
отводами |
равна |
7"д |
(периоду |
|||||||||||||||||
следования |
|
|
импульсов |
|
дальности), |
т. е. |
соответствует |
||||||||||||||||||
разрешающей |
способности |
Р Л С |
по дальности. |
Выходное |
|||||||||||||||||||||
напряжение |
|
от |
среднего |
|
отвода |
ЕІ |
|
сравнивается |
в ком |
||||||||||||||||
параторах |
Кь |
|
Кг, |
Кз, |
К/* с |
|
напряжениями |
|
|
предшест |
|||||||||||||||
|
|
+ ь |
|
£ , + 2 |
И |
|
|
|
от -2. |
отводов. |
|||||||||||||||
вующих £ , - |
|
|
|
|
|
|
ПОСЛеДуЮЩИХ |
Еі-і, |
|
Ei |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
175
|
|
АД |
|
|
Л |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ui-! |
|
|
4-і |
|
|
|
|
|
Z , |
\1 |
2\ \t |
|
|
|
|
г |
|
|
||
|
|
|
|
|
|
|
|
|
F |
|
|
|
|
|
|
Рис. |
2.35. Модифицированный знаковый обнаружитель. |
|
|||||||||||
Выходное |
напряжение |
каждого |
|
компаратора равно |
|
|||||||||
если Ег превышает |
напряжение на |
|
|
входе |
компа |
|||||||||
|
|
|
I , |
|||||||||||
ратора. Выходные |
напряжения |
компараторов |
|
совместно |
||||||||||
|
|
втором |
|
|
|
|
|
|||||||
с импульсами дальности |
поступают |
на |
|
|
|
схемы |
И. |
|||||||
|
1выходе фчрмирователя F при |
наличии |
|
единиц |
на |
|||||||||
На |
|
входы |
|
|
|
|||||||||
всех |
входах схемы |
И |
образуется |
импульс стандартной |
длительности и амплитуды. Дальнейшая обработка по следовательности единиц и нулей осуществляется, одним из методов бинарного обнаружения, рассмотренных ра нее в этой главе [например, методом весового бинарного накопления (см. п. 2.3.4)].
Для оценки эффективности данной схемы определим вероятности появления единицы под действием одного шума pN и сигнала с шумом pS w- Для одного шума на пряжения Еі+2, ЕІ+І, Ей Еі-і, ЕІ-2 одинаково распределе
ны и независимы |
друг от друга. Последнее утверждение |
|||
справедливо, если |
TKp>\f&f, |
где |
А/ — полоса пропуска |
|
ния УПЧ (см. п. |
2.3.1). |
|
Fi(u)=p(E<u)—интег |
|
|
|
Пусть |
ральная функция распределения Е. При постоянном зна
чении |
Ei = U |
ВерОЯТНОСТЬ ТОГО, ЧТО |
Еі+2, |
ЕІ+І, |
ЕІ, |
ЕІ-І, |
||||
Ei-Zi |
. . . |
|
|
Для |
определения |
|||||
|
|
меньше чем и, равна |
[Fi(u)]m. |
|
||||||
вероятности |
pN необходимо yqpeAHHTb |
[Fi(u)]m |
по |
всем |
||||||
и, т. е. |
|
|
|
|
|
|
|
|
|
|
|
|
|
оо |
|
со |
|
|
|
|
|
|
р ы = |
|
j [ ^ ( " ) l m t M « W " = |
^\F,(u)\^dF(a): |
т+ |
1 ' |
||||
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
(2.-186) |
176
где Wi(u) =dFi(u)/du\ |
m — число |
отводов от линии за |
|
держки. |
|
|
|
Різ (2.186) следует, что вероятность появления шумо |
|||
вой |
единицы рк не зависит от мощности шума и его за |
||
кона |
распределения, |
а заданная |
вероятность рк обеспе |
чивается выбором т. |
Тем самым |
при дальнейшей обра |
ботке обеспечивается постоянство вероятности ложной
тревоги. Пусть теперь в і-м |
элементе разрешения |
при |
|||||||||
сутствует сигнал |
ЕІ (импульс от цели), |
а в т—-1 |
сосед |
||||||||
н и х — только |
шум. Напряжение ЕІ имеет функцию |
рас |
|||||||||
пределения |
F'i |
такую, что |
F\(u)<Fi(u) |
|
для |
всех |
Е. |
||||
Тогда, интегрируя |
по частям, |
получим |
|
|
|
|
|
||||
|
|
|
|
со |
|
|
|
|
|
|
|
|
|
Рм= |
j l F > ) ] W ' i ( " ) |
= |
|
|
|
|
|||
|
|
|
|
—со |
|
|
|
|
|
|
|
|
|
оо |
|
|
|
|
|
|
|
|
|
= 1 - / и |
\jF'1(u)[F1(u)]^-4Fl(u)>-^~. |
|
|
(2.187) |
|||||||
|
|
—со |
|
|
|
|
|
|
|
|
|
При отсутствии |
флуктуации |
цели для |
получения |
pSN |
|||||||
в (2.187) следует подставить (2.42), (2.43), |
а при шумо- |
||||||||||
подобных |
флуктуациях — (2.33), |
(2.34) |
с |
подстановкой |
|||||||
a- = as-\~aN . |
Как |
показывают соответствующие |
расчеты |
[83], модифицированный знаковый обнаружитель проиг рывает всего около 1 дБ обнаружителю движущегося окна и около 2 дБ идеальному аналоговому видеоинтег ратору. Однако постоянство вероятности ложной тревоги оправдывает практическое применение такого устройст ва. Недостатком рассмотренного обнаружителя является подавление слабых сигналов: наличие сильного сигнала в пределах ±(/п—1)/2 элементов разрешения относи тельно слабого сигнала снижает вероятность обнаруже ния последнего. Уменьшение этого эффекта достигается за счет усложнения схемы обработки [83].
2.6. СИНТЕЗ Ц И Ф Р О В Ы Х О Б Н А Р У Ж И Т Е Л Е Й ПОСЛЕ А Н А Л О Г О - Ц И Ф Р О В О Г О
ПР Е О Б Р А З О В А Н И Я
2.6.1. Вводные замечания. В задачах радиолокацион ного обнаружения принятый сигнал может быть пол ностью известным (детерминированным) 5 ( t ) , содержать
несколько неизвестных случайных параметров c c = ( a i , . . .
12—1410 |
177 |
а , „ ) т , не зависящих от времени S(a, |
t) |
(квазпдетер- |
-> |
|
|
минированным) или случайным S(a(\t), |
t). |
Комбинация |
сигнала, шума и помехи чаще всего представляет их ад дитивную смесь при гипотезах Н0 и Ну.
при |
#о |
« ( 0 =p('t) |
+ti(t), |
при |
Ну |
|
u{t)=S{t)+p{t)+n{t). |
Зная плотности вероятностей для сигнала, помехи и шума, можно определить плотность вероятности аддитив ной смеси u(t), используя теорему свертки [1].
Пусть аналитических выражений для плотностей wn(u\\Qo)
и wn{u[Qy) получить не удается и они записаны в виде интегралов свертки. После аналого-цифрового преобра зования процесса-смеси u(t) последовательности дискрет
ных |
переменных |
щ, |
/ = 1 , 2, ... , .принимающих |
значения |
ха, |
а = 1 , г можно |
аппроксимировать цепями |
Маркова |
|
с г |
состояниями |
(см. п. 1.2.1). Приводимые ниже алго |
ритмы обнаружения в рамках такой аппроксимации име ют функциональную форму, не зависящую от распреде лений сигналов и помех на входе аналого-цифрового пре
образователя. Структурная |
инвариантность |
алгоритмов |
обнаружения имеет большое практическое значение. |
||
На практике перед аналого-цифровым |
преобразова |
|
нием процесса ('1.127) u(t, |
X(t), a[t)) обычно произво |
дится предварительная его обработка, связанная с уст ранением избыточной информации и облегчающая после дующие аналого-цифровое преобразование и цифровую обработку. Предварительная обработка может быть сле дующей: преобразование частоты, частотная фильтрация, «идеальное» ограничение, центрирование, нормирование, инвариантные преобразования [44], приводящие к доста точным статистикам.
Врезультате таких операций происходит редукция
выборочного /Угф/У и параметрического 0 ф 6 пространств ( Я £ Ф , а(=_Л). Дальнейшая редукция выборочного про
странства U осуществляется квантизатором.
Оценка потерь информации при таких преобразова ниях и определение оптимальных параметров аналогоцифрового преобразования могут быть выполнены с по мощью информационных критериев, рассмотренных в п. 1.7.4.
178
Оценка изменения среднего риска при редукции пара метрического пространства дается соотношением (1.129),
где w — ш(м,Л, a); P = w(u, Я ) = |
|"ОУ(«, Я, a) da; |
А — про- |
||||
странство иеизмеряемых |
параметров. |
|
|
|||
Пусть в результате редукции параметрического 9 и |
||||||
выборочного |
пространств |
U |
мы |
имеем |
распр^делени; |
|
|
|
—• |
|
|
|
|
дискретных |
переменных Р(и\1), |
обеспечивающ е |
допусти |
|||
мое увеличение средн?го |
риска |
(1.129). В |
соответствии с |
1.119) решающая функция у(«) должна быть построена на редуцированном выборочном пространстве О. При проверке гипотез при функции штрафов (1.128) решаю-
щая функция есть коэффициент правдоподобия у(й) =
= Л(гТ). В задаче оценки параметров при функции штра фов (1.131), равномерном распределении w(k) и симмет-
- >
ричном апостериорном распределении Р(Х\й) решающая функция определяется из уравнения максимального правдоподобия (3.1):
Р(и\Х) = тахР(и\ї).
Таким образом, в соответствии с уравнением (1.119)
мы приходим к задаче синтеза решающей функции у (її) на редуцированном выборочном пространстве П, т. е. к синтезу обнаружителей и измерителей после аналогоцифрового преобразования.
2.6.2. Синтез обнаружителей и анализ эффективности при независимых отсчетах. В результате аналого-цифро вого преобразования узкополосного процесса или неко
торых |
его параметров |
(В.1) — (В.5) |
получается |
последо |
||
вательность |
отсчетов |
дискретных |
переменных |
ЇЇі, |
t = l , |
|
2, ... , |
каждая из которых принимает значения |
xi, |
х%,... |
|||
...,хг. |
Если |
отсчеты статистически |
независимы, |
то |
веро |
|
ятность ра |
попадания |
случайной величины щ в интервал |
квантования Диа определяется в соответствии с (2.46). Независимость отсчетов соответствует ситуации, когда интервал временной дискретизации Та превышает время корреляции процесса u(t) или его параметров: TR^xit.
12* |
179 |