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

Формула изобретения

1. Способ определения координат источника радиоизлучения (ИРИ), основанный на приеме сигналов ИРИ на три антенны, отличающийся тем, что располагают три антенны в вершинах треугольника АВС, измеряют разности времен t_AC и t_BC приема сигнала ИРИ антеннами, образующими измерительные базы {А,С} и {В,С}, вычисляют значения r_АС, r_BC и r_АВ разностей дальностей от ИРИ до пар антенн, вычисляют значение азимута ИРИ с использованием выражения

=arctg(a( r_AC+ r_BC)-x₃( r_AC- r_BC))/(у₃( r_BC- r_АС)), где х₃=(b²-с ²)/а; у₃= 4b²-(а+x₃)²; а=|АВ|/2; b=|AC|/2; с=|ВС|/2, вычисляют координаты {x_F, y _F) точки F, принадлежащей линии пеленга ИРИ, с использованием выражений:

x_F=(2а r_АС-(а+х₃) r_AB)(a²+ r_AC( r_AC- r_AB)-g²) r_AB/К,

y_F=-у₃(4а ²- r² _АВ)(a²-g²+ r_АС( r_АС- r_АВ))/К,

К=4а²(у₃ ²+ r_AC( r_AB- r_AC))- r² _АВ(а²+g²)+2ах₃ (2- r² _АВ);

g= х₃ ²+у₃ ²,

вычисляют значение угла места ИРИ с использованием выражения =arctg(l/а₁), где l= (b₁ ²cos² -а₁ ²sin² ); а₁= r_АВ/2; b₁= а²- r² _AB/2, измеряют плотности потока мощности P ₁ и Р₂ сигнала ИРИ в точках А и В соответственно, вычисляют значение дальности L до ИРИ от точки J, которая является центром сферы с радиусом L, являющейся поверхностью положения ИРИ, вычисляют также координаты точки J с использованием выражений:

L=2a P₁P₂/|P₁-P₂ |, x_J=2aP₂/(P₁-P₂ )-a, y_J=0, z_J=0, вычисляют значение дальности L' до ИРИ от точки F с использованием выражения L'=( х·cos( )-у_F·sin( )·cos( )+ L²-( x·sin( )-y_F-cos( ))²-( x-cos( )-y_F·sin( ))²·sin²( ), где x=x_J-x_F, и отображают полученные результаты.

2. Устройство для определения координат источника радиоизлучения (ИРИ), содержащее три приемные антенны, два измерителя разности времен приема сигнала ИРИ, вычислительный блок и блок индикации, причем выходы первой и второй антенн подключены соответственно к первым входам первого и второго измерителей разности времен приема сигнала ИРИ, вторые входы которых соединены с выходом третьей антенны, а выходы подключены соответственно к первому и второму входам вычислительного блока, выход которого подключен к входу блока индикации, отличающееся тем, что введены два измерителя плотности потока мощности сигнала ИРИ, предназначенные для измерения плотности потока мощности сигнала ИРИ в точках размещения первой и второй антенн, входы первого и второго измерителей плотности потока мощности сигнала ИРИ соединены с выходами первой и второй антенн соответственно, а выходы подключены соответственно к третьему и четвертому входам вычислительного блока, предназначенного для вычисления азимута , угла места ИРИ, координат точки F, принадлежащей линии пеленга ИРИ, дальности L до ИРИ от точки J, координат точки J, дальности L' до ИРИ от точки F, при этом точка J является центром сферы с радиусом L, являющейся поверхностью положения ИРИ.

Описание изобретения к патенту

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

Координатами объекта называется упорядоченный набор линейных и(или) угловых числовых величин, определяющих положение объекта на поверхности или в пространстве в соответствующей системе координат [1].

Поэтому устройства для определения координат разделяются на две группы:

1) системы определения местоположения (ОМП) - устройства для определения трех координат объекта;

2) пеленгаторы - устройства для определения одной или двух угловых координат объекта.

В современных пассивных системах определения координат широко используется разностно-дальномерный (РДМ) принцип.

Так известен [2] РДМ способ определения трех линейных координат источника радиоизлучения (ИРИ), предполагающий выполнение следующего:

