RU2069866C1 - Фазовый пеленгатор - Google Patents

Фазовый пеленгатор Download PDF

Info

Publication number
RU2069866C1
RU2069866C1 SU4738175A RU2069866C1 RU 2069866 C1 RU2069866 C1 RU 2069866C1 SU 4738175 A SU4738175 A SU 4738175A RU 2069866 C1 RU2069866 C1 RU 2069866C1
Authority
RU
Russia
Prior art keywords
phase
outputs
inputs
antenna elements
ambiguity
Prior art date
Application number
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 SU4738175 priority Critical patent/RU2069866C1/ru
Application granted granted Critical
Publication of RU2069866C1 publication Critical patent/RU2069866C1/ru

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

Относится к пассивной радиолокации, а именно к оценке углового положения источника электромагнитного излучения в двух ортогональных плоскостях.
Фазовый пеленгатор состоит из N расположенных в одной плоскости антенных элементов, N приемников, входы которых соединены с соответствующими антенными элементами, N-1 фазометров, первый и второй входы каждого из которых подключены к выходам приемников соответственно i и i+1,
Figure 00000001
, блока устранения неоднозначности, входы которого соединены с выходами соответствующих фазометров, двух весовых сумматоров, подключенных к выходам блока устранения неоднозначности, весовые коэффициенты которых определяются выражениями :
Figure 00000002

где B -1 Φ - матрица, обратная корреляционной матрице фазовых ошибок,
Figure 00000003
- векторы масштабных коэффициентов, зависящие от расположения антенных элементов в плоскости,
Figure 00000004
- квадратичная матрица. 5 ил.

Description

Устройство относится к области пассивной радиолокации, а именно, к оценке углового положения источника электромагнитного излучения в двух ортогональных плоскостях, и может быть использовано для измерения траекторий движущихся объектов в радиолокационных и радионавигационных системах.
Угловое положение источника электромагнитного излучения определяется двумя угловыми координатами α (азимут Аz) и b (угол места) или парой направляющих косинусов к осям прямоугольной системы координат, расположенной на поверхности Земли. Для совместного определения Аz и El фазовыми методами необходимо наличие двух пар антенн, расположенных на перпендикулярных прямых и образующих базы интерферометра (фиг. 1). Разность фаз v12 и Φ34 определяются следующими выражениями:
Figure 00000011
(1)
где l12, l34 расстояние, на которое разнесены соответственно антенны 1, 2, 3, 4;
β угол места, угол между векторами ОС' и ОС;
a азимут, угол между вектором ОС' и осью Х.
Решая эту систему уравнений, находим оценки величин a и β:
Figure 00000012
(2)
Ошибки, возникающие при оценке α и β, могут характеризоваться условными дисперсиями этих величин и определяются следующими выражениями:
Figure 00000013
;
Figure 00000014
(3)
при условии, что l12=l34=l
Если база между антеннами превышает λ/2, то возникает неоднозначность отсчета в секторе [-π; π] или сужается сектор однозначности. Для увеличения точности оценок углового положения необходимо увеличивать базу между антеннами. Для решения этой проблемы вводят дополнительные антенны, которые необходимы для устранения неоднозначности. Базы между дополнительными антеннами меньше, чем база между антеннами для точного определения угловых координат, причем по одной из дополнительных пар антенн должен быть получен однозначный отсчет в секторе [-π; π] по каждой координате.
Известен фазовый пеленгатор системы миниТРЭК, предназначенный для определения траектории искусственных спутников Земли. Известный фазовый пеленгатор построен по принципам, изложенным выше и содержит восемь антенн точного определения угловых координат источника сигнала и пять антенн для устранения неоднозначности по методу уточнений.
Недостатками такого пеленгатора являются следующие. Во-первых, в процессе измерений с фазометров снимается n=m+k разностей фаз, m во одному плечу, к по другому плечу, перпендикулярному первому, а в процессе оценки угловых координат α и β используются только две, измеренные между антеннами, база между которыми максимальна. Информация об угловых координатах, содержащаяся в остальных разностях фаз утрачивается. Во-вторых, на точность оценки угловых координат влияют только антенны, расстояние между которыми максимально, а не все, участвующие в измерениях. В-третьих, для улучшения точности оценок необходимо увеличивать расстояние между антеннами, но при этом уменьшается сектор однозначности для этой базы и, следовательно, уменьшается вероятность правильного устранения неоднозначности, то есть достоверность оценок.
Наиболее близким по технической сущности является известный фазовый пеленгатор, содержащий пять слабонаправленных антенных элемента, располагающихся на двух перпендикулярных плечах, образующих две пары на плоскости, где один из антенных элементов расположен в точке пересечения плеч, пять приемно-усилительных тракта и четыре фазометра, два из которых подключены между антеннами одного плеча, а два других между антеннами другого плеча, причем выходы первых двух фазометров подключены в входам первого устройства устранения неоднозначности, а выходы вторых двух фазометров к входам второго устройства устранения неоднозначности, выходы устройств устранения неоднозначности подключены к входам вычислителя координат углового положения источника сигнала.
Недостатком этого фазового пеленгатора является то, что в процессе оценки утрачивается информация об угловых координатах, содержащаяся в разности фаз между антеннами, имеющих меньшую базу, и наличие этих антенн не влияет на точность оценок, которые зависят только от антенн, имеющих большую базу.
Применение метода максимального правдоподобия позволит улучшить качественные характеристики оценок и избавиться от недостатков, изложенных выше.
Рассмотрим две произвольные точки на плоскости ХОУ (фиг. 2). Положение точек определяется их координатами, а именно: точка А с координатами (ХA, YA) и (XB, YB). Пусть источник сигнала расположен в точке С, удаление которой от точек А и В значительно больше отрезка АВ. Тогда разность фаз между точками А и В, создаваемая источником сигнала С равна:
Figure 00000015

