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

книги / Радиолокационные измерители дальности и скорости. [Т.1]

.pdf
Скачиваний:
13
Добавлен:
12.11.2023
Размер:
21.62 Mб
Скачать

Два замечания по пользованию этой и последующими табли­ цами.

1 . Если отношение сигнал/шум ЦдБ задано в децибелах, то при

его подстановке в формулы таблицы, следует перейти к q = 10 20

2 . Для получения абсолютных значений погрешностей необхо­ димо умножить величину овд на половину длительности строба Тстр/2 или половину дискрета дальности ЛД/2 .

Поскольку формулы для вычисления СКО процессов на выхо­ де ВД при квадратичном детектировании справедливы для любого значения q, то, положив q=0 , получим выражение для вычисле­ ния максимального значения СКО

(3.3.34)

Наличие такой предельной погрешности наблюдается и на флуктуационных характеристиках (рис. 3.3.7, 3.3.8) при больших рассогласованиях |тв|>2. Объяснение такого свойства характери­ стик ВД состоит в том, что с ростом шума увеличивается СКО, но одновременно падает величина коэффициента передачи ВД.

Аналогичная особенность проявляется и для ВД с линейным детектором. К сожалению, выражение, аналогичное (3.3.34), нель­ зя получить непосредственно из формул табл. 3.3.1 для линейного детектора, т.к. они справедливы лишь для дл^ 1 б дБ. Однако ис­ следования показали, что для системы обработки с линейным де­ тектором предельная погрешность при для= 0 описывается соотно­ шением

(3.3.36)

Сами по себе выражения (3.3.34) и (3.3.36) не имеют практи­ ческой ценности, но устанавливают факт, что вне зависимости от числа п предельная погрешность в системе с квадратичным детек­ тором в два раза превышает аналогичную величину во ВД с ли­ нейным детектором.

Втабл. 3.3.2 представлены формулы относительных СКО вре­ менных дискриминаторов с минимальным числом стробов (п= 2 и п=3).

Вкачестве иллюстрации применения формул табл. 3.3.1 и 3.3.2, на рис. 3.3.9 и 3.3.10 представлены графики зависимости

Тип детектора

 

п= 2

 

п=3

Линейный

 

4

_

4 (

4 - х )

 

<Гвд= —

вд_

q

г

 

 

 

 

 

 

 

-7

= + Vw

 

 

_Vo,5q2+2

 

л/2

 

 

 

 

2 л/2

Квадратичный

вд

0,25q2 + 2

° вд

0 6q2 + 3

 

 

 

 

 

CKO авд, как функции от числа п стробов в составе ВД, при отно­ шении сигнал/шум 0дв=20 дБ (рис. 3.3.9), и зависимости ствд от от­ ношения сигнал/шум при минимальном числе стробов п=2 и п=3 (рис. 3.3.10).

Для удобства пользования графиками, в табл. 3.3.3 даны но­ мера вариантов построения ВД.

Число стробов

Тип детектора

Таблица 3.3,3.

Номер варианта

Четное

Линейный

1

Квадратичный

2

 

Нечетное

Линейный

3

Квадратичный

4

 

Анализ приведенных зависимостей показывает, что мини­ мальные флуктуационные погрешности даёт вариант 4 (нечетное число стробов при квадратичном детектировании). Причина состо­ ит в том, что для этого варианта в области та«0 , для которой вы­ числялась СКО, дискриминационная характеристика имеет пло­ скую часть (рис. 3.3.4,б), т.е. практически нулевой коэффициент передачи, что приводит к подавлению шума. Отрицательным мо­ ментом такой формы ДХ будет затягивание переходных процессов при малых рассогласованиях.

При квадратичном детектировании и четном числе стробов (вариант 2 ) погрешности существенно возрастают, ввиду высокой крутизны дискриминационной характеристики в области та*0 (рис. 3.3.6). При линейном детекторе (варианты 1 и 3) разброс между погрешностями для четного и нечетного значений чисел стробов в ВД незначителен. Для малых величин п флюктуационные погрешности этих вариантов занимают промежуточные значе­ ния между погрешностями вариантов 2 и 4, а рост авд с увеличе­ нием числа п происходит практически по линейному закону.

