RU2599257C1 - Способ пространственной обработки радиосигналов - Google Patents

Способ пространственной обработки радиосигналов Download PDF

Info

Publication number
RU2599257C1
RU2599257C1 RU2015151075/07A RU2015151075A RU2599257C1 RU 2599257 C1 RU2599257 C1 RU 2599257C1 RU 2015151075/07 A RU2015151075/07 A RU 2015151075/07A RU 2015151075 A RU2015151075 A RU 2015151075A RU 2599257 C1 RU2599257 C1 RU 2599257C1
Authority
RU
Russia
Prior art keywords
sin
aperture
signal
channel
channels
Prior art date
Application number
RU2015151075/07A
Other languages
English (en)
Inventor
Борис Николаевич Горевич
Original Assignee
Борис Николаевич Горевич
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Борис Николаевич Горевич filed Critical Борис Николаевич Горевич
Priority to RU2015151075/07A priority Critical patent/RU2599257C1/ru
Application granted granted Critical
Publication of RU2599257C1 publication Critical patent/RU2599257C1/ru

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14Systems for determining direction or deviation from predetermined direction
    • G01S3/46Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/292Extracting wanted echo-signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

Изобретение относится к радиотехнике и может быть использовано в системах пеленгования узкополосных сигналов с известными несущей частотой, в том числе в радиолокации, радионавигации, связи. Достигаемый технический результат - сужение основного лепестка характеристики угловой направленности, формируемой на выходе системы пеленгования, повышение углового разрешения принимаемых сигналов. Указанный результат достигается за счет того, что способ пространственной обработки сигналов заключается в приеме с помощью фазированной антенной решетки с плоским прямоугольным раскрывом, имеющим размеры X×Y, сигнала с плоским волновым фронтом, приходящего с направления (αсс), определяемого углами пеленгования α, β, отсчитываемыми в прямоугольной системе координат с осями, ориентированными вдоль сторон раскрыва, от нормали к раскрыву, восстановленной из опорной точки, находящейся в одном из углов раскрыва, его согласованной пространственной фильтрации и квадратурного детектирования суммы напряжений всех приемных элементов после их фильтрации, отличающийся тем, что для обработки напряжения, снимаемого с каждого приемного элемента решетки, используют по четыре канала с выполнением в них одинаковых для данного приемного элемента операций согласованной пространственной фильтрации, причем в трех каналах в обрабатываемые в них напряжения дополнительно вводят фазовые сдвиги, равные: kX(sin(α)-sin(αc)) - для одного канала, kY(sin(β)-sin(βc)) - для другого канала, kX(sin(α)-sin(αc))+kY(sin(β)-sin(βc)) - для третьего канала, где k=2π/λ, λ - длина волны принимаемого сигнала. 9 ил.

Description

