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

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

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

где {k} и {/} — множества, состоящие

из одного

индекса

/г и / соответственно, /, /гєМ.

 

 

Обозначим

 

 

 

 

 

Л

= 2

S

(gilSh)(l

-gilg*

(2.86)

в

* =

П

gjfgb)'

 

(2.87)

Перепишем формулу (2.85) в виде

 

 

 

л/2

і

 

 

 

да> (У)

 

 

 

(2.88)

 

 

 

 

Формула (2.88) выражает зависимость искомой плотно­ сти вероятности от величины накопленного напряжения в виде комбинации ^-распределений с двумя и четырьмя степенями свободы. Для получения зависимости вероят­ ности ложной тревоги от величины порогового уровня проинтегрируем выражение (2.88):

 

 

 

я/2

 

 

 

 

 

F = J и», (у) dy ••

 

• Л 1 1 + с / 2 ^ - Л ] , ( 2 . 8 9 )

где

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

Аи,

Bh

определены в

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

с (2.86)

и (2.87). Пользуясь

формулой (2.89), можно рас­

считать

зависимость

F{c)

при заданных п и

G2 ,

/г = 1,

2

п.

 

 

 

 

 

 

 

Практически целесообразным оказалось

вести

расче­

ты не по формуле (2.89), а по разложению плотности ве­

роятности

в ряд Эджворта

(см. приложение 1, алгоритм

3). Для этого определяются

кумулянты

плотности

(2.88)

по соотношению

(2.70). Имея в виду соотношение

(2.84),

получаем

 

 

 

 

 

 

 

 

л/2

п/2

 

 

«/2

 

 

и, =

4 £ gk,

K 2 = 8 S

gl,

к3 =

32 J

g3h,...,

 

 

fe=I

fe=l

 

 

A=.l

 

 

 

KJ = 2 ( / - 1 ) I 2 » X ;

^ .

 

(2.90)

При использовании алгоритма 3 приложения 1 и со­ отношения (2.90) была вычислена зависимость F(c). Ре­ зультаты расчетов приведены на рис. 2.15. Для оценки 120

погрешности разложения в ряд Эджворта отдельно вы­ числялся последний учитываемый в данной задаче член разложения

•В табл. 2.2 приведена относительная погрешность 6 = =R/F разложения в ряд при я = 20, 50. Кроме того, оцен­ ка погрешности разложения в ряд производилась 'путем

сравнения

с

зависимостью

F(c),

полученной

методом

Монте-Карло

иа Э Ц В М при и =6 0 и 7V = 2 • 105

повторе-

 

 

 

 

 

 

 

 

 

 

 

Т А Б Л И Ц А 2.2

 

 

 

 

 

 

 

 

F

 

 

 

 

11

 

10"»

 

 

10-з

ю-«

 

ю-»

ю-°

20

 

—5,76-10-»

 

—6,76-10"»

2,77-10-°

 

 

50

 

4,53-10"»

 

—1,86-10"»

1,4.10"°

 

2,1.10-'

2,65-10-'

ниях

опыта

 

(пунктирная кривая

на рис. 2.15). Из сопо­

ставления рис. 2.15

и табл. '2.2

видно,

что погрешность

разложения в ряд Эдж­

 

 

 

 

 

ворта

 

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

 

 

 

 

 

(2.88),

оцениваемая

мето­

 

 

 

 

 

дом Монте-Карло

при за­

 

 

 

 

 

ведомо

большом

 

числе

 

 

 

 

 

повторений

 

опыта,

не­

 

 

 

 

 

сколько выше, чем давае­

 

 

 

 

 

мая

таблицей. (Получение

 

 

 

 

 

точных

аналитических со­

 

 

 

 

 

отношений

для

вероятно­

 

 

 

 

 

сти

правильного

обнару­

 

 

 

 

 

жения

D (с, qz)

в

данной

Рис.

2.15.

Зависимость F(c)

задаче

затруднительно.

 

 

 

 

 

Поэтому в дальнейших рассмотрениях эта зависимость будет получена методом Монте-Карло по методике, опи­ санной в п. 6 приложения 1.

На рис. 2.16 приведены характеристики обнаружения весового и равновесного накопителей при бинарном кван­ товании и идеальной аналоговой обработке. Нефлуктуирующей цели соответствуют кривые /, шумоподобным

флуктуациям — 2,

совместным

флуктуациям— 3.

В табл. 2.3 приведены потери в пороговой мощности

обнаружителей с

равновесным

(2.61) и весовым (2.68)

121

суммированием по сравнению с обнаружителем (2.57). При этом было взято л=-50, F = 10- 6 , /Э = 0,5. Потеря в обоих случаях не превышают 1 дБ.

 

Т А Б Л И Ц А

2.3

Tim флуктуации

Потери в пороговой мощности,

дБ

Отсутствие

флуктуации

0,8

Шумоподобные

0,7

Совместные

 

1,0

2.4.Ц И Ф Р О В Ы Е Э К В И В А Л Е Н Т Ы

КО Г Е Р Е Н Т Н Ы Х О Б Н А Р У Ж И Т Е Л Е Й

2.4.1.Статистические свойства пассивных помех и эхосигналов. Процесс на выходе линейной части радиолока­ ционного лриемника описывается квазигармоническом функцией (В . 1); он может быть представлен совокуп­

ностью

значений

огибающей

E(t),

En(t),

фазы т|}(£),