Помимо рассмотренных погрешностей, шум вызывает дефор­ мацию дискриминационной характеристики, которая сводится к снижению крутизны ДХ и уменьшению ее линейного участка. Да­ лее, для краткости, эти проявления будем называть сминанием ДХ шумом. Эффект сминания по-разному проявляется при линей­ ном и квадратичном детектировании, что видно из рис. 3.3.5 и 3.3.6.

Получение количественных закономерностей, определяющих проявление эффекта сминания ДХ, начнем с простейшего примера для п=3. На рис. 3.3.11 изображены ШДХ при п=3 для четырёх значений отношения сигнал/шум: длк=15 дБ, 20 дБ, 25 дБ и для (]дБ=ао. Графики, отображающие ШДХ для линейного и квадра­ тичного детекторов представлены соответственно на рис. 3.3.11,а, и 3.3.11,6. На рисунках введены следующие обозначения кривых: 1-<1дБ=15дБ, 2 -< ^ = 2 0 дБ, 3 - <3двв 25 дБ, 4 - <Здв=оо.

Проявления эффекта сминания будем оценивать количествен­ но по степени уменьшения максимума кривой ШДХ с увеличени­ ем шума. В приведенном примере такой максимум соответствует относительной расстройке та=2. Это значение та приходится на на­ чало плоского участка ДХ, когда при формировании дискримина­ ционной характеристики, остается всего один отсчет, который совпадает с максимумом временной свертки, равным UN. Здесь, как и ранее, U - амплитуда принимаемого импульса на входе сис-

1

темы обработки, а N - максимальное число накопленных отсчётов в предварительном сумматоре.

На основании формулы (3.3.7) при Рд=10 и соотношения (3.3.29) запишем выражение для ШДХ при п=3 и та=2

(п\ _

2 т 6(2 ) - 2 ш 4(2 )

(3.3.36)

ш4(2) + т б(2) + т в(2)'

В этой формуле при сдвиге сигнала на та=2 только отсчет Шв(2) содержит сигнал в остальных отсчетах присутствуют только шу­ мы. Используя (3.3.32) и (3.3.33), найдем

me(2) * UN,

т 4(2) = т б(2) = ауД 5^ .

(3.3.37)

Тогда из (3.3.36) получим

 

2 ( Ш

- стуД б я )

(3.3.38)

т шх(2) =

 

UN + 2ауЛ/0,5я

 

Соотношение (3.3.38) раскрывает механизм сминания дискри­ минационной характеристики шумом. Оно происходит за счет то­ го, что математическое ожидание релеевского шума стуЛ/0,5я

уменьшает числитель и одновременно увеличивает знаменатель в уравнении ДХ. Если шум отсутствует сгу=0, то шшх(2)=2, что и следует из рис. 3.3.11,а, т.е. сминания ДХ не происходит.

Последовательно увеличивая число п в (3.3.29), нетрудно обобщить соотношение (3.3.37) на п стробов

(3.3.39)

С учётом (3.3.21), выражение (3.3.39) преобразуется к виду

(3.3.40)

В силу условия (3.3.32) формула (3.3.40) справедлива лишь для q>15 дБ.

Для системы обработки сигналов с квадратичным детектором формулы, аналогичные (3.3.32) и (3.3.33), получаем непосредст­ венно из соотношений (3.3.24), которые пригодны для любых зна­ чений щ(та). Поэтому для (3.3.37), учитывая (3.3.24), запишем

me(2) = 2ст* + (UN)2, т 4(2) = т в(2) = 2ст*.

(3.3.41)

Следовательно,

(3.3.42)

Сопоставление (3.3.38) с (3.3.42) позволяет на элементарном уровне установить причину лучшего сохранения формы ШДХ в системе обработки с квадратичным детектором. В частности, здесь отсутствует математическое ожидание шума в числителе формулы для ШДХ, что обусловлено различием статистических характери­ стик сигнала и шума на выходах линейного и квадратичного де­ текторов. Обобщение формулы на п стробов дает

(3.3.43)

Количественную меру проявления эффекта сминаемости ДХ при сопоставлении двух рассмотренных вариантов построения ВД назовём коэффициентом несминаемости, который определим как

(3.3.44)

где тал " рассогласование на пределе линейного участка дискрими­ национной характеристики. Развернутое выражение для Кнс полу­ чим из (3.3.40) и (3.3.43)

q 2(n - l) + ^jOjbnq

^2n + q2)(Vo^5rcq - l)

