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

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

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

Комплексный

 

вектор принимает значения

 

и; =

(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

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