Методика расчёта антенных структур многобазовых фазовых пеленгаторов

УДК 621.317.08

Г.Г. Порубов, В.П. Денисов

Методика расчёта антенных структур многобазовых фазовых пеленгаторов

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

Проектирование многобазовых фазовых пеленгаторов (МФП) начинается с выбора структуры антенной решётки (АР), т.е. с выбора числа фазометрических баз, их относительных размеров и положения в АР, обеспечивающих выполнение заданных на пеленгатор технических требований.

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

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

Расчёт начинается с выбора наибольшей базы пеленгатора, которая полностью определяет точность пеленгования при построении пеленгатора по «методу уточнений» и в значительной степени при использовании метода максимального правдоподобия [2]. Обозначим наибольшую базу пеленгатора символом /i, пусть /1» X , где X - длина волны.

Полная разность фаз сигналов на базе /1

Ф1 = 2п—sinа = ф1 + 2пк ,

X

где ф1 - измеренная разность фаз в радианах п <ф1 <п; к - полное число периодов разности фаз на базе /1, утраченное при измерениях; а - угол прихода плоской волны, отсчитанный от нормали к антенной системе.

Устранение неоднозначности измерений на базе /1 заключается в отыскании к по результатам измерений на меньших базах.

Обозначим — = %1, sina = v, а разность фаз будем измерять в рад/2п, так что 0,5 <ф1 <0,5 .

X

Формулу (1) перепишем в виде

Ф1 = nx1v - к .

При изменении v = sin а в пределах от -1 до 1 полная разность фаз на базе /1 ф1 + к изменяется

/1 /1 D « X в пределах от — до —. В указанных пределах ось v можно разделить на интервалы длиной —, в

XX /1

пределах которых разность фаз изменяется на 2п радиан. Каждому из этих интервалов, называемых интервалами однозначного измерения на базе /1, соответствует определенное значение к. На

26 ЭЛЕКТРОНИКА, ИЗМЕРИТЕЛЬНАЯ ТЕХНИКА, РАДИОТЕХНИКА И СВЯЗЬ

рис. 1, а показана ось V , разбитая на интервалы однозначности в пределах V <1. Принято, для примера, — = 5 . Вертикальными стрелками в каждом интервале однозначности показано рассчитанное

X

по формуле (1) значение V , соответствующее некоторой разности фаз ф1 и определённому к.

-1ж -5я -Зя -ж

-1

-К

Ч-

0

я 3 ж 5ж 1ж 9ж <р+2пк -*Ч-*Ч-*Ч-*Ч-!-

а

ж

Зж %+2пк б

Н-*-!-

Истинный пеленг 1

Рис. 1. Геометрическая интерпретация неоднозначных фазовых измерений

Допустим, мы хотим найти истинное значение V , используя измерение на другой базе /2, также неоднозначной. На рис. 1, б показаны измеренные пеленги на оси V. Для построения принято

— = 2, считается, что погрешности фазовых измерений отсутствуют. Тогда в точке истинного пе-

X

ленга отсчёты по двум шкалам совпадают, а в других интервалах однозначности они расходятся. Совпадение отсчётов является критерием истинности пеленга.

При наличии фазовых погрешностей результаты измерений пеленга на различных базах не совпадают ни в одной точке V . За результат измерений принимается тот отсчёт по базе /1, где отсчёты по базам /1 и /2 наиболее близки.

/1 /2

При построении рис. 1 предполагалось, что относительные базы —, — - взаимно простые це-

X X

лые числа. Из рис. 1 видно, что совокупность отсчётов разностей фаз по базам обладает периодичностью: если при некотором V результаты измерений равны ф1 и ф2 , то они точно повторяются в точках V ± 1.

Данный результат легко проверить аналитически. Представим совокупность относительных баз вектор-столбцом пх :

( ñх1) \ñx2/'

1 ф :

(Ф).

\Ф2/

а совокупность результатов измерении - вектором ф :

ф =

Тогда вектор полных разностей фаз

Ф = ф + к = ñxv.

Если v изменяется на ±1, то вектор) полных разностей фаз

Ф = Пх (v ±1) = ñxV ± Пх .

На каждой базе полная разность фаз изменилась на целое число периодов, поэтому измеренная разность фаз не изменилась ни на одной из баз.