Изобретение относится к радиотехнике и может быть использовано в системах пеленгования узкополосных сигналов с известной несущей частотой, в том числе в радиолокации, радионавигации, связи.
Наиболее близким аналогом заявляемого способа является классический способ согласованной пространственной фильтрации, используемый для обработки узкополосных сигналов с плоским волновым фронтом, принимаемых фазированными антенными решетками (ФАР). Реализация способа включает сдвиг фаз гармонических напряжений, снимаемых с выхода каждого приемного элемента ФАР и взвешенное суммирование этих напряжений. При этом фазовые сдвиги, придаваемые напряжениям, снимаемым с различных точек раскрыва антенны, для их когерентного суммирования, линейно зависят от координаты данной точки раскрыва и определяются лишь направлением на источник сигнала, а веса суммируемых напряжений определяются распределением амплитуды сигнала по раскрыву [Пространственно-временная обработка сигналов / И.Я. Кремер, А.И. Кремер, В.М. Петров и др.: Под ред. И.Я. Кремера. - М.: Радио и связь, 1984, с. 55-57; Чердынцев В.А. Радиотехнические системы: Учеб. пособие для вузов. - Мн.: Выш. шк., 1988, с. 124-125]. Для приемного элемента (точки раскрыва) с координатой x∈[0, Х] линейной антенной решетки для фронта волны, приходящего на раскрыв антенны под углом αс, отсчитываемым от нормали к краю раскрыва, придаваемый сдвиг фазы равен -kxsin(α), где k=2π/λ - волновое число λ - длина волны. При угле пеленгации α=αc происходит синфазное суммирование напряжений на выходах каналов приемных элементов, при этом суммарный сигнал максимален [Григорьев Л.Н. Цифровое формирование диаграммы направленности в фазированных антенных решетках. - М.: Радиотехника, 2010, с. 13, 19; Теоретические основы радиолокации / Под ред. В.Е. Дулевича. - М.: Сов. радио, 1964, с. 68; Многофункциональные радиолокационные системы: учеб. пособие для вузов / П.И. Дудник, А.Р. Ильчук, Б.Г. Татарский; под ред. Б.Г. Татарского. - М.: Дрофа, 2007, с. 43].
Зависимость суммы напряжений на выходах каналов антенной решетки от угла пеленгации α характеризуется диаграммой направленности антенны. Диаграмма направленности антенны определяет угловые разрешающие возможности системы пеленгования сигналов, использующей ФАР.
При реализации известного способа пространственной обработки сигналов в системе пеленгования с линейной антенной решеткой с раскрывом размером X, для случая непрерывной интерпретации дискретного раскрыва решетки, что возможно, когда расстояние между приемными элементами решетки меньше длины полуволны, при приеме сигнала с равномерным амплитудным распределением по раскрыву, диаграмма направленности антенны описывается формулой [Коростылев А.А. Пространственно-временная теория радиосистем: Учеб. пособие для вузов. - М.: Радио и связь, 1987, с. 83; Справочник по радиолокации / Под ред. М.И. Сколника. пер с англ. под общей ред. В.С. Вербы. В 2 книгах. Книга 1. Москва: Техносфера, 2014, с. 609-611]
Figure 00000001
Формула (1) описывает многолепестковую функцию угла пеленгации α с одним главным лепестком. Угловая разрешающая способность системы пеленгования, реализующей способ пространственной обработки сигнала, определяется шириной главного лепестка диаграммы направленности. Как следует из формулы (1), расстояние между нулями главного лепестка диаграммы направленности по координате sin(α) при реализации известного способа пространственной обработки сигнала, равно 2λ/X. Эта величина характеризует предельную (потенциальную) разрешающую способность известного способа пространственной обработки [Коростылев А.А. Пространственно-временная теория радиосистем: Учеб. пособие для вузов. - М.: Радио и связь, 1987. - С. 66].
Недостатком известного способа пространственной обработки сигналов является недостаточно высокая угловая разрешающая способность системы пеленгации, потенциально ограниченная, для случая прямоугольного раскрыва с размерами X×Y, величинами 2λ/X, 2λ/Y в соответствующих угловых плоскостях (по нулевым значениям главного лепестка диаграммы направленности).
Изобретение направлено на решение задачи повышения углового разрешения принимаемых системой пеленгации сигналов.
Техническим результатом использования предлагаемого способа является сужение основного лепестка характеристики угловой направленности, формируемой на выходе системы пеленгования.
Достижение заявленного технического результата основывается на использовании дополнительных каналов обработки напряжений, снимаемых с приемных элементов решетки, с введением в них специальных фазовых сдвигов.
Схема, поясняющая прием волны антенной с линейным раскрывом, приведена на фиг. 1. Здесь используются обозначения:
0x - система координат раскрыва;
α - угол пеленгации сигнала, отсчитываемый по часовой стрелке от нормали No, восстановленной в плоскости 0xNo из начала раскрыва;
αс∈[0, ±π/2] - угол, характеризующий направление на источник сигнала;
kxsin(αc) - фазовая задержка фронта волны в x-й, x∈[0, X], точке раскрыва антенны, отсчитываемой относительно опорной точки.
Принципиальным для предлагаемого способа является выбор опорной точки для отсчета фаз фронта падающей на раскрыв волны. В предлагаемом способе в качестве точки опоры используется крайняя точка раскрыва (применительно к фиг. 1 - нулевая точка).
Пространственная интерпретация предлагаемого способа заключается в формировании дополнительного виртуального раскрыва антенны и совместной обработке напряжений, снимаемых с реального и виртуального раскрывов с общей опорной точкой, находящейся в центре объединенного раскрыва. Формируемый виртуальный раскрыв антенны на фиг. 1 показан полужирной штриховой линией.
Для раскрытия изобретения рассмотрим предлагаемый порядок обработки напряжений, снимаемых с приемных элементов линейного раскрыва. Схема построения приемных каналов ФАР и выполнения в них обработки напряжений представлена на фиг. 2. Здесь обозначено:
1 - раскрыв ФАР;
2 - управляемый фазовращатель;
3 - выключатель «P(α)»-«F(α)» дополнительных каналов;
4 - сумматор;
5 - смеситель;
6 - интегратор;
7 - фазосдвигающий на π/2 элемент;
8 - устройство детектирования.
Для каждого приемного элемента решетки используется по два канала обработки напряжений - основной, в котором реализуется известный способ фильтрации снимаемого с приемного элемента гармонического напряжения, и дополнительный, в котором кроме известного способа фильтрации применяется дополнительная фильтрация напряжения.
Создаваемое падающей на раскрыв волной напряжение в произвольной точке раскрыва x, отсчитываемой от опорной точки, описывается формулой
Figure 00000002
.
где ϕоп - начальная фаза сигнала в опорной точке;
Figure 00000003
- несущая частота сигнала.
Напряжения u(t,x,ϕоп) подвергаются фильтрации в каналах приемных элементов известным способом - за счет введения фазового сдвига -kxsin(α), компенсирующего имеющуюся фазовую задержку фронта волны относительно опорной точки.
Компенсирующий фазовый сдвиг -kxsin(α) может вводиться одним из двух способов - фильтровым, с помощью управляемого фазовращателя, или с помощью пространственного коррелятора, включающего смеситель принятого и гетеродинного напряжения, в которое введен указанный канальный фазовый сдвиг [Теоретические основы радиолокации: Учеб. пособие для вузов / А.А. Коростылев, Н.Ф. Клюев, Ю.А. Мельник и др., под ред. В.Е. Дулевича - М.: Сов. радио, 1978. - С. 85].
На фиг. 2 в качестве примера использован фильтровой способ введения фазовых сдвигов в каналах приемных элементов ФАР. В результате введения компенсирующего фазового сдвига напряжение в основном приемном канале, соответствующем точке раскрыва x, описывается формулой
Figure 00000004
Напряжение (2) основного канала со скомпенсированным фазовым сдвигом, согласно известному способу обработки, поступает непосредственно в сумматор 4. Одновременно, в отличие от известного способа, через выключатель 3 это напряжение поступает в дополнительный канал обработки напряжения, обеспечивающий формирование виртуального раскрыва антенны.
В качестве опорной точки для обработки принятой виртуальным раскрывом волны используется край, противоположный по отношению к точке опоры реального раскрыва. Для фиг. 1 - это точка X. Набег фазы в этой точке равен Ф(αс)=kXsin(αc). Напряжения в точках x виртуального раскрыва получают из соответствующих напряжений реального раскрыва за счет введения в них сдвига фазы Ф(α)-Ф(αс), где Ф(α)=kXsin(α). В итоге напряжение на выходе дополнительного канала фильтрации сигнала, после введение фазовых сдвигов -Ф(αс) и Ф(α) имеет вид:
Figure 00000005
Подстановка значений Ф(αс), Ф(α) в (3) дает:
Figure 00000006
Отметим, что формально, учитывая симметричность функции косинуса, в формуле (4) важны не знаки фазовых добавок Ф(αс), Ф(α), а разность этих добавок.
Из анализа фазовых соотношений напряжений (2) и (4) видно, что фронт волны, падая на раскрыв антенны под углом αс, пробегает координаты 0…X реального раскрыва, и одновременно координаты -Х…0 виртуального раскрыва.
Следующая операция пространственной обработки - суммирование напряжений основных и дополнительных каналов всех приемных элементов раскрыва и квадратурное детектирование суммарного напряжения выполняется известным способом.
Суммирование выходных напряжений основного и дополнительного каналов всех x приемных элементов раскрыва для случая непрерывного раскрыва описывается интегралом:
Figure 00000007
После подстановки (2) и (4) в (5) и выполнения интегрирования получим выходное напряжение сумматора, описываемое формулой
Figure 00000008
.
Полученное выходное напряжение сумматора зависит от высокочастотной составляющей сигнала
Figure 00000009
и значения неизвестной фазы ϕоп. Для устранения этой зависимости на практике при выполнении пространственно-временной обработки сигнала выполняют, соответственно, временную фильтрацию напряжения, и его квадратурную обработку.
Для узкополосных сигналов операции пространственной и временной фильтрации являются независимыми и могут выполняться в произвольном порядке, исходя из предпочтений их технической реализации [Пространственно-временная обработка сигналов / И.Я. Кремер, А.И. Кремер, В.М. Петров и др.: Под ред. И.Я. Кремера. - М.: Радио и связь, 1984, с. 56], в том числе как до сумматора (в каналах приемных элементов), так и после сумматора.
Выполнение временной фильтрации сигнала осуществляется умножением напряжения в смесителе на гетеродинное напряжение и интегрированием на интервале обработки сигнала Tобр. Гетеродинное напряжение имеет вид
Figure 00000010
где ϕг - начальная фаза гетеродинного напряжения.
Временная фильтрация сигнала после сумматора для рассматриваемого случая, представленного схемой на фиг. 2, заключается в вычислении интеграла
Figure 00000011
где Δϕ=ϕопг - неизвестная разность фаз входного и гетеродинного напряжений.
После фильтрации по формуле (7) суммарное напряжение имеет вид
Figure 00000012
Устранение зависимости обрабатываемого сигнала от значения неизвестной фазы Δϕ достигается выполнением операции квадратурного детектирования сигнала [Коростылев А.А. Пространственно-временная теория радиосистем: Учеб. пособие для вузов. - М.: Радио и связь, 1987. - С. 36].
Квадратурное детектирование суммарного напряжения заключается в вычислении корня квадратного из квадратов ортогональных по фазе реализаций суммарного напряжения:
Figure 00000013
Выполнение вычислений по формуле (9) с учетом (8) дает функцию пеленгации, описывающую зависимость напряжения на выходе устройства пространственной обработки сигнала от угла пеленгации:
Figure 00000014
Отметим, что при отключении выключателями 3 всех дополнительных каналов, при выполнении операций согласно формулам (5), (7), (9), на выходе устройства пеленгации сигнала формируется напряжение, описываемое функцией F(α) в соответствии с формулой (1).
Полученная функция пеленгации P(α) имеет такой же вид, как и функция диаграммы направленности F(α) - она имеет по координате sin(α) симметричную многолепестковую структуру с одним главным лепестком, направленным под углом αс. Однако ширина главного лепестка функции P(α) по нулевому уровню по координате sin(α) равна λ/X, то есть в два раза меньше ширины главного лепестка функции F(α).
Таким образом, предлагаемый способ пространственной обработки сигнала обеспечивает повышенное по сравнению с известным способом угловое разрешение сигнала.
Рассмотрим особенности предлагаемого способа пространственной обработки сигнала для более общего двумерного случая, когда плоская волна, несущая сигнал, падает на плоский прямоугольный раскрыв антенны. Рисунок, характеризующий приход волны на раскрыв антенны, представлен на фиг. 3.
На фиг. 3 используются обозначения: x, y - координатные оси, привязанные к угловой точке раскрыва, выбранной в качестве точки опоры; X, Y - линейные размеры раскрыва; α, β - углы отсчета направления прихода сигнала относительно нормали No к раскрыву, соответственно в плоскостях 0xNo, 0yNo; (αс, βс) - направление на источник сигнала.
Создаваемое волной напряжение в произвольной точке раскрыва антенны с координатами (x, y), отсчитываемыми от опорной точки, описывается формулой
Figure 00000015
где ϕоп - неизвестная фаза волны в точке опоры с координатами (0, 0).
Предлагаемый способ обработки сигнала для двумерного случая заключается в том, что для координаты (x, y) раскрыва формируются основной и три дополнительных фильтровых канала обработки (по числу угловых точек раскрыва).
В основном канале реализуется известный способ обработки, заключающийся в том, что фильтровым или корреляционным способом в соответствии с значениями координаты (x, y) раскрыва в сигнал (11) вводятся компенсирующие фазовые сдвиги -kxsin(α) и -kysin(β), отсчитываемые относительно опорной точки.
В дополнительных каналах, реализующих пространственную фильтрацию сигнала, принимаемого тремя виртуальными раскрывами антенны, в качестве опорных используют другие угловые приемные элементы решетки, а именно, угловые элементы с координатами (0, X), (0, Y), (Y, X).
В дополнительных каналах помимо общих для координаты раскрыва (x, y) компенсирующих фазовых сдвигов -kxsin(α), -kysin(β), обеспечивающих реализацию известного способа пространственной обработки сигнала, вводят (фильтровым или корреляционным способом) следующие сдвиги фаз:
- для одного канала: -kXsin(αc)+kXsin(α);
- для другого канала: -kYsin(βc)+kYsin(β);
- для третьего канала: -kXsin(αc)+kXsin(α)-kYsin(βc)+kYsin(β).
Так, например, при использовании фильтрового способа введения фазовых сдвигов, напряжения на выходах каналов приемного элемента решетки с координатами (x, y) описываются следующими формулами:
на выходе основного канала:
Figure 00000016
,
на выходах дополнительных каналов, условно обозначенных «X», «Y», «XY», соответственно:
Figure 00000017
,
Figure 00000018
,
Figure 00000019
.
Далее, как и в случае с линейным раскрывом антенны, согласно известному способу обработки, выполняют суммирование напряжений всех каналов со всех приемных элементов раскрыва.
После суммирования приведенных выше напряжений на выходах основного и дополнительных каналов и интегрирования по координатам x∈[0,X], y∈[0,Y], на выходе сумматора получим
Figure 00000020
После временной фильтрации суммарного напряжения (12) с использованием гетеродинного напряжения (6) и квадратурного детектирования в соответствии с формулой (9), получим напряжение на выходе устройства пространственной обработки в виде функции пеленгации:
Figure 00000021
.
Как следует из приведенной формулы, ширина главного лепестка функции пеленгации P(α,β) по нулевому уровню по координатам sin(α) и sin(β) равна, соответственно, λ/X и λ/Y, то есть в два раза меньше ширины главного лепестка функции диаграммы направленности F(α,β).
Таким образом, предлагаемый способ обработки сигнала и для двумерного случая обеспечивает повышенное по сравнению с известным способом угловое разрешение сигнала.
Заметим, что формула (10) одномерного случая является частным случаем формулы (12) двумерного случая.
Далее изобретение раскрывается на примерах технической реализации предлагаемого способа обработки сигнала.
Пример 1. Рассматривается устройство пространственной обработки сигнала с использованием линейной N-элементной антенной решетки размером X. Для демонстрации вариативности способа в отличие от схемы, приведенной на фиг. 2, выбран смешанный корреляционно-фильтровой вариант введения компенсирующих фазовых сдвигов, причем временная фильтрация сигнала совмещена с корреляционным способом введения фазовых сдвигов Ф(α) за счет использования общего смесителя.
Структурная схема устройства пространственной обработки сигнала приведена на фиг. 4. В уточнение и дополнение к обозначениям фиг. 2 здесь используются обозначения:
1 - антенная решетка из N приемных элементов, имеющих нумерацию i=0, 1, …, N-1;
9 - блок управления фазовращателями.
В качестве опорной точки выбран крайний (нулевой) приемный элемент раскрыва решетки.
Принцип функционирования устройства пространственной обработки состоит в следующем. На приемные элементы антенной решетки 1, размещенные с шагом X/N, с направления αс приходит волна. Неизвестная начальная фаза волны, принимаемая опорным элементом решетки, равна ϕоп. Напряжения на выходах i-х элементов антенной решетки описываются формулой
Figure 00000022
.
С помощью управляемых фазовращателей 2 выполняют операции согласованной пространственной обработки - в напряжения ui(t,ϕоп) вводят компенсирующие фазовые сдвиги
Figure 00000023
, соответствующие текущему углу пеленгации α. Фазовые сдвиги δϕi(α) для каждого i-го канала вырабатываются блоком управления фазовращателями 9. С учетом введенных фазовых сдвигов напряжение на выходе фазовращателя 2, присоединенного к i-му приемному элементу антенной решетки, имеет вид
Figure 00000024
.
Через включенные выключатели 3 напряжения uосн.i(α,t,ϕоп) поступают в дополнительные каналы обработки. Здесь в них, с целью формирования виртуального раскрыва, с помощью фазовращателей 2 вводят общий для всех дополнительных каналов фазовый сдвиг -Ф(αс)=-kXsin(αc), обеспечивающий смещение точки опоры на другой край раскрыва (в точку X):
Figure 00000025
.
Далее напряжения дополнительных каналов на смесителях 5 перемножают с гетеродинным напряжением (6), в которое предварительно с помощью смесителя 2 введен фазовый сдвиг -Ф(α)=-kXsin(α):
Figure 00000026
.
Выходные напряжения дополнительных каналов имеют вид
uг.доп.i(α,t,Δϕ)=uдоп.i(α,t,ϕоп)uг.сдв.(t).
Гетеродинное напряжение (6) без фазового сдвига подается на смесители 5 основных каналов с целью выполнения временной фильтрации напряжений. Напряжения на выходах смесителей 5 основных каналов имеют вид:
uг.осн.i(α,t,Δϕ)=uосн.i(α,t,ϕоп)uг(t).
Выходные напряжения основных и дополнительных каналов всех No приемных элементов решетки суммируются в сумматоре 4:
Figure 00000027
.
Сигнал с выхода сумматора подается в интегратор 6 для интегрирования на интервале повторения сигнала T в соответствии с формулой (7):
Figure 00000028
.
Выходное напряжение интегратора U(α,Δϕ) по двум каналам, в одном из которых с помощью фазосдвигающего элемента 7 вводится сдвиг фазы на π/2, поступает в устройство детектирования 8, в котором согласно формуле (9) вычисляется результирующее напряжение в виде функции пеленгации P(α).
В случае отключенных выключателями 3 каналов дополнительной обработки напряжений реализуется известный способ согласованной пространственной фильтрации, и напряжение на выходе устройства детектирования описывается функцией диаграммы направленности F(α).
Для оценки работоспособности предложенной схемы устройства пространственной обработки выполнено численное моделирование в соответствии с описанным порядком действий с сигналом, снимаемым с приемных элементов антенной решетки, для следующих исходных данных: αс=-0.2 рад, N=20, X=0.3 м, λ=0.03 м. Неизвестные значения фаз ϕг, ϕоп могут быть выбраны любыми, так как они не оказывают влияния на результат вычислений.
Результаты моделирования в виде функций P(α), F(α) приведены на фиг. 5.
Вычисление функций F(α) и P(α) по теоретическим формулам (1) и (10) соответственно, дало их полное совпадение с представленными на фиг. 5 результатами моделирования.
Из результатов моделирования, представленных на фиг. 5 видно, что ширина главного лепестка пеленгационной функции P(α) по нулевым значениям составляет ~0.1 рад, в то время, как соответствующая ширина главного лепестка функции F(α) составляет ~0.2 рад, то есть разрешение сигнала с использованием функции P(α) в два раза выше, чем с использованием функции F(α). Таким образом, для рассмотренного случая требуемый технический результат использования предлагаемого способа подтвержден.
Дополнительно тестирование приведенной на фиг. 4 схемы устройства пространственной обработки проведено для случая прихода на раскрыв антенны сигналов от двух источников, находящихся на направлениях, заданных, соответственно, углами αc1=-0.2 рад и αc2=-0.3 рад. Угловое расстояние между источниками выбрано таким образом, что оно оказывается меньше потенциальных разрешающих возможностей диаграммы направленности F(α).
Напряжения на выходах i-х элементов антенной решетки для этого случая описываются формулой
Figure 00000029
.
Значения фаз в опорной точке для сигналов произвольно приняты равными ϕоп1=0.6π и ϕоп2=2π соответственно. Значение фазы ϕоп1 оказывает влияние на результат, так как оно не компенсируется при фильтрации напряжения в дополнительном канале, настроенном на прием сигнала с направления αc2, однако это влияние незначительно.
Для рассматриваемого случая принято, что направление на первый источник известно, в связи с чем с целью формирования виртуального раскрыва с помощью фазовращателей 2 в дополнительные каналы устройства обработки вводят общий фазовый сдвиг -Ф(αc1)=-kXsin(αc1). В связи с этим напряжения дополнительных каналов после ввода данного фазового сдвига имеют вид:
Figure 00000030
.
В остальном параметры схемы и порядок расчетов остаются прежними.
Результаты моделирования пространственной обработки суммарного сигнала от двух источников в виде функции пеленгации P(α) и функции диаграммы направленности F(α), полученных на выходе детектора 8, приведены на фиг. 6. Здесь же показаны функции диаграмм направленности F1(α), F2(α) для каждого источника сигнала в отдельности. Функция пеленгации P(α) формируется при подключенных выключателями 3 дополнительных каналах устройства обработки, а функция диаграммы направленности F(α) - при их отключении. Функции F1(α), F2(α) рассчитаны по формуле (1).
Расчеты подтверждают повышенные разрешающие возможности предлагаемого способа обработки сигнала - для заданных условий он обеспечил полное разрешение источника сигнала с направления αc1 по отношению к источнику сигнала с направления αc2, что отражает функция пеленгации P(α), в то время как известный способ пространственной обработки, основанный на построении функции диаграммы направленности F(α), не позволяет разрешить сигналы.
Пример 2. Рассматривается устройство пространственной обработки сигнала с использованием плоской прямоугольной антенной решетки с размерами раскрыва в двух плоскостях X×Y и с количеством элементов в этих плоскостях N×K соответственно. С целью упрощения представления схемы выбран фильтровой вариант введения фазовых сдвигов. Структурная схема устройства пространственной обработки сигнала приведена на фиг. 7, где используются те же обозначения элементов, что и для фиг. 4, за исключением следующих обозначений: 1 - антенная решетка из N×K приемных элементов, имеющих нумерацию i=0, 1, …, N-1, j=0, 1, …, N-1; дополнительные каналы условно обозначены «X», «Y», «XY».
Принцип функционирования устройства соответствует описанному в примере 1, со следующими отличиями: направление прихода сигнала задается двумя углами αс и βс; в качестве опорного выбран приемный элемент раскрыва решетки с координатами (0,0).
Напряжения на выходах элементов антенной решетки описываются формулой
Figure 00000031
.
С помощью управляемых фазовращателей 2 в напряжения ui,j(t,ϕоп) вводят компенсирующие фазовые сдвиги
Figure 00000032
,
Figure 00000033
. С учетом этого напряжение на выходе основного канала, присоединенного к приемному элементу антенной решетки с номером i,j, описывается формулой
Figure 00000034
.
Напряжения на выходах дополнительных каналов «X», «Y», «XY», присоединенных к приемному элементу решетки с номером i,j, с учетом введенных фазовых сдвигов -Ф(αс)=-kXsin(αc), -Ф(βс)=-kYsin(βc), Ф(α)=kXsin(α), Ф(β)=kYsin(β) описываются формулами:
Figure 00000035
,
Figure 00000036
,
Figure 00000037
.
Выходные напряжения основных и дополнительных каналов всех N×K приемных элементов решетки суммируются в сумматоре 4:
Figure 00000038
.
Сигнал с выхода сумматора, пройдя через смеситель 5, где он умножается на гетеродинное напряжение (6), интегрируется на интервале повторения сигнала T в интеграторе 7:
Figure 00000039
.
Выходное напряжение интегратора U(α,β,Δϕ) подвергается квадратурному детектированию:
Figure 00000040
.
При отключенных выключателями 3 дополнительных каналах реализуется известный способ пространственной обработки принятого сигнала и на выходе детектора 8 формируется напряжение, которое описывается функцией диаграммы направленности F(α,β).
Для оценки работоспособности предложенной схемы устройства пространственной обработки выполнено численное моделирование в соответствии с описанным порядком действий с сигналом, снимаемым с приемных элементов антенной решетки, для следующих исходных данных: αс=-0.2 рад, βс=-0.2 рад, N=20, K=30, X=0.3 м, Y=0.4 м, λ=0.03 м. Неизвестные значения фаз ϕг, ϕоп, как и в примере 1, выбирались произвольно и не влияли на результат вычислений.
Результаты моделирования в виде функций P(α,β), F(α,β) приведены на фиг. 8 и фиг. 9 соответственно. Как следует из анализа функций P(α,β), F(α,β), главный лепесток функции P(α,β) в различных угловых сечениях в два раза уже главного лепестка функции F(α,β).
Проведенные расчеты для двумерного случая также подтверждают достижение технического результата использования предлагаемого способа обработки сигнала.
Таким образом, использование предложенного способа пространственной обработки обеспечивает повышение углового разрешения принимаемых радиосигналов.
В связи с тем, что двумерный случай является более общим по отношению к одномерному, формула изобретения дается применительно к двумерному случаю.

