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

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

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

Проведем анализ свойств ФН ПДКЧС и ДСЧЧМ в сечениях при нулевой задержке = 0) и нулевом доплеровском сдвиге отраженного сигнала = 0).

Для составных ДКЧС амплитуды элементарных импульсов а5 Гпо­ стоянны и приняты равными единице, следовательно, в формуле (3.2)

Z=l/(LN).

В сечении при нулевой задержке, х = 0 и, следовательно, I, к и х тоже равны нулю. В этом случае, в соответствии с (3.2), выражение для ФН СДКС не зависит от вида т]5Г и, учитывая, что ^ 2(0,у) = 0 , для

ПДКЧС и ДСЧЧМ при х = 0 принимает вид:

^^°’^ =7^ 2 2 ~ ^ ^ exP[j^(25l)+2r+1)] •

(36)

LN 1=0 г=0

ХУ

 

Используя соотношение (2.9), формула (3.6) сводится к виду

Модуль полученного выражения:

 

 

sin (nLDy) sin (nNy)

sin(flrLDy)

sin {nNy)

|* ( М 1 = Lsin(xDy) TuNy

Lsm(7cDy)

(3.7)

xNy

Проведя анализ формулы (3.7), можно показать, что для достижения компромисса между разрешающей способностью по частоте и УБЛ се­ чения |jf(0,^)| целесообразно выбирать значение параметра D=N. В

этом случае сечение ФН при х = 0 для ПДКЧС и ДСЧЧМ соответствует ДКЧС размерности LN, т. е. (0,>>)| = |sinc(?rLNy)| [5,7].

Таким образом, использование составных сигналов на основе ДКЧС размерности N при выборе D - N позволяет в L раз повысить раз­ решающую способность по частоте, так как общая размерность состав­ ного сигнала равна LN, по сравнению с применением отдельного ДКЧС размерности N.

Рассмотрим сечение ФН ПДКЧС и ДСЧЧМ при нулевом допле­ ровском сдвиге частоты = 0). При у = 0 выражение для ФН СДКС (3.2) в случаях ПДКЧС и ДСЧЧМ имеет различный вид, вследствие от­ личия нормированных частот т]5у.

И * . ° ) |= т ^ к (*’° ) + *2 (*>0)1>

<3-8)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L^-\- l

 

A47V 4‘^S^1(;r( ^ 'r

Vl+sMr}^)

*))

*i(* >' E {('-*) Z

- » (, „ - w

„ ) ( , - j )

 

 

5=0

 

 

 

 

 

 

 

 

 

 

X e x p ^ j ^ ( ^ j r

T/l*s,k+r

{jh ,r

+rh * s ,k * r )^ ) + j ( ^ s.r

 

9?/+ i,* + r )] +

+*

У

— M