Таким образом, если относительные базы - взаимно простые целые числа, интервал однозначного измерения величины v = sin а равен 1, что соответствует угловому сектору ±30° относительно нормали к антенной системе. Фазовые соотношения на базах достаточно рассмотреть на указанном интервале.

При работе пеленгатора в диапазоне частот условие целочисленности относительных баз — не

X

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

ñx =

Пх А^оди = 4 ,

где ех - вектор взаимно простых целых чисел, а величина Ауоди имеет смысл интервала однозначного пеленгования системой баз пх . Действительно, подставляя в

пх =

имеем

Х АУодн

Ф = еху—1—

А^дн

Пусть V изменится на Ауодн . Тогда Ф = ех ( ±Ауодн)—1— = еху—1— ± ех - на каждой базе

АЦздн Аvодн

полная разность фаз изменилась на целое число периодов, так что измеряемые разности фаз ф, остались прежними.

Таким образом, рассмотрение фазовых соотношений на интервале V < 0,5 при целочисленных базах равносильно рассмотрению этих соотношений на интервале А^дн , если относительные базы пх, не целочисленные, но относятся между собой как взаимно простые целые числа. Будем считать, что это условие выполняется, и рассматривать фазовые соотношения на интервале А^дн .

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

Пусть истинный пеленг задается некоторой величиной V . В этой точке на базе ¡2 измеренная разность фаз

Ф2 = Пх2'^—У,

где к = 0, ±1, ±2,..., ех2.

Если изменить пеленг на у интервалов однозначности на большой базе ¡1 А^дн =-, то из-

пх1

меряемая разность фаз на базе ¡2 будет

1

Ф2у = Пх2\ V+у— | — (У + /) = Пх2V+у—-(к + /). I Пх1) Пх1

Разность ф2 (5) и ф2 у (6) характеризует величину отклонения измеренной разности фаз на базе ¡2 с направления V+уА^д^ относительно измеренной с направления V

А2 у =Ф2 у —Ф2 = —/.

Пх1

Величина А2у определяет вероятность исключения из возможных результатов измерения отсчёта пеленга, отстоящего от истинного на у интервалов однозначного измерения по большой базе ¡1. Задаваясь величиной А2у, найдём из (7) отношение баз

¡2 А2 , + /

¡2 =_2^ , (8) ¡1 у

где / = 0,1,2,... - целые числа; у - количество интервалов однозначного измерения по большой базе относительно истинного пеленга.

¡2

С другой стороны, при выбранном отношении баз — из (8) можно найти А2 у для любого ин-

¡1

тервала однозначного измерения по базе ¡1

Д 2 у =-360° = 180°

Д 2 у =-360° = 120°.

Д2у = А - Г ,

11

где у - номер интервала однозначности на базе Д; величина Д2 j выражена в рад/2п; / = 0,1,2,... -целые числа.

Величина Д2у не может быть больше 180°: на базе /2 соседние отсчёты отстоят друг от друга на 360° (см. рис. 1), поэтому если расстояние от разности фаз ф2у до разности фаз ф2, соответствующей истинному пеленгу, больше 180°, оно меньше 180° в соседнем интервале однозначности. В формуле (9) это соответствует изменению целого числа / на единицу.

В приводимых ниже примерах величина Д2 у принималась равной следующим значениям:

1

2

Д 2 у = у360° = 154,29°

Д 2 у = -2-360° = 144°,

1

3

Минимальная величина Д2у определяется погрешностями фазовых измерений. Если, к примеру, фазовые погрешности распределены нормально, имеют нулевые средние значения и дисперсии стф , можно принять Д2у = (4 6)стф .

Порядок действий при выборе баз следующий.

По базе /1 , определённой исходя из требуемой точности пеленгования, и принятой величине Д2у , для первой зоны однозначности измерений на базе /1 (у = 1) по формуле (8) находится база /2 . Для других интервалов однозначного измерения на базе /1 (т.е. для других у) находятся фактические значения Д2у, которые определяются выбранной базой /2. Находится минимальный номер интервала однозначного измерения у на базе /1 , для которого выбранная база /2 не обеспечивает желаемого значения Д2 у. Далее находится база /3 с помощью формулы, аналогичной (8). Запишем её в виде

,

/1 ]