Claims (1)

  1. Способ пространственной обработки сигнала, заключающийся в приеме с помощью фазированной антенной решетки с плоским прямоугольным раскрывом, имеющим размеры X×Y, узкополосного сигнала с плоским волновым фронтом, приходящего с направления (αсс), определяемого углами пеленгования α, β, отсчитываемыми в прямоугольной системе координат с осями, ориентированными вдоль сторон раскрыва, от нормали к раскрыву, восстановленной из опорной точки, находящейся в одном из углов раскрыва, его согласованной пространственной фильтрации и квадратурного детектирования суммы напряжений всех приемных элементов после их фильтрации, отличающийся тем, что для обработки напряжения, снимаемого с каждого приемного элемента решетки, используют по четыре канала с выполнением в них одинаковых для данного приемного элемента операций согласованной пространственной фильтрации, причем в трех каналах в обрабатываемые в них напряжения дополнительно вводят фазовые сдвиги, равные: kX(sin(α)-sin(αc)) - для одного канала, kY(sin(β)-sin(βc)) - для другого канала, kX(sin(α)-sin(αc))+kY(sin(β)-sin(βc)) - для третьего канала, где k=2π/λ, λ - длина волны принимаемого сигнала.
RU2015151075/07A 2015-11-30 2015-11-30 Способ пространственной обработки радиосигналов RU2599257C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2015151075/07A RU2599257C1 (ru) 2015-11-30 2015-11-30 Способ пространственной обработки радиосигналов

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2015151075/07A RU2599257C1 (ru) 2015-11-30 2015-11-30 Способ пространственной обработки радиосигналов

