СПОСОБ ОПРЕДЕЛЕНИЯ КООРДИНАТ ИСТОЧНИКА РАДИОИЗЛУЧЕНИЙ С БОРТА ЛЕТАТЕЛЬНОГО АППАРАТА С ИСПОЛЬЗОВАНИЕМ ТРИОРТОГОНАЛЬНОЙ ВИБРАТОРНОЙ АНТЕННОЙ СИСТЕМЫ Российский патент 2024 года по МПК G01S1/08 G01S5/02 

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

Известен способ определения местоположения пользовательского терминала с использованием двух спутников-ретрансляторов (CP) [1]. Указанный способ заключается в том, что на основе измерений временных задержек и частотных сдвигов между переданными и принятыми тестовыми сигналами системы, с учетом известных координат первого и второго спутников-ретрансляторов CP₁ СР₂, векторов их скоростей координат узла межсетевого сопряжения (УМС) и предварительно заданных частот конвертирования рабочих частот CP₁, СР₂, вычисляют широту ϕ_ПТ и долготу λ_ПТ пользовательского терминала (ПТ).

Для реализации указанного известного способа выполняют следующие этапы:

определяют расстояние между первым CP₁ и ПТ;

определяют расстояние между вторым СР₂ и ПТ;

измеряют модуль азимут α_пт вектора скорости пользовательского терминала и его высоту h_ПТ относительно земной поверхности;

вычисляют в УМС доплеровские сдвиги частот первого и второго узкополосных тестовых сигналов, обусловленные радиальными скоростями ПТ относительно CP₁ и СР₂, для чего предварительно определяютвероятные местоположения ПТ с учетом известных координат CP₁, СР₂ и определенных параметров определяют, по меньшей мере, один из параметров: радиальную скорость перемещения первого CP₁ относительно ПТ, и/или радиальную скорость перемещения второго СР₂ относительно ПТ, с учетом доплеровских сдвигов частот

вычисляют широту ϕ_ПТ и долготу λ_ПТ ПТ.

При таком способе обеспечивается высокая точность определения координат пользовательского терминала в широком диапазоне скоростей его перемещения за счет определения частотных сдвигов сигналов системы.

Недостатками аналога являются: длительное время определения координат ПТ, связанное с необходимостью проведения дополнительных измерений модуля, азимута вектора скорости ПТ и его высоты относительно земной поверхности; необходимость ответной передачи тестовых сигналов с ПТ, чьи координаты необходимо определить в УМС через СР.

Известен способ определения координат ИРИ с использованием летательного аппарата [2]. Указанный способ заключается в том, что:

выбирают 3-мерную декартову систему координат (ДСК) для производства измерений и расчетов;

размещают измеритель на ЛА;

перемещают указанный Л А в пространстве;

принимают радиосигнал ИРИ в заданной полосе частот ΔF;

измеряют и запоминают первичные координатно-информативные параметры принятого радиосигнала, в качестве которых используют амплитуды напряженности электрического поля (АНЭП), создаваемые ИРИ в точках приема;

измеряют и запоминают вторичные параметры (ВП) - координаты местоположения ЛА;

многократно измеряют и запоминают совокупности АНЭП E_n (где n=1…N) и ВП - координаты ЛА в 3-х мерной ДСК x_n y_n и z_n в процессе перемещения ЛА;

вычисляют N-1 коэффициентов окружностей Апполония, как отношений АНЭП, перемноженных на обратные отношения высот ЛА в соответствующих точках измерения;

формируют N-1 сферических поверхностей положения (СПП) ИРИ путем построения и последующего вращения окружностей Апполония вокруг осей, соединяющих соответствующие фокусы;

в качестве координат ИРИ в пространстве принимают координаты точки пересечения N - 1 указанных СПП ИРИ.

В указанном способе компенсируется ошибка определения координат ИРИ на плоскости на основе использования СПП ИРИ, формируемых вращением окружностей Апполония вокруг осей, соединяющих соответствующие фокусы.