- размещают четыре антенны в разнесенных точках A_i, i=1, 2, 3, 4 таким образом, что объем ;

- одновременно принимают сигнал ИРИ на все антенны;

- измеряют три независимых разности времен t_j1=t_j-t₁, j=2, 3, 4 приема сигнала парами антенн, образующих измерительные антенные базы (где t_j - время приема сигнала в точке А_j );

- составляют систему нелинейных уравнений, связывающих искомые координаты x,у,z ИРИ с координатами х_i,у_i,z_i точек приема А_i, i=1, 2, 3, 4 и измеренными разностями дальностей r_j1= t_j1/с (где с - скорость распространения электромагнитных волн):

- вводят дополнительную переменную r₁ и преобразуют систему нелинейных уравнений (1) к виду системы линейных уравнений:

где

, j=2, 3, 4;

- подставляют решения системы линейных уравнений (2) вида x(r₁), у(r₁ ), z(r₁) в уравнение (3) и решают квадратное уравнение относительно неизвестного значения r₁;

- находят два набора корней системы уравнений (2) для двух значений r₁;

- выбирают из двух наборов корней один, соответствующий искомым координатам ИРИ.

Недостатками данного способа является необходимость априорной информации о местоположении объекта для выполнения последней операции (неоднозначность решения) и необходимость использования 4-х антенн для определения трех координат ИРИ.

Известен [3] РДМ способ пеленгования, который предполагает выполнение следующих операций:

- располагают три антенны в вершинах треугольника АВС;

- принимают сигнал ИРИ на все три антенны;

- измеряют разности времен t_AC и t_BC приема сигнала ИРИ антеннами, образующими измерительные базы {A,C} и {В,С};

- вычисляют значения разностей дальностей от ИРИ до пар антенн с использованием выражений:

r_АС= t_АС· _ЭМВ, r_ВС= t_ВС· _ЭMB, r_АВ= r_АС- r_ВС,

где _ЭМБ - скорость распространения электромагнитной волны,

- вычисляют значение азимута ИРИ с использованием выражения

где x₃=(b²-c ²)/a; ;

а=|АВ|/2; b=|AС|/2; с=|BС|/2;

- вычисляют координаты {х_F,у_F} точки F, принадлежащей линии пеленга ИРИ, с использованием выражений:

где ;

- отображают полученные результаты.

Недостаток данного способа состоит в том, что определяется только одна угловая координаты ИРИ - азимут (см. фиг.1), тогда как для однозначного определения координат ИРИ необходимо определить также угол места и дальность до ИРИ.

Данный способ выбран в качестве прототипа.

Целью изобретения является обеспечение возможности определения трех координат ИРИ.

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

Предлагаемый способ предполагает выполнение следующих операций:

- располагают три антенны в вершинах треугольника АВС, при этом антенны 1, 2 и 3 располагают в точках А, В и С соответственно;

- принимают сигнал ИРИ на все три антенны;

- измеряют разности времен t_AC и t_BC приема сигнала ИРИ антеннами, образующими измерительные базы {А,С} и {В,С},

- измеряют плотности потока мощности Р₁ и Р₂ сигнала в точках размещения антенн 1 и 2;

- вычисляют значения разностей дальностей от ИРИ до пар антенн с использованием выражений:

r_АС= t_АС· _ЭMB, r_ВС= t_ВС· _ЭМВ, r_АВ= r_АС- r_ВС,

где _ЭМВ - скорость распространения электромагнитной волны,

- вычисляют значение азимута ИРИ с использованием выражения

где x₃=(b²-с ²]/a;

;

а=|АВ|/2;

b=|АС|/2;

с=|ВС|/2;

- вычисляют координаты {х _F,y_F} точки F, принадлежащей линии пеленга ИРИ, с использованием выражений:

где ;

- вычисляют значение угла места ИРИ с использованием выражения

=arctg(l/a₁),

где ,

a₁= r_AB/2; ;

- вычисляют значение дальности L до ИРИ от точки J и координаты точки J с использованием выражений:

,

x_J=2aP₂/(P₁ -P₂)-a,