--------------- ^ - ^ e x p [j^ (2 ^ ir - ( ^ , r +/7/+J>+r+i)xJ +

 

r=0

^ \ls * r

*1l+s,k+r+\ ) *

 

 

 

 

+j {<Ps,r-<Pl+5.Л+Г+1

 

 

 

 

 

 

 

 

 

*.(*«)- Z

* I

-

 

X

--------

)

i<

 

 

 

5=0

[

r = 0

 

^

^

5tD -Jfc-l + r

V /+ 5+ l , r / *

 

 

 

X exp[j Я- (2T]S'D- k - l+r -

{isJJ-k-U -r + 'ft+ J+ l.r) * ) + j ( % .D - * - l+r - ? W l . r ) ] +

/ .

A \

W+X 7 > 1 s i n ( ^ (

? »

j -

h r

~ 7 / + j + i . r ) 0 ~ * ) )

 

 

 

x e x p [ j ^ r ( ^ 0 _t+r -

4 u s+\ j - { n s.D-k+r + '7 / « +i . r ) * )

+ ){<Ps,D-k+r - W + » i .r ) ] } ;

Tjs>rдля ПДКЧС определяется по формуле (3.4), а для ДСЧЧМ - по фор­ муле (3.5).

Для удобства анализа разрешающей способности и УБЛ в сечении |/(х,0)| формулу (3.8) целесообразно разбить на две, соответствующие

окрестности центрального максимума (|^ ц (х,0)|, 0 < х < 1 ) и области

пьедестала ФН.

Для анализа разрешающей способности ПДКЧС и ДСЧЧМ по за­ держке рассмотрим сечение ФН (3.8) при .у = 0 в области центрального

максимума, т. е. при / = к = 0. При D>N х (*,0)

имеет вид:

 

- L-1 ДМ

X ^ 1 ,^ 2 sin( х ( т ^ - 7 ^ ) х )

*.(**°)=T ^Z Z exp[-J2,t^ i]+^ Z Z

 

 

5=0 Г=0

5=0 r=0

 

X e x p [ j x

(r}s r + n*.r+\) * ) + J {<Ps.r-<Ps.r+1 ) ] .

(3.9)

а при D = N\

х

(ЗЛО)

Основной вклад в значение х п(*,0) при х, близких к нулю, в (3.9) и

(ЗЛО) вносят первые слагаемые. Так как эти слагаемые имеют одинако­ вый вид, то разрешающие способности по задержке для СДКЧС с D > N

ис D = N будут примерно равны. Тогда

-L -1 JV-1

ехр[ ■ 2

, * « i

(3.11)

 

5=0 г=0

 

 

Подставляя в (3.11) значение rjsr из выражения (3.4), используя со­ отношение (2.9), и принимая во внимание, что код частоты ПДКЧС {6sr} для каждого 5-го ДКЧС представляет собой перестановку из //элементов и, следовательно, содержит все целые числа от 1 до N, получаем:

(3.12)

Формула (3.12) соответствует выражению, полученному ранее для разрешающей способности по задержке для ДКЧС размерности N. Из этого следует, что разрешающая способность по задержке ПДКЧС не зависит от числа ДКЧС (L) в системе и определяется их размерностью N. При этом, для достижения компромисса между разрешающей спо­ собностью по задержке и УБЛ значение М целесообразно выбирать близким к единице.

Для оценки разрешающей способности по задержке ДСЧЧМ опре­ делим его сечение |^ц(х,0)|. Подставляя rjsr в соответствии с (3.5) в

(3.11), и учитывая, что код частоты исходного ДКЧС {£} представляет собой перестановку из L элементов от 1 до Z,, а код частоты производя-

щей ПДКЧС {0s,r} для каждого s-го ДКЧС представляет собой переста­ новку из N элементов от 1 до N, с учетом соотношения (2.9), имеем

sin (xNMx) sin (KLNMMnx)

Ы * ° ) И 1~И ) N sm{xMx) Zsin (nNMMJc)

sine (xNMx) smo^TcLNMMnx)

(3.13)

sinc(*rM£) sinc(nNMMnx)

Проводя анализ полученного выражения, можно сделать вывод о том, что для достижения компромисса между разрешающей способно­

стью по задержке и УБЛ в сечении ФН |z n(*,0)| для ДСЧЧМ параметры

М и МПцелесообразно выбирать равными единице [8]. В соответствии с (3.13), при Мп- 1 ДСЧЧМ обеспечивает разрешающую способность по за­ держке на уровне ДКЧС размерности LN, т. е. в L раз выше, чем ПДКЧС.

Проанализируем теперь в общем случае ФН в области пьедестала для ПДКЧС и ДСЧЧМ. Так как оба этих сигнала образуются на основе ДКЧС, то для обеспечения низкого УБЛ ФН составных сигналов целе­ сообразно использовать ДКЧС Костаса с различающимися кодами частоты {вг}.

Как уже отмечалось, выбор периода следования ДКЧС в составных сигналах D =N позволяет добиться разрешающей способности по частоте и УБЛ в сечении при нулевой задержке соответствующего ДКЧС размерно­ сти LN. При этом на разрешающую способность по задержке и УБЛ

\хц (*»0)| в окрестности центрального максимума ФН выбор значения D

влияния не оказывает. Однако от выбора значения параметра D и кодовых последовательностей {6s>r} зависит УБЛ ФН в области пьедестала.

При выборе D = N в случае ПДКЧС для обеспечения малого УБЛ ФН необходимо использовать кодовые последовательности {в5>г}, обра­ зующие систему ДКЧС с минимальными значениями взаимокорреляционных функций [3]. Учитывая, что общей аналитической конструкции для матриц Костаса пока не существует, построение подобной системы кодовых последовательностей является сложной задачей, которая в слу­ чае небольшой размерности N может быть решена поиском на ЭВМ.