где a′ угол между отрезком АВ и осью ОХ,
α азимут,
b угол места.
Figure 00000016
длина отрезка АВ.
Рассмотрим направляющие косинусы вектора ОС:
Figure 00000017

Figure 00000018

где С (Хo, Yo, Zo) точка, расположенная на векторе ОС.
Рассмотрим сомножитель cos(α-α′) в (4):
cos(α-α′)=cosαcosα′+sinαsinα′;
Figure 00000019

Figure 00000020

Подставив это в (4), получим:
Figure 00000021

при этом учтено, что
Figure 00000022

Определим, чему равны cosα, cosβ, sinα, исходя из фиг. 2.
Figure 00000023

Подставив это в уравнение (6) с учетом уравнений (5), получим:
Figure 00000024

Далее заменив
Figure 00000025

получим более компактный вид формулы:
FAB=2π[nxv+nyu] (7)
Пусть антенная система состоит из N антенн. При наличии погрешностей в виде фазовых ошибок, распределенных нормально с нулевыми средними значениями и известной корреляционной матрицей получим условную плотность распределения вероятностей фазовых ошибок, при фиксированных V и U:
Figure 00000026

где
Figure 00000027
вектор полных разностей фаз в единицах 2π,
K коэффициент пропорциональности,
nx, ny векторы масштабных коэффициентов,
B -1 Φ матрица, обратная корреляционной матрице фазовых ошибок.
Определим оценки направляющих косинусов U и V по максимуму их условной плотности распределения:
Figure 00000028

и получим:
Figure 00000029

где Н квадратная матрица размером (N-1)(N-1):
Figure 00000030

Figure 00000031
векторы весовых коэффициентов.
Так как получены максимально правдоподобные оценки, то они имеют нулевые смещения и минимальные значения условных дисперсий.
Условные оценки дисперсий получаются из формул для U* и V* по правилам теории вероятностей.
Figure 00000032

Figure 00000033

где М оператор усреднения.
Если учесть, что М
Figure 00000034
.
где В ковариационная матрица фазовых ошибок, то получим:
Figure 00000035
(8)
В случае, если B=σ 2 Φ •BΦ, то
Figure 00000036

Для получения оценок необходимо иметь полный вектор фазовых измерений
Figure 00000037
который представляет собой
Figure 00000038
где
Figure 00000039
вектор фазовых измерений,
Figure 00000040
- вектор целых значений разностей фаз, утраченных при измерении (вектор неоднозначности).
Для отыскании вектора
Figure 00000041
необходимо найти, при каком значении вектора из совокупности возможных значений достигается максимумов функции правдоподобия или нижняя граница следующей формы:
Figure 00000042

где R квадратная матрица размерности (N-1)(N-1), имеющая следующий вид:
Figure 00000043

