
Литература / Гришин Ю.П., Ипатов В.П., Казаринов Ю.М. Радиотехнические системы (1990)
.pdfется в более сильных от ражениях от близких посто ронних объектов (построек, морской поверхности, эле ментов рельефа, специально устанавливаемых ложных целей); в дальней радиона вигации используемый для местоопределения слабый сигнал поверхностной вол ны маскируется накладыва ющимся на него более силь ным, отраженным ионосфе рой, и т. д. Поэтому так актуальны попытки избе
жать осложнений, связанных с наличием у ФН Т(т) ЛЧМ-импульсов больших боковых лепестков. Определен ного снижения уровня последних можно добиться, приме няя ЛЧМ-сигналы с непрямоугольными, гладкими (типа колоколообразной) огибающими, переходя к нелинейной частотной модуляции либо заменяя согласованную фильт рацию специальной весовой обработкой.
Дискретные и фазоманипулированные сигналы. Обра тимся к дискретным сигналам, т. е. к последовательностям регулярно повторяющихся радиоимпульсов одинаковой формы и центральной частоты, отличающихся друг от друга лишь значениями комплексных амплитуд. Дискрет ные сигналы, содержащие конечное число импульсов N, уже анализировались в § 3.4, где были названы пакетами. Модель ^импульсного дискретного сигнала может быть
эоггттг>аттп |
т> тэттгта Гг>гч |
. |
г» |
('Х |
-Jfу |
('Х |
1 |
|
|
|
V |
|
|
) J |
|
||
"N- 1 |
|
|
|
|
x -1 |
|
|
|
s(z) = Re |
E d,j0(z —z'Tn) |
|
= Re |
x df50(z-zTn}exp(/’27t/oO |
||||
N- 1 |
‘■ = 0 |
|
|
|
|
i = 0 |
|
|
|
|
|
|
|
|
|
|
|
= E l«i|50(r-«rn)cos[2K/0Z + Y0(z-zrn) + (pi], |
(6.10) |
|||||||
1=0 |
|
|
|
|
|
|
|
|
где .y0 (z) = S0 (z)exp (y2n/oz)— аналитический сигнал, реаль ная часть которого ,y0(z) = Re [j0(z)] описывает отдельный
(одиночный) импульс пакета; 50(z) = 50(z)exp [jYo(z)] — комплексная огибающая одиночного импульса s0 (t), учиты вающая форму его действительной огибающей So (z),
а также закон |
внутриимпульсной угловой модуляции Yov)i |
|
Т„— период |
повторения импульсов в |
пакете; а; = |
= |«,|ехр(у<р,) — комплексная амплитуда |
Z-ro импуль- |
160
са, модуль |а;| и аргумент <p,= arga; которой задают
действительную |
амплитуду |
|
||
и начальную фазу z-ro им |
|
|||
пульса |
пакета. |
Набор |
|
|
i = 0, 1, |
N— 1, |
устанавли |
|
|
вающий |
закон |
изменения |
Рис. 6.7 |
|
амплитуд и фаз |
дискретного |
|||
|
||||
сигнала |
от импульса к импульсу, называют кодовой |
последовательностью или кодом. Пример дискретного сиг нала длиной (числом импульсов) N = 4 для прямоугольного одиночного импульса и кода zi0=l, a1=2j, л2=—2, а3= — j приведен на рис. 6.7.
Вычислим корреляционную функцию ф(т) комплексной огибающей
5(0 = Nf1«i5o('-/rn) |
(6.11) |
i = 0
дискретного сигнала (6.10). Подставив выражение (6.11)
в (6.4), |
найдем |
|
|
= |
i. 1 = 0 |
f 5o('-'Tn)5o(^-T-/rn)dr, |
(6.12) |
|
-оо |
|
|
гцеЕ = |
N -1 |
1 00 |
• |
|л;12-£0—энергия всего пакета; Ео=- f |
|S0(r)|2dr — |
энергия одиночного импульса .v0(z). Введя корреляционную
функцию ф0 (т) комплексной огибающей |
50 (г) одиночного |
||||
импульса |
|
|
|
|
|
• |
1 |
® |
• |
• |
(6.13) |
|
= |
f |
S0(t)S'0(t-T)dt, |
||
перепишем |
(6.12) |
в виде |
|
||
• |
Е |
х |
|
• |
(6.14) |
'l/(T) = “F |
Z |
aia^0(T:-(i-l)T„\, |
Е i.l- -да
условившись считать нулевыми те а;, а1 в последней
сумме, |
индексы которых отрицательны или превышают |
|||||
N— 1. |
Заменив индекс суммирования / |
в (6.14) на m = i—l, |
||||
получим |
Ф(т)=-^ |
£ |
£ |
а,т^-тфо(т-шТп) = |
||
ту |
оо |
оо |
I— — оо |
т— — оо |
|
|
|
|
|
|
|||
= -~ |
Ё |
Ё а,а--тф0(т-теТп). |
|
|||
т — — гл i-~ со |
|
|
|
[т) кодовой после |
||
Введем корреляционную |
функцию |
|||||
довательности й0, |
...» |
|
|
6 Заказ 3173 |
161 |
|
Е |
00 |
|
|
|
|
(6.15) |
|
= |
X “^-т, |
|
|
|||
|
|
i= -00 |
|
|
|
|
|
где по-прежнему все d„ |
d;_„, с индексами |
вне диапазона |
|||||
О, 1, ..., N— 1 |
считаются |
нулевыми. Тогда |
соотношение |
||||
(6.14) |
можно |
записать в |
виде |
|
|
||
vi'(T)= |
f |
|
|
тТ„). |
|
(6.16) |
|
|
m = — со |
|
|
|
|
|
|
При m^O d,d^-m = O |
для всех i<m и |
V^-N, поэтому |
|||||
подробная |
расшифровка |
записи (6.15) такова: |
|
||||
|
|
|
г n-i |
|
|
|
|
|
|
|
/ t |
О. i — |
|
|
|
|
|
|
E i = m |
|
|
|
|
|
|
|
|
|
|
|
(6.17) |
Последнее из соотношений (6.17), обобщающее верхнее |
|||||||
на случай произвольного целого т, следует из |
(6.15) |
||||||
после замены индекса суммирования i на |
i + m. |
пока |
|||||
Название |
(ш) |
отражает смысл этой величины, |
зывающей, насколько близки, коррелированы между собой кодовая последовательность d0, а1г ..., aN_ х и ее копия, сдвинутая на т позиций. Для вычисления, например, фа(2) нужно под исходной кодовой последовательностью запи сать ее сдвинутую на две позиции вправо и комплексно-со пряженную копию:
«о, ах, аг, а3, ..., d.v-i do, dj, ..., *yd_ 3
и, перемножив попарно стоящие в каждом столбце этой записи элементы, просуммировать полученные произведе ния, после чего умножить результат на нормирующую константу Ео/Е. При прямом вычислении (— 2) сдвинутая копия кода должна быть смещена относительно исходной также на две позиции, но влево:
«О, «1, •••, Hn-3> aN-2> aN-l |
|
||
i,*d *аз, |
..., d‘N_ i. |
|
|
Ясно, |
что |
фа(га) = О при |
m^N и m^-N [см. (6.17)J |
потому, что |
сдвинутая на |
такое число позиций кодовая |
162
последовательность не перекроется с исходной. Окон
чательно |
для |
корреляционной функции |
(6.16) получим |
Ф(т) = |
N-1 |
Фа("»)Фо(т — тТп). |
|
X |
(6.18) |
т= -(N- 1)
Отсюда видно, что корреляционная функция комплексной огибающей дискретного сигнала (6.10) есть сумма повто ряющихся с интервалом Т„ корреляционных функций ком плексной огибающей одиночного импульса, взвешенных кор реляционной функцией кодовой последовательности
Так как ф0 (т) воспроизводит комплексную огибающую одиночного импульса 50 (z), пропущенного через согласован ный с ним фильтр, то длительность т10 ф0(т) не превышает удвоенной длительности одиночного импульса ти:тк0^
^2ти^2Гп. Таким образом, корреляционная функция ф(т) (6.18) комплексной огибающей пакета (6.10) представляет собой пакет из 2(V— 1 повторяющихся с интервалом Тп импульсов вида ф0(т) с длительностью не более 2ГП, причем комплексная амплитуда ш-го импульса этого нового пакета равна значению корреляционной функции кодовой последовательности дискретного сигнала (6.10) при m-м сдвиге.
Нетрудно понять, какое необходимое и достаточное условие должно быть выполнено для того, чтобы дискрет ный сигнал (6.10) сжимался в СФ. Если бы кодовая последовательность была некоррелирована со всеми своими
сдвигами |
[фа('») = 0, ш^О], то в |
сумме |
(6.18) |
осталось |
||
бы единственное ненулевое слагаемое с |
т = 0 |
[фо(0) = |
||||
=(Е0/Е) |
1ф|2 = Е/Е=1 |
и корреляционная функция |
ф(т) |
|||
|
/ = о |
сигнала |
(6.10) |
имела |
бы |
вид |
комплексной огибающей |
одиночного импульса ф0 (т) длительности 2 Гп. Отсчиты вая эту длительность по некоторому условному ненулевому уровню, можно в первом приближении положить ткй7'п. В то же время длительность сигнала (6.10) Tc~NTn, так
что при выполнении перечисленных требований к \|/a(w) произошло бы сжатие сигнала в N раз. Так как ширина спектра сжатого сигнала, а следовательно, и самого пакета не меньше 1/тк«1/Гп, то база дискретного сигнала (6.10) длительности TcxNTn B = AfcTc^ N, т. е. не меньше длины кодовой последовательности N.
Так как первый и последний импульсы в (6.10) по определению имеют ненулевые амплитуды do/0, aN^1^0, то 1))^0 и обращение
6* |
163 |
в нуль фа(га) при всех ненулевых т заведомо невозможно. Поэтому функция (6.18) кроме основного пика будет иметь и боковые, уровень которых определяется значениями фа(/и) при т^=0. Для того чтобы сделать максимальный боковой лепесток ф(т) приемлемо малым, необходимо добиться минимума максимального по всем ненулевым т уровня корреляции фамакс кодовой последовательности со своими сдвигами: фамакс = тах|ф0(ш)| = т1п.
*0т
На практике элементы а, кодовых последовательностей выбирают не произвольно, а из некоторого заранее оговоренного множества (алфавита кода) по возможности небольшого объема. Подобное ограничение, как правило, отражает желание максимально упростить формирование и обработку дискретных сигналов. По этим соображениям предпочтение нередко отдается двоичным фазоманипулированным сигналам, т. е. дискретным сигналам с действи тельными элементами кода, принадлежащими двоичному алфавиту {— 1, + 1}. При а, = ± 1 амплитудная модуляция импульсов в (6.10) отсутствует, так что энергия двоичного ФМ-сигнала равномерно рассредоточена в пределах его длительности, что обычно и требуется от сложного сигнала (см. § 6.3). В то же время начальные фазы радиоимпульсов такого сигнала могут принимать только два значения: 0 и л. что, в свою очередь, заметно упрощает и удешевляет
соответствующие схемотехнические |
средства. |
|
||
|
|
N-1 |
|
|
При |
|д,.| = 1 |
EjE0= £ |а;|2 = М |
фа(АГ-1) = |
|
= аолл_X/Ar= ± 1/jV, |
i = 0 |
|
|
|
вследствие чего максимальный боко |
||||
вой лепесток |
фамакс |
функции ф0(га) |
[наибольший уровень |
фа(ш) при /н#0] не может быть меньше 1/N. Таким образом, среди двоичных ^элементных кодовых после
довательностей а0, , ..., |
aN_ j |
с d;=±l, f= 0,1, ..., — 1, |
лучшими следует считать |
те, |
для которых фамакс = 1/N. |
Двоичные последовательности с указанным свойством, известные под названием кодов Баркера, существуют только
при N=2, |
3, 4, |
5, 7, 11, 13. Рассмотрим в качестве |
примера |
код |
Баркера длины W=7. Для него |
<70 = а, = а2 = as = 1, |
а3 — а4 = а6 = — 1, т. е. дискретный сигнал, |
манипулированный таким кодом, содержит семь радиоим пульсов с начальными фазами (ро = <р1=<р2 = (р5 = 0, <р3 = <р4 = <р6 = л. Такой сигнал и его комплексная огибающая
показаны на рис. 6.8, а, б, причем |
на рис. 6.8,6 знаки « + » |
и « —» отвечают комплексным |
амплитудам a, = ej0=l |
и = е-"' = — 1. Следуя данному ранее рецепту вычисления
164
Рис. 6.8
фа(ш), составим таблицу, в которой под кодом Баркера
запишем его |
сдвиги на одну, две, три, четыре, пять, |
шесть позиций |
(габл. 6.1). |
|
Таблица 6.1 |
Перемножая элементы исходного кода (строки с ш = 0) со стоящими под ними элементами т-го сдвига кода и складывая произведения, получим числа, стоящие в по
следнем |
столбце таблицы. |
После деления |
на |
Л'= 7 |
это |
||
и |
будут |
значения 'ta(w). |
Обращаясь |
к |
(6.18) |
и |
имея |
в |
виду, |
что корреляционная функция |
^0(т) одиночного |
прямоугольного видеоимпульса длительности Та есть рав
нобедренный треугольник с основанием 2Т„, приходим
165
к выводу, что для двоичного ФМ-сигнала на основе кода Баркера комплексная огибающая сжатого сигнала ф(т) имеет вид основного треугольного пика и «отрицательных» боковых лепестков той же формы, в семь раз меньшей высоты с вершинами в точках т=+2Гп, ±4ГП, +6ГП (рис. 6.8,в).
Одна из возможных структур СФ для двоичного ФМ-сигнала показана на рис. 6.9, а (для примера взят прежний семиэлементный код Баркера). Сигнал проходит через N~ 1 последовательных элементов задержки на время Т„, выходы которых подключены к сумматору с весами
+1 или —1, взятыми в порядке, зеркальном по отношению
кпорядку следования dt в коде. В результате на N входах
166
сумматора при поступлении на вход СФ сигнала, показан
ного |
на рис. 6.9, б, появляются |
его сдвиги, |
причем те из |
них, |
которые умножаются на |
коэффициент |
—1, меняют |
начальные фазы импульсов на л (рис. 6.9, в, где нумерация эпюр соответствует номерам входов сумматора на рис. 6.9, а). В результате с выхода сумматора снимается колебание (рис. 6.9, г), каждый прямоугольный импульс которого после обработки в согласованном с одиночным импульсом фильтре СФОИ принимает вид треугольника (рис. 6.9, д). Получаемый в итоге сигнал имеет огибающую
вида Р(т)* |
= |ф(т)| (ср. |
с рис. 6.8,в). |
|
Несуществование |
кодов Баркера для |
14 стимули |
ровало поиски двоичных ФМ-сигналов, обладающих доста точно малым (хотя и большим 1/7V) уровнем бокового ле пестка при больших длинах N. Таких сигналов известно сейчас довольно много, хотя вопрос об их оптимальности и существовании лучших кодов не решен. К числу после довательностей с особенно интересными корреляционными свойствами относятся двоичные М-последовательности, существующие для любых длин вида N=2"—1 (л — нату ральное число). Чтобы проиллюстрировать достоинства этих последовательностей, введем еще одно общее понятие.
на |
Если |
любую |
последовательность |
длины N сдвигать |
||
т позиций не |
так, |
как |
прежде, а |
циклически, т. е. |
||
/?„,и 5 |
••П.1 ’ |
Z4-’ |
,> |
-• /V — 1 |
’ |
|
-я*2> |
*й-1,. |
|
..., |
*а-з, |
|
|
то корреляция кода с его т-м циклическим сдвигом будет определяться несколько иначе, чем в (6.17): в соответст вующей сумме при любом т всегда будет точно N слага емых. Корреляционные свойства последовательностей при циклических сдвигах характеризуют периодической корре ляционной функцией
Е N~1 |
(6.19) |
|||
Фап(™) = -^ |
X |
|||
£ |
| = о |
|||
|
|
|||
где обозначение |
((i — т)) символизирует вычитание |
по |
||
модулю N, |
т. е. |
взятие остатка от деления i—m иа |
N. |
На практике дискретные сигналы иногда передаются не изолированными пакетами, а повторяются с интервалом NTn, образуя бесконечную периодическую последователь ность манипулированных радиоимпульсов: спустя Гп секунд
167
после импульса с комплексной амплитудой следует импульс с комплексной амплитудой а0 и сигнал (6.10) повторяется. В этих случаях важно, чтобы малый уровень
боковых лепестков имела |
функция фоп(м), а не ф0(/и). |
Так как функция фоп(/и) |
периодична по т с периодом |
N [((/-m))=((z-m±2V))], |
то фап(^) = фап(О)= 1 и жела |
емое свойство фап(ги) можно сформулировать как равен ство нулю ее значений при всех т, не кратных N [под
боковыми лепестками ф„п(га), |
и понимают |
ее значения |
|
при всех т, не кратных |
N ]. |
|
|
Многочисленными исследованиями установлено, что для |
|||
двоичных последовательностей |
с элементами |
+1, — 1 это |
|
требование невыполнимо |
(за |
исключением |
тривиального |
случая N=4). Поэтому М-последовательности (наряду с не которыми другими) являются оптимальными среди двоич ных, так как их боковые лепестки [фап(ш) = — 1 /W при всех не кратных N значениях т ] имеют минимально возможный при нечетных длинах уровень. Что же касается «Пакетных» (не периодических) свойств таких кодов, характеризуемых кор реляционной функцией (6.17), то максимум бокового лепес тка фама1С корреляционной функции фа(т) для М-последова-
тельностей имеет значение, близкое к |
что считается |
|
вполне приемлемым, ибо |
пока никаких |
двоичных кодов |
с существенно меньшим |
уровнем фомакс |
неизвестно. |
Отобранные здесь в качестве примеров виды сложных сигналов не должны рассматриваться как самые типичные для современных РТС. Повсеместное упоминание именно этих сигналов в учебной литературе связано с их большей мето дической наглядностью по сравнению с другими. В реальных же системах помимо рассмотренных находят применение и другие виды сложных сигналов. Так, для получения нулевых боковых лепестков фоп(ш) вместо двоичной часто применяют р-ичную фазовую манипуляцию либо двоичную ФМ дополняют введением пассивных пауз. Ценными для ряда приложений свойствами обладают сигналы с манипу ляцией частоты, импульсно-временным кодированием и т. п.
§ 6.5. РАЗРЕШЕНИЕ ПО ВРЕМЕНИ ЗАПАЗДЫВАНИЯ
И ЧАСТОТЕ. ЧАСТОТНО-ВРЕМЕННАЯ ФУНКЦИЯ НЕОПРЕДЕЛЕННОСТИ СИГНАЛА
В локационных системах интерферирующие на входе приемника сигналы, отраженные различными целями, обыч но отличаются друт от друта не только временем запаз-
168
<F(T,F)
дывания, но и доплеровскими сдвигами. В таких случаях, характерных и для других приложений (радионавигация, связь по многолучевым трассам и др.), приходится гово
рить о разрешении |
сигналов по времени запаздывания т |
и по частоте F, т. е. |
по двумерному векторному параметру |
X, компонентами которого служат т и F:1 = (t, F)t. Качест во разрешения при этом определяется видом частотно-вре
менной ФН, |
введенной в § 5.2 |
[см. (5.16)], т. е. |
|
|||
1 |
” . |
|
|
|
(6.20) |
|
»Р(т, F) = — f |
(/5(?)5* -т)ехр(-у2лГг)А |
|
||||
2F _ж |
|
|
|
|
|
|
Геометрически 'Р(т, F) задает некоторую поверхность |
||||||
над координатной |
плоскостью |
т, F, причем |
в |
начале |
||
координат т = 0, F=0 высота этой поверхности фиксирована |
||||||
и равна единице: |
Т (0, 0)=1 |
(рис. 6.10, а). |
Две |
копии |
сигнала, сдвинутые друг относительно друга по времени запаздывания на т и по частоте на F, разрешить тем легче, чем ниже уровень ФН (6.20) при данных т, F, т. е.
чем меньше высота поверхности на рис. 6.10, а |
в точке |
|
с координатами т, F. Таким образом, желательно, чтобы |
||
ФН Т(т, F) как можно быстрее спадала по мере удаления |
||
точки т, F от начала координат. Нужно отметить, что |
||
сечение Т(т, 0) ФН Т(т, F) плоскостью |
F=0, как |
следует |
из (6.20), есть ФН по запаздыванию |
Т(т) [см. |
(6.3)] и, |
следовательно, протяженность Т(т, 0) по оси т харак теризует достижимую для данного сигнала разрешающую способность только по времени запаздывания (т. с. воз можность разрешения двух копий сигнала, имеющих разное время запаздывания, но одинаковые частоты). Аналогич
но, |
протяженность |
вдоль |
оси |
F сечения Т(0, F) |
|||
ФН |
(6.20) |
плоскостью |
т = 0 |
определяет |
разрешающую |
||
способность |
только |
ио |
частоте |
(когда |
разрешаемые |
169