Использование D 2 2N позволяет разделить по времени частотно­ временные матрицы ДКЧС Костаса, входящие в ПДКЧС, на интервалы, большие или равные длительности ДКЧС (NT). Тогда УБЛ ФН стремит­ ся к уровню УБЛ ФН ДКЧС Костаса размерности N.

Для ДСЧЧМ, построенного с использованием исходного ДКЧС Костаса и производящего ПДКЧС на основе ДКЧС Костаса при исполь­ зовании D = N частотно-временные матрицы ДКЧС, образующие сиг-

нал, разделены по частоте в соответствии с частотно-временной матри­ цей Костаса исходного сигнала, обеспечивая при этом УБЛ ФН не хуже УБЛ ФН ДКЧС Костаса размерности N. Для обеспечения более низкого УБЛ возможно использование D > 2 N , тогда УБЛ ФН ДСЧЧМ стре­ мится к уровню УБЛ ФН ДКЧС Костаса размерности LN.

Следует отметить, что увеличение D>N, кроме ухудшения разре­ шающей способности по частоте и УБЛ сечения |;г(0,у)|, приводит к

росту длительности сигнала, как для ПДКЧС, так и для ДСЧЧМ, а следо­ вательно возрастает минимальная дальность обнаружения импульсной РЛС и становится проблематичным использование данных сигналов в системах селекции движущихся целей (СДЦ) с квазинепрерывным излу­ чением. Тем не менее проблемы, связанные с ростом длительности зон­ дирующего сигнала, могут быть решены на основе компромисса между сложностью технической реализации устройств формирования и обра­ ботки за счет уменьшения длительности элементарного импульса Т.

Таким образом, при построен™ ПДКЧС и ДСЧЧМ с использовани­ ем ДКЧС Костаса параметры М и Мп рекомендуется принимать равными единице, a D - равным N или 2N в зависимости от выбранных кодовых последовательностей Костаса {в5зГ} и требований, предъявляемых к РЛС.

3.3. Системы фазоманипулированных сигналов Система фазоманипулированных сигналов (ФМС) в общем случае пред­

ставляет собой набор из L ФМ-сигналов на основе различных кодовых последовательностей одинаковой размерности N.

Рассмотрим корреляционные свойства кодовых последовательно­ стей ФМС. С этой целью представим совокупность L кодовых последо­ вательностей системы ФМС и соответствующую им последователь­ ность весовых коэффициентов L фильтров сжатия в виде комплексных квадратных кодирующей S и весовой Н матриц размерности L.

Представим матрицы S и Н по строкам и столбцам соответственно в виде

s = № 0№ (у).. А ( у)Г; н = [ Ш Ш - Ш Ъ

где sjj) - кодовая последовательность /-го ФМ-сигнала; hJJ) - весовые коэффициенты соответствующего фильтра сжатия, символ «т» обозна­ чает операцию транспонирования матрицы, i= l...L - число кодовых по­ следовательностей системы ФМС.

Представим процесс сжатия /-го сигнала из системы ФМС в виде произведения весового вектора-столбца на кодирующую вектор-строку

R,-= *,(/)$,СО

(3.14)

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

Аналогично (3.14) представим процесс обработки системы ФМС с учетом когерентного суммирования L сжатых сигналов

R , = 1=1

(315)

Согласно теореме о блочном разбиении матриц, выражение (3.15)

представляет собой матричное произведение R

весовой матрицы на ко­

дирующую матрицу

 

R = H S .

(3.16)

Для формирования квазиидеальной формы суммарной функции неопределенности (СФН) матрица R должна быть диагональной. Тогда сумма элементов каждой боковой диагонали, параллельной главной, тождественно равна нулю, а значит равен нулю каждый боковой лепе­ сток СФН.

Рассмотрим условия, которым должны удовлетворять матрицы S и

Ндля того, чтобы матрица R была диагональной.

Вобщем случае кодирующая матрица S должна быть унитарной, тогда весовая матрица Н является эрмитово-сопряженной кодирующей

матрицей H = S*T Данный случай соответствует ФМС, состоящей из набора L многофазных сигналов. В случае бинарно-фазовых сигналов кодирующая матрица S должна быть ортогональной, тогда весовая мат­

