книги / Широкополосные дискретно-кодированные сигналы в радиотехнике и радиолокации
..pdfТаким образом, сравнительный анализ ФН ДКС и ее сечений дает возможность отдать предпочтение наиболее перспективным сложным сигналам, позволяющим обеспечить высокие энергетические возможно сти РЛС и требуемое совместное разрешение целей по дальности и ско рости при малом УБЛ.
При выборе зондирующего сигнала необходимо учитывать еще одно важное свойство, связанное со скрытностью работы РЛС. В этом отношении ДКЧС Костаса дают возможность обеспечить высокую скрытность работы РЛС, так как общее число перестановочных матриц для выбранного массива N равно N1, и при увеличении размерности ко да существенно повышается помехозащищенность РЛС.
Принимая во внимание эти привлекательные возможности и учи тывая, что ДКЧС не нашли по сравнению с другими сложными сигна лами достаточного отражения в литературе, следует рассмотреть их свойства более подробно.
2.3. Свойства функции неопределенности дискретно-кодированных по частоте сигналов
Дискретно-кодированным по частоте сигналом будем называть дискрет ный частотный сигнал первого порядка [3] при следующих условиях:
амплитуда апэлементарных импульсов в формуле (2.4) постоянна и равна единице, а начальные фазы срп - постоянны или равны нулю;
дискрет частоты /„ в каждом элементарном импульсе длительно
стью Т равен /„ = Bnb f , где d f - значение дискрета сетки частот сиг
нала, 9п - элемент частотно-временной матрицы сигнала, определяю щий способ кодирования частоты;
код частоты {&„} ДКЧС размерности N должен содержать все зна
чения от 1 до N. |
|
Базу ДКЧС размерностью N представим в виде |
B = FNT, где |
F - полоса частот, занимаемая сигналом около несущей, равная Д/Г, а |
|
NT - длительность сигнала. |
|
Тогда значение базы можно определить как |
|
В = AJNNT = N2AJT = N 2M , |
(2.5) |
где М = А/Т - масштабный коэффициент.
Повышение помехозащищенности и скрытности работы РЛС на прямую связано с увеличением базы зондирующего сигнала и количест вом возможных вариантов его кодирования.
Рассмотрим свойства ФН ДКЧС при различных значениях базы
сигнала В = N2M
Для этого воспользуемся результатами расчета на ЭВМ трехмер ного тела ФН, используя формулу (2.4). Для ДКЧС код амплитуды {аг} элементарных импульсов постоянный и равен единице, поэтому пара метр Z в (2.4) равен 1IN. Наличие обратно пропорциональной зависимо сти от N приводит к снижению УБЛ ФН с ростом N, зона пьедестала расширяется и, таким образом, соблюдается требование постоянства объема трехмерного тела ФН. Кроме того, необходимо выявить зависи мость ФН от значения параметра М.
На рис. 2.11, 2.12 и 2.13 представлены результаты расчета ФН ДКЧС с размерностью N= 11 дляМ= 1 (база сигнала В = N2), М - 0,8 и М = 1,2 соответственно. Выбор сравнительно небольшой размерности N объясняется ограниченной разрешающей способностью графиков и не накладывает ограничений на общность результатов. Как видно из этих рисунков, при М <1 происходит сжатие пьедестала ФН по оси допле ровских частот, что приводит к увеличению УБЛ в сечениях, парал лельных оси задержки около нулевых доплеровских частот, а приЛ/>1 зона пьедестала по оси частот расширяется, что приводит к уменьше нию УБЛ в этой области.
Для анализа разрешающей способности ФН ДКЧС по задержке и частоте рассмотрим вид ФН в сечениях при нулевой задержке (х=0) и ну левом доплеровском сдвиге отраженного сигнала (у=0). Для этого преоб разуем выражения Х\ и Х2 в формуле (2.4) к виду, более удобному для аналитического исследования. Подставляя F] и F2 в экспоненциальных частях X] иХ2соответственно и приводя подобные члены, получаем:
s in Q f i( l - i) ]
X
|
irFx(\- x ) |
|
|
xexp{[jxly(/c + r) + x+l) + (0r(l-£)-6>k+r(l + x))M] + j(<pr -<pi+r)}; |
(2.6) |
||
X |
sinH~2*].. |
|
|
2 |
7Lp2^ |
|
|
Xexp{jЯ-[y(2(k + r +1) + x) + (2вг - (0r + Ok+r, x) x ) M |
?>t+r+1)}. |
(2.7) |
Тогда, учитывая, что при х=0 к и х равны нулю и используя (2.6), выражение, определяющее ФН ДКЧС при нулевой задержке, примет ввд:
(2.8)
Из (2.8) следует, что %(0,у) , в общем случае, зависит только от
размерности ДКЧС N и кода амплитуды. Для случая ДКЧС код ампли-
туды {аг} постоянный и равен единице, а код фазы {(рг) - равен нулю. Тогда, используя соотношение
YVmr |
iW-IW2 Sin(AW2) |
|
(2.9) |
|
sin(m/2) ’ |
|
|
|
|
|
|
форму лу (2.8) можно представить в виде |
|
||
|
sin\жЫу\ |
|
|
х(о,у) = |
exp(j^A^). |
|
(2. 10) |
жЫу |
|
||
Полученное выражение для |
>у) |
показывает, что для ДКЧС |
размерностью N разрешающая способность по частоте и УБЛ в сечении при х = 0 соответствует прямоугольному высокочастотному импульсу длительностью NT, и чем больше NT, тем выше разрешающая способ ность по частоте.
Рассмотрим сечение ФН вдоль оси нормированной задержки х = г/Т. Подставив (2.6). (2.7) и выражения для F\ и F2 в (2.4), для ДКЧС при
у = 0, получим: |
|
N -\-k sin[a~(flr -в к+г)М (l-x )] |
|
0 - i ) £ |
X |
r=0 |
|
xexp {j *■ [в, (1 - x) - вк+г (l ♦ x)] M } - |
|
Ы -2 -к |
|
r=0 |
(2. 11) |
|
Из (2.11) видно, что сечение %(х,0) зависит и от размерности сиг
нала N, и от М, составляющих базу В ДКЧС, а также от конкретного ко да частоты {вг}. Для удобства расчета и анализа разрешающей способ ности по задержке и УБЛ в сечении, формулу (2.11) целесообразно раз бить на две, соответствующие окрестности центрального максимума
( x t (*.0), 0 й х < 1) и области пьедестала ( х 2(*.0), \й х< Ы ) функции
неопределенности.
Для области центрального максимума 0 й х <1 к = 0 и ,v = ,v, под ставляя к в (2.11). получаем выражение для % (х,0)
Х х (*.0) = ^ £ e x p { - j • 2п 9 гх М } +
(2.12)
х |
sin Г7t(О —0 .лlA'/vl |
, г |
. |
. п v |
Код частоты {0Г} в ДКЧС представляет собой перестановку из N элементов и, следовательно, содержит все целые числа от 1 до N. Тогда, заменяя Эгна ((0r -1) + 1) в экспоненциальной части первого слагаемого
в (2.12) и используя (2.9), для ;^(х,0) получаем:
*. (*-0) ■ (' |
|
|
+1)>- |
|
|
jc v? sinT^rf#г - в г+1)Мх~\ |
, |
г |
. |
ч |
ч |
Рассмотрим теперь область пьедестала ФН: 1 < х <>N |
Формула для |
||||
Х2(х,0) соответствует (2.11), но для |
1^х< У , &=1,2,...,Л/-1. Наличие в |
выражении (2.11) обратно пропорциональной зависимости * 2(х,0) от
N и присуствие в обоих слагаемых функций sinc(...), аргумент которых зависит от М9указывает на то, что с ростом В значение ФН х 2(*,0)
стремится к нулю, и УБЛ ДКЧС будет также уменьшаться. Оба слагае мых в (2.11) зависят от кода частоты {Вг}, поэтому вид пьедестала (бо ковых лепестков) для каждого конкретного ДКЧС будет различным.
Для удовлетворения требования уменьшения УБЛ ФН можно вос пользоваться ДКЧС, синтезированным на основе матриц Костаса [4]. Среди возможных вариантов частотно-временных кодировок, матрицы Костаса позволяют построить сигнал с ФН «кнопочного» вида и УБЛ, стремящимся к 1IN. Как было показано, частотно-временная матрица сигнала Костаса представляет собой перестановочную матрицу размер ностью N*N. В общем случае существует М вариантов перестановоч ной матрицы, но только те из них, которые обеспечивают не более од ного совпадения элементов при сдвиге матрицы по координатным осям, являются матрицами Костаса. Экспериментальное определение на ЭВМ матриц Костаса при значительном росте N становится трудно выполни мой задачей, поэтому очень важны аналитические формы синтеза мат риц Костаса. Общей формы алгебраической записи всех возможных ва риантов матриц Костаса любой размерности пока не найдено. Но суще
ствуют алгебраические формы для некоторых размерностей А/, в частно сти конструкция Уэлча [6].
Используя полученные выражения для ^ 2(х,0), #(х,0) и # (х ,0 )
(формулы (2.11) при Л=1,2,...,ЛМ, (2.10) и (2.13)), проведены вычисле ния разрешающей способности по задержке (DX), доплеровской частоте (DY) и уровня боковых лепестков в сечениях при х и у, равных нулю (УБЛу и УБЛх соответственно). В табл. 2.2 приведены результаты рас
чета на ЭВМ параметров ФН ДКЧС Костаса, построенных по конструк ции Уэлча для следующих размерностей N: 11 <>=13, g=2), 52 (р=53, g=3), 100 0=101, g=3), 198 0=199, g=3) и 400 (р=401, g=3), где g - при митивный элемент конечного поля GF(p). Параметр М для всех сигна лов принят равным единице.
Таблица 2.2 |
|
|
|
|
|
|
|
|
N |
В |
DX |
DY |
УШЬадБ |
УБЛу, |
Хувл |
YУБЛ |
|
|
|
|
|
М<1 |
\<\x\<N |
дБ |
||
И |
|
0,110520 |
|
|
|
|||
121 |
0,109702 |
-6,87 |
-8,92 |
-6,63 |
0,131501 |
0,130027 |
||
52 |
2704 |
0,023244 |
0,023206 |
-6,70 |
-14,72 |
-6,63 |
0,027552 |
0,027506 |
100 |
10000 |
0,012078 |
0,012068 |
-6,67 |
-15,68 |
-6,63 |
0,014314 |
0,014303 |
198 |
39204 |
0,006098 |
0,006096 |
-6,65 |
-19,07 |
-6,63 |
0,007227 |
0,007224 |
400 |
160000 |
0,003018 |
0,003018 |
-6,64 |
-20,23 |
-6,63 |
0,003576 |
0,003576 |
Параметры Хуш и Уувл в табл. 2.2 соответствуют положению мак симального бокового лепестка в сечениях.
Результаты расчета показывают, что с ростом N совместная разре шающая способность ДКЧС Костаса улучшается (DX и DY уменьшаются), и снижается УБЛх ФН в области пьедестала. Уровень максимального бокового лепестка в сечении при нулевой задержке (УБЛу) постоянен и равен уровню первого бокового лепестка модуля функции sincfaNy) в (2.10), но его положение (Уувл) стремится к нулю при увеличении раз мерности ДКЧС. Максимальный боковой лепесток в сечении ФН ДКЧС Костаса при нулевой доплеровской частоте (УБЛх) сравним со значе нием УБЛу, слабо зависит от У и расположен в области центрального максимума ФН, т. е. при |*|<1. Причем с ростом размерности ДКЧС его положение Хувл становится много меньше единицы.
Было показано [5], что с ростом размерности ДКЧС вклад второго
слагаемого в значение \хх(*,0)| (2.13) при х « 1 уменьшается, поэтому
ФН для этого интервала задержки можно представить следующим при ближенным выражением:
sin (nNMx) |
sin (nNMx) |
,Х « 1 . |
(2.14) |
M H 1-W N sin (яМг) |
nNMx |
||
При М =1 полученное приближенное равенство стремится к точ |
|||
ному, соответствующему |^(х,0)| |
для у - х, |
с ростом |
размерности |
ДКЧС, что подтверждается полученными в табл. 2.2 результатами для разрешающей способности и УБЛ ФН.
На основе полученных в [5] аналитических выражений показано, что для достижения компромисса между обеспечением высокой разре-
тающей способности по задержке и низким УБЛ ФН параметр ДКЧС M-tsfT целесообразно выбирать равным единице. В этом случае база ДКЧС B - N 2, следовательно для повышения помехозащищенности и улучшения основных характеристик ФН необходимо увеличивать раз мерность ДКЧС. Однако это сопряжено с ростом сложности техниче ской реализации устройств формирования и обработки зондирующего сигнала, поэтому представляет интерес использование составных ДКЧС для достижения компромисса между сложностью технической реализа цией и характеристиками РЛС [7].
Вопросы для самоконтроля
1.Что характеризует функция неопределенности сложных дискретнокодированных сигналов?
2.Какими свойствами обладает функция неопределенности?
3.Чем достигается оптимизация формы трехмерного тела неопределен
ности?
4.Отчего зависит разрешающая способность по задержке и частоте?
5.Чем определяется уровень боковых лепестков функции неопределен
ности дискретно-кодированных сигналов?
6. Из каких условий выбираются дискретные значения частоты при ис пользовании дискретно-кодированных по частоте сигналов?
3. СИСТЕМЫ ДИСКРЕТНО-КОДИРОВАННЫХ СИГНАЛОВ
3.1. Функция неопределенности составных дискретно-кодированных сигналов
Одним из способов оптимизации трехмерного тела ФН, позволяющего сформировать ФН с заданными свойствами и в частности более низкий УБЛ, является использование систем ДКС.
В этом случае для получения суммарной функции неопределенно сти (СФН) совокупность ДКС можно рассматривать как составной сиг нал, включающий в себя несколько различных ДКС.
Математическое описание составного дискретно-кодированного сигнала (СДКС), представляющего собой периодическую последова
тельность L различных ДКС одинаковой размерности N, имеет вид: L-1 ЛМ
пРт.п (0 exPfj • 2JT(/o + f m,n)t + j P m,„ l
m -0 n=О |
|
/е[0,(1-1)Г п+ЛТ], |
(3.1) |
1, t e[m7j, +пТ,тТ„ + (и + 1)г];
гДе Pm_n(t) =
[о, /й[отГп+и7’,/я7’ + (и + 1)7’];
/о - несущая частота; Т„ - период следования ДКС в СДКС; Г - длитель ность элементарного импульса с амплитудой ат „ , частотой заполнения
fo+fn,„ и начальной фазой |
. |
|
Комплексная огибающая СДКС (3.1), необходимая для вычисле |
||
ния автокорреляционного интеграла, имеет следующий вид: |
||
I - 1 N- 1 |
|
|
и (0 = |
а”."ехР[j • 2* /а д ('-■тТп-.пТ) + j cpm_„] . |
|
m=0п=О |
|
|
Для оптимизации времени вычисления ФН СДКС представим г в |
||
виде суммы |
дискретных |
и непрерывной составляющих (рис. 3.1): |
т= 1Тп + кТ+8, где 1 = 0, |
; к = 0, 1.... Z>-1; 8 е[0,Т] ; D равно |
целой части отношения Тп ГГ [7].
Для упрощения аналитических выражений и удобства вычисле ний будем считать Тп кратным Т и введем нормированную задержку сигнала х = т/Т и нормированную частоту y - v T [5,7]. Тогда:
x = (lTn+kT + 5)IT = lD + k + x,vj& D = Tn/T , а £ = £ /7 \т .е . хе[0,1].
В соответствии с введенными обозначениями выражение для ФН СДКС, после проведения ряда преобразований, можно представить в следующем, удобном для анализа, виде:
И * .°)Н * 1 (*’°) + Хг (*>°)| ’ |
(3,2) |
||
где |
|
|
|
L-1-/ |
N -\-k |
sin {я{ъ .г-4и*м г+у )(1- xj) |
|
X,{x,0) = Z ]Г U 1- * ) |
X ai+sM ra *.r |
||
x {r!s,r-’7usMr+y)(l - x ) |
|||
5=0 |
r = 0 |
Xexp[j * (t]s,r - |
m+t jk*r - (^ .r + n u s M r ) x + y(2(l + s)D + 2(k + r) + \ + x ) ) + |
||||||
|
|
|
N -2 -к |
|
sin (x{ns,r -4 b sMrA +y)*) |
|
|
J{<Ps,r |
^/+J,jfc+r)] |
|
|
|
|||
X |
al+sM r*\a s,r |
K (77*.r - ^l+s,k+r+\ + -у) ^ |
|
||||
|
|
|
|
|
|
|
|
Xexp[ jж^2tjs r - (rjS'F+ 7jl+sMr+l )x+y(2(/ + s)D + 2(A: + r + l) + x)) + |
|
||||||
+j(%,r-W«.*+r+i)] }•> |
|
|
|
|
|||
, M , v ' l * v D |
|
|
М х( ч ^ ^ - ч ^ +у)*). |
||||
Xexp[ |
jК (2T]s.D-k-u r - { ls j> - k - b r + th+s+\,r) *+ у (2 (/ + Л+1) D + 2r + x)) + |
||||||
./ |
|
|
м |
|
|
sinU {ns,r-m ^M r*\+y)x) |
|
+ J( |
W |
M |
« |
I |
W |
V я (7, г - 7;+а+г+1+у)х X |
|
xexp [jж{lt]s r - (r)s r + T}US'k+r+1) * + |
|
|
|||||
+{T}sj>-ic+r+ih+s*i,r)x+y{2(l+s +\)D+2r+\ +x))+j(<pS'D-k+r ft+JH-l.r)] |
| > |
||||||
/'M -l AM |
|
|
|
|
|
|
|
z = ,2 Z <4 |
|
|
|
|
|
|
|
V 5=0 r=0 |
|
|
|
|
|
|
|
7lsr = fs,fT ~ |
нормированная |
частота |
заполнения элементарного |
им |
пульса с индексами s, г.
Полученное соотношение (3.2) является достаточно общим и по зволяет осуществлять расчет и анализ ФН СДКС с кодированием как одного, так и нескольких параметров сигнала.
3.2. Анализ свойств функции неопределенности составных дискретно-кодированных по частоте сигналов
В качестве составных дискретно-кодированных по частоте сигналов (СДКЧС) рассмотрим последовательности дискретно-кодированных по частоте сигналов (ПДКЧС) и дискретные составные частотные сигналы с частотной манипуляцией (ДСЧЧМ) [3,7].
Последовательность дискретно-кодированных по частоте сигналов представляет собой периодическую последовательность ДКЧС, дискрет частоты f s>rкоторого в каждом элементарном импульсе с индексами s,r равен:
(з.з)
где А/ - шаг сетки частот, который для ПДКЧС можно рассматривать как масштабный коэффициент полосы частот сигнала, а Qs,r - элемент частотно-временной матрицы, определяющий правило кодирования частоты 5-го ДКЧС, входящего в составной сигнал. Тогда rjStr в (3.2) для ПДКЧС имеет вид
ъ .г = А ? = е ЖшГцт = 0Я9Гм , |
(3.4) |
где М = Д/Г - масштабный коэффициент полосы ПДКЧС относительно
длительности элементарного импульса Г [5].
Дискретный составной частотный сигнал с частотной манипуляци ей строится на основе исходного (ДКЧС размерности L) и производяще го (ПДКЧС из L ДКЧС размерности N) сигналов [3]. В соответствии с этим дискрет частоты f5tr в каждом элементарном импульсе с индексами s.r в ДСЧЧМ можно представить в виде f s r = /„ , + f ljpsr, где/„, иf nps,r-
дискреты частот исходного и производящего сигналов соответственно. Дискрет частоты исходного сигнала (/^), представляет собой под
несущую частоту для 5-го ДКЧС производящего сигнала f HS=<5sAfu, где
& - элемент частотно-временной матрицы исходного сигнала, а А/н - частотный интервал между соседними значениями поднесущих ДСЧЧМ Afn = NAf -Мп, где Мп - положительное действительное число, являю
щееся масштабным коэффициентом поднесущих частот ДСЧЧМ. Дискрет частоты производящей ПДКЧС ( f np s>r), в соответствии с
(3.3), равен: / пр1,г = ^ ,гА/ |
|
Таким образом, rjsr в (3.2) для ДСЧЧМ имеет вид |
|
ч,.г= А гТ = ( / „ , + и , . г ) т + 6s r) Д/Г={%SNMn +esr)M |
(3.5) |