ф п ( 0 'или квадратурных

составляющих A(t),

C(t).

Как

известно

(43,

57],

эффективная отражающая

поверхность цели

(ЭОП)

есть

случайная

величина

сэкспоненциальной плотностью вероятности. В этом

случае для пачки из п импульсов 2п-мерная плрт122

ность вероятности является .нормальной и имеет вид

% (z с ) = (2i.)"

 

dzr' M X

 

4

Ї ] 2

4

 

2nS

Х е х р

 

 

 

ГДЄ Z S = (Z b 2 2 , . . ., 2 2 п ) т = («'і, «Л,

u'%, ll"z, . . ., U'n, м " п ) т

2/г-мерный вектор-столбец с компонентами и'А и и"л, соот­ ветствующими А (і) и C(i/) для Л-го импульса; Мгий — корреляционная матрица порядка (2пХ2/г) с элементами

 

'"aft - 1 .

2 Z - i = Z 2

f c _ 1

Z 2

j _ , =

«'/l«'j, r2fc, гг

 

=u"ku"i,

->

raft-1. ai =

u'ku"i,

r 2 f t , 2u1=u"ku'i,

k,

l=\,n;

(2.92)

(a,, a3 ,

 

a „ n ) T

— вектор,

равный

 

математическому

a =

 

 

 

 

 

 

 

->

 

 

 

 

 

 

 

ожиданию

вектора

z 5 ;

det M 2 n S — определитель

матрицы

Показатель

экспоненты — положительно

определенная

квадратичная

форма

 

 

 

 

 

 

 

 

 

 

Й - а ] ) Л Г - ' 5

( z s - a s )

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= Е Е ГЦ (% - a i ) (zs, ~ ai) > °.

і 2 - 9 3 )

 

 

 

І=І

/=і

 

 

 

 

 

 

 

где

A f 2 / l S — невырожденная матрица ( d e t M 9 ( [

3 ^

0), имею­

щая

обратную

 

 

 

 

 

 

 

 

 

 

 

 

 

М^=Ш,

 

 

M2nSM^s

=

I;

'

(2.94)

/ — единичная (2/гХ2/г) матрица. Конкретный вид ма­ трицы M2ns определяется энергетическим спектром сиг­ нала и величиной допплеровского смещения средней ча­ стоты зондирующего сигнала. Можно показать (57], что элементы матрицы определяются из соотношений

rS (2Й - 1, 21—I) = r s (2k,2l)= C 0 S t^OS ~~ ^ X

 

X

J F s

(f) cos [2« (f-

f0)

(k

- I) Ta]

df,

(2.95)

 

 

о

 

 

 

 

 

 

 

 

 

r S

(2fc - l,

2/)

~ r S

(2ft, 2/-1)

=

S l 1 1

i^DS

~

I)] X

 

 

X]ps(f)

 

cos [2* ( f - f 0 )

( k - l )

Ta)

df,

(2.96)

 

 

о

 

 

 

 

 

 

 

 

 

где f D S

— & 0 3

Т п

=

4п^-Та

допплеровский сдвиг

фазы

123

за период повторения; vr—радиальная

составляющая

скорости цели; X — длина волны; Fs(f)

—энергетический

спектр сигнала, симметричный относительно промежуточ­ ной частоты fo-

Часто на практике ширина энергетического спектра эхо-сигналов составляет единицы герц. Поэтому за вре­ мя облучения цели пачка эхо-сигналов практически по­

стоянна,

т. е. полезный сигнал является неслучайной

функцией

времени

(нефлуктуирующая цель):

 

 

S(t)=Es(i)

COS_(MO + Q B S K + <PO].

(2.96а)

Пассивная і(коррелированная) помеха представляет собой суперпозицию большого числа независимых сла­ гаемых, т. е. является случайным гауссовым процессом. В теоретических исследованиях обычно предполагают, что она стационарна ![42]. Многомерная плотность вероятно-

—>

сти помехи W2n(zP), а также элементы корреляционной матрицы описываются соотношениями, аналогичными

(2.91),

(2.95), (2.96).

 

 

 

В радиолокационном приемнике на помеху наклады­

ваются

его собственные

шумы, поэтому

корреляционная

матрица смеси помехи

и шума

равна

