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