у_J=0, z_J =0,

- вычисляют значение дальности L' до ИРИ от точки F с использованием выражения

где x=x_J-x_F;

- отображают полученные результаты.

Устройство, реализующее данный способ, состоит из трех функционально взаимосвязанных элементов (фиг.2):

- антенной системы, содержащей три антенны 1, 2 и 3;

- системы измерения, содержащей блоки 4 и 5, предназначенные для измерения разностей времен приема сигнала ИРИ парами антенн {1;3} и {2;3}, и блоки 6 и 7, предназначенные для измерения плотности потока мощности принимаемых сигналов в точках размещения антенн 1 и 2;

- системы обработки и отображения, содержащей вычислительный блок 8 и блок 9, осуществляющий визуализацию результатов.

Принцип работы устройства, реализующего предлагаемый способ, состоит в следующем:

антенны 1, 2 и 3 располагают в трех точках трехмерного пространства А, В, С, имеющих координаты x_A,у_A,z_A , х_B,у_B,z_B и x_С,у_С,z_С соответственно.

Для удобства и наглядности дальнейшего изложения предположим, что точка расположения ИРИ совпадает с некоторой точкой D, имеющей неизвестные координаты x,у,z .

Любое радиоизлучение во временной области может быть представлено в виде потока бесконечно коротких импульсов [4]. Представим излучение ПРИ в виде последовательности импульсов, задаваемых функцией

,

;

t_i - момент времени появления i-го импульса;

U_i - амплитуда i-го импульса.

Тогда излучение ИРИ можно характеризовать мощностью P₀ и временем t₀ излучения i-го импульса.

Местоположение ИРИ относительно опорных точек (ОТ) характеризуется дальностью r. Плотность потока мощности Р_П, создаваемая i-м импульсом излучения в ОТ, и время t_П прихода i-го импульса излучения в ОТ содержат информацию о дальности r. При этом плотность потока мощности излучения в ОТ и излученная мощность связаны со значением дальности мультипликативным соотношением

P_П=P₀·k ₁{r),

а время прихода в ОТ и время излучения i-го импульса - аддитивным

t_П=t ₀+k₂(r),

где k₁(r), k₂(r) - функции дальности.

При распространении в изотропной среде для точечного ИРИ

k₁ (r)=1/(4 r²), k₂(r)=r/c,

где с - скорость распространения электромагнитного излучения.

Измерив значения Р_П и t_П в ОТ, при известных значениях P₀ и t₀, значение дальности r может быть определено с использованием формул:

,

r=c(t_П-t₀),

где Ф - коэффициент пропорциональности, учитывающий форму диаграммы направленности передающей антенны.

Методы измерения дальности, основанные на таком подходе, называются соответственно амплитудными и временными. Они широко применяются в радиолокации и радионавигации. Однако в других практических случаях, в частности при пассивном определении координат ИРИ, значения P₀ , Ф и t₀ неизвестны.

При наличии двух ОТ неизвестные значения P₀ и t₀ могут быть исключены из рассмотрения. Если ИРИ является изотропным или квазиизотропным (ИРИ излучает во все стороны с одинаковой интенсивностью в одной плоскости, что свойственно для ИРИ УКВ диапазона, использующих вертикальные штыревые антенны), значение Ф также исключается из рассмотрения.

Из системы уравнений

где Р_П1 и Р_П2 - значения плотности потока мощности излучения в двух пространственно разнесенных ОТ О₁ и О₂,

можно составить пропорцию

.

Выбрав систему координат (СК) таким образом, чтобы выполнялись условия:

, ,

где R - расстояние между ОТ (фиг.3),

можно записать

,

откуда следует, что

х ²Р_П1+у²Р_П1+z² Р_П1-x²Р_П2+2xRP_П2-R ²Р_П2-у²Р_П2-z² Р_П2=0.

Данное уравнение можно переписать в виде

x²[Р_П1-P_П2 )+у²(P_П1-Р_П2)+z²(Р _П1-Р_П2)+2xRР_П2-R²РП ₂=0.

Поскольку

,

то можно записать

.