Publications (1)

Publication Number Publication Date
RU2599257C1 true RU2599257C1 (ru) 2016-10-10

Family

ID=57127643

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2015151075/07A RU2599257C1 (ru) 2015-11-30 2015-11-30 Способ пространственной обработки радиосигналов

Country Status (1)

Country Link
RU (1) RU2599257C1 (ru)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108008346A (zh) * 2017-10-10 2018-05-08 上海交通大学 一种基于两单元时间调制阵列的雷达系统
RU2682165C1 (ru) * 2018-02-14 2019-03-15 Акционерное общество "Центральное конструкторское бюро автоматики" Фазовый пеленгатор

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4626859A (en) * 1983-10-07 1986-12-02 Racal Research Limited Direction finding systems
US6606058B1 (en) * 1999-03-26 2003-08-12 Nokia Networks Oy Beamforming method and device
RU2226704C2 (ru) * 2001-03-11 2004-04-10 Открытое акционерное общество "Научно-исследовательский институт измерительных приборов" Способ пространственной обработки сигнала (варианты)
WO2006087783A1 (ja) * 2005-02-17 2006-08-24 Fujitsu Limited 到来方向推定装置
SU1840570A1 (ru) * 1983-10-24 2007-08-27 Нижегородский научно-исследовательский институт радиотехники Многолучевая адаптивнная антенная решетка
RU139809U1 (ru) * 2013-10-31 2014-04-27 Федеральное государственное автономное образовательное учреждение высшего профессионального образования "Санкт-Петербургский государственный университет аэрокосмического приборостроения", ГУАП Устройство пространственной обработки сигналов

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4626859A (en) * 1983-10-07 1986-12-02 Racal Research Limited Direction finding systems
SU1840570A1 (ru) * 1983-10-24 2007-08-27 Нижегородский научно-исследовательский институт радиотехники Многолучевая адаптивнная антенная решетка
US6606058B1 (en) * 1999-03-26 2003-08-12 Nokia Networks Oy Beamforming method and device
RU2226704C2 (ru) * 2001-03-11 2004-04-10 Открытое акционерное общество "Научно-исследовательский институт измерительных приборов" Способ пространственной обработки сигнала (варианты)
WO2006087783A1 (ja) * 2005-02-17 2006-08-24 Fujitsu Limited 到来方向推定装置
RU139809U1 (ru) * 2013-10-31 2014-04-27 Федеральное государственное автономное образовательное учреждение высшего профессионального образования "Санкт-Петербургский государственный университет аэрокосмического приборостроения", ГУАП Устройство пространственной обработки сигналов

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Пространственно-временная обработка сигналов. Под ред. КРЕМЕРА И.Я., Москва, Радио и связь, 1984, с.55-57. *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108008346A (zh) * 2017-10-10 2018-05-08 上海交通大学 一种基于两单元时间调制阵列的雷达系统
RU2682165C1 (ru) * 2018-02-14 2019-03-15 Акционерное общество "Центральное конструкторское бюро автоматики" Фазовый пеленгатор

