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

книги / Широкополосные дискретно-кодированные сигналы в радиотехнике и радиолокации

..pdf
Скачиваний:
12
Добавлен:
20.11.2023
Размер:
3.29 Mб
Скачать

Таким образом, сравнительный анализ ФН ДКС и ее сечений дает возможность отдать предпочтение наиболее перспективным сложным сигналам, позволяющим обеспечить высокие энергетические возможно­ сти РЛС и требуемое совместное разрешение целей по дальности и ско­ рости при малом УБЛ.

При выборе зондирующего сигнала необходимо учитывать еще одно важное свойство, связанное со скрытностью работы РЛС. В этом отношении ДКЧС Костаса дают возможность обеспечить высокую скрытность работы РЛС, так как общее число перестановочных матриц для выбранного массива 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)