Недостатком аналога являются длительное время определения координат ИРИ, связанное с необходимостью измерения N≥5 раз совокупности АНЭП и ВП в процессе перемещения ЛА.

Из известных способов наиболее близким аналогом (прототипом) предлагаемого способа по своей технической сущности является способ определения координат источника радиоизлучения с борта летательного аппарата с использованием триортогональной вибраторной антенной системы (ТОВАС) [3] заключающийся в том, что: измеряют с помощью ТОВАС ортогональные компоненты Е_х1, Е_у1, E_z1 и Е_х2, Е_у2, Ez_{2 г} векторов напряженности электрического поляс интервалом равным четверти периода электромагнитной волны измеряемого радиосигнала в моменты времени t₁ и t2, определяют ориентации векторовв пространстве, осуществляют построение вспомогательных плоскостей Ω_E1 и Ω_E2, осуществляют построение линии положения ИРИ l, выбирают поверхность Земли как поверхность положенияИРИ, вычисляют координаты ИРИ в точке пересечения линии положения ИРИ l с поверхностью Земли.

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

Техническим результатом изобретения является повышение точности определения координат ИРИ на основе использования перемещаемого в пространстве Л А, оснащенного ТОВАС.

Поставленная цель достигается тем, что в известном способе определения координат ИРИ с борта ЛА с использованием ТОВАС, состоящей из трех АЭ в виде несимметричных вибраторов штыревого типа, заключающемся в том, что выбирают топоцентрическую декартову систему координат (ДСК) O₃X₃Y₃Z₃ и ДСК O₁X₁Y₁Z₁, связанную с ЛА, для производства измерений и расчетов, размещают измерительную аппаратуру, оснащенную ТОВАС, на ЛА, измеряют координаты центра ТОВАС O₁ и углы ориентации ТОВАС в моменты времени t₁ и t₂, измеряют в моменты времени t₁ и t₂ с помощью ТОВАС ортогональные компоненты Е_х1, Е_у1, E_z1 и Е_х2, E_y2, E_z2 векторов напряженности электрического поля принятого аналогового радиосигнала, определяют ориентацию векторов напряженности электрического поляв ДСК O₁X₁Y₁Z₁ путем векторного сложения ортогональных компонент Е_х1, Е_у1, E_z1 и Е_х2, E_y2, E_z2 соответственно, строят вспомогательные плоскости Ω_EI и Ω_E2, так, что бы они были перпендикулярны векторам и соответственно и проходили через центр ТОВАС О₁, строят линию положения ИРИ l₁, соответствующую пересечению плоскостей Ω_E1 и Ω_E2, вычисляют координаты ИРИ с использованием линии положения ИРИ l₁, дополнительно перемещают в пространстве ЛА.

Измеряют координаты центра ТОВАС О₂. Выбирают ДСК O₂X₂Y₂Z₂ для производства измерений и расчетов. Измеряют углы ориентации ТОВАС в моменты времени t₃ и t₄.

Измеряют в моменты времени t₃ и t₄ с помощью ТОВАС ортогональные компоненты Е_х3, Е_у3, E_z3 и Е_х4, E_y4, E_z4 векторов напряженности электрического поля принятого аналогового радиосигнала. Определяют ориентацию векторов напряженности электрического поляв ДСК O₂X₂Y₂Z₂ путем векторного сложения ортогональных компонент Е_х3, Е_у3, Е₂з и Е_х4, E_y4, E_z4 соответственно.

Строят вспомогательные плоскости Ω_E3 иΩ_E4, так, что бы они были перпендикулярны векторам соответственно и проходили через центр ТОВАС О₂. Строят линию положения ИРИ l₂, соответствующую пересечению плоскостей Ω_E3 иΩ_E4. Вычисляют координаты ИРИ находящегося в точке пересечения линий положения ИРИ l₁ и l₂.