Similar Documents

Publication Publication Date Title
US5019793A (en) Digitally implemented variable phase shifter and amplitude weighting device
RU2374724C1 (ru) Периметрическая антенная решетка радара
US3202990A (en) Intermediate frequency side-lobe canceller
US20080100510A1 (en) Method and apparatus for microwave and millimeter-wave imaging
RU2596018C1 (ru) Способ амплитудной пеленгации источников радиосигналов
RU2495447C2 (ru) Способ формирования диаграммы направленности
US3460145A (en) Electronic scanning system for wave energy beam forming and steering with receptor arrays
RU2599257C1 (ru) Способ пространственной обработки радиосигналов
Welsby et al. Multiplicative receiving arrays
US3603992A (en) Dual harmonic frequency synthetic aperture radar
US3339199A (en) Single-channel signal-processing network and monopulse receiver systems embodying the same
RU2460087C2 (ru) Радиолокационная станция с широкополосным непрерывным линейно частотно-модулированным излучением
Ryle et al. The synthesis of large radio telescopes by the use of radio interferometers
JP6245797B2 (ja) レーダ装置
US5430453A (en) Cylindrical phased array antenna system to produce wide-open coverage of a wide angular sector with high directive gain and moderate capability to resolve multiple signals
US3500303A (en) Signal generator for producing a set of signals of common frequency and adjustable phase slope
US5339284A (en) Signal processor for elimination of sidelobe responses and generation of error signals
Feng et al. Two-way pattern grating lobe control for distributed digital subarray antennas
US3618093A (en) Obstacle detection radar system
Nowok et al. 3D imaging system based on a MIMO approach at 360GHz for security screening
Shcherbyna et al. Prospect for using low-element adaptive antenna systems for radio monitoring stations
US4521893A (en) Clock distribution circuit for active aperture antenna array
RU2553270C1 (ru) Способ и устройство определения угловой ориентации летательных аппаратов
Kim et al. Digital beamforming technique with high resolution digital phase shifter and digital phase calibration using SDR
RU2514197C1 (ru) Способ и устройство определения угловой ориентации летательных аппаратов

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20171201