Разделив все слагаемые на (Р _П1-Р_П2), окончательно можно получить, что

.

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

Таким образом, для случая наличия двух точек приема и при наличии информации о значениях плотности потока мощности сигнала в этих точках, априорная неопределенность пространственного положения ИРИ снижается до принадлежности поверхности положения в виде сферы с параметрами, определяемыми из приведенных выражений.

Аналогично, из системы уравнений

где t_П1 и t_П2 - значения моментов времени приема излучения в двух пространственно разнесенных ОТ O₁ и О₂,

можно составить уравнение

t_П1-t _П2=k₂(r₁)-k₂(r₂ ).

С учетом того, что

k₂ (r_i)=r_i/c, i=1, 2.

можно записать

t_П1-t_П2=(r ₁-r₂)/c.

Выбрав СК таким образом, чтобы выполнялись условия (фиг.4)

, ,

можно получить

.

Возведя в квадрат правую и левую части уравнения, получим

,

а, следовательно,

Если раскрыть скобки в левой части и произвести упрощения, то уравнение (4) можно представить в виде уравнения двуполостного гиперболоида вращения

,

где a₁= r/2; .

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

Обозначим разность расстояний от точки ИРИ до антенн 1 и 3 через r_AC, до антенн 2 и 3 - через r_BC, до антенн 1 и 2 - через r_AB, плотность потока мощности сигнала на входе антенны 1 - через Р₁, на входе антенны 2 - через P ₂.

На основе обработки двух значений разностей дальностей r_AB и r_AC вычисляют азимут , угол места ИРИ и координаты {х_F,у_F) точки F, через которую проходит линия положения ИРИ (фиг.5) с использованием выражений [3]:

где x₃=(b²-c ²)/a; ;

а=|АВ|/2; b=|АС|/2; с=|ВС|/2;

; ;

a₁= r_AB/2; .

На основе обработке двух значений плотности потока мощности Р₁ и Р₂ вычисляют значение дальности L до ИРИ от точки J и координаты точки J (фиг.5) с использованием выражений:

,

x_J=2aP₂/{P ₁-P₂)-a,

у_J=0, z _J=0.

Поскольку линия пеленга проходит через точку F, то для практического удобства целесообразно определить дальность до ИРИ не от точки J, а от точки F.

Линия пеленга ИРИ, проходящая через точку F, и заданная углами и , является линией положения ИРИ. Сфера, заданная центром в точке J и радиусом L, является поверхностью положения ИРИ. Следовательно, ИРИ расположен в точке пересечения линии пеленга с этой сферой.

Поскольку ИРИ принадлежит линии пеленга, следовательно, принадлежит и плоскости , перпендикулярной плоскости хОу и проходящей через точку F под углом к оси Ох (поскольку эта плоскость содержит линию пеленга). Линия пересечения плоскости и сферы, соответствующей дальности L, представляет собой окружность радиуса |NS| (фиг.6).

Координата точки М по оси х может быть вычислена с использованием выражения

x_М=x_F-у_Fctg( ).

Следовательно

|JN|=|JM|ctg( )=(x_J-x_М)ctg( ).

Следовательно, можно вычислить

.

Поскольку L'=|FD|, то с учетом фиг.7 можно записать

.

В итоге

,

где x=x_J-x_F.

Таким образом, на заключительном этапе вычисляют значение дальности L' до ИРИ от точки F с использованием выражения

.

и отображают полученные результаты.

В состав заявляемого устройства входят (фиг.2):

1 - антенна;

2 - антенна;

3 - антенна;

4 - измеритель разности времен;

5 - измеритель разности времен;

6 - измеритель плотности потока мощности;

7 - измеритель плотности потока мощности;

8 - вычислительный блок;

9 - блок индикации.

Выходы антенн 1 и 2 соединены с первыми входами измерителей разности времен 4 и 5 соответственно, на вторые входы которых подается сигнал с выхода антенны 3. Выходы антенн 1 и 2 также соединены со входами измерителей плотности потока мощности 6 и 7 соответственно. Выходы измерителей разности времен 4 и 5 подключены к первому и второму входам вычислительного блока 8 соответственно. Выходы измерителей плотности потока мощности 6 и 7 подключены к третьему и четвертому входам вычислительного блока 8 соответственно. Выход вычислительного блока 8 подключен ко входу блока индикации 9.