Это выражение, как и (3.3.40), справедливо лишь для q£l6 дБ. Удобно также ввести более надежный показатель устойчиво­

сти к сминаемости, определив его как

Дне = (Кнс-1)100%.

(3.3.46)

Величина Днс, выраженная в процентах, характеризует выигрыш, который даст применение в системе обработки квадратичного де­ тектора вместо линейного, в соответствии с принятым критерием.

На рис 3.3.12 приведены зависимости Днс от отношения сигнал/шум для ВД с двумя п=2, тремя п=3 и четырьмя п=4 строба­ ми.

Рис. 3.3.12.

Выигрыш возрастает при уменьшении сигнал/шум и увеличении числа стробов, используемых при программировании временного дискриминатора.

3,3.4. П р о гр а м м и ру е м ы й в рем ен н о й д и с к р и м и н а т о р п р и

ис п о л ьзо в а н и и СЛОЖНЫХ РАДИОЛОКАЦИОННЫХ с и гн а л о в

ВБРЛС с низкими и средними частотами повторения импуль­ сов повышение энергии сигнала, а, следовательно, и дальности действия, достигается применением сложных сигналов значитель­ ной протяженности. Характерными примерами таких сигналов будут фазоманипулированные (ФМ) последовательности, построен­ ные по коду Баркера (КБ) [17]. Описанные в п. 3.3.2 общие прин­ ципы построения программируемых временных дискриминаторов сохраняются и для сложных сигналов. Необходимо только дать ряд дополнительных пояснений к схеме рис. 3.3.1. Аппаратная часть схемы до ППС остается без изменений. Однако содержание выполняемых здесь операций несколько меняется. Если при про­ стом сигнале результатом суммирования в предварительном сум­ маторе будет сам накопленный сигнал, который передается в ППС, то при использовании сигнала с КБ в ППС передаются на­ копленные отсчеты элементов кода Баркера. В программируемом процессоре сигналов эти данные аппаратно или программно обра­ батываются в фильтре, согласованным с применяемым кодом. Па­ ра согласованных фильтров обеспечивают обработку вещественной

имнимой частей комплексного сигнала. Модуль свернутого сигна­ ла КБ содержит один основной лепесток, амплитуда которого в NB раз (Na - число элементов КБ) больше исходного сигнала, и (N0-l) боковых лепестков, относительная амплитуда которых составляет 1/Na, т.е. убывает с ростом числа элементов в коде. При таком уровне боковых лепестков возникают значительные искажения в дискриминационной характеристике БД, если число стробов, на которых строится БД, больше двух.

Проблема снижения уровня боковых лепестков (БЛ) путем их подавления с помощью специальных фильтров является актуаль­ ной. В работе [36] обсуждается эта проблема и предлагается не­ сколько вариантов таких фильтров.

Известно [60], что нормированный спектр свернутого сигнала кода Баркера записывается как

(3.3.47)

где S(co) - спектр центрального лепестка свертки КБ, a F(o) - спектр так называемой кодовой последовательности, порождаю­ щей боковые лепестки.

В [60], где приводится вывод развернутого выражения для F(©), показано, что

F(<D) = N3 -1 + sin(aN3TK6)^

(3.3.48)

sin(e>TK6)

 

где Ткб - протяженность элемента КБ.

Из (3.3.47) следует, что, если создать импульсный (цифровой) фильтр с амплитудно-частотной характеристикой (АЧХ)

(3.3.49)

N 1 |sin(g>N3TK6) *

э яп(юТкб)

то при прохождении сигнала свертки КБ через такой фильтр боко­ вые лепестки будут устранены полностью. Задача состоит в том, чтобы реализовать соотношение (3.3.49) в виде программы цифро­ вого или нерекурсивного (трансверсального) фильтра.

Поскольку функция Н(со) периодична, первый шаг состоит в разложении ее в ряд Фурье

H (f) = А 0 +

оо

Oirkf

(3.3.50)

Z A k - c o s = ^ ,

 

k=i

£г

 

где fT= ------ -

тактовая

частота,

с которой следуют отсчеты ре-

2 Ткб

 

 

 

