Интерферометрический интерфейс системы определения относительных координат радиоизлучающих объектов Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

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

радиоизлучающих объектов

1 12 1

В.Г. Манжула , С.Г.Крутчинский , А.В.Савенко , В.В.Воронин

!ФГБОУ ВПО «ЮРГУЭС», г. Шахты Ростовской обл.

2 ФГБОУ ВПО «ПВГУС», г. Тольятти Самарской обл.

ВВЕДЕНИЕ

Использование активных методов вычисления относительных координат цели в силу критичной загрузки радиодиапазона, очевидно, не позволит развивать эти методы для целей безопасности движения [1]. На эту проблему в первую очередь обратили внимание отечественные специалисты из ВВИА им Н.Е. Жуковского [2]. Этот вопрос детально обсуждался за «круглым столом» на 5 международном конгрессе «Авиационные технологии 21 века» (г. Жуковский, 1999 г.). Детальное исследование этого вопроса выполнено в РАН в рамках федеральной программы под руководством академика А.А. Красовского [3-5]. В настоящей работе показана техническая осуществимость предложенного А.А. Красовским принципа пассивной локации и вычисления относительных координат цели в рамках компьютерных интегрированных мобильных систем.

Содержательная сторона технического предложения связана с развитием принципа специализированного интерферометрического интерфейса [6], основными нестандартными узлами которого являются прецизионные широкополосные умножители аналоговых сигналов и специальные ограничители спектра, обеспечивающие выделение и необходимое прецизионное преобразование в реальном масштабе времени информации об относительных координатах цели.

1. Постановка задачи

Интерферометрический метод [1] мало чувствителен к передаточным коэффициентам приемных каналов. Поэтому можно предполагать, что калибровка, выполняемая после монтажа комплекса достаточна для его нормального функционирования.

Для интерферометрического определения разностей дальностей [2] в любой плоскости сигналы каналов, соответствующие приемным точкам, попарно переключаются при наличии сдвига (задержки) по времени одного из умноженных сигналов. При отсутствии помех, пренебрегая зависимостью амплитуд от дальности и указанных выше условиях, определяется взаимная корреляционная функция первой пары сигналов

Я12 (¿, Т12 ) = М[в! ^>2 (г - Т12 )] =

где М [*] - математическое ожидание,

^(0, £2(^-т12) - входные, принимаемые системой напряжения (напряжения в приемных антеннах),

т12 - временной сдвиг (время задержки),

иш - шумовое напряжение в приемных антеннах,

с - скорость света,

А Г ~ Г2

Ат12 - Т!2 ,

С

и - частота излучения,

г - относительная дальность,

т12 - задержка одного из сигналов во времени,

Лш (Ат12 ) - корреляционная функция огибающей, полагаемой стационарным случайным процессом.

По условию [1] время корреляции огибающей тш много больше периода несущей, так что тши >> 1. Для истинной корреляционной функции (1) необходимо использовать специальный фильтр, выделяющий квазипостоянную составляющую

Лг - 0 5 Лш (Ат12 )со8(иАТ12 ). (2)

Поиск главного экстремума корреляционной функции составляет суть интерферометрического метода. Практически всегда вместо истинной корреляционной

функции Л12 (¿1зт12) приходится пользоваться ее оценкой Я\2 (¿1зт12), вычисленной на некотором скользящем интервале времени Ти, существенно превышающем тш [3]. Вычисление может производиться на базе различных аппаратных средств. Однако необходимо отметить, что для чисто цифрового варианта частота считывания,

преобразования и фиксации сигналов должна определяться не несущей частотой и, а т Ш1. Именно это обуславливает требования к цифровым процессорам системы. Так, при считывании 10 значений входящих сигналов за время тш — 10 мкс и вычисления Лш(Ат12) на интервале Тш —100 тш —10 мс частота считывания будет составлять 100 кГц, а необходимое быстродействие, даже без использования специального ЗУ, будет соответствовать 2 х105 оп/с.

2 Структурная схема интерфейса системы

Приведенная оценка показывает, что для реализации в реальном масштабе времени микрокомпьютерной системы пассивной локации необходим специализированный интерфейс, обеспечивающий выделение квазипостоянной составляющей (2) с дополнительными функциями:

- предварительная частотная селекция целей;

- вычисление корреляционной функции и выделение огибающей входного сигнала;

- точная автоматическая настройка смежных приёмных каналов.

Базовая структура такого интерфейса приведена на рис. 1.

УВЧ1, УВЧ2 -усилители высокой частоты; Г - управляемый генератор; А1, А2 - антенны; АП1 -аналоговый перемножитель; ФНЧ - фильтр низких частот; Б1, Б2 - ключи

Рис. 1 Структура интерферометрического интерфейса

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

После того, как будет определён диапазон частот, в котором излучение цели максимально, система должна переключиться в режим определения параметров движения цели.

Рассмотрим это на примере двух приемных точек. От двух приёмных антенн А1 и А2 через ключи Б1 и Б2, работающие в противофазе, сигнал поступает на умножитель АП1. При этом из-за того, что антенны разнесены в пространстве, перемножаемые сигналы сдвинуты по фазе. В результате на выходе перемножителя АП1 появляется сигнал, содержащий низкочастотную функцию от угла сдвига фаз и высокочастотные гармоники.

Пусть на входы умножителя поступают следующие сигналы:

и - А с°з(ш^ + ф1), ^

и - А СОБ(ш^ + Ф2).

Сам перемножитель реализует следующую функцию:

ип - (и • ^ )х(и • ^) + и • к + иг • к2. (4)

где к1з к2 - коэффициенты, определяющие неидеальность АП;

_Р1з ¥2 - коэффициенты преобразования каналов АП.

Тогда сигнал на его выходе будет иметь следующий вид:

ип -1 А1А2^1со8(ф1 -ф2) +1А1А2Е2соб^ю? + ф1 +ф2)+к^соБ^Ш + ф1)+к2А2оо$(2Ш + ф2). (5)