рица Н является транспонированной кодирующей матрицей H = ST Практический интерес представляют кодирующие матрицы, содержа­ щие только элементы +1 и -1, т.е. матрицы Адамара. В принципе могут быть использованы и другие кодирующие матрицы, например матрицы функций Хаара, что соответствует системе амплитудно-фазоманипу- лированных сигналов. Но применение матриц Адамара обеспечивает максимальное использование энергетического потенциала РЛС, поэто­ му использование матриц Адамара предпочтительнее.

Обработка системы ФМ-сигналов заключается в сжатии каждого из L сигналов в соответствующем фильтре и последующем когерентном суммировании всех L сжатых сигналов.

Вслучае бинарно-фазовых сигналов Moryij использоваться систе­ мы кодовых последовательностей Баркера.

Вчастности, суммарная автокорреляционная функция (АКФ) двух различных кодов Баркера-ипределяется выражением

*(»») = *,(»»)+ *,(»»),

(3.17)

где %1(т)= £ а„ап_т,

%г(т)

Ь

л . .

-

соответственно АКФ

п=т+1

 

п=т+1

 

 

кодов Баркера {ап} и {Ьп} с числом элементов N\ и N2, имеющих разно­

полярные боковые лепестки.

 

 

'

 

С целью получения максимальной области компенсации боковых

лепестков необходимо

обеспечить

выполнение

условия |N 1- N 2| = 2.

Данное требование удовлетворяется для N-разрядных кодов Баркера с числом элементов N]=5 и N2=3, N\=l и N2=5, Ni=13 и N2=11. Модуль суммарной АКФ в этом случае описывается выражением

'NX+N2 при т - 0;

|*(»)| = 0

при w = ±l,±2,...,±(N1-2);

(3.18)

1

при т = ±(NX- 1).

 

Как показали расчеты, при использовании пары кодов Баркера с числом элементов М=13 и N2= ll УБЛ суммарной АКФ составляет -27,6 дБ, а с учетом весовой обработки - -55,2 дБ [8].

3.4. Формирование функции неопределенности дискретнокодированных сигналов с заданными свойствами

Перед современной радиолокацией стоит важная задача обнаружения малозаметных объектов. Трудность ее решения заключается в том, что подобные объекты являются слабоотражающими и поэтому требуется обнаруживать слабые сигналы при наличии шумов и помех.

В этих условиях одним из путей повышения надежности обнару­ жения является использование широкополосных зондирующих сигна­ лов. При этом функция неопределенности таких сигналов должна удов­ летворять определенным требованиям и в частности обеспечения УБЛ.

Однако оптимизировать ФН по УБЛ только за счет выбора пара­ метров сложного зондирующего сигнала, учитывая особенности меха­ низма формирования тела неопределенности ДКС в полной мере не удается. Дальнейшее совершенствование ФН по снижению УБЛ может осуществляться на основе весовой обработки [1].

Методы весовой обработки с целью уменьшения УБЛ могут быть основаны на амплитудном взвешивании огибающей передаваемого сиг­ нала или частотного отклика согласованного фильтра приемника [1].

Рассмотрим возможности метода взвешивания по амплитуде пар­ циальных радиоимпульсов ДКЧС для снижения УБЛ.

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

37

Т, осуществлять взвешивание дискретных компонент по определенному закону. Воспользуемся вторым способом, наиболее приемлемым для практического применения. В соответствии с формулой (2.4) свойства ФН взвешенного по амплитуде ДКЧС определяются двумя массивами: частотно-кодовым Nfn и весовым ап.