Интервалы между моментами времени t₁ и t₂ и между моментами времени t₃ и t₄ выбирают равным четверти периода электромагнитной волны измеряемого радиосигнала.

Благодаря указанной новой совокупности существенных признаков, в том числе за счет размещении ТОВАС на перемещаемом в пространстве ЛА, выбора интервалов между моментами времени t₁ и t₂ и между моментами времени t₃ и t₄ равным четверти периода электромагнитной волны измеряемого радиосигнала, измерении с помощью ТОВАС ортогональных компонент Е_х1, Е_у1, E_z1 Е_х2, Е_у2, E_z2, Е_х3, Е_у3, E_z3 и Е_х4, Е_у4, E_z4 векторов напряженности электрического поля в моменты времени t₁, t₂, t₃ и t₄, построения линий положения ИРИ l₁ и l₂, соответствующих пересечению плоскостей Ω_E1, и Ω_E2, Ω_E3, Ω_E4, вычисления координат ИРИ, находящегося в точке пересечения линий положения ИРИ l₁ и l₂, достигается технический результат изобретения: повышение точности определения координат ИРИ на основе использования перемещаемого в пространстве ЛА, оснащенного ТОВАС.

Заявленное изобретение поясняется чертежами, на которых показаны: на фиг. 1 положение вспомогательной плоскости Ω_E1 включающей вектор и перпендикулярной вектору в момент времени t₁;

на фиг. 2 положение вспомогательной плоскости Ω_E2, включающей вектор и перпендикулярной вектору в момент времени t₂;

на фиг. 3 построение линии положения ИРИ l₁ как линии пересечения вспомогательных плоскостей Ω_E1 и Ω_E2;

на фиг. 4 конфигурация ТОВАС состоящая из трех АЭ в виде несимметричных вибраторов штыревого типа в ДСК O₁X₁Y₁Z₁;

на фиг. 5 ортогональные компоненты Е_х1, Е_у2, E_z3 вектора напряженности электрического поля в момент времени t₁ в ДСК O₁X₁Y₁Z₁;

на фиг. 6 ортогональные компоненты Е_х2, E_y2, E_z2 вектора напряженности электрического поля t₂ в момент времени в ДСК O₁X₁ Y₁Z₁;

на фиг. 7 временные диаграммы ортогональных компонент Е_х, Е_у, E_z, принятых на АЭ 3, 4 и 5 ТОВАС, а также их значения Е_х1, Е_у1, Е₂₁ и Е_х2, E_y2, E_z2, измеренные в моменты времени t₁ и t₂ соответственно;

на фиг. 8 положение вспомогательной плоскости Ω_E3, включающей вектор и перпендикулярной вектору в момент времени t₃;

на фиг. 9 положение вспомогательной плоскости Ω_E4, включающей вектор и перпендикулярной вектору в момент времени t₄;

на фиг. 10 построение линии положения ИРИ l₂, как линии пересечения вспомогательных плоскостей Ω_E3 и Ω_E4;

на фиг. 11 графическое представление определения координат ИРИ как пересечения линий положения ИРИ l₁ и l₂ в ДСК O₃X₃Y₃Z₃.

Определение координат ИРИ является важным составным элементом мониторинга сигналов. Достоинством системы ОМП ИРИ является скрытность при определении координат вследствие отсутствия активного излучения. Размещение технических средств системы ОМП на ЛА и в том числе на беспилотных ЛА позволяет существенно расширить зону мониторинга с возможностью осуществлять обнаружение и определение координат ИРИ в труднодоступных районах.

Использование ЛА в качестве платформы для развертывания средств радиомониторинга приводит к возникновению ряда проблем, основными из которых являются:

увеличение уровня помех и связанное с ним снижение отношения сигнал/шум на входе бортового радиоприемного устройства;