Из выражения (5) видно, что в результате перемножения двух сдвинутых по фазе гармонических сигналов в спектре произведения присутствует постоянная составляющая, являющаяся функцией разности фаз, гармоники с частотой исходного сигнала и с удвоенной частотой. Так как информацию о параметрах движения объекта несёт разность фаз, то для её выделения необходимо использовать фильтр нижних частот (ФНЧ). Кроме ограничения спектра такой фильтр должен характеризоваться относительно низкой неравномерностью его амплитудно-частотной характеристики и предельно низким дрейфом нуля. Действительно, сравнением соотношений (2) и (5) несложно установить, что дрейф нуля будет непосредственно определять точность вычисления разности фаз и, следовательно, относительной дальности.

В качестве аппроксимирующей функции ФНЧ будем рассматривать полиномы Баттерворта, как имеющие минимальную неравномерность в полосе пропускания, и полиномы Чебышева, как имеющие наибольшую крутизну характеристики в полосе задерживания. Оценим порядок фильтра, необходимый для обеспечения следующих его характеристик:

- неравномерность АЧХ в полосе пропускания - <0,1 дБ;

- подавление сигнала на частоте 1 МГц - >40 дБ;

- частота среза - 500 кГц.

Для аппроксимации по Баттерворту порядок фильтра связан с требуемыми характеристиками следующим соотношением:

n =

1об

f 1QQ1Amin 1qO. 1Amax

-1 - 1

Y

J

2 log

v®, J

(6)

где ^min, Атах - минимальное и максимальное затухания сигнала;

шс, - граничные частоты полосы пропускания и полосы заграждения фильтра.

Для аппроксимации по Чебышеву порядок фильтра связан с требуемыми характеристиками следующим соотношением:

n =

Arch

Arch

v®c J

(7)

Подставляя численные значения, получим, что минимальный порядок для Баттервортовой аппроксимации равен 5, а для Чебышевской - 4. Реализовать устойчивую систему с низким дрейфом нуля лестничного типа четвёртого порядка на D-элементах [ 10] невозможно. Поэтому будем рассматривать оба типа фильтра 5-го порядка. При этом нормированные полиномы для знаменателей передаточной функции будут иметь вид:

- для Баттервортовой аппроксимации

Р5Б =(1 + р)( + 0.618 р + р2 )( +1.618 р + р2 ); (8)

- для Чебышевской аппроксимации

Р5Ч =(0.539 + р)(1.195 + 0.333р + р2 )0.636 + 0.872р + р2 ). (9)

Настройку каналов интерферометра (рис. 1) предлагается производить следующим образом: одно плечо является неизменным и образцовым, а второе может изменять ёмкость колебательного контура. На оба плеча подаётся гармонический сигнал от образцового генератора. Если сдвиг фазы точно равен 90°, то в результате перемножения должно получиться нулевое напряжение. Появление нуля на выходе умножителя при изменении параметров колебательного контура и будет говорить о том, что оба плеча интерферометра стали идентичными. При этом вовсе не обязательно, чтобы частоты настройки точно совпадали с несущей сигнала.