Вид формы (10) получается при отыскании максимального значения функции правдоподобия и отбрасывания членов, не зависящих от
Figure 00000044
.
Таким образом, устранение неоднозначности фазовых измерений заключается в подборе такого вектора
Figure 00000045
, который давал бы нижнюю границу формы
Figure 00000046
.
На фиг. 1 изображена геометрическая интерпретация фазового пеленгатора с двумя перпендикулярными плечами; на фиг. 2 пояснение к алгоритму работы фазового пеленгатора с плоской антенной решеткой; на фиг. 3 функциональная схема фазового пеленгатора; на фиг. 4 функциональная схема устройства устранения неоднозначности; на фиг. 5 функциональная схема весового сумматора.
Фазовый пеленгатор содержит N антенных элементов 1, приемников 2, (N-1) фазометров 3, устройство устранения неоднозначности 4, два весовых сумматора 5.
Устройство устранения неоднозначности 4 (фиг.4) содержит m(N-1) первых весовых сумматоров 6, m запоминающих устройств 7, m вторых весовых сумматоров 8, третьих сумматоров 9, устройство выбора вектора неоднозначности 10, сумматор 11.
Весовой сумматор 5 (фиг. 5) является входовым весовым сумматором, который включает m умножителей 12, запоминающее устройство 13, m-входовый сумматор 14.
Фазовый пеленгатор работает следующим образом.
N приемников 2 подключены к (N-1) фазометрам 3, причем первый и второй входы каждого i-го фазометра 3 подключены к выходам соответственно первого и (i+1) приемников 2,
Figure 00000047
а входы приемников 2 подключены к соответствующим антенным элементам 1. Антенные элементы 1, приемники 2 и фазометры 3 подключаются в соответствии с принципом линейной независимости, например, одна из антенн является опорной, фазометры подключены между (N-1) остальных антенн и опорной. Выходы фазометров 3 подключены к (N-1) входу устройства устранения неоднозначности 4, к выходам которого подключены два весовых сумматора 5 с весовыми коэффициентами, определяемыми выражениями:
Figure 00000048

где B -1 Φ матрица, обратная корреляционной матрице фазовых ошибок,
Figure 00000049
векторы масштабных коэффициентов, зависящие от расположения антенных элементов в плоскости,
Figure 00000050
квадратичная матрица.
Выходы фазометров 3 соединены с входами первых весовых сумматоров 6, выходы которых соединены с входами вторых весовых сумматоров 8, входы третьих сумматоров 9 подключены к выходам вторых весовых сумматоров 8 и запоминающих устройств 7, а выходы сумматоров 9 подключены к входам устройства выбора вектора неоднозначности 10. Входы сумматора 11 соединены с (N-1) выходами устройства выбора вектора неоднозначности 10 и (N-1) выходами фазометров 3.
На плоскую антенную решетку падает плоский волновой фронт, и, так как антенные элементы 1 разнесены на плоскости, то сигналы, поступающие на входы приемников 2, имеют разные фазы. В приемниках 2 эти сигналы усиливается и подаются на входы фазометров 3, которые измеряют разность фаз между сигналами, но при этом, так как область измерения фазометров лежит в пределах [-π; π], теряется целое число периодов разностей при измерении. Измеренная разность фаз подвергается дополнительной обработке в соответствии с алгоритмом (10). Информация о разности фаз поступает на вход устройства устранения неоднозначности 4. Коды, содержащие информацию о разности фаз Φ1, ..., ΦN-1 от фазометров 3 поступают на входы первых весовых сумматоров 6, где формируется первое слагаемое в выражении (10), где вектор весов
Figure 00000051
. На вторых весовых сумматорах 8 формируется сумма
Figure 00000052
.
Затем с выхода вторых весовых сумматоров 8 поступают на третьи сумматоры 9, где к ним добавляется кодовое расстояние каждого вектора неоднозначности, которое считывается с запоминающего устройства 7:
Figure 00000053
.
Устройство выбора вектора неоднозначности 10 определяет вектор
Figure 00000054
, которому соответствует минимальное значение
Figure 00000055
. Полученные коды dj образуют форму, соответствующую выражению (10). Все m значений фазы поступают на устройство выбора вектора неоднозначности 10, или схемы сравнения, где необходимый вектор выбирается из общей совокупности из векторов по минимуму формы
Figure 00000056
. Полученный вектор неоднозначности
Figure 00000057
складывается с измеренным вектором разностей фаз
Figure 00000058
и таким образом получается полная разность принятых сигналов. Оценки V* и U* направляющих косинусов получаются на весовых сумматорах 5.
При работе весового сумматора 5 на входы умножителей 12 поступают коды чисел Ф1, Ф2.Ф полных разностей фаз с устройства устранения неоднозначности 4, которые необходимо просуммировать с весовыми коэффициентами Q1, Q2,Qm. Коды весовых коэффициентов поступают на входы умножителей 12 с запоминающего устройства 13. С выхода умножителей 12 результаты произведений поступают на вход сумматора 14, выход которого является выходом весового сумматора 5. Весовые сумматоры устройства устранения неоднозначности могут быть реализованы по этому же принципу.
Устройство выбора вектора неоднозначности может быть реализовано различными способами, один из них на запоминающих элементах. Так как известны все возможные коды разностей фаз, то можно получить все возможные коды формы
Figure 00000059
, которые являются минимальными для данного вектора. Таким образом, в запоминающем устройстве записываются коды векторов, а коды значений формы являются адресами кодов векторов. Другой заключается в том, что процесс поиска минимального значения формы сводится к показанному сравнению значений формы с помощью схем сравнения, образующих пирамидальную структуру. В сумматоре 11 формируется вектор полных разностей фаз добавлением целой части p разности и ее дробной части.
Фазовый пеленгатор может быть реализован как в аналоговом, так и в цифровом виде.