ограничение массогабаритных показателей полезной нагрузки на ЛА, которые не позволяют разместить на нем эффективные антенные системы и многоканальные радиоприемные устройства;

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

Большинство способов определения координат ИРИ основано на пеленговании радиосигналов несколькими измерителями, либо одним, перемещающимся в пространстве измерителем. Повышение точности определения координат ИРИ в большинстве случаев достигается увеличением базы пеленгования. В частности, в способе прототипе в ходе определения координат ИРИ предусмотрено перемещение ЛА на расстояние, соизмеримое с расстоянием от ЛА до ИРИ.

В заявленном способе координаты ИРИ определяют с помощью сосредоточенной ТОВАС, способной определить поляризацию радиосигнала в точке приема.

Поляризация электромагнитной волны является ее пространственно-временной характеристикой и определяется видом траектории, описываемой концом вектора напряженности магнитного или электрического поля в фиксированной точке пространства.

На фиг. 1 и фиг. 2 отображены ИРИ 1 и измеритель 2 параметров электромагнитной волны в моменты времени t₁ и t₂.

Вектор Пойтинга совпадает с направлением распространения электромагнитной волны и является результатом векторного произведения векторовнапряженности электрического и магнитного полей, то есть образует вместе с ними правую тройку векторов.

На фиг. 1 и фиг. 2 отображены положения вектора Пойтинга , а также векторов напряженности электрического и магнитного полей в моменты времени t₁ и t₂ соответственно. Причем временной интервал между моментами времени t₁ и t₂ выбирают равным четверти периода электромагнитной волны измеряемого радиосигнала. Перемещением в пространстве ЛА за указанный временной интервал можно пренебречь.

