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

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

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

Формулы (3.3.70), (3.3.73), (3.3.75), (3.3.77), (3.3.79) и (3.3.80), описывающие элементы дискриминационной характери­ стики в различных зонах, позволяют, в конечном итоге, воспроиз­ вести полную ДХ программируемого временного дискриминатора рассматриваемого типа. При выводе перечисленных формул удоб­ но было относительную ДХ нормировать к длительности строба тст. Применительно к рассматриваемому варианту построения ВД эле­ ментами относительной ДХ будут безразмерные множители в пра­ вых частях упомянутых формул. Однако, во всём предыдущем ма­ териале по временным дискриминаторам и последующем описании программируемых частотных дискриминаторов фигурируют отно­ сительные дискриминационные характеристики, нормированные к половине длительности строба (временного или частотного). Для сохранения единства стиля изложения в приводимой ниже табл. 3.3.4 даны элементы относительной ДХ, нормированной к полови­ не длительности строба, т.е. к 0,бтст. Для их получения во всех перечисленных выше формулах безразмерные множители увели­ чиваются вдвое.

1-я зона

2-я зона

 

3-я зона

Таблица 3.3.4.

 

4-я зона

А,_ _

Ai, Aj?

 

й

й

- 6А, +7А2

 

Ат, А , As

Ат, А , Ая, Ai

и,Д1=-3

 

—4,6А^ +1,6А2 —Ag

Ai +Ag - А2 - Aj

 

Um2“ 2(AJ- А2)

и,д3_

2А, +А2 4

*** Aj +A2 +A3 + A4

 

 

 

Таблица 3.3.4, продолжение.

 

5-я зона

 

6-я зона

7-я зона

 

Ай, Ая, Ai

 

Ая, Ai

A.

 

A2 - 1,6AS + 4,5A4

5A„-7A3

UM7=3

 

A2 +Ag +2A4

j

вдв 2(A4 —Ag)

 

При переходе к размерной дискриминационной характеристи­ ке содержащиеся в табл. 3.3.4 элементы относительной ДХ умно­ жаются на 0,51^,. В результате получаем те же соотношения, ко­ торые содержатся в перечисленных выше формулах. На рис. 3.3.19 изображена относительная дискриминационная характери­ стика ВД при четырёх стробах, перекрывающихся на половину своей длины, составленная из элементов, алгоритмы которых за­ писаны в табл. 3.3.4. Для перехода к абсолютным отсчётам шкала оси ординат умножается на 0,51^. Строго линейный участок ДХ составляет ±1,25^ (протяжённость - 2,5тст). Плоские части ДХ возникают там, где имеется только один отсчёт в отклике. Здесь

г"

,4 I----------------------

 

----------- ----------------------

---------------------

1— 1

-4

-3

- 2 - 1 0

1

2 т а = ^ 4

Рис. 3.3.19.

образуется методическая погрешность в измерениях. Как следует из рис. 3.3.19, максимальное значение такой погрешности состав­ ляет 0,25тст (на краях плоских частей ДХ).

Линеаризация ДХ в рассматриваемом дискриминаторе дости­ гается усложнением программы его работы. Она должна содер­ жать анализатор ситуаций, который выявляет набор отсчётов на каждом такте формирования ДХ. Состав таких наборов показан во второй строке табл. 3.3.4. По результатам анализа логическая схема выбирает и запускает один из алгоритмов, записанных в третьей строке табл. 3.3.4. Естественно, что сбои в работе анали­ затора приведут к появлению аномальных ошибок в измерениях временного рассогласования. Это та цена, которую приходится платить за высококачественную работу такого ВД в штатном ре­ жиме.

3.4.ПРОГРАММИРУЕМЫЕ ЧАСТОТНЫЕ ДИСКРИМИНАТОРЫ

3.4.1.П ри н ц и п ы п о с т р о е н и я и а л го р и т м ы ф у н к ц и о н и р о в а н и я

ПРОГРАММИРУЕМЫХ ЧАСТОТНЫХ ДИСКРИМИНАТОРОВ

Общая схема предварительной обработки радиолокационных сигналов (рис. 3.3.1) является достаточно универсальной для со­ временных БРЛС с цифровой обработкой. Её общая структура со­ храняется для многих режимов от обнаружения цели, до её авто­ матического сопровождения с оцениванием фазовых координат. После каждого зондирования, отражённый от цели сигнал, пройдя

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

Применительно к системе автоматического сопровождения це­ ли по скорости процедура М-точечного БПФ трактуется как алго­ ритм многофильтрового спектроанализатора, содержащего М+1 частотных каналов, которые для краткости назовем фильтрами БПФ.