где у - номер интервала однозначного измерения на базе /1 , в котором база /2 не обеспечивает желаемой величины Д2 у ; Д3 у - отклонение измеренной разности фаз на базе /3 с направления V+у'Дуодн1 относительно измеренной с направления V . Интервалы Д2у, Д3у не обязательно должны быть одинаковыми.

Далее по найденной базе /3 находятся величины Д3 у в других секторах однозначности базы /1 . Находится минимальный номер у сектора однозначности, в котором ни база /2 , ни база /3 не обеспечивают выбранных величин Д2у, Д3у . Для этого сектора по изложенной методике находится дополнительная база /4 . Процесс выбора дополнительных баз продолжается, пока во всех секторах неоднозначности базы /1 не будет разности фаз Ду, равной принятой величине.

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

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

Результаты расчётов ряда структур антенных решёток фазовых пеленгаторов с числом баз от двух до пяти приведены в таблице, где они представлены соотношением баз.

В графах таблицы указаны величины Дгу, для которых производились расчёты. Из таблицы

видно, что существует целый ряд антенных структур, для которых Дгу не больше заданных величин.

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

Сравним между собой по вероятности правильного устранения неоднозначности два пеленгатора. Пусть один из них имеет антенную решетку, оптимизированную по критерию максимального правдоподобия, а другой имеет антенную решетку, найденную предложенным в статье методом, но не являющуюся оптимальной по указанному критерию. Предположим, что обработка совокупности результатов фазовых измерений выполняется в соответствии с принципом максимального правдоподобия и все необходимые для этого условия выполняются. Естественно, габаритные размеры антенных решеток и число их элементов должно быть одинаковым. Примем, что число элементов решетки N =4, а количество фазометрических баз п = N -1=3. Будем считать, что в первом случае (оптимизированная решетка) фазометры подключены к приемным трактам по схеме «с опорной антенной». Структурная схема такой системы приведена на рис. 2, а. Видно, что наибольшая база определяет размеры всей антенной системы.

Ао V

Аз V

А 2

• ф2.

А 1 V

•Г

А о V

А 3 V

А 2 V

• <р2.

А1 V

.ш.

а б

Рис. 2. Схема организации фазометрических баз на антенной решетке: а - пеленгатор с опорной антенной; б - «свернутая» схема организации баз

Среди приведенных в таблице трехбазовых пеленгаторов оптимальными по критерию максимального правдоподобия являются системы со структурами антенной системы, задаваемой цело-

Т т т

численными векторами ех =(6,3,2) , ех =(6,4,1) , ех =(9,6,4) . Примем для сравнения пеленгатор

Т

со структурой антенной системы ех =(9,6,4) . Для удобства читателя заметим, что на рис. 2, а изменена закономерность обозначений баз, принятая в монографиях [1, 2], которые являются основой для анализа: в монографиях индексом 1 обозначена наименьшая база, а остальные индексы увеличиваются в порядке увеличения размера баз; на рисунке, как это принято в статье, индексом 1 обозначена наибольшая база.

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

рис. 2, б, ех = (9,6,4) . Такой вектор имеется в таблице, где он представлен соотношением баз.

Структуры антенных решеток многобазовых пеленгаторов

Вари- Число Соотношение баз: еХ1, еХ2, еХ3, еХ4, еХ5

ант баз ех1 еХ2 еХ3 еХ4 еХ5

1 2 2 1

2 3 1, 2

3 4 2

4 6 3

5 3 6 2, 4 1, 3, 5

6 9 3, 6

7 8 4

8 12 6

9 12 4, 8

10 12 6

11 4 18 9 3, 6, 12, 15 1, 2, 4, 5, 7, 8, 10, 11,

12 18 6, 12 3, 9, 15 1, 2, 4, 5, 7, 8, 10, 11,

13 18 6, 12 2, 4, 8, 10, 14, 1, 3, 5, 7, 9, 11, 13, 15,

9, 18 (120°) 3, 6, 12, 15, 1, 2, 4, 5, 7, 8, 10, 11,

14 27 21, 24 13, 14, 16, 17, 19, 20,

15 5 16 8

16 24 12 4, 8, 16, 20 2, 6, 10, 14, 18, 22 1, 3, 5, 7, 9, 11, 13, 15, 17, 19,

17 24 8, 16 4, 12, 20 2, 6, 10, 14, 18, 22 1, 3, 5, 7, 9, 11, 13, 15, 17, 19,

5

