RU2787952C1 - Способ определения направления прихода радиосигнала - Google Patents
Способ определения направления прихода радиосигнала Download PDFInfo
- Publication number
- RU2787952C1 RU2787952C1 RU2021130285A RU2021130285A RU2787952C1 RU 2787952 C1 RU2787952 C1 RU 2787952C1 RU 2021130285 A RU2021130285 A RU 2021130285A RU 2021130285 A RU2021130285 A RU 2021130285A RU 2787952 C1 RU2787952 C1 RU 2787952C1
- Authority
- RU
- Russia
- Prior art keywords
- arrival
- signals
- elements
- antenna array
- vector
- Prior art date
Links
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 230000014509 gene expression Effects 0.000 claims description 9
- 230000000694 effects Effects 0.000 abstract description 3
- 239000000126 substance Substances 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 11
- 230000005670 electromagnetic radiation Effects 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 6
- 230000000875 corresponding Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 239000005433 ionosphere Substances 0.000 description 2
- 238000000034 method Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000000051 modifying Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000002194 synthesizing Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Abstract
Изобретение относится к радиотехнике и может быть использовано для определения направления прихода радиосигналов. Требуемый технический результат, заключающийся в повышении быстродействия и обеспечении равномерной во всем диапазоне пространственных углов точности определения направления прихода радиосигнала, а также в расширении арсенала технических средств, достигается в способе, основанном на регистрации сигналов элементами приемной антенной решетки, отличающийся тем, что используют сигналы N=4 ненаправленных или одинаково направленных элементов антенной решетки, фазовые центры которых лежат в вершинах правильного тетраэдра на расстоянии R от центра тетраэдра, где R - радиус описанной сферы, и формируют измеренный вектор полных фаз на N антенных элементах, трехмерные координаты которых заданы матрицей координат, определяют трехмерный волновой вектор k прихода плоской электромагнитной волны , откуда пеленг θ и угол места β определяют из соотношений: θ=arctg(kx/ky), β=arcsin(kz/|k|). 1 ил.
Description
Изобретение относится к радиотехнике и может быть использовано для определения направления прихода радиосигналов.
Известен способ [RU 2739486, C1, G01S 3/10, 24.12.20], основанный на приеме сигнала с помощью ненаправленных антенн, образующих антенную решетку, и измерении разности фаз между сигналами, принятыми в парах антенн, перпендикулярно ориентированных и попарно подключенных к входам двухканальных приемных трактов, выходы которых подключены к фазометрам, при этом измеряют разность фаз между сигналами, принятыми парой диаметрально расположенных антенн в виртуальной круговой эквидистантной решетке, образованной методом перемещения в пространстве по кругу в параллельной земной поверхности плоскости пеленгования беспилотными летательными аппаратами с установленными на них бортовыми ненаправленными антеннами бортовых радиопеленгаторов с двухканальными приемными трактами, соединенными входами каналов с ненаправленной антенной, фазовыми измерителями, блоками приема-передачи пеленгуемых радиосигналов и сигналов управления, бортовыми навигационными системами, устройствами вычисления пеленга и наземным модулем управления и индикации пеленга, причем настройка на частоту бортовых радиопеленгаторов на пеленгуемый источник радиосигнала производится наземным модулем управления и индикации пеленга от поискового командного радиоприемного устройства через блоки приема-передачи пеленгуемых радиосигналов и сигналов управления, при этом пеленгуемым радиосигналом второго из каналов каждого двухканального приемного тракта после прохождения преселектора в двухканальном приемном тракте радиосигнал перенаправляется для модуляции повышенной частоты в блок приема-передачи пеленгуемых радиосигналов и сигналов управления летательного аппарата приема, ретранслируется в блок приема-передачи пеленгуемых радиосигналов и сигналов управления другого летательного аппарата, где демодулируется и подается во второй канал, двухканального приемного тракта вместо перенаправленного радиосигнала для последующей обработки и измерения фазовыми измерителями, с исключением набега фаз из-за ретрансляции радиосигналов между беспилотными летательными аппаратами, определения результатов устройствами вычисления пеленга, при учете данных бортовых навигационных систем на каждом беспилотном летательном аппарате и передачей на наземный модуль управления и индикации для обработки результатов пеленгования, индикаций и регистрации.
Недостатком этого технического решения является относительно узкая область применения, поскольку для его реализации требуется использование беспилотных летательных аппаратов.
Кроме того, известен способ [RU 2625094, C1, G01S 3/890, 11.07.2017], заключающийся в том, что регистрируют время прихода электромагнитного излучения (ЭМИ) на однопозиционный пункт наблюдения с двумя точками регистрации инфразвука, а также время прихода инфразвука на две точки регистрации и определяют для каждой точки регистрации разность времени прихода ЭМИ и инфразвука, причем дополнительно, до прихода инфразвука на две точки регистрации, регистрируют магнитные компоненты сигнала ЭМИ двумя взаимно перпендикулярными в горизонтальной плоскости магнитными антеннами, ориентированными максимумами диаграмм направленности, соответственно, на Север-Юг и Запад-Восток, определяют по соотношению сигналов магнитных антенн азимут, и определяют приближенную дальность до источника сигналов по изменению спектра сигнала ЭМИ в зависимости от пройденного ЭМИ расстояния, для чего определяют суммарный спектр сигналов двух магнитных антенн, верхнюю и нижнюю частоты спектра и амплитуды сигналов на этих частотах, отношение верхней частоты к нижней частоте спектра и отношение амплитуды нижней частоты к амплитуде верхней частоты спектра; по полученным отношениям, нижней частоте спектра, скорости света и азимуту определяют приближенные дальность до источника электромагнитного излучения и его местоположение, по приближенному местоположению для каждой точки регистрации определяют угол прихода сигнала между направлением на источник сигнала и прямой, соединяющей точки регистрации; определяют приближенное расстояние до источника сигнала и по заданной скорости инфразвука определяют ожидаемый интервал времени прихода инфразвука для каждой точки регистрации с учетом погрешности заданной скорости инфразвука и определения приближенной дальности, и прекращают анализ сигналов до наступления ожидаемых интервалов времени прихода инфразвука, а в течение ожидаемых интервалов времени после прихода инфразвука и определения разности времени прихода ЭМИ и инфразвука для каждой точки регистрации, по углам прихода сигналов, известному расстоянию между точками регистрации и разностям времени прихода ЭМИ и инфразвука для каждой точки регистрации уточняют скорость инфразвука во время прохождения сигналов, по разностям времени прихода ЭМИ и инфразвука для каждой точки регистрации и уточненной скорости инфразвука во время прохождения сигналов уточняют значение дальности до источника сигналов, а по азимуту и уточненному значению дальности уточняют местоположение источника сигналов.
Недостатком этого технического решения также является относительно узкая область применения, поскольку оно может быть использовано для инфразвуковых длин волн.
Наиболее близким по технической сущности к предложенному является способ [RU 2379709, C1, G01S 13/95, 20.01.2010], основанный на регистрации сигналов приемными станциями, образующими протяженную решетку, с последующим формированием временных рядов полного электронного содержания и их фильтрацией в диапазоне периодов колебаний, соответствующих отклику ионосферы на воздействие источника ионосферного возмущения, причем используют протяженную приемную решетку и последовательно проверяют гипотезу о значениях направления прихода и скорости распространения плоского фронта ионосферного возмущения путем формирования диаграммы направленности приемной решетки и ее сканирования в заданном секторе обзора пространства «направление прихода - скорость распространения ионосферного возмущения» за счет синтеза выходного сигнала приемной решетки при синфазном суммировании рядов вариаций полного электронного содержания отдельных элементов решетки с временными сдвигами, рассчитанными исходя из проверяемых значений направления прихода и скорости распространения ионосферного возмущения и расстояний, пройденных фронтом ионосферного возмущения между элементами приемной решетки в проверяемом направлении внутри сферического слоя ионосферы Земли, решение о правильности проверяемой гипотезы и обнаружении ионосферного возмущения принимается при превышении суммарным сигналом заданного порогового уровня, соответствующие значения направления прихода и фазовой скорости распространения ионосферного возмущения считаются оценочными значениями.
Недостатком способа является относительно узкая область применения, поскольку способ предназначен, преимущественно, для определения направления прихода ионосферного возмущения. Кроме того, способ обладает относительно низким быстродействием, обусловленным необходимостью сканирования диаграммы направленности приемной решетки в заданном секторе обзора пространства. Все это сужает арсенал технических средств, которые могут быть использованы для определения направления прихода радиосигнала.
Задача, которая решается в изобретении, направлена на создание способа определения направления радиосигнала, обладающего повышенным быстродействием и обеспечивающим равномерную во всем диапазоне пространственных углов точность углометрии, и расширении на этой основе арсенала технических средств, которые могут быть использованы для определения направления прихода радиосигнала.
Требуемый технический результат заключается в повышении быстродействия и обеспечения равномерной во всем диапазоне пространственных углов точности определения направления прихода радиосигнала и расширении арсенала технических средств, которые могут быть использованы для решения такой задачи.
Поставленная задача решается, а требуемый технический результат достигается тем, что в способе, основанном на анализе данных регистрируемых приемной антенной решеткой, с последующим формированием временных рядов полного электронного содержания, согласно изобретению, используют сигналы N=4 ненаправленных или одинаково направленных элементов антенной решетки, фазовые центры которых лежат в вершинах правильного тетраэдра на расстоянии R от центра тетраэдра, где R - радиус описанной сферы, и формируют измеренный вектор полных фаз
на элементах антенной решетки, состоящей из N антенных элементов, трехмерные координаты которых заданы матрицей координат
в системе координат с началом в геометрическом центре антенной решетки
и определяют трехмерный волновой вектор к прихода плоской электромагнитной волны
где λ - длина волны в метрах, λ=300/F, F - частота сигнала в МГц, который связан с углами прихода сигнала пеленгом θ и углом места β выражением
при условии, что ось X системы координат направлена на восток, ось Y - на север, ось Z - вертикально вверх, откуда пеленг θ и угол места β определяют из соотношений выражениями:
θ=arctg(kx/ky), β=arcsin(kz/|k|).
На чертеже представлена схема размещения N=4 ненаправленных или одинаково направленных элементов антенной решетки, фазовые центры которых лежат в вершинах правильного тетраэдра на расстоянии R от центра тетраэдра, где R - радиус описанной сферы.
Способ определения направления прихода радиосигнала реализуется следующим образом.
Предварительно проведем теоретическое обоснование способа.
Математическая постановка задачи оценки направления прихода сигнала сводится к составлению уравнения измерения и его решению относительно искомых углов прихода сигнала, т.е. по измеренному вектору полных фаз ϕ' на элементах антенной решетки (АР), состоящей из N антенных элементов (АЭ), трехмерные координаты которых заданы матрицей координат
в системе координат с началом в геометрическом центре АР
необходимо оценить трехмерный волновой вектор k прихода плоской электромагнитной волны
где λ - длина волны в метрах, λ=300/F, F - частота сигнала в МГц, который связан с углами прихода сигнала пеленгом θ и углом места β выражением
при условии, что ось X системы координат направлена на восток, ось Υ - на север, ось Z - вертикально вверх.
Пеленг θ и угол места связаны с волновым вектором k=(kx,ky,kz) выражениями:
Можно показать, что оптимальной линейной оценкой фазы ϕ0 в точке, совпадающей с началом выбранной системы координат, является фаза, усредненная по элементам АР
Тогда для вектора полных фаз относительно фазы сигнала в центре АР
и искомого волнового вектора к плоской электромагнитной волны, падающей на решетку А, уравнение измерения может быть записано в следующем виде
где ε - вектор ошибок измерения фаз, относительно которого полагаем выполненными условия
где Е(⋅) - математическое ожидание; (⋅)T - знак транспонирования; I - единичная матрица размерности N×N; - среднеквадратическая ошибка фазовых измерений.
Для решения сформулированной задачи применим метод наименьших квадратов (МНК)
где Ф(k) - функционал квадратичной невязки фазовых измерений
Ф(k)=(Ak-ϕ)T(Ak-ϕ).
В качестве обоснования применения МНК укажем на тот факт, что к выражению (4) сводится получение оценки волнового вектора по принципу максимума правдоподобия при дополнительном предположении о нормальности закона распределения ошибок фазовых измерений.
Для трехмерной АР задача является существенно нелинейной. Путем сведения ее к алгебраическому уравнению с одним неизвестным оценим степень этой нелинейности и предложим способ ее преодоления.
Функционал невязки Ф(k) может быть представлен в следующем виде
где В=АТА - симметричная положительно определенная характеристическая матрица антенной решетки или матрица пространственной ориентации АР, - линейная МНК оценка волнового вектора по вектору фаз, доставляющая безусловный минимум функционалу невязки Ф(k)
Действительно, в силу положительной определенности В
Приходим к задаче условной оптимизации с нелинейным ограничением:
найти
при условии
Дадим геометрическую интерпретацию задачи (5, 6). В трехмерном пространстве волновых векторов уравнение
для C>0 задает концентрическое семейство подобных эллипсоидов с общим центром в точке . Условие (6) задает сферу с центром в начале координат. Тогда в геометрической формулировке задача (5, 6) означает следующее: из семейства (7) необходимо выбрать эллипсоид минимального размера, имеющий общую точку со сферой (6), то есть касающийся сферы. Точка касания и будет искомым решением .
Очевидно, что если случайный вектор попадает на сферу (6), то он и будет решением задачи, . В дальнейшем рассмотрении этот тривиальный случай исключаем и считаем, что .
Вводя множитель Лагранжа L, сводим задачу к безусловной оптимизации поиска минимума функционала Дифференцируя по k и приравнивая производную к 0, получим:
Расписывая покоординатно это равенство в системе координат, оси которой - собственные вектора матрицы В, получим систему уравнений для координат вектора k=(k1,k2,k3)T и параметра L:
где b1, b2, b3 - неотрицательные собственные значения матрицы В.
Подставляя выражения для kl, k2, k3 из первых трех уравнений в четвертое, получаем
или в виде алгебраического полинома относительно L
уравнение 6-й степени с одним неизвестным L.
Исходя из приведенной выше геометрической интерпретации задачи, уравнение (10) имеет либо 2, либо 4 вещественных корня, и решением задачи оптимизации (5, 6) является корень, ближайший к 0.
Можно показать, что в случае осевой симметрии антенной решетки два из трех собственных значений матрицы В равны (b2=b3), две из трех осей эллипсоида тоже равны и он становится эллипсоидом вращения, и задача сводится к уравнению 4-й степени, а в случае полной симметрии АР (b1=b2=b3) - к тривиальной задаче касания сфер, т.е. к линейному уравнению относительно L и к линейной оценке
К линейной оценке (11) приходим и в случае, когда вектор совпадает с одним из собственных векторов матрицы В.
В общем случае для решения уравнения (10) можно применить быстро сходящийся итерационный метод Ньютона:
где ƒ'(L) - производная функции ƒ(L):
После нахождения ближайшего к 0 корня уравнения (10), подставляя его в систему (9), находим координаты волнового вектора в системе координат, оси которой совпадают по направлению с собственными векторами матрицы В. Переходя в исходную систему координат, из (1) получаем оценки углов прихода сигнала - пеленга и угла места.
В предложенном устройстве используется практический случай применения координат четырех (N=4) элементов АР. Фазовые центры элементов находятся в вершинах правильного тетраэдра на расстоянии R от центра тетраэдра, где R - радиус описанной сферы. Пусть для определенности в системе координат с началом в центре тетраэдра три вершины находятся в горизонтальной плоскости, первая - на оси Y (в направлении на север), две других - по часовой стрелке, если смотреть сверху, четвертая - на вертикальной оси Z в положительном (вверх) направлении (см. фиг.).
Тогда матрица координат АР имеет вид
Матрица вторых моментов координат имеет вид
где I- единичная матрица размерности 3×3.
МНК-оценка волнового вектора по измеренному на элементах АР вектору фаз имеет вид
Пеленг и угол места связаны с волновым вектором k=(kx,ky,kz) и вычисляются на основе выражений (2). Из этих соотношений видно, что пространственные углы прихода волны не зависят от модуля волнового вектора, и так как и отличаются лишь скалярным множителем и оба коллинеарны вектору q=ATφ, то для исключения избыточных действий при оценке углов прихода в качестве вектора к в выражениях (2) достаточно взять k=q.
где σо - среднеквадратическая ошибка фазовых измерений.
Откуда относительная ошибка оценки волнового вектора, она же угловая ошибка оценки направления прихода в радианах
При этом ошибка одинакова и не зависит от направления прихода радиосигнала, что является замечательным свойством симметрии АР.
Таким образом, предложенный способ оценки направления прихода сигнала по измеренному распределению фаз на элементах трехмерной АР не требует процедур сканирования и двумерных итераций, что упрощает способ и исключает возникновением погрешностей этими процедурами. Определение направления прихода ЭМ волны производится по четырехканальным фазовым измерениям, что упрощает реализацию соответствующих устройств и повышает их быстродействие. По предварительным оценкам, основанным на подсчете количества выполняемых операций, эффект от использования предложенного способа позволяет повысить оперативность определения направления прихода в 5-7 раз. Исследования показали, что при среднеквадратической ошибке (СКО) фазовых измерений σо=5°, СКО определения направления прихода излучения σ<2° при удалении фазовых центров элементов антенной решетки от ее центра на расстояние R=01λmax и σ≤1° при R=0,3λmax, где λmax - длина волны радиосигнала, соответствующая нижней границе диапазона рабочих частот.
Claims (13)
- Способ определения направления прихода радиосигнала, основанный на регистрации сигналов элементами приемной антенной решетки, отличающийся тем, что используют сигналы N=4 ненаправленных или одинаково направленных элементов антенной решетки, фазовые центры которых лежат в вершинах правильного тетраэдра на расстоянии R от центра тетраэдра, где R - радиус описанной сферы, и формируют измеренный вектор полных фаз
- на N антенных элементах, трехмерные координаты которых заданы матрицей координат
- в системе координат с началом в геометрическом центре антенной решетки
- и определяют трехмерный волновой вектор k прихода плоской электромагнитной волны
- где λ - длина волны в метрах, λ=300/F, F - частота сигнала в МГц,
- который связан с углами прихода сигнала пеленгом θ и углом места β выражением
- при условии, что ось X системы координат направлена на восток, ось Y - на север, ось Z - вертикально вверх, откуда пеленг θ и угол места β определяют из соотношений:
- θ=arctg(kx/ky), β=arcsin(kz/|k|).
Publications (1)
Publication Number | Publication Date |
---|---|
RU2787952C1 true RU2787952C1 (ru) | 2023-01-13 |
Family
ID=
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5030959A (en) * | 1986-06-03 | 1991-07-09 | Southwest Research Institute | Enhanced automatic bearing indicator system |
RU2379709C1 (ru) * | 2008-06-25 | 2010-01-20 | Государственное образовательное учреждение высшего профессионального образования "Иркутское высшее военное авиационное инженерное училище (военный институт)" Министерства обороны Российской Федерации | Способ определения скорости распространения и направления прихода ионосферного возмущения |
RU2739486C1 (ru) * | 2020-04-09 | 2020-12-24 | Владимир Иванович Милкин | Способ пеленгования источника радиосигнала |
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5030959A (en) * | 1986-06-03 | 1991-07-09 | Southwest Research Institute | Enhanced automatic bearing indicator system |
RU2379709C1 (ru) * | 2008-06-25 | 2010-01-20 | Государственное образовательное учреждение высшего профессионального образования "Иркутское высшее военное авиационное инженерное училище (военный институт)" Министерства обороны Российской Федерации | Способ определения скорости распространения и направления прихода ионосферного возмущения |
RU2739486C1 (ru) * | 2020-04-09 | 2020-12-24 | Владимир Иванович Милкин | Способ пеленгования источника радиосигнала |
Non-Patent Citations (1)
Title |
---|
Ю.Б. Нечаев и др., Алгоритмы диаграммообразования адаптивных антенных решеток в условиях многолучевого распространения радиоволн, Научные ведомости, выпуск 21/1, 2012, номер 1 (120), СТР. 193-202. Гайдученко А.Б., АЛГОРИТМ ОПРЕДЕЛЕНИЯ УГЛА ПРИХОДА СИГНАЛА С ПОМОЩЬЮ НЕЭКВИДИСТАНТНОЙ АНТЕННОЙ РЕШЕТКИ, ЖУРНАЛ РАДИОЭЛЕКТРОНИКИ, номер 1, 2002 г. * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7187327B2 (en) | Method and system for determining the position of an object | |
RU2258242C2 (ru) | Разностно-дальномерный способ пеленгования источника радиоизлучения и реализующее его устройство | |
WO2005119288A9 (en) | Method and system for determining the position of an object | |
US10324160B2 (en) | Geolocation of beyond LOS HF emitters | |
CN104569625B (zh) | 一种基于可转动辅助天线的大型天线方向图测量方法 | |
CN110764059B (zh) | 一种收发垂直波束三坐标相控阵雷达方法 | |
CN111381215A (zh) | 相位校正方法以及流星位置获取方法 | |
CN103487798A (zh) | 一种相控阵雷达测高方法 | |
RU2560094C2 (ru) | Способ определения скорости распространения и направления прихода ионосферного возмущения | |
US4387376A (en) | Phase linear interferometer system and method | |
Grabbe et al. | Geo-location using direction finding angles | |
Jain et al. | Efficient time domain HF geolocation using multiple distributed receivers | |
US6583755B2 (en) | Method and apparatus for locating a terrestrial transmitter from a satellite | |
Kim et al. | Precise phase calibration of a controlled reception pattern GPS antenna for JPALS | |
RU2787952C1 (ru) | Способ определения направления прихода радиосигнала | |
Tsai et al. | HF Radio Angle-of-Arrival Measurements and Ionosonde Positioning. | |
RU2527943C1 (ru) | Способ обработки сигналов при многоканальной фазовой пеленгации источников радиоизлучений коротковолнового диапазона | |
Sadeghi et al. | Positioning of geostationary satellite by radio interferometry | |
CN110208741A (zh) | 一种基于多圆阵测相的超视距单目标直接定位方法 | |
RU2758979C1 (ru) | Способ автоматического измерения параметров диаграммы направленности антенны в дальней зоне методом облета с помощью бла | |
CN109856597B (zh) | 一种新体制超视距短波定位***及定位方法 | |
RU2583450C1 (ru) | Способ определения местоположения наземного источника радиоизлучения системы спутниковой связи | |
RU2325666C2 (ru) | Разностно-дальномерный способ пеленгования источника радиоизлучения | |
RU2686113C1 (ru) | Способ амплитудного двухмерного пеленгования | |
Aleshechkin | Algorithm of GNSS-based attitude determination |