Исходными данными для алгоритмов частотных дискримина­ нт

торов (ЧД) будут комплексные амплитуды UFP, сформированные

в результате выполнения операции БПФ. Здесь U - амплитуда

*

комплексного цифрового сигнала, a Fp - комплексная амплитуд­ но-частотная характеристика р-го фильтра БПФ, входящего в со­

став частотного дискриминатора. Величина р принимает дискрет-

______

*

ные значения р=1,М + 1. Далее вычисляются модули

UFP этих

(

* V

амплитуд, либо их квадраты I UFP I Систему обработки первого

вида будем называть линейным

детектором, а второго - квадра­

тичным.

Фильтры БПФ располагаются в диапазоне частот от 0 до Fn, где Fn= l/T n - частота повторения зондирующих импульсов. Часто переходную частоту f0, т.е. частоту входного сигнала, при которой дискриминационная характеристика ЧД обращается в нуль, рас­ полагают вблизи Fn/2. Это достигается соответствующим выбором частоты опорного сигнала, который формируется цифровым гете­ родином, входящим в ППС. Точное значение fo задается соотно­ шением

fo =

M -e v

(3.4.1)

2M Fn,

где переменная ev=l при чётном числе фильтров, образующих ЧД, и ev= 0 - при нечётном.

Номер фильтра p(i), в зависимости от их числа г в алгоритме

ЧД и количества М точек БПФ, определяется как

 

М- r(ev - 1)

+ i .

(3.4.2)

2

 

 

Здесь i = 1, г - текущий номер фильтра в их гребенке, входящий в состав ЧД.

При традиционном выполнении ЧД на расстроенных конту­ рах, число фильтров г=2. Тогда, для М=32 и линейном детекторе алгоритм ЧД принимает вид

F17-F ie

(3.4.3)

ичд=^У-

+ F17

Fie

 

Введение нормировки за счет суммарного сигнала устраняет влия­ ние амплитуды цифрового сигнала U на форму ДХ.

Расширение апертуры ДХ достигается увеличением числа фильтров в составе ЧД. Продуктивным оказался метод весового суммирования, когда весовые коэффициенты при увеличива­ ются с ростом номера фильтра p(i).

Детально процедура весового суммирования при получении алгоритма временного дискриминатора рассмотрена в п. 3.3.2. По­ этому приведем конечную формулу алгоритма (3.3.15), адаптиро­ ванную к обозначениям, связанным с частотным дискриминатором

- ------ ;-----------------

(3.4.4)

Сохраняется и общее правило формирования весовых коэффи­ циентов: при четном числе г фильтров, входящих в состав ЧД, ко­ эффициенты составляют последовательность положительных и от­ рицательных нечетных чисел, а при г нечетном - последователь­ ность четных (включая нуль). Соответственно, в ЧД с квадратич­ ным детектором в формуле (3.4.4) F ^ заменяется на

При плавном изменении частоты f входного сигнала перемен­ ные Fp^)(f) и F ^ (f) отображают соответственно модули амлитуд-

но-частотных характеристик фильтров БПФ и квадраты модулей, а зависимость u ^ f) в (3.4.4) является аналитическим выражением дискриминационной характеристики ЧД.

Для реально используемых величин М и г справедлива сле­ дующая форма записи набора АЧХ фильтров БПФ

з т | ( 2я£Тп - ^ ( р { 1 ) - 1 ) ]

0,5^fT n-^ (p (i)-l)]

где f и p(i) изменяются соответственно в пределах f = 0,Fn и

p(i) = 1»N+1 .

Дискриминационную характеристику ЧД принято изображать

в расстройках Af относительно переходной частоты f0, т.е.

 

Af = f-fo = f - ^ p F n.

(3.4.6)

Введем относительную расстройку fa, разделив величину Af на по­ ловину полосы пропускания фильтра БПФ AF/2 = Fn / 2М, т.е.

A f

(3.4.7)

L = - = - = — M -M + ev.

a AF/2 Fn

 

Используя (3.4.7), выразим текущую частоту f через относитель­ ную расстройку fa

4 + M - e v

 

(3.4.8)

"*

 

Проведя замену в (3.4.5), получим

 

Msin£(fa +r + l-2i)

(3.4.9)

______ а________________

| (4 + r + l-2i)

АЧХ достигают максимальных значений, равных М, при рас­ стройках fao* удовлетворяющих уравнению fa0+ r + l - 2 i = 0 , т.е.