Сама процедура подстройки может осуществляться как с использованием алгоритма двоичного приближения [2], так и линейным изменением варьируемого заведомо смещённого параметра. Второй алгоритм более предпочтителен, так как делает более плавным переходной процесс.

3 Динамические характеристики и параметры интерфейса

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

Для повышения вероятности точного определения параметров движения объекта предлагается проводить подстройку интерферометра в каждом цикле измерения.

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

Поэтому найдём вид и время переходного процесса в каждом линейном звене, а затем на основании этих данных оценим быстродействие всей системы в целом.

Пусть на вход системы поступает гармонический сигнал

^ ) = зт(о/), (10)

где а - обобщенная круговая частота.

Найдём сигнал, появляющийся на выходе полосового фильтра (резонансного контура), воспользовавшись преобразованием Лапласа:

U1 (р ) =

Ар + В1

■ +

Qp + D1 2 2 1 2 2 * р + pdю р + ш р р + а

Здесь параметры определяются следующим образом:

' ,);

(її)

Iа -юр;

A1 = М 1 (а В1 = -M1;

С1 = М1 (-а2 + юр — = М1,

);

(12)

где параметр М1 = ■

а^р ю р

а2шР +(а2 “ЮР )2

Переходя обратно во временную область, получим, что сигнал состоит из двух составляющих: затухающей и незатухающей.

Затухающая часть имеет вид

dрюр

U13 (t) = A1e 2 cos

f dP dn ю ^

ю р^І 1 - • t + Bx - A1 • -р-^ 1 1 2

4

V V

1

dp ш р

ю,

1

sin

1 - -í

4

f -р 1

ю р^І 1 —р • t

4

V

(13)

Незатухающая часть имеет вид

U1h (t) = C1 cos(at) + — sin(at)

а

(14)

Здесь параметры А1, В1, С1, П1 описаны выше.

Эти выражения справедливы для случая, когда центральная частота полосового фильтра шр не совпадает с частотой а принимаемого сигнала. Более частный случай, когда резонансный контур точно настроен на принимаемый сигнал (ш ^ а), дает

следующие выражения для затухающей и незатухающей составляющих:

U13 (t)= -а "

1

dpшp

ю,

1

sin

f dp 1

ю р]і 1 —р • t

4

V

(15)

1 - £

4

U,„ (t) = sin(at). (16)

Для второго (не перестраиваемого) плеча интерферометра аналогичным образом находим составляющие выходного сигнала

I------ \

d„ ю„ 1

d„ ш„

U2з (t) = A2 • e 2

cos

f 2 1 f

ю р^ V 1 —р • t +

4 У V

В - A- р р

x

t

2

X

e

2

e

X

1

-п ю„

ш,

2

81П

1 - -р 4

ю

В,

у

и2н (t) = С2 cos(at) н—2 s1п(at).

а

(17)

(18)

Здесь параметры формулы (17) определяются следующим образом:

А2 = -М 2а2 -р ю р;

В2 = -М2 (шр -а2 }

■).

^2 = М а dp ш р;

Вт = -М2ар2 (юр - а2

(19)

где параметр М2 =

dР ш р

а2^юр +(а2 -юр )2

При точной настройке колебательного контура (ш ^ а) переходной процесс будет иметь следующий вид:

и2, ^ )=—-1

dn

• е

dpа

+ ■

1

cos

—„ ю„

Г dp 1

ш р^1 1 р • t

4

V У

+

2а2

ш,

V

■ • е

2

sm

Г dp 1

ш р]1 1 - -р • t

4

V У

(20)

и2н (t) = cos(at). (21)

Как видно из приведённых выражений, в обоих плечах интерферометра

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

примерно

^ер =, (22)

dp ш р

где п - порядок аппроксимирующей передаточной функции ФНЧ.

Тогда при п = 3 получаем точность 3 %, а при п = 5 точность установления сигнала составляет 1 %.

При прохождении сигнала через ФНЧ, включенный на выходе аналогового перемножителя (рис. 1), естественным образом возникает некоторый переходной процесс, характер которого обусловлен только видом корней характеристического полинома, а длительность - их численными значениями.

Как видно из выражений (8) и (9), для обоих типов фильтров переходный процесс будет колебательным, так как характеристические полиномы содержат комплексно-сопряжённые пары корней, а его длительность можно оценить по формуле (22). Для типовых случаев систем охранной сигнализации допустимо использовать граничную частоту ФНЧ в 300 кГц, поэтому в приведенных выше соотношениях частота полюса будет равна шр = 2л • 500000. Отсюда можно найти время установления сигнала при 1 % точности выходного сигнала