шётчатой функции свертки. Для коэффициентов ряда Фурье име­ ем

А0 =

; Ak = Y H (f)co s^ ^ d f

(3.3.51)

1т о

*т 0

 

_ f

В целях сокращения записи введем относительную частоту f =—.

ч

Представим косинус в (3.3.50) как

cos(27ikf) = -^cos(2Tckf) + jsin(2Tckf) + соз(2тгк?) - jsin(2xkf)j=

2

= -(z+k+z'k),

(3.3.52)

2 '

'

 

где z - аргумент функции, подвергшейся Z-преобразованию. Тогда (3.3.50) можно представить в виде соотношения

H(z) = А0 + l4 *-(zk + z_k),

(3.3.53)

k=i 2 '

7

 

в котором H(z) - дискретная передаточная функция некоторого трансверсального фильтра, устраняющего боковые лепестки в свертке КБ.

Реализация (3.3.53) требует значительного объема вычисле­ ний, поэтому выполнение такого фильтра в практических разра­ ботках затруднительно. Следует найти способы упростить это вы­ ражение даже за счет некоторых потерь в качестве фильтрации. Необходимо отметить, что предлагаемые здесь методы построения фильтров подавления боковых лепестков пригодны лишь для ко­ дов Баркера с положительными по знаку боковыми лепестками. Это значит, что число Na составляет 2, 5, 13 или составные коды из этих элементов.

Анализ функции H(f) для N8=13, представленной на рис. 3.3.13,

Рис. 3.3.13.

показывает, что простейшим аппроксимирующим выражением для H(f) будет

1

______ 1______ sin(27tN3fTK6)

(3.3.54)

Ho(f) = N .- 1

(N8 - 1)(2N3 -1) sin(2*fTK6)

Оно найдено полуэмпирическим путем на основе приближенного равенства 1/(1+х)«1-х, при х « 1 . Выражение (3.3.54) сводится к виду:

1

1

%1 eJ2*(N. - 1- 2^

(3.3.55)

НоМ = N , - 1 (N3 - 1)(2N3 -1)

kt 0

 

Обратное преобразование Фурье дает импульсную характеристику фильтра

hJt)

= — 6(t)---------------------

Nf

6ft- (N8 -1 - 2kVy. (3.3.56)

W

N . - 1 ^

(Na-1X2N3- I ) k^

^

Окончательное выражение для весовой характеристики фильтра подавления боковых лепестков имеет вид

(3.3.57)

Смысл выполняемых операций состоит в следующем. Сумма невзвешенных сигналов на входах сдвигового регистра умножает­ ся на коэффициент 1/(2N8-1)=0,04. Поскольку при k=(Na-l)/2 ар­ гументом второй дельта-функции становится t, вычитание данных осуществляется с выхода центрального отвода сдвигового регистра. Применение рассмотренного фильтра снижает боковые лепестки 13-и элементного кода Баркера с -22,5 дБ до -34 дБ, т.е. на 12 дБ. Более эффективное подавление боковых лепестков обеспечивают приведенные в [35] усложненные варианты фильтров подавления.

З.З.б. Особенности построения программируемого временного

ДИСКРИМИНАТОРА ПРИ НИЗКОЙ СКВАЖНОСТИ СИГНАЛОВ

Низкая скважность сигналов характерна для режима ВЧПИ. Здесь величина скважности достигает Q=5 и менее. В этой ситуа­ ции проблема расширения апертуры дискриминационной характе­ ристики БД простым увеличением числа стробов в составе ВД не имеет удовлетворительного решения. Так, при Q=5 в зоне приема нельзя расположить впритык более 3-х стробов. В силу того, что половину зоны приема займут плоские части ДХ, то на линейный участок ДХ останется 50% зоны приема.

Положение можно улучшить, если располагать стробы с пере­ крытием. Такой вариант при аналоговом выполнении ВД пред­ ставлен в [33]. Его можно реализовать и при цифровой обработке сигналов, если воспользоваться описанным в [38] методом частич­ ного суммирования. Суть метода состоит в том, что виртуальный строб, образуемый на выходе аналого-цифрового преобразователя гребенкой импульсов начала суммирования, дробится на более мелкие элементы, внутри которых образуются частичные суммы. Затем из этих частичных сумм формируются отсчеты самих вир-