Сопоставление формул (3.4.4) и (3.4.10) показывает, что вели­ чины fao и весовые коэффициенты ДХ равны при одинаковых i, поэтому при такой нормировке расстроек линейный участок дис­ криминационной характеристики имеет единичную крутизну.

Относительные (безразмерные) ДХ формируются, если Fp^fa) подставить в формулу (3.4.4). Для перехода к абсолютным дис­ криминационным характеристикам вместо fa вводится абсолютная расстройка Af (3.4.6), а ординаты относительной ДХ умножаются на AF/2 = Fn /2M .

В качестве примера на рис. 3.4.1 приведены несколько вари­ антов относительных ДХ, рассчитанных по приведенным выше формулам.

 

 

л

 

 

 

r v

 

 

 

Г=7

 

-10 -6

0

5

10

а)

б)

 

 

Сделаем несколько предварительных замечаний. Следует иметь ввиду, что для снижения уровня боковых лепестков в АЧХ фильтров БПФ в тракт обработки сигналов перед операцией БПФ вводятся весовые окна (Хемминга, Кайзера и т.д.). Второе замеча­ ние касается условности термина «линейный участок» ДХ при

многофильтровых ЧД. В дальнейшем линейным участком ДХ бу­ дем называть ту ее часть, которая располагается в диапазоне час­ тот, занимаемом главными лепестками АЧХ фильтров БПФ. По­ яснения к рисункам: на рис. 3.4.1,а отображена ДХ для линейно­ го детектора без весового окна при чётном (г^б) числе фильтров в составе ЧД; на рис. 3.4.1,б - ДХ для линейного детектора без ве­ сового окна при нечётном (г=7) числе фильтров; на рис. 3.4.1,в - ДХ для квадратичного детектора без весового окна при нечётном (г®7) числе фильтров; на рис. 3.4.1,г - ДХ линейного детектора с весовым окном Хемминга при чётном (г=6) числе фильтров.

Для рис. 3.4.1,а характерна абсолютная линейность ДХ на участке fa±l. Здесь выполняется равенство u^=fa. Эта линейность нарушается при нечетном числе фильтров. Особенно заметны ее проявления для ЧД с квадратичным детектором (рис. 3.4.1,в). Не­ линейность ДХ в области, близкой к fa= 0 приводит к заметному затягиванию времени переходных процессов в следящей системы с таким ЧД. При больших рассогласованиях на рис. 3.4.1,а, б и в видны аномалии, выражающиеся в отклонениях от линейной за­ висимости 1Цд от fa. Равенство u^^dfa (где d=const) выполнятся лишь для расстроек fa0, удовлетворяющих условию (3.4.10). Эти аномалии вызваны боковыми лепестками в АЧХ фильтров. Они не оказывают существенного влияния на качество работы следящей системы, но вызывают большие погрешности в измерениях часто­ ты прямоотсчётным методом. Аномалии эффективно устраняются при включении в программу ЧД какого-либо весового окна, на­ пример, окна Хемминга (рис. 3.4.1,г), но при этом неизбежно снижение отношения сигнал/шум приблизительно на 1,37 дБ.

3.4.2. Ст а т и с т и ч е с к и е х а р а к т е р и с т и к и и о ц е н к а то ч н о с ти ч а с т о т н ы х д и с к р и м и н а т о р о в

Сопоставляя принципы построения программируемых вре­ менных и частотных дискриминаторов, можно отметить, что наи­ большие различия содержатся лишь в формировании исходных данных. Так, при построении БД это образование временной свертки в пределах строба, а для функционирования ЧД формиру­ ется АЧХ фильтров БПФ. Если эти элементы получены, то ос­ тальная технология образования алгоритмов функционирования дискриминаторов содержит по существу однотипные процедуры.

Такая особенность учитывается и при получении статистиче­ ских характеристик дискриминаторов и оценки их погрешностей.

Поэтому в материалах этого пункта будут довольно частые ссылки на п.3.3.3, чтобы избежать повторов одинаковых, по существу, формул.

Отношение сигнал/шум после когерентного накопления в фильтрах БПФ определяется соотношением q=A/ay, где A=UM, U - уровень сигнала на выходе предварительного сумматора (см. рис. 3.3.1), М - число когерентных накоплений (число точек БПФ), ау - СКО шума на выходе системы обработки.

Для произвольной расстройки по частоте fa введем относи­ тельный сигнал

(3.4.11)

В дальнейшем будем полагать, что СКО шума на выходах различ­ ных фильтров БПФ равны, т.е. tfyp(i)=ayp(2)=...=ayp(r)=Gy.