г

2

х

е

V

2

1

2

2п

2 • 5

перФНЧБ

d р т1пШ р

0.618 • 3141500 Аналогично, для Чебышевского фильтра

^ерФНЧЧФ = 2.7 х 10 с.

Таким образом, существует свобода выбора пропускания при построении ограничителя спектра

= 5 х 10-6с .

(23)

(24)

граничной частоты полосы по второму варианту, что существенно снижает требования к активным элементам схемы. Так, в случае, если граничная частота фильтра равна 300 кГц, длительность переходных процессов будет определяться значением (23).

Таким образом, зная задержки в каждом линейном звене этой нелинейной системы, можно предположить, что минимальное время для определения параметров движения цели будет не меньше суммы времени максимальной перестройки колебательного контура и времени установление сигнала на выходе ФНЧ. Учитывая, что частота настройки входных контуров значительно больше граничной частоты ФНЧ, время настройки интерфейса будет определяться только полученными соотношениями.

4 Алгоритм функционирования интерфейса

Рассмотренная структурная схема интерферометрического интерфейса является цифроуправляемой периодически перестраиваемой системой.

Ниже приведён алгоритм функционирования, включающий в себя функции первичного обнаружения несущей и вычисления разности фаз принимаемых сигналов (рис. 2).

Рис. 2 Алгоритм функционирования интерфейса

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

Основные выводы

Полученные результаты показывают, что для успешного решения общей задачи в рамках предложенной структуры интерферометрического интерфейса оказывается достаточным реализация прецизионного аналогового перемножителя (АП) и прецизионного бездрейфового фильтра нижних частот (ФНЧ). Эти базовые узлы интерфейса являются составной частью многочисленных устройств РЭА, и их схемотехника в настоящее время достаточно хорошо отработана. В классе аналоговых перемножителей достаточно отметить работы [8,9], где предложены прецизионные структуры этих устройств и рассмотрено их применение в сложно-функциональных блоках. Что касается измерительных фильтров нижних частот, то их реализация предпочтительна на базе D-элементов с расширенным диапазоном рабочих частот [10].

Литература:

1.Красовский А.А. Определение относительных координат радиоизлучающих объектов в пространстве интерферационным методом // Изв. РАН. Теория и системы управления. 1997, № 3.

2.Красовский А.А., Сизых В.Н. Алгоритмическое обеспечение пассивных интерферометрических бортовых комплексов радиолокации // Теория и системы управления. 2000, №2.

3.Красовский А.А., Белоглазов И.Н., Чигин Г.П. Теория корреляционно -экстремальных навигационных систем. М.: Наука. 1979.

4.Красовский А.А., Наумов А.И. Интерферометрическая радиопеленгация как пассивная радиолокация // Изв. РАН. Теория и системы управления. 1998, № 1.

5.Красовский А.А. Развитие теории дальнего пасивного мониторинга, навигации и резервного управления движением // Изв. РАН. Теория и системы управления. 1999, № 2.

С. 77 - 83.

6.Красовский А.А. Пассивная макроволновая радиолокация, мониторинг, навигация и резервное управление движением // Изв. РАН. Теория и системы управления. 1998, № 3. С. 156 - 163.

7.Крутчинский С.Г. Нестандартные аппаратные средства адаптивных регуляторов и систем пассивной локации // Изв. РАН. АиТ, 2001, № 7. С. 138 - 147.

8.Starchnko E. I., Gavlicky A. I. Krutchinsky S.G. Analogous voltage multiplier based on bipolar transistors and MOSFET. Proceeding ICCSC'08, Bucharest, Romania, 2008 p.p. 155160.

9.Старченко Е.И., Гавлицкий А.И., Крутчинский С.Г., Малышев И.В. Квадратурные модуляторы для технологического процесса SGB25VD. Опыт практической разработки // Проблемы разработки перспективных микро - и наноэлектронных систем : сб. науч. трудов всерос. науч.- практ. конф. М.: ИППМ РАН, 2008, С. 313-317.

10. Крутчинский С.Г., Щекин ДА. Структурный синтез D-элементов с расширенным частотным диапазоном. // Международный научно-технический журнал “Электроника и связь”, №8. Т.2, 2001 С.83-87.

Статья подготовлена при выполнения НИР по теме «Разработка и исследование аналоговой электронной компонентной базы нового поколения для систем связи, радиоэлектроники и технической кибернетики» в рамках федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 годы»