18 24 12 6, 18 2, 4, 8, 10, 14, 16, 20, 22 1, 3, 5, 7, 9, 11, 13, 15, 17, 19,

19 24 12 6, 18 3, 9, 15, 21 1, 2, 4, 5, 7, 8, 10, 11, 13, 14, 16,

20 36 18 6, 12, 24, 30 3, 9, 15, 21, 27, 33 1, 2, 4, 5, 7, 8, 10, 11, 13, 14, 16, 17, 19, 20, 22, 23, 25, 26, 28, 29, 31, 32, 34, 35

21 5 36 12, 24 6, 18, 30 3, 9, 15, 21, 27, 33 1, 2, 4, 5, 7, 8, 10, 11, 13, 14, 16, 17, 19, 20, 22, 23, 25, 26, 28, 29, 31, 32, 34, 35

18 (180°) 9, 27 (180°) 3, 6, 12, 15, 21, 24, 30, 1, 2, 4, 5, 7, 8, 10, 11, 13, 14, 16,

22 36 33

23 4 7 3, 4 2, 5 1, 6

24 14 6, 8 4, 10 2, 12 1, 3, 5, 7, 9, 11, 13

25 5 21 9, 12 6, 15 3, 18 1, 2, 4, 5, 7, 8, 10, 11, 13, 14, 16,

26 3 5 2, 3 1, 4

27 4 10 4, 6 2, 8 1, 3, 5, 7, 9

На схеме рис. 2, б имеет место полная система линейно независимых фазометрических баз. При выполнении этого условия вероятность правильного устранения неоднозначности /0 не зависит от того, как подключены фазометры к приемно-измерительным трактам, а зависит только от

структуры решетки [2]. Перейдем к подключению фазометров по схеме «с опорной антенной» и

Т

получим вектор баз еХ =(9,6,5) .

Методика расчетов вероятности / для пеленгаторов, реализующих метод максимального правдоподобия, приведена в монографии [2] и в трехбазовом случае сводится к вычислению двумерного нормального интеграла. На рис. 3 приведена рассчитанная по этой методике зависимость вероятности правильного устранения неоднозначности / от СКО фазовых погрешностей для опТ

тимизированного пеленгатора со структурой антенной системы еХ =(9,6,4) , а также со структурой

Т

еХ = (9,6,5) , полученной на основании предложенного метода.

0, 0,98 0,97 0,96 0,95 0,94 0,93 0,92 0,91 0,9

Рис.

21 24 27

СКО разности фаз, град 3. Зависимость вероятности правильного устранения неоднозначности от СКО фазовых погрешностей

Как видно из графиков, при малых фазовых погрешностях (СКО не больше 20 град), вероятность р практически одинакова для сравниваемых структур, но становится существенно различной при дальнейшем увеличении СКО: при СКО, равном 36 град, разница достигает 0,05.

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

Идея метода принадлежит Г.Г. Порубову. В.П. Денисов по его просьбе подготовил представленные материалы к опубликованию в журнале. Нетрудно видеть, что метод применим не только к пеленгаторам, но и к другим видам многошкальных фазовых измерительных систем.

Литература

1. Белов В.И. Теория фазовых измерительных систем. - Томск: Изд-во ТУСУРа, 1994. - 102 с.

2. Денисов В.П. Фазовые радиопеленгаторы / В.П. Денисов, Д.В. Дубинин. - Томск: Изд-во ТУСУРа, 2002. - 252 с.

Порубов Геннадий Гаврилович

Инженер, ОАО «НИИАП», г. Новосибирск

Тел.: 8-961-871-47-25

Эл. почта: porub27@mail.ru

Денисов Вадим Прокопьевич

Д-р техн. наук, профессор каф. радиотехнических систем ТУСУРа

Тел.: 8 (382-2) 41-36-70

Эл. почта: dvp15@sibmail.com

Porubov G.G., Denisov V.P.

Method to calculate the structures of antenna arrays for multibase phase direction finders

The paper considers the method that can be used to calculate the structures of antenna arrays for multibase phase direction finders. As the criterion for calculating the relative size of bases was taken the condition, that ensures obtaining of the specified amount of deviation of the expected value in the ambiguity areas of the maximum base-efficient sign is also presented. The paper gives the results of mathematical modeling of phase direction finders.

Keywords: direction finder, phase difference, ambiguity resolution.