Для ЧД с линейным детектором математические ожидания mp<i)(fa) и среднеквадратические отклонения Op(i)(fa) процессов на выходе системы обработки (после вычисления модуля) описывают­ ся формулами (3.3.23), если в них аргумент та заменить на fa, а переменная Up(i)(fa) берется из (3.4.11). В варианте квадратичного детектора те же характеристики описываются формулами (3.3.24), если в них также произвести замену та на fa и использовать отно­ шение (3.4.11).

Статистическими характеристиками для частотного дискри­ минатора также будут: среднее значение сигнала на выходе дис­ криминатора как функция рассогласования fa (шумовая дискри­ минационная характеристика ШДХ), и флуктуационная характе­ ристика, определяемая как зависимость среднеквадратического отклонения выходного отсчета дискриминатора от относительной расстройки fa. Первая определяет систематическую погрешность, а вторая - случайную.

Расчет статистических характеристик ЧД выполняется по формулам (3.3.29)-(3.3.31), если в них всюду аргумент та заменить на fa. В качестве примера на рис. 3.4.2 приведены дискриминаци­ онная и флуктуационная характеристики для системы с линейным детектором и числом фильтров БПФ г=6 . Графики построены в относительных величинах по осям абсцисс и ординат. При перехо­ де к абсолютным значениям отсчеты по указанным осям умножа­ ются на половину полосы пропускания фильтра БПФ AF/2 = Fn / 2М.

Как следует из рис. 3.4.2,а, действие шума приводит к дефор­ мациям ДХ, что следует из сравнения графиков с исходной бесшумовой кривой (1 ). Деформации приводят к уменьшению кру­ тизны ДХ на ее линейном участке и сокращению размера этого участка (эффект сминания). Флуктуационная характеристика име­ ет приблизительно постоянную величину в области линейности ДХ. Вне этой области СКО резко возрастает. Можно показать, что при любой конфигурации ЧД с уменьшением отношения сигнал/шум до нуля Стфх стремится к некоторой постоянной величине софтах* Величина сТфвддо определяется лишь числом фильтров г. Так, для варианта ЧД с линейным детектором эта зависимость имеет вид

_

IC4 - ^

2- 1)

^ф х max у

З л г

Для г= 6 имеем афх max

®.73, что согласуется с рис. 3.4.2,б. Объ­

яснение такого свойства ЧД состоит в том, что с ростом шума уве­ личивается СКО, но одновременно падает величина коэффициента передачи.

В п. 3.3.3 было показано, что достаточно содержательной ха­ рактеристикой, определяющей погрешность дискриминатора, яв­ ляется среднеквадратическое отклонение которое находится как значение флуктуационной характеристики при нулевой рас­ стройке, когда аЧд=афх(0). При линейном детекторе в системе об­ работки для вычисления о ^ используются формулы (3.3.23), (3.3.30), (3.3.31) в которых аргумент та= 0 (fa=0). При этом соот­ ношение (3.3.23) заменяется двумя другими, пригодными для двух предельных значений q - большом и малом. Для частотного дискриминатора эти соотношения получаются из (3.3.32) и (3.3.33) при замене в них аргумента та на fa, а относительная ам­ плитуда Up(i)(fa) определяется формулой (3.4.11).

Аналогичная методика применяется при вычислении для системы обработки с квадратичным детектором. Однако здесь применяются модифицированные для ЧД формулы (3.3.24), кото­ рые справедливы для любого значения q.

Окончательные формулы расчета СКО для варианта линейного детектора представлены в табл. 3.4.1.

Здесь использованы следующие обозначения: г - число фильт­ ров БПФ, образующих ЧД; q - отношение сигнал/шум на выходе ЧД; bi - параметр, учитывающий долю, вносимую в а^д боковыми лепестками АЧХ фильтров; Кпо=0,46 - относительный уровень пересечения АЧХ двух фильтров, смежных с центральным, fa=0; КП1=0,83 - относительный уровень пресечения АЧХ двух цен­ тральных фильтров при fa=0 .

При наличии в системе обработки квадратичного детектора формулы для расчёта относительных СКО ЧД приведены в табл.

3.4.2. К ранее указанным обозначениям здесь введено: 1^=0,153

-относительный уровень пересечения АЧХ двух фильтров, смеж­ ных с центральным, при fa=0 . Если в формулах табл. 3.4.2 поло­ жить q=0, то получим предельную величину СКО

<*чд тах= Vo*2 —1) / Sr .

(3.4.13)