На фиг. 1 и фиг. 2 показана часть фазового фронта волны Ω_EH, определяемого как поверхность одинаковых фаз векторов поля перпендикулярная направлению распространения электромагнитной волны (т.е. перпендикулярная вектору Пойтинга

Вспомогательные плоскости Ω_E1 и Ω_E2, ортогональны векторам напряженности электрического поля в моменты времени t₁ и t₂ соответственно и проходят через центр ТОВАС O_1.

Пересечение плоскостей Ω_E1 и Ω_E2, соответствует линии положения ИРИ l₁ (см. фиг. 3). Поскольку временной интервал между моментами времени t₁ и t₂ выбирают равным четверти периода электромагнитной волны измеряемого радиосигнала, угол между плоскостями Ω_E1 и Ω_E2, равен π/2.

Для измерения ортогональных компонент векторов напряженности электрического поля в заявленном способе используют ТОВАС, состоящую из трех АЭ 3, 4 и 5 в виде несимметричных вибраторов штыревого типа (см. фиг. 4).

В заявленном способе для моментов времени t₁ и t₂ используют ДСК O₁X₁Y₁Z₁, в которой центр координат O₁ совмещен с центром ТОВАС, оси O₁X₁, O₁Y₁ и O₁Z₁ направлены вдоль АЭ 3, 4 и 5 (см. фиг. 4).

Векторная сумма ортогональных компонент напряженностей электрического поля Е_х1, Е_у1, E_z1, измеренных на АЭ 3, 4 и 5 ТОВАС соответственнов момент времени t₁ составляет вектор напряженности электрического поля (см. фиг. 5).

Аналогично, векторная сумма ортогональных компонент напряженностей электрического поля Е_х2, Е_у2, E_z2, измеренных на АЭ 3, 4 и 5 ТОВАС в момент времени t₂ составляет вектор напряженности электрического поля (см. фиг. 6).

На фиг. 7 представлены временные диаграммы ортогональных компонент Е_х, Е_у, E_z, принятого в общем случае эллиптически поляризованного аналогового радиосигнала на АЭ 3, 4 и 5 ТОВАС соответственно. В моменты времени t₁ и t₂ измеряют значения компонент Е_х1, Е_у1, E_z1 и Е_х2, Е_у2, E_z2 с помощью АЭ 3, 4 и 5 ТОВАС соответственно.

На фиг. 8 и фиг. 9 отображены ИРИ 1 и измеритель 6 параметров электромагнитной волны после перемещения его в новую позицию O₂ в моменты времени t₃ и t_4.

На фиг.8 и фиг.9 отображены положения вектора Пойтинга а также векторов напряженности электрического и магнитного полей в моменты времени t₃ и t₄ соответственно. Причем временной интервал между моментами времени t₃ и t₄ выбирают равным четверти периода электромагнитной волны измеряемого радиосигнала. Перемещением в пространстве ЛА за указанный временной интервал можно пренебречь.

Вспомогательные плоскости Ω_E3 и Ω_E4, ортогональны векторам напряженности электрического поля в моменты времени t₃ и t₄ соответственно и проходят через центр ТОВАС 0₂.

Пересечение плоскостей Ω_E3 и Ω_E4, соответствует линии положения ИРИ l₂ (см. фиг. 10). Поскольку временной интервал между моментами времени t₃ и t₄ выбирают равным четверти периода электромагнитной волны измеряемого радиосигнала, угол между плоскостями Ω_E1 и Ω_E2, равен π/2.

Для измерения ортогональных компонент векторов напряженности электрического поля в заявленном способе используют ТОВАС, состоящуюиз трех АЭ 3, 4 и 5 в виде несимметричных вибраторов штыревого типа.

В заявленном способе для моментов времени t₃ и t₄ используют ДСК O₂X₂Y₂Z₂, в которой центр координат О₂ совмещен с центром ТОВАС, оси О₂Х2, O₂Y₂ и O₂Z₂ направлены вдоль АЭ 3, 4 и 5 (см. фиг. 11).

Векторная сумма ортогональных компонент напряженностей электрического поля Е₃, Е_у3, E_z3, измеренных на АЭ 3, 4 и 5 ТОВАС соответственно в момент времени t₃ составляет вектор напряженности электрического поля .

Аналогично, векторная сумма ортогональных компонент напряженностей электрического поля Е_х4, E_y4, E_z4, измеренных на АЭ 3, 4 и 5 ТОВАС в момент времени t₄ составляет вектор напряженности электрического поля .

Строят вспомогательные плоскости Ω_E1 и Ω_E2, ортогональные векторам напряженности электрического поля соответственно и проходящие через центр ТОВАС О₂ в моменты времени t₁ и t₂, который совмещен в свою очередь с началом ДСК O₁X₁Y₁Z₁ (см. фиг. 1 и фиг. 2). Строят линию положения ИРИ l₁ на пересечении плоскостей Ω_E1,Ω_E2 (см. фиг. 3).

Строят вспомогательные плоскости Ω_E3 и Ω_E4, ортогональные векторам напряженности электрического поля соответственно и проходящие через центр ТОВАС О₂ в моменты времени t₃ и t₄ который совмещен в свою очередь с началом ДСК O₂X₂Y₂Z₂ (см. фиг. 8 и фиг. 9). Строят линию положения ИРИ l₂ на пересечении плоскостей Ω_E3, Ω_E4 (см. фиг. 10).

Вычисляют координаты ИРИ, находящегося в точке точки пересечения линий положения ИРИ l₁ и l₂ (см. фиг. 10).

Вариант алгоритма вычисления координат ИРИ, находящегося в точке пересечения линий положения ИРИ l1 и l₂ по заявляемому способу представлен в приложении А.

Реализация заявленного способа преимущественно целесообразна при размещении ТОВАС на подвижном объекте, в частности на ЛА. При этомнеобходимо с высокой точностью определять координаты ЛА и углы ориентации ЛА.

Имитационное моделирование заявленного способа определения местоположения ИРИ в специальном программном обеспечении Matlab [5] показало повышение точности измерения координат ИРИ по сравнению со способом прототипом на 10…20% (в зависимости от топологии - взаимного расположения ИРИ и ЛА), что указывает на возможность достижения указанного технического результата.

ИСТОЧНИКИ ИНФОРМАЦИИ:

1. Волков Р.В., Саяпин В.Н., Севидов В.В. Способ определения местоположения пользовательского терминала с использованием двух спутников-ретрансляторов. Патент РФ на изобретение №2605457 от 18.09.2015. Опубл. 20.12.2016. Бюл. 35.

2. Агиевич С.Н., Дворников СВ., Земсков Д.С., Севидов В.В., Федоренко И.В. Способ определения координат источника радиоизлучения с использованием летательного аппарата. Патент РФ на изобретение №2644580 от 31.05.17. Опубл. 13.02.2018. Бюл. 5.

3. Григорьев В.В., Куликов М.В., Севидов В.В., Симонов А.Н., Фокин Г.А. Способ определения координат источника радиоизлучения с борта летательного аппарата с использованием триортогональной вибраторной антенной системы. Патент №2741072. МПК G01S 5/04 (2006.01). Бюл. №3 от 22.01.21. Заявка №2020125594 от 23.07.2020.

4. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. - М.: Наука, 1973. - 832 с.

5. Севидов В.В. Программа оценки точности угломерной системы координатометрии, реализующей итерационный способ Ньютона-Рафсона. Свидетельство о регистрации программы для ЭВМ №2022669254 от 19.10.2022. Заявка №2022668612 от 12.10.2022.

Приложение А

Алгоритм вычисления координат ИРИ, находящегося в точке пересечения линий положения ИРИ l₁ и l₂.

Векторы напряженности электрического поля из ДСК O₁X₁Y₁Z₁ преобразуют в за счет последовательного перемножения каждого из них на три матрицы поворота соответствующие измеренным углам поворота ДСК O₁X₁Y₁Z₁ относительно топоцентрической ДСК O₃X₃Y₃Z₃, называемыми углами Эйлера [4]:

где

Тогда, вспомогательные плоскости Ω_E1 и Ω_E2 в топоцентрической ДСК O₃X₃Y₃Z₃ будут описываться уравнениями:

где х₁, у₁, z₁ - координаты точки O₁, являющейся началом ДСК O₁X₁Y₁Z₁ в топоцентрической ДСК O₃X₃Y₃Z₃.

Векторы напряженности электрического поля из ДСК O₂X₂Y₂Z₂ преобразуют в за счет последовательного перемножения каждого из них на три матрицы поворота соответствующие измеренным углам поворота ДСК O₂X₂Y₂Z2 относительно топоцентрической ДСК O₃X₃Y₃Z3, называемыми углами Эйлера [4]:

где

Тогда, вспомогательные плоскости Ω_E3 и Ω_E4 в топоцентрической ДСК O₃X₃Y₃Z₃ будут описываться уравнениями:

где х₂, у₂, z₂ - координаты точки О₂, являющейся началом ДСК O₂X₂Y₂Z₂

Строят линию положения ИРИ l₁ на пересечении плоскостей Ω_E1 и Ω_E2. Алгебраически такое построение соответствует решению системы уравнений (1) и (3):

Аналогично строят линию положения ИРИ l₂ на пересечении плоскостей Ω_E3 и Ω_E4 Алгебраически такое построение соответствует решению системы уравнении (3) и (4):

Решают систему, составленную из уравнений (1), (2), (3) и (4):

Решением данной системы (5) будут координаты точки пересечения линий положения ИРИ l₁ и l₂ в топоцентрической ДСК O₃X₃Y₃Z₃. Полученные координаты принимают в качестве координат ИРИ.

Изобретение относится к радиотехнике и предназначено для использования в системах радиомониторинга при решении задачи скрытого определения координат источника радиоизлучения (ИРИ) с борта летательного аппарата (ЛА). Техническим результатом изобретения является повышение точности определения координат ИРИ на основе использования перемещаемого в пространстве ЛА, оснащенного триортогональной вибраторной антенной системой (ТОВАС). Технический результат достигается посредством размещения ТОВАС на перемещаемом в пространстве ЛА, измерения с помощью ТОВАС ортогональных компонент Е_х1, Е_у1, E_z1, Ех₂, E_y2, E_z2, Е_х3, Е_у3, E_z3 и Е_х4, E_y4, E_z4 векторов напряженности электрического поля в моменты времени t₁, t₂, t₃ и t₄, определения ориентации векторов в пространстве, построения вспомогательных плоскостей Ω_E1 , Ω_E2 , Ω_E3 и Ω_E4, построения линий положения ИРИ l₁ и l₂ и вычисления координат ИРИ в точке пересечения линий положения ИРИ l₁ и l₂. 1 з.п. ф-лы, 11 ил., 1 прил.

1. Способ определения координат источника радиоизлучения (ИРИ) с борта летательного аппарата (ЛА) с использованием триортогональной вибраторной антенной системы (ТОВАС), состоящей из трех антенных элементов в виде несимметричных вибраторов штыревого типа, заключающийся в том, что выбирают топоцентрическую декартову систему координат (ДСК) O₃X₃Y₃Z₃ и ДСК O₁X₁Y₁Z₁ связанную с ЛА для производства измерений и расчетов, размещают измерительную аппаратуру, оснащенную ТОВАС, на ЛА, измеряют координаты центра ТОВАС O₁ и углы ориентации ТОВАС в моменты времени t₁ и t₂, измеряют в моменты времени t₁ и t₂ с помощью ТОВАС ортогональные компоненты Е_х1, Е_у1, E_z1 и Е_х2, E_y2, E_z2 векторов напряженности электрического поля принятого аналогового радиосигнала, определяют ориентацию векторов напряженности электрического поля в ДСК O₁X₁Y₁Z₁ путем векторного сложения ортогональных компонент Е_х1, Е_у1, E_z1 и Е_х2, E_y2, E_z2 соответственно, строят вспомогательные плоскости Ω_E1 и Ω_E2, так, чтобы они были перпендикулярны векторам соответственно и проходили через центр ТОВАС O₁, строят линию положения ИРИ l₁, соответствующую пересечению плоскостей Ω_E1 и Ω_E2, вычисляют координаты ИРИ с использованием линии положения ИРИ l₁, отличающийся тем, что дополнительно перемещают в пространстве ЛА, измеряют координаты центра ТОВАС О₂, выбирают ДСК O₂X₂Y₂Z₂ для производства измерений и расчетов, измеряют углы ориентации ТОВАС в моменты времени t₃ и t₄, измеряют в моменты времени t₃ и t₄ с помощью ТОВАС ортогональные компоненты Е_х3, Е_у3, E_z3 и Е_х4, E_y4, E_z4 векторов напряженности электрического поля принятого аналогового радиосигнала, определяют ориентацию векторов напряженности электрического поля в ДCK O₂X₂Y₂Z₂ путем векторного сложения ортогональных компонент Е_х3, Е_у3, E_z3 и Е_х4, E_y4, E_z4 соответственно, строят вспомогательные плоскости Ω_E3 и Ω_E4, так, чтобы они были перпендикулярны векторам соответственно и проходили через центр ТОВАС O₂, строят линию положения ИРИ l₂, соответствующую пересечению плоскостей Ω_E3 и Ω_E4 вычисляют координаты ИРИ, находящегося в точке пересечения линий положения ИРИ l₁ и l₂.

2. Способ по п. 1, отличающийся тем, что интервалы между моментами времени t₁ и t₂ и между моментами времени t₃ и t₄ выбирают равным четверти периода электромагнитной волны измеряемого радиосигнала.