Claims (1)

  1. Фазовый пеленгатор, включающий N расположенных в одной плоскости антенных элементов, N приемников, входы которых соединены с соответствующими антенными элементами, N-1 фазометров, первый и второй входы каждого i-го фазометра подключены к выходам соответственно первого и (i+1)-го приемников,
    Figure 00000060
    блок устранения неоднозначности, входы которого соединены с выходами соответствующих фазометров, отличающийся тем, что к выходам блока устранения неоднозначности подключены два весовых сумматора с весовыми коэффициентами, определяемыми выражениями
    Figure 00000061
    Figure 00000062

    где B -1 Φ матрица, обратная корреляционной матрице фазовых ошибок;
    Figure 00000063
    векторы масштабных коэффициентов, зависящие от расположения антенных элементов в плоскости;
    Figure 00000064
    квадратичная матрица.
SU4738175 1989-09-18 1989-09-18 Фазовый пеленгатор RU2069866C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
SU4738175 RU2069866C1 (ru) 1989-09-18 1989-09-18 Фазовый пеленгатор

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
SU4738175 RU2069866C1 (ru) 1989-09-18 1989-09-18 Фазовый пеленгатор

Publications (1)

Publication Number Publication Date
RU2069866C1 true RU2069866C1 (ru) 1996-11-27

Family

ID=21470116

Family Applications (1)

Application Number Title Priority Date Filing Date
SU4738175 RU2069866C1 (ru) 1989-09-18 1989-09-18 Фазовый пеленгатор

Country Status (1)

Country Link
RU (1) RU2069866C1 (ru)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2631422C1 (ru) * 2016-06-02 2017-09-22 Акционерное общество "Особое конструкторское бюро Московского энергетического института" Корреляционно-фазовый пеленгатор

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Jacobs E., Ralson E.W., Ambignibg Resolution in Jnterferometry, JEEE Transactions on Aerospacc and Electronic Sustems, vol. AES-17, N 6, 1981. *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2631422C1 (ru) * 2016-06-02 2017-09-22 Акционерное общество "Особое конструкторское бюро Московского энергетического института" Корреляционно-фазовый пеленгатор

Similar Documents

Publication Publication Date Title
US7233285B2 (en) Correlation interferometer geolocation
US5451964A (en) Method and system for resolving double difference GPS carrier phase integer ambiguity utilizing decentralized Kalman filters
US7268728B1 (en) Moving transmitter correlation interferometer geolocation
US7187327B2 (en) Method and system for determining the position of an object
WO2005119288A9 (en) Method and system for determining the position of an object
CA2535538A1 (en) A method for using signal frequency change to differentially resolve long baseline interferometer measurements
US9213100B1 (en) Bearing-only tracking for horizontal linear arrays with rapid, accurate initiation and a robust track accuracy threshold
Park et al. Integer ambiguity resolution for GPS based attitude determination system
Luebbers et al. GTD terrain reflection model applied to ILS glide scope
Xue et al. Single-differenced models for GNSS-acoustic seafloor point positioning
RU2069866C1 (ru) Фазовый пеленгатор
US20020126043A1 (en) Method and apparatus for locating a terrestrial transmitter from a satellite
RU2052832C1 (ru) Фазовый пеленгатор
RU2787319C1 (ru) Способ фазовой пеленгации
RU2444740C1 (ru) Способ определения местоположения и мощностей источников излучения
RU2066458C1 (ru) Способ определения координат источника электромагнитного излучения
JP2570110B2 (ja) 水中音源位置推定システム
Kozlov et al. An analytic approach to the relation between GPS attitude determination accuracy and antenna configuration geometry
CN112230194A (zh) 一种基于平移阵列的解模糊方法、设备及存储介质
RU2298804C2 (ru) Подвижный пеленгатор
CN116500657B (zh) 提升接收机的定位精度的方法、装置及存储介质
RU2143123C1 (ru) Устройство определения высокоточного относительного местоположения движущегося объекта по сигналам спутниковых радионавигационных систем
RU2786495C1 (ru) Способ пеленгации излучающего объекта в фазовых многошкальных угломерных системах
Ivanov et al. Comparative Analysis of Search Methods for the Source of a Magnetic Field of a Dipole Model
RU2029963C1 (ru) Многошкальное фазовое устройство измерения m параметров