Антенны 1, 2 и 3 располагают в вершинах треугольника АВС соответственно.

Сигналы с выходов антенн 1 и 3 поступают на первый и второй входы измерителя разности времен 4 соответственно, аналогично сигналы с выходов антенн 2 и 3 поступают на первый и второй входы измерителя разности времен 5 соответственно. Измерители разности времен 4 и 5 осуществляют операцию измерения разностей времен t_AC и t_BC прихода сигнала ИРИ на пары антенн {1,3} и {2,3}.

Измерители разности времен 4 и 5 реализуют один из известных [например, 5] способов измерения разности времен.

С выходов измерителей разностей времен 4 и 5 измеренные значения t_AC и t_BC поступают на первый и второй входы вычислительного блока 8 соответственно.

Сигналы с выходов антенн 1 и 2 также поступают на входы измерителей плотности потока мощности 6 и 7 соответственно. Измерители плотности потока мощности 6 и 7 осуществляют операцию измерения плотности потока мощности P₁ и Р₂ сигнала ИРИ на входах антенн 1 и 2 соответственно.

Измерители плотности потока мощности 6 и 7 реализуют один из известных способов измерения плотности потока мощности. Из уровня техники известны измерители плотности потока мощности ПО-1, М3-10, П3-9, П3-13, П3-18 и др.

С выходов измерителей плотности потока мощности 6 и 7 измеренные значения Р₁ и Р₂ поступают на третий и четвертый входы вычислительного блока 8 соответственно.

Вычислительный блок 8 представляет собой специализированное вычислительное устройство, в котором последовательно выполняются следующие операции:

- вычисляются значения разностей дальностей r_АC, r_BC и r_AB с использованием выражений:

r_АС= t_АС· _ЭMB, r_ВС= t_ВС· _ЭМВ, r_АВ= r_АС- r_ВС.

- вычисляется значение азимута ИРИ с использованием выражения

- вычисляются значения x_F , у_F координат точки F, принадлежащей линии положения ИРИ, с использованием выражений:

где ,

- вычисляется значение угла места ИРИ с использованием выражения

,

- вычисляется значение дальности L до ИРИ от точки J и координаты точки J с использованием выражений:

,

х_J=2аР₂/(Р ₁-Р₂)-а,

- вычисляется значение дальности L' до ИРИ от точки F с использованием выражения

,

где x=x_J-x_F.

Необходимые для проведения вычислений априорно известные значения:

- _ЭМВ - скорости распространения электромагнитной волны;

- а - половины расстояния между антеннами 1 и 2;

- x₃, у₃ - координат антенны 3 в системе координат Oxyz хранятся в памяти вычислительного блока 8.

Вычисленные значения , , x_F, у_F, L' с выхода вычислительного блока 8 поступают в блок индикации 9, который предназначен для визуализации результатов предлагаемого способа.

Таким образом, предлагаемые способ и устройство для его реализации, по сравнению с прототипом, обеспечивают возможность определения трех координат ИРИ за счет определения дальности до ИРИ.

СПИСОК ЛИТЕРАТУРЫ

1. Шебшаевич B.C. Введение в теорию космической навигации. - М.: Сов. радио, 1971. - 296 с.

2. Дулевич В.Е., Коростелев А.А., Мельник Ю.А. и др. Теоретические основы радиолокации / Под ред. В.Е.Дулевича. - М.: Сов. радио, 1964. - 732 с.

3. Сайбель А.Г. Разностно-дальномерный способ пеленгования источника радиоизлучения и реализующее его устройство. Патент РФ № 2258242 от 31.05.2005.

4. Седякин Н.М. Элементы теории случайных импульсных потоков. - М.: Сов. радио, 1965. - 263 с.

5. Wuu Chenn, Pearson Allan E. On time deley estimation involving received signals // IЕЕЕ Trans. Acount, Speech and Signal Process., 1984, 32, N 4, C.828-835.