CC BY
138
УДК 621.3.095.22
А.В. Новиков, В.А. Хлусов
Способ формирования диаграммы направленности антенной решетки на прием с низким уровнем боковых лепестков
Предложен способ формирования диаграммы направленности (ДН) линейной решетки на прием с низким уровнем боковых лепестков (УБЛ) без расширения основного. Способ основан на современном методе спектрального оценивания Spatially Variant Apodization (SVA). Математическое моделирование показывает работоспособность способа при наличии случайных ошибок в амплитудно-фазовом распределении (АФР) поля. Ключевые слова: линейная решетка, приемная антенна, адаптивный фильтр, пространственный фильтр. doi: 10.21293/1818-0442-2017-20-2-5-9
В настоящее время в современных речных и морских радиолокаторах широко используются антенны в виде линейных решеток. Это, например, радары «Raymarine» и «Река» АО НПФ «Микран». В последнем реализована активная линейная антенная решетка из N элементов с N малошумящими усилителями, N смесителями и одним аналоговым сумматором [1].
Угломестная ширина диаграммы направленности (ДН) таких решеток составляет 25...30 градусов для обеспечения достаточного обзора пространства. Однако ширина по азимуту задается на уровне одно-го-двух градусов для достаточного разрешения наблюдаемых объектов по углу; эта ширина зависит от отношения длины волны поля к длине решетки (к ее апертуре).
Линейная решетка состоит из N антенных элементов, разнесенных обычно на половину длины волны поля. Выходной сигнал антенны в режиме приема равен взвешенной сумме сигналов всех антенных элементов. Весовая функция обычно задается так, что ее значения спадают к краям апертуры антенны, что позволяет снизить уровень боковых лепестков (УБЛ) по сравнению со случаем «прямоугольной» весовой функции (-13 дБ). Отметим, что способ весового суммирования работает как на прием, так и на передачу.
Весовое суммирование сигналов с целью снижения УБЛ называется аподизацией, причем линейной.
Существует множество весовых окон [2]: Блэк-мана, Кайзера, Хеннинга, Хемминга и т.д. Они отличаются, например, УБЛ, шириной основного лепестка и скоростью спадания боковых лепестков. В целом, чем ниже УБЛ, тем шире основной лепесток. Таким образом, весовое суммирование ограничено тем, что чем ниже требуется УБЛ, тем длиннее должна быть антенна для получения заданного разрешения по углу.
В данной статье предлагается способ снижения УБЛ ДН решетки на прием, сохраняющий ширину основного лепестка. Способ основан на нелинейной многооконной аподизации, которая известна в зарубежной периодике как Spatially Variant Apodization (SVA) [3-8]. SVA широко используется при синтезе
апертуры [4, 9]. Предложенный способ позволяет снять ограничение весового суммирования в плане расширения основного лепестка.
Весовое суммирование с последующим вычислением квадрата модуля равносильно оценке мощности нулевой пространственной частоты. Так оценивается уровень сигнала, приходящего с нормального к апертуре решетки азимутального направления. Предложенный способ требует дополнительной оценки уровней сигналов, которые приходят с направлений, соответствующих первой и минус первой пространственным частотам; при этом весовая функция равна константе.
8УА - это один из способов спектрального оценивания сигналов, позволяющий оценивать амплитуды всех частот в спектре, давая низкий УБЛ и не расширяя основной лепесток. Предложенный способ требует вычисления амплитуд всего лишь трех пространственных частот для текущего азимутального положения антенны. По трем комплексным амплитудам с помощью параметрического фильтра формируется амплитуда, соответствующая оценке уровня нулевой пространственной частоты. Физическое вращение антенны даст оценку всего пространственного спектра, который пропорционален ее ДН.
Подготовительные формулы (8УЛ)
ДН и АФР антенны для малых углов связаны между собой преобразованием Фурье [10]. Для линейных решеток используется одномерное дискретное преобразование Фурье (ДПФ)
N-1
^=1 sne ^ ы;, (1)
п=0
где - сигнал с выхода п-го антенного элемента; к - индекс пространственной гармоники, -Ы/2< к < N/2.
Индекс к соответствует углу прихода плоской волны [10]
а = агс5ш(кХ/Ь),
где Х - длина волны падающего поля; Ь - размер апертуры антенны (длина).
Если рассмотреть нормальное падение плоской волны на антенную решетку (рис. 1), только нулевая
гармоника будет иметь ненулевую мощность, т.е.
Ы2 *0.
ф-const
\L
YYYYYYYY
Антенная • решетка
Рис. 1. Нормальное падение плоской волны на антенную решетку
Если рассмотреть падение плоской волны на решетку под любым углом из счетного множества (рис. 2)
6 = агс5ш(Х/I), дей ,
то только гармоника с индексом д будет иметь нену-I |2
левую мощность, т.е. Ф0.
YYYYYYYY
I Личинная # реш&гка Рис. 2. Падение плоской волны на антенную решетку под углом 6
Сказанное справедливо при выполнении условия полуволнового размещения антенных элементов.
В общем случае пространственные гармоники будут иметь ненулевую мощность, так как угол падения плоской волны не обязан принадлежать счетному множеству. При этом возникает эффект растекания спектра (spectral leakage) и появляются ложные сигналы («боковые лепестки»), хотя изначально на решетку падала одна плоская волна, приходящая с одного направления. Дополнительно сама решетка имеет произвольное азимутальное положение.
Весовое суммирование вместо (1) предполагает взвешенное ДПФ
N-1
i2nknl
Vk = Z wnsn e
I N I
- N/2< k < N/2, (2)
n=0
где - отсчеты оконной (взвешивающей) функции. Данный способ приводит к расширению основного лепестка в спектре, в целом тем большему, чем ниже УБЛ.
Известно, что (2) может быть записано в спектральной области через свертку
Ук = 1 ^_, (3)
гей
где - ДПФ оконной функции м>„ подобно (1).
Вычисление (3) потребует минимум ресурсов, если подобрать такую оконную функцию, ДПФ от которой содержит минимальное число частот. Как показано в [3], в качестве «легкой» функции подходит ряд по косинусам кратных частот
K / l4k (2%kn wn = Z(-1) akcosl——-
k=0 ^ N
0 < n < N .
(4)
Тогда можно записать явное выражение для W¡
^ =11 НЧ (8г_к +§г+к), (5)
2 к=0
где 5г - символ Кронекера.
Ограничимся в (5) окном первого порядка, т.е. К = 1. Это позволит минимизировать вычислительные затраты, не теряя в практически реализуемом УБЛ, который ограничен точностью настройки антенных элементов. Подставляя (5) в (3), получим
Ук = а0Як _0,5а1 ((_1 + рк+1).
Как показано в [3], коэффициенты а0 и а1 могут быть выбраны и удобно обозначены следующим образом: а0 =1, а1 = а , 0 < а <1. Таким образом, если а = 0, взвешивающая функция переходит в прямоугольное окно; при а = 1 - в окно Хеннинга. В этом и кроется суть параметричности 8УЛ. Подчеркнем еще раз, что спектральные амплитуды рк вычислены без взвешивания согласно (1).
Итак, для текущего индекса гармоники к параметр а выбирается так, чтобы обеспечить минимум квадрата модуля спектральной амплитуды, но при ограничении 0 < а < 1. Эта оптимизация была проделана в [3] для К = 1,2,3 . Ограничение на а позволяет сохранять основной лепесток; критерий минимума может дать такое а , которое выйдет за рамки ограничения и окно при этом потеряет свою функцию ослабления сигнала на краях.
В данной работе К = 1, поэтому искомый параметр а может быть вычислен как
a = Re
Sk
(6)
|0,5( _1 + Як+1)1 '
В соответствии со способом 8УЛ оценка спектра делается по алгоритму
Як _0,5а(Як_1 + Як+1), если 0 <а <1;
Ук = Ьк, если а <0; (7)
Як _ 0,5 (Як_1 + Як+1), если а > 1.
Делить числитель на знаменатель в (6) не требуется: достаточно сравнить числитель и знаменатель между собой и вынести одно из трех решений (7). В реальности в комплексных амплитудах будет присутствовать шум, поэтому деление на ноль исключено (в знаменателе после преобразования к вещественным квадратурам будет квадрат модуля).
Аппаратная функция анализатора спектра (7) имеет такую же ширину основного лепестка, как и
для анализатора (1) с прямоугольной взвешивающей функцией, зато УБЛ значительно ниже (по крайней мере -31,5 дБ против -13 дБ). Более того, скорость спадания боковых лепестков аппаратной функции также будет выше (18 дБ/октаву против 6 дБ/октаву). Предлагаемый способ формирования ДН Предположим, что на линейную антенную решетку нормально падает плоская волна (см. рис. 1). Тогда суммирование (1) для к = 0 даст максимальный уровень выходного сигнала. Известно, что если антенна при этом вращается в азимутальной плоскости, то такое суммирование даст форму ДН решетки
по азимуту. Величина |^о|2, зависящая от азимутального угла ß, по форме будет совпадать с функ-2
цией (sinx/x) , где x = C sinß ~ Cß для малых углов, C - константа.
Способ SVA (7) предполагает вычисление амплитуд всех пространственных гармоник. Мы предположили, что хорошие свойства SVA сохранятся,
если вычислять лишь j0|2 для произвольного азимутального положения решетки. Математическое моделирование подтвердило наше предположение. Результаты моделирования Рассмотрим случай падения двух плоских волн на вращающуюся линейную антенную решетку (рис. 3). Согласно рисунку можно записать комплексную амплитуду узкополосного сигнала на выходе n-го антенного элемента
sn = A1e j™smß+ A2e j™sin(e+ß) (8)
где ß - азимут антенны против часовой стрелки; e - угол между фронтами двух плоских волн по часовой стрелке.
со -
Рис. 3. Падение двух плоских волн на вращающуюся антенную решетку
На рис. 3 углы р и 6 положительны. Пусть амплитуды плоских волн равны единице, т.е.
А1 = А2 =1.
Для некоторых углов р и 6 (в момент «фотосъемки») принятый суммарный сигнал (8) преобразуется в спектральную область
.[2лп] , . [2яп
N-1 Я
N-1
S0 = Z sn , S-1 = Z sn e
n=0 n=0
N
N-1 -í\- Ar Si = Z Sn e I N
n=0
после чего по (6) и (7) вычисляются параметр a и амплитуда выходного сигнала антенны Vq соответственно.
Введем в каждый антенный элемент амплитудно-фазовые ошибки с помощью коэффициента передачи
Kn = wn10°'05AAn ejAф" , где wn - коэффициент передачи идеального n-го элемента, пусть wn = 1; AAn - амплитудная ошибка n-го элемента, дБ; Афп - фазовая ошибка, радиан.
Предположим, что амплитудная ошибка распределена нормально с заданным средним квадратичным отклонением, а фазовая ошибка - равномерно на интервале (-фтах...фтах). Примем также, что
ошибки в антенных элементах независимы.
Современные цифровые фазовращатели, например продукции Analog Devices, HMC642A, дают q,5 дБ амплитудной ошибки и интервал фазовой (-10°...10°).
Пусть две падающие плоские волны разделены по углу на координату первого бокового лепестка, тогда пространственная частотаp будет равна
sin0 1,5 . 0 1,5А, 3 (10)
p=-=—^ sin0=—--=—. (10)
А L (А/2)N N
Рассмотрим пример 128-элементной антенной решетки, т.е. N = 128.
Рассчитанные ДН покажем в области малых углов, так как Фурье-связь между ДН и АФР справедлива лишь для малых углов [10].
На рис. 4 показаны десять реализаций ДН антенной решетки на прием, полученных предложенным способом и способом простого суммирования (1).
Предлагаемый способ формирования ДН
i г
М| Г Ii МП Д i НШит Т !i 1/Ш ДО. i|¿ ш пк ÉÉt ш лк mi iff
|ГТР)|2, дБ
о -10 -20 -30 -40 -50 -60
|S0(ß)f, ДБ Способ простого суммирования
0 -10 -20 -30 -40 -50
3 4
ß, град
-60
! ! ! [
V
1 lí
i í i i i
-1
-з
-2
-1
(9)
и 1 2 3 4
в, град
Рис. 4. Десять реализаций ДН 128-элементной антенной решетки на прием
Из рис. 4 следует, что простое суммирование дает УБЛ - 13 дБ, малочувствительный к ошибкам АФР. Ширина основного лепестка совпадает с расчетной
© = 2агсзт(2/ N ) = 2а1гат(2/128)«1,8°.
Предлагаемый способ сохраняет ширину основного лепестка и его амплитуду и дает УБЛ не хуже -30 дБ. Разброс УБЛ в этом случае заметно больше, чем разброс для простого суммирования.
Рисунок 5 аналогичен рис. 4; отличие лишь в наличии второго сигнала, принимаемого по первому боковому лепестку (10).
|ГТР)|2, дБ
Предлагаемый способ формирования ДН
|50(Р)|% дБ Способ простого суммирования
1 6
в, град
-L
-2
4 6
в, град
Рис. 5. Десять реализаций двулучевых ДН 128-элементной антенной решетки на прием
Из рис. 5 следует, что предлагаемый способ, по сравнению с простым суммированием, так же разрешает два сигнала (провал -7,5 дБ), один из которых принимается по основному, а другой - по первому боковому лепестку. При этом средний УБЛ возрастает на 3...4 дБ по сравнению со случаем одного сигнала.
Если использовать взвешенное суммирование с окном Хеннинга (по (4) при K = 1 и ao =1), то два рассматриваемых сигнала не разрешаются (провал -0,5 дБ), но обеспечивается УБЛ не хуже -25 дБ.
Если рассмотреть вариант приема двух сигналов, разделенных по углу на координату первого нуля, т. е.
sin6 1 _ X 2
X L
>sin 6=-
(X/2) N N
то максимумы, соответствующие оценкам угла прихода принятых сигналов, опустятся на 3,5 дБ и совпадут с результатом расчета для окна Хеннинга, причем вплоть до провала между этими максимума-
ми (-22 дБ). Расчет с прямоугольным окном, наоборот, даст подъем максимумов на 0,7 дБ (провал между максимумами -18 дБ). Это позволяет сделать вывод о том, что в предлагаемом способе комбинируются оценки спектра по многим оконным функциям из заданного класса и выбирается та функция, кото -рая дает минимум оцениваемой мощности. Причем минимальным угловым разрешением для предлагаемого способа является (10), а разрешение (11) подходит лишь для прямоугольной оконной функции, которая дает максимальный из возможных УБЛ (-13 дБ).
Выводы
Предложен способ формирования ДН приемной АР, основанный на применении современных методов спектрального оценивания (SVA) и позволяющий существенно снизить УБЛ без расширения основного лепестка ДН АР. Необходимо отметить, что предлагаемый способ неприменим для передающих антенн, поскольку является параметрическим способом формирования ДН на основе информации о принимаемом сигнале, полученной по нескольким направлениям приема.
Численное моделирование показывает, что для 128-элементной решетки достижим УБЛ не хуже -30 дБ при наличии случайных амплитудных (0,5 дБ) и фазовых (±10°) ошибок; при этом сохраняется ширина основного лепестка. Предложенный способ обеспечивает одинаково (по сравнению с простым суммированием) угловое разрешение двух сигналов, разнесенных на координату первого бокового лепестка.
Исследование способа SVA в двумерной области [4] позволяет расширить предложенный способ формирования ДН на плоские решетки.
Работа выполнена при финансовой поддержке Министерства образования и науки РФ, соглашение № 14.577.21.0188 от 27.10.2015 г., идентификатор RFMEFI57715X0188.
Литература
1. Доценко В.В. Повышение энергетического потенциала РЛС с непрерывным ЛЧМ-сигналом / В.В. Доценко, М.В. Осипов, В. А. Хлусов // Доклады ТУСУРа. - 2011. -№ 1(23). - С. 29-33.
2. Harris F.J. On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform // Proc. IEEE. -1978. - Vol. 66. - P. 51-83.
3. Pat. 6 298 363 US. Adaptive windowing of FFT data for increased resolution and sidelobe rejection / R.J. Ian-nuzzelli (US). - No 09/191 665; submitted 13.11.98; issued 02.10.01.
4. Stankwitz H.C. Nonlinear Apodization for Sidelobe Control in SAR Imagery / H.C. Stankwitz, R.J. Dallaire, J.R. Fienup // IEEE Transaction on Aerospace and Electronic Systems. - 1995. - Vol. 31, No. 1. - P. 267-279.
5. Smith B.H. Generalization of spatially variant apodization to noninteger Nyquist sampling rates // IEEE Transactions on Image Processing. - 2000. - Vol. 9, No. 6. - P. 10881093.
6. Lee J.A.C. Effectiveness of spatially-variant apodization / J.A.C. Lee, D.C. Munson // Proc. IEEE ICIP. - 1995. -Vol. 1. - P. 147-150.
7. Lee J.A.C. Spatially Variant Apodization for Image Reconstruction from Partial Fourier Data / J.A.C. Lee, D.C. Munson // IEEE Transactions on Image Processing. -2000. - Vol. 9, No. 11 - P. 1914-1925.
8. Stankwitz H.C. Advances in non-linear apodization / H.C. Stankwitz, S.P. Taylor // IEEE Transactions on Aerospace and Electronic Systems. - 2006. - Vol. 21, No. 1. - P. 3-8.
9. Resolution enhancement for LASAR 3D imaging via 11 regularization and S VA / G. Xiang, X. Zhang, J. Shi, S. Wei // Progress In Electromagnetic Research M. - 2015. -Vol. 41. - P. 95-104.
10. Levy M. Parabolic equation methods for electromagnetic wave propagation. - London: The Institution of Electrical Engineers, 2000. - 336 p.
Новиков Анатолий Викторович
Канд. техн. наук,
доцент каф. радиотехнических систем (РТС) ТУСУРа
Тел.: +7-952-155-96-99
Эл. почта: naw1983@yandex.ru
9
Хлусов Валерий Александрович
Д-р техн. наук, вед. науч. сотр. НИИ РТС ТУСУРа Тел.: +7 (382-2) 41-35-62 Эл. почта: hva@micran.ru
Novikov A.V., Khlusov V.A.
Receiving Linear Array Beamforming Using Spatially Variant Apodization
The authors present a novel beamformer to receive linear array based on Spatially Variant Apodization. This technique reduces the side-lobe level compared with a uniform linear array, at the same time keeping the main lobe beam-width. For example, the side-lobe level of a 128-element linear array would be better than -30 dB in the presence of phase and amplitude errors based on the proposed technique (5.8 degrees r.m.s. and 0.5 dB r.m.s. respectively).
Keywords: Linear antenna arrays, receiving antennas, adaptive filters, spatial filters.