M2npN = M2np +

+ M2nN,

а элементы матрицы шума

 

 

 

 

Гз2 при k = l,

 

^ ( 2 * - 1 . 2 / - 1 , = > , 2 * , 2 / ) - | 0

п р и

 

 

'N{2k—\,2l)

'N (2ft, 2/—1)

и <

 

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

^ . н - W , (2.966)

где

р2 =

а2 2 ,

а суммарная

мощность

о2 э Л ,°=а2 -(-о2 .

 

Выражение

для R N P ( T T I )

получается

с учетом того,

что

для

функций корреляции

имеем

соотношения:

 

 

 

ВР^(ТП)=ВР-(Ти)+В(Тп),

 

 

 

 

 

BN (Ta) = fRN а),

Bv

а) = o2PRV

а)

124

где о2 +стр2 мощность смеси помехи и шума. Статисти­ ческие свойства реализаций смеси сигнала с помехой также описываются распределением типа (2.91) с матри­ цей MznSy

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

соотношения [42, 43, 67]:

 

 

 

для сигнала — резонансный

спектр

 

 

+ 2

f>-:f.

—і

(2.97)

 

для помехи — гауссов спектр

 

 

/Ч/)г = е х р [ - 2 , 8 ( ^ - )

 

(2.98)

Соответствующие коэффициенты корреляции имеют вид:

 

« ( х ) р = ехр [ — иД/1х/

],

(2.99)

 

едг

= е х р [ - ^ ( Д Н 2

] .

(2.100)

При

этом элементы

корреляционных

матриц помехи

и сигнала с помехой с учетом

(2.95),

(2.96) имеют вид:

для резонансного спектра

 

 

 

 

=cos[fD(k~l)exp{—nkfTa

|

k—l\},

 

 

 

 

 

 

(2.101)

 

= sin [<?D

(k -

I)] exp { -

«bfTa\k~l\j,

для

гауссова спектра

 

 

 

=

cos [<fD( k -

I)} exp { - ^

[AfГп (A -

/)]*},

 

 

 

 

 

 

(2.102)

=

sin [<?D (k -

/)] exp I

[ A f T a

(k -

/)]2}.

125

2.4.2. Цифровые обнаружители нефлуктуирующих ко­ герентных сигналов. Синтез структуры оптимального или достаточного обнаружителя пачки когерентных эхо-сиг налов импульсной РЛС обзора на фоне пассивных помех осуществляется в соответствии с методикой .п. 2.2.1, ана­ лиз эффективности — п. 2.2.2, построение многоканальной

системы — п.

2.2.3.

Коэффициент

правдоподобия

для

пачки нефлуктуирующих когерентных сигналов в

смеси

с пассивной

помехой

записывается

с учетом распределе­

ний типа (2.91).

 

 

 

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

где u'j, u"j, / = 1 , /г, — последовательности выборочных значений на выходах фазовых детекторов двух квадра­ турных каналов:

п

K't

=

£

^iESi C 0 S iQDS (ІТ*)

+

Vol

 

 

 

;=i

 

 

 

 

K

-

= t

?i}Eslsm[QDS(iTa)

+

<t0]

(2.104)

— последовательности весовых коэффициентов, опреде­ ляемых как свойствами помехи (элементами гц матрицы М~'р ), так и радиальной составляющей скорости цели

(2.96а). Как следует из (2.103), достаточный приемник является двухканальным. В каждом из каналов осуще­ ствляется линейная обработка.

•Ближе к практике является обнаружение пачки коге­

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

неизвестной начальной

фазой.

В этом случае образуется

статистика

 

y =

U'* +U"2

(2.105)

126

где

пп

(2.106)

Г л

гп

(2.107)

h , , h определяются из (2.104) при подстановке <р0 =0.

Достаточный приемник (2.106), (2.107), осуществляю­ щий фильтрацию, является четырехканальным.

Техническая реализация алгоритма (2.105) с учетом (2.106) и (2.107) сложна как в аналоговом, так и в циф­ ровом исполнении из-за многоканальности системы по дальности и радиальной скорости движения цели. Воз­ можен переход к многоканальной по частоте Допплера системе, когда каждый канал настроен на определенное значение неизвестной частоты. При этом в импульсной когерентной Р Л С обзора с детерминированным периодом повторения T^ — const структура оптимального приемника определяется в соответствии с (2.4), так как однознач­ ное измерение радиальной составляющей скорости цели невозможно. При псевдослучайной модуляции периода повторения 7 \ = 7 " O ± A 7 V ( C M . п. 2.2.3) Р Л С обладает раз­ решением как по дальности, так и по радиальной со­ ставляющей скорости и структура оптимального прием­ ника определяется из (2.6). Пример квазиоптимальной системы с разрешением по дальности и скорости рас­ смотрен в п. 2.2.3.

Вследствие сложности технической реализации алго­ ритм, определяемый соотношениями '(2.105) — (2.107), на практике конкретизируется на случай сильно коррелиро­ ванной помехи ( ^ р ( Г п ) ~ 1 ) , а радиальная составляющая скорости цели принимается равной оптимальной vr=

— Vrovt, т. е. допплеровский сдвиг фазы за период повто­ рения

Ф д 5 = ± 2 ( 6 — 1 ) Я , fe = 0, 1, 2, . . .

(2.108)

При этом считается, что скорость помехи, если она дви­ жется, скомпенсирована системой компенсации скорости ветра, т. е.

фср = 0:

(2.109)

127

В этом случае алгоритм оптимальной обработки сво­ дится к многократному череспериодному вычитанию (ЧПВ) 1[42, 69]

 

 

S ^ ( - 1 ) V ,

V +

+ І

f s ^ t - D V

' ,

V = і ;

(2.1Ю)

/=/( + 1

\v=0

/

/ = К + 1

 

где К — кратность вычитания; С"к — число сочетаний из К элементов по v. В сущности величины (—1)V C^ представ­ ляют собой знакопеременные биномиальные коэффициенты.

Достаточный приемник, выполняющий фунхияи по­ давления (режекцнн) помехи (2.110), хорошо известен

 

Sip <р

п

 

 

 

 

ОЗУ

Блок обра­

 

 

ботки.

 

 

 

u(t)

\sinco0t

гА Синхро­

 

ЦК

 

\C0StDst

низация

 

 

 

 

 

 

п

 

 

 

 

ОЗУ

Блок

обра.- _

 

 

ботка

 

 

 

 

COStp

 

 

 

Рис. 2.17.

Цифровая система многократного

череспериодного вычи-

 

 

'танип.

 

 

и в аналоговом исполнении реализуется с помощью уль­ тразвуковых линий задержки (УЛЗ) или вычитающих потенциалоскопов [43, 63, 69, 70]. Многократное ЧПВ воз­ можно как на промежуточной, так и на видеочастоте. Б последнем случае снижаются требования на стабиль­ ность УЛЗ, но система содержит теперь два квадратур­ ных канала.

Часто перед фазовыми детекторами квадратурных каналов ставится «идеальный» амплитудный ограничи­ тель с целью устранения амплитудной модуляции пачки эхо-сигналов диаграммой направленности антенны при сканировании последней и уменьшения динамического диапазона сигналов. Система с ограничителем, в сущно128

сти, являющаяся фазовой (амплитудная информация устранена ограничителем), имеет и преимущества и не­ достатки перед системой без ограничителя.

Рассмотрим возможную техническую реализацию цифрового эквивалента алгоритма (2.110) на видеоча­ стоте (рис. 2.17), когда кратность вычитания К—п—1. Для системы без ограничителя выходное напряжение усилителя промежуточной частоты поступает на входы двух фазовых детекторов, опорные напряжения которых есть sin mt и cos ©о4 Видеонапряжения обоих фазовых детекторов подаются на входы преобразователей напря­ ж е н и е — код, принципы построения которых рассмотрены в монографии '[2]. В системе с ограничителем в качестве преобразователя может использоваться дискретный фазо­ метр совместно с блоками, вычисляющими sin ip и cos ip. Последние могут быть выполнены на основе постоянных запоминающих устройств. При бинарном квантовании фазы можно использовать преобразователь, аналогичный преобразователю, представленному на рис. 3.23,а.

Период следования импульсов временной дискрети­ зации 7д, подаваемых в преобразователь, равен длитель­

ности

зондирующего импульса

to. В каждом

квадратур­

ном

канале используется

оперативное

запоминающее

устройство

(ОЗУ) . Параметры

ОЗУ должны

быть

таки­

ми: время

цикла Т ц ^ Г ц ,

число ячеек равно

числу

даль-

номерных

каналов

M = TJT}1,

число двоичных разрядов

в ячейке — пт, где

п — число

импульсов

в

пачке,

т —

число разрядов в выходном коде аналого-цифрового пре­ образователя.

Квадратурные каналы построены одинаковым обра­ зом, поэтому рассмотрим лишь один. В ОЗУ хранится и постоянно обновляется информация, соответствующая всем М дальномерным каналам. Обновление информации в каждом канале (ячейке ОЗУ) сводится к ее сдвигу на одну из п позиций в сторону возрастания их номеров и записи нового /п-разрядного кода в освободившуюся первую позицию (режим «движущегося окна»). Обраще­ ние к ячейкам ОЗУ производится в порядке возрастания их естественных номеров, с первой по М-ю.

Для обработки информации каждого дальномерного ка­ нала используется одно и то же устройство, реализую-

л—1

 

щее операцию свертки У] Киі—,< г д е j=K~\-l,

п. Оно

«=о

9—1410

129

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