Как известно, получение выходного отклика фильтра сжатия сложного сигнала с низким УБЛ должно сопровождаться коррекцией спектра сигнала [1]. Для модулированных по частоте сигналов коррекция спектра может осуществляться методом весовой обработки во временной области, поскольку временной и частотный параметры в данном случае взаимосвязаны. Следовательно, первоначально выбирается частотная весовая функция, распространяющаяся на область, занятую сеткой кодовых частот, и определяющая частичное снижение уровня тех компонент, которые расположены на краях указанной области. Формирование весового вектора [ап] можно упростить и сделать универсальным, если весовую функцию Щх) определить на нормированном интервале х = [-1, 1]. Алгоритм получения [ап] будет состоять в разбиении отрезка х = [-1, 1] на (N- 1) элементарных интервалов величиной Ах = xi+l- x { = 2/(N -1), значе­

ния весовой функции на краях которых будут представлять собой веса

дискретных компонент ДКЧС: ап = РГ(хп), где хп =-1+///пДх; Nfn =Nfn -1.

В результате анализа и расчетов была выбрана весовая функция общего вида

^ (* ) = /> + (l-/> )c o s " (j* j,

(3.19)

г д е = 0,..., 1 - уровень пьедестала.

Таким образом, механизм формирования весового массива заключа­ ется в том, что значения частотной весовой функции определяются в соот­ ветствии с частотно-временным кодом. На рис. 3.2 представлены значения весов для ДКЧС с размерностью кода N = 10 и степени косинуса п =1.

Проведена проверка влия­ |* м | ния весовых функций с раз­ личивши значениями степени п на свойства ФН ДКЧС.

На рис. 3.3 представлено сечение |^(г,0)| тела ФН по

оси задержки (F= 0), рассчитан­

 

 

 

ное при помощи ЭВМ с исполь­

 

 

 

зованием соотношения (2.4).

01 ог

аз аз о*

0 7 о * 0 9 1 rjj

При этом были взяты зна­

 

 

 

чения параметра р = 0,3; 0.5; 1;

 

Рис. 3.3

 

2; 4 весовой функции (3.19) и матрица Костаса размерностью Af= 10, кото­ рая имеет последовательность чисел Nfn =2,4,8,5,10,9,7,3,6,1. Тогда

rffn =N/n -1=1,3,7,4.9,8,6,2,5,0, Ах=х,+1- JC,. = - ^ = | ; хп= - ] Щ ь .

В результате расчетов было показано, что наиболее предпочти­

тельной является весовая функция со степенью косинуса п = 1:

 

fV(x) = р + (1 - /?)cos

(3.20)

Использование в (3.19) более высоких степеней косинуса является нежелательным, поскольку нулевые значения производной таких функций в точках х = ±1 способствует медленному изменению амплитуд дискрет­ ных компонент на краях частотного интервала, занятого сеткой частот, и существенному искажению ФН, а когда п <1, то УБЛ ФН увеличивается.

Рассмотрим сечение |^(г,0)| с использованием весовой функции

(3.20), матрицы размерности N= 10 и различных уровней пьедесталар. На рис. 3.4,а представлено сечение |^(г}0)| с уровнями пьедестала: р = \\р =0.6;

р = 0,3 (значение р = 1 соответствует отсутствию весовой обработки).

Как видно из рис. 3.4данный метод подавления боковых лепест­ ков ФН является весьма эффективным в пределах одного дискрета Т по задержке. Снижение УБЛ происходит со снижением уровня пьедестала р весовой функции. Однако при этом снижение УБЛ сопровождается расширением главного лепестка, т.е. ухудшением разрешающей спо­ собности сигнала, что является характерным в процессе осуществления весовой обработки сигнала [1]. На рис. 3.4,6 показаны удаленные боко­ вые лепестки, соответствующие задержке Т< тйЪТ и значениям пара­ метров уровня пьедестала весовой функции: р = 0,8; р = 0,3.

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

Рассмотрим сечение тела неопределенности по оси доплеровских частот (т = 0), построенное с использованием соотношения (2.4). Учи­ тывая, что к = 0 и при г = 0, выражение, определяющее ФН примет вид:

\X{0,F)\ = Z

sin(TiP) N-1

 

 

тф

Х а»ехр ОяР(2и + 1)} •

 

 

п=0

 

 

 

 

 

 

На рис. 3.5 представлено сечение |^(0,F )|

со значениями /7 = 1 ;

 

 

р = 0,6; р = 0,3 весовой функ­

 

 

ции (3.20) и матрицей N = 10.

 

 

Отметим, что с понижением

 

 

уровня

пьедестала р

повыша­

 

 

ется УБЛ в пределах одного

 

 

дискрета Af по частоте.

 

 

На рис. 3.6л и б показаны

 

 

результаты расчета значения ко­

 

 

эффициента расширения глав­

 

 

ного лепестка Кр(т) и уровня

 

 

первого

бокового

лепестка

 

 

Упбл(т) по задержке,

коэффи­

 

 

циента

расширения

главного

лепестка KP(F) и уровня первого бокового лепестка по частоте УПбл(^) при значениях р = 0,1... 1.

Показано, что уровень первого бокового лепестка по задержке УпвяСО снижается с -14 дБ до -32 дБ с понижением р - с 1 до 0,3. Когда р < 0,3, то Упбл(т) повышается. Так, значение уровня пьедестала р = 0,3

40