Способ определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата Российский патент 2017 года по МПК G01S1/08 G01S5/02 

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

Задача определения координат наземных источников электромагнитного излучения в радиочастотном диапазоне является важной составляющей радиомониторинга. Достоинством радиопеленгования в системе местоопределения ИРИ является скрытность при определении координат ИРИ вследствие отсутствия необходимости генерации и излучения зондирующих радиосигналов из точки радиомониторинга. Размещение радиопеленгационных средств на ЛА, и в том числе на беспилотных летательных аппаратах [1. Василии Н.Я. Беспилотные летательные аппараты. - Мн.: ООО «Попури», 2003.], позволяет существенно расширить зону радиомониторинга и осуществлять превентивное обнаружение и определение координат ИРИ, находящихся на поверхности Земли. При этом ввиду ограничений на массо-габаритные характеристики радиопеленгационных средств, размещаемых на ЛА, в качестве бортовых пеленгаторных антенн (БПА) в абсолютном большинстве практически важных случаев используются плоские антенные решетки, обеспечивающие возможность радиопеленгования с требуемым качеством в одной (азимутальной) плоскости [2. Кукес И.С., Старик М.Е. Основы радиопеленгации. - М.: Сов. радио. 1964. - 640 с.].

Известны способ и устройство определения координат наземного ИРИ [3. Патент Российской Федерации №2419106, МПК G01S 13/46, опубл. 20.06.2011], в котором на ЛА размещена бортовая пеленгаторная антенна, приемоиндикатор, обеспечивающий измерение радиопеленга с использованием бортовой пеленгаторной антенны, измеритель угловой ориентации ЛА, измеритель местоположения ЛА, преобразователь угловых координат линии радиопеленга в координаты наземного ИРИ.

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

Известен способ определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата [4. Патент Российской Федерации №2510618, МПК G01S 1/08, опубл. 10.04.2014], включающий радиоприем радиосигналов бортовой пеленгаторной антенной, частотную селекцию радиосигналов, определение линий радиопеленгов в азимутальной плоскости бортовой пеленгаторной антенны, регистрацию полученных данных периодически отсчетами, формирование по результатам регистрации полученных данных и с использованием соответствующей отсчетам априорно известной текущей информации о координатах фазового центра бортовой пеленгаторной антенны и угловой ориентации продольной, поперечной и нормальной осевых линий бортовой пеленгаторной антенны в нормальной земной прямоугольной системе координат не менее двух полуплоскостей положения наземного источника радиоизлучения, ортогональных азимутальной плоскости бортовой пеленгаторной антенны, проходящих через каждую полученную линию радиопеленга и ограниченных нормальной осевой линией бортовой пеленгаторной антенны, определение линий положения наземного источника радиоизлучения как линий пересечения каждой полуплоскости положения наземного источника радиоизлучения с поверхностью Земли, усреднение полученных данных в случае, если сформировано более двух вышеупомянутых полуплоскостей положения наземного источника радиоизлучения и текущее значение полученных данных отличается от ранее полученного среднего значения на величину, не превышающую заданное априорно известное значение, и определение координат наземного источника радиоизлучения как точки пересечения линий положения наземного источника радиоизлучения.

Недостатками известного способа определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата являются низкая точность определения координат наземных источников радиоизлучения и большие вычислительные затраты при радиопеленговании с борта летательного аппарата, которые обусловлены следующими факторами.

Во-первых, линии положения наземного источника радиоизлучения определяются как линии пересечения каждой полуплоскости положения наземного источника радиоизлучения с поверхностью Земли, что при отклонении поверхности Земли от плоскости, изрезанности рельефа земной поверхности в зоне радиопеленгования и заданной априорно известной кусочно-постоянной аппроксимации поверхности Земли приводит к погрешностям определения линий положения и высоким вычислительным затратам при их определении.

Во-вторых, определение координат наземного источника радиоизлучения осуществляется путем определения координат точки пересечения пар вышеупомянутых линий положения наземного источника радиоизлучения, что при наличии вышеупомянутых погрешностей линий положения, отклонении поверхности Земли от плоскости, изрезанности рельефа земной поверхности в зоне радиопеленгации и заданной априорно известной кусочно-постоянной аппроксимации поверхности Земли приводит к погрешностям определения координат наземного источника радиоизлучения и высоким вычислительным затратам при их определении.

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

Наиболее близким по технической сущности к заявляемому способу определения координат наземного источника радиоизлучения с борта летательного аппарата и принятым за прототип является способ определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата [5. Патент Российской Федерации №2432580, МПК G01S 1/08, опубл. 27.10.2011], включающем радиоприем радиосигналов бортовой пеленгаторной антенной, частотную селекцию радиосигналов, определение линий радиопеленгов в азимутальной плоскости бортовой пеленгаторной антенны, регистрацию полученных данных периодически отсчетами, формирование по результатам регистрации полученных данных и с использованием соответствующей отсчетам априорно известной текущей информации о координатах фазового центра бортовой пеленгаторной антенны и угловой ориентации продольной, поперечной и нормальной осевых линий бортовой пеленгаторной антенны в нормальной земной прямоугольной системе координат не менее одной независимой пары пересекающихся полуплоскостей положения наземного источника радиоизлучения, ортогональных азимутальной плоскости бортовой пеленгаторной антенны, проходящих через каждую полученную линию радиопеленга и ограниченных нормальной осевой линией бортовой пеленгаторной антенны, определение пересекающихся линий положения наземного источника радиоизлучения как линий пересечения каждой полуплоскости положения наземного источника радиоизлучения с поверхностью Земли, причем при формировании двух и более независимых пар пересекающихся полуплоскостей положения наземного источника радиоизлучения, и, соответственно, при определении двух и более независимых пар пересекающихся линий положения наземного источника радиоизлучения, выбирают пары независимых отсчетов данных, соответствующие наименьшим расстояниям от фазового центра бортовой пеленгаторной антенны до наземного источника радиоизлучения и значениям модулей разностей между парами радиопеленгов, наиболее близким к 90 градусам, весовую обработку полученных данных в случае, если сформировано две и более независимые пары пересекающихся полуплоскостей положения наземного источника радиоизлучения, которая зависит от значений радиопеленгов, соответствующих парам независимых отсчетов данных, и увеличивает веса пар независимых отсчетов данных, соответствующих наименьшим расстояниям от фазового центра бортовой пеленгаторной антенны до наземного источника радиоизлучения и значениям модулей разностей между парами радиопеленгов, наиболее близким к 90 градусам, и определение координат наземного источника радиоизлучения как точки пересечения линий положения наземного источника радиоизлучения.

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

Во-первых, линии положения наземного источника радиоизлучения определяются как линии пересечения каждой полуплоскости положения наземного источника радиоизлучения с поверхностью Земли, что при отклонении поверхности Земли от плоскости, изрезанности рельефа земной поверхности в зоне радиопеленгования и заданной априорно известной кусочно-постоянной аппроксимации поверхности Земли приводит к погрешностям определения линий положения и высоким вычислительным затратам при их определении.

Во-вторых, определение координат наземного источника радиоизлучения осуществляется путем определения координат точки пересечения пар вышеупомянутых пересекающихся линий положения наземного источника радиоизлучения, что при наличии вышеупомянутых погрешностей пересекающихся линий положения, отклонении поверхности Земли от плоскости, изрезанности рельефа земной поверхности в зоне радиопеленгации и заданной априорно известной кусочно-постоянной аппроксимации поверхности Земли приводит к погрешностям определения координат наземного источника радиоизлучения и высоким вычислительным затратам при их определении.

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

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

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

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

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

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

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

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

Решение поставленной задачи с достижением технического результата обусловлено следующим.

В предложенном способе определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата, в отличие от известного:

во-первых, линии положения наземного источника радиоизлучения определяются в аналитическом виде как линии пересечения соответствующих пар пересекающихся полуплоскостей положения наземного источника радиоизлучения, параметры которых определяются из условия ортогональности к нормалям вышеупомянутых пар пересекающихся полуплоскостей положения наземного источника радиоизлучения, и, вследствие этого, не зависят от параметров поверхности Земли, за счет чего повышается точность их определения при одновременном снижении вычислительных затрат;

во-вторых, определение координат наземного источника радиоизлучения осуществляется путем определения координат точки пересечения линии положения наземного источника радиоизлучения, параметры которой заданы в аналитическом виде, с поверхностью Земли, что при заданной априорно известной кусочно-постоянной аппроксимации поверхности Земли позволяет повысить точность определения координат наземного источника радиоизлучения при одновременном снижении вычислительных затрат;

в-третьих, при формировании двух и более независимых пар пересекающихся полуплоскостей положения наземного источника радиоизлучения выбирают пары независимых отсчетов данных, которые дополнительно соответствуют не значениям модуля разности между парой радиопеленгов, наиболее близким к 90 градусам, а дополнительно соответствуют независимым парам пересекающихся полуплоскостей положения наземного источника радиоизлучения, наиболее близким к взаимно ортогональным полуплоскостям, обеспечивающим углы пересечения линий положения наземного источника радиоизлучения с поверхностью Земли, наиболее близкие к нормали к поверхности Земли в зоне радиопеленгования, имеющим нормали, расположенные в плоскостях, наиболее близких к горизонтальной плоскости, и при которых модули углов возвышения наземного источника радиоизлучения характеризуются максимальными значениями, а модули углов места наземного источника радиоизлучения характеризуются минимальными значениями, что позволяет из общего числа несовпадающих пар пересекающихся полуплоскостей положения наземного ИРИ выбрать совокупность независимых пар пересекающихся полуплоскостей положения наземного ИРИ, для которых дисперсии круговых (радиальных) ошибок единичных оценок координат наземного ИРИ имеют минимальные значения, и за счет этого повысить точность определения координат при многократном радиопеленговании с разных точек траектории движения ЛА;

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

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

Указанные преимущества, а также особенности настоящего изобретения поясняются лучшим вариантом его осуществления со ссылками на прилагаемые фигуры.

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

Фиг. 2 - системы координат при радиопеленговании с борта летательного аппарата.

Фиг. 3 - процесс построения линии положения наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата.

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

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

Устройство определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата (фиг. 1), реализующее предлагаемый способ, содержит бортовую пеленгаторную антенну 1, включающую в своем составе Н антенных элементов 1.1, 1.2, …, Н, обеспечивающих возможность формирования в ее азимутальной плоскости Н пересекающихся диаграмм направленности (амплитудных или фазовых), с числом Н антенных элементов, не меньшем трех, позволяющим в круговом азимутальном секторе 360° однозначно измерять азимут наземного ИРИ и, соответственно, однозначно формировать линию радиопеленга наземного ИРИ в азимутальной плоскости бортовой пеленгаторной антенны 1. Устройство, кроме того, имеет радиоприемный блок 2, выполненный Н-канальным с идентичными каналами и возможностью частотной селекции и измерения параметров (амплитуд или фаз) сигналов в каждом канале, вычислитель азимута 3, выполненный с возможностью определения в пределах от 0 до 360° азимута наземного ИРИ по измеренным значениям параметров (амплитуд или фаз) сигналов в каждом канале и априорно заданным параметрам вычислений, блок 4 определения угловой ориентации бортовой пеленгаторной антенны, выполненный с возможностью определения угловой ориентации продольной, поперечной и нормальной осевых линий бортовой пеленгаторной антенны в нормальной земной системе координат в виде углов курса, тангажа и крена, блок 5 определения координат фазового центра бортовой пеленгаторной антенны, выполненный с возможностью определения абсциссы, ординаты и аппликаты фазового центра бортовой пеленгаторной антенны в нормальной земной системе координат, блок 6 определения полуплоскости положения и формирования нормали к полуплоскости положения, выполненный с возможностью определения координат единичного вектора азимутальной линии полуплоскости положения и нормали к полуплоскости положения в нормальной земной подвижной системе координат, запоминающее устройство 7 координат единичных векторов азимутальных линий полуплоскостей положения и координат нормалей к полуплоскостям положения, блок 8 вычисления параметров линии положения, выполненный с возможностью вычисления параметров линии положения как параметров линии пересечения несовпадающей пары пересекающихся полуплоскостей положения наземного ИРИ, определяемых из условия ортогональности к нормалям несовпадающей пары пересекающихся полуплоскостей положения наземного ИРИ, в виде координат вектора, находящегося на линии положения, в нормальной земной подвижной системе координат, и с возможностью вычисления координат точки начального приближения координат наземного ИРИ в нормальной земной системе координат, запоминающее устройство 9 параметров линий положения, вычислитель 10 единичной оценки координат наземного ИРИ, выполненный с возможностью вычисления координат наземного ИРИ в нормальной земной системе координат как координат точки пересечения линии положения наземного источника радиоизлучения с поверхностью Земли, которые определяются с использованием итерационной процедуры их поиска, шаг сходимости и минимальное число шагов которой выбирают зависящими от изрезанности рельефа поверхности Земли в зоне радиопеленгования, априорно известной заданной точности определения координат наземного ИРИ и угла между линией положения наземного ИРИ и местной вертикалью поверхности Земли, вычислитель 11 дисперсий единичной оценки координат наземного ИРИ, выполненный с возможностью вычислений дисперсии круговой (радиальной) ошибки единичной оценки координат наземного ИРИ и дисперсий составляющих линейных ошибок единичной оценки координат наземного ИРИ по осям абсцисс и ординат нормальной земной системы координат, запоминающее устройство 12 единичных оценок координат наземного ИРИ, запоминающее устройство 13 дисперсий единичных оценок координат наземного ИРИ, вычислитель 14 координат наземного ИРИ, выполненный с возможностью выбора линий положения, соответствующих независимым парам пересекающихся полуплоскостей положения и характеризуемых наименьшими дисперсиями круговых (радиальных) ошибок единичных оценок координат наземного ИРИ, вычисления весовых коэффициентов и вычисления координат наземного ИРИ с учетом полученных значений весовых коэффициентов, модуль 15 картографирования и индикации и модуль 16 программного управления, синхронизирующий работу радиоприемного блока 2, вычислителя азимута 3, блока 4 определения угловой ориентации бортовой пеленгаторной антенны, блока 5 определения координат фазового центра бортовой пеленгаторной антенны, блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения, запоминающего устройства 7 параметров единичных векторов азимутальных линий полуплоскостей положения и нормалей к полуплоскостям положения, блока 8 вычисления параметров линии положения, запоминающего устройства 9 параметров линий положения, вычислителя 10 единичной оценки координат наземного ИРИ, вычислителя 11 дисперсий единичной оценки координат наземного ИРИ, запоминающего устройства 12 единичных оценок координат наземного ИРИ, запоминающего устройства 13 дисперсий единичных оценок координат наземного ИРИ, вычислителя 14 координат наземного ИРИ и модуля 15 картографирования и индикации.

При этом выходы антенных элементов 1.1, 1.2, …, 1.Н, служащие выходами бортовой пеленгаторной антенны 1, подсоединены к входам соответствующих каналов радиоприемного блока 2. Выходы каналов радиоприемного блока 2 подсоединены к соответствующим входам вычислителя азимута 3. Выход вычислителя азимута 3 подсоединен к объединенным первым входам блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения и вычислителя 11 дисперсий единичной оценки координат наземного ИРИ. Первый, второй и третий выходы блока 4 определения угловой ориентации бортовой пеленгаторной антенны подсоединены соответственно к объединенным второму входу блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения и четвертому входу вычислителя 11 дисперсий единичной оценки координат наземного ИРИ, объединенным третьему входу блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения и пятому входу вычислителя 11 дисперсий единичной оценки координат наземного ИРИ и объединенным четвертому входу блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения и шестому входу вычислителя 11 дисперсий единичной оценки координат наземного ИРИ. Первый и второй выходы блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения подсоединены соответственно к первому и второму входам запоминающего устройства 7 параметров единичных векторов азимутальных линий полуплоскостей положения и нормалей к полуплоскостям положения. Первый и второй выходы запоминающего устройства 7 параметров единичных векторов азимутальных линий полуплоскостей положения и нормалей к полуплоскостям положения подсоединены соответственно к первому входу блока 8 вычисления параметров линии положения и объединенным второму входу блока 8 вычисления параметров линии положения и третьему входу вычислителя 11 дисперсий единичной оценки координат наземного ИРИ. Выход блока 5 определения координат фазового центра бортовой пеленгаторной антенны подсоединен к объединенным третьему входу блока 8 вычисления параметров линии положения и седьмому входу вычислителя 11 дисперсий единичной оценки координат наземного ИРИ. Первый и второй выходы блока 8 вычисления параметров линии положения подсоединены соответственно к первому и второму входам запоминающего устройства 9 параметров линий положения, первый и второй выходы которого подсоединены соответственно к первому и второму входам вычислителя 10 единичной оценки координат наземного ИРИ, выход которого подсоединен к объединенным второму входу вычислителя 11 дисперсий единичной оценки координат наземного ИРИ и входу запоминающего устройства 12 единичных оценок координат наземного ИРИ. Первый, второй и третий выходы вычислителя 11 дисперсий единичной оценки координат наземного ИРИ подсоединены соответственно к первому, второму и третьему входам запоминающего устройства 13 дисперсий единичных оценок координат наземного ИРИ. Выход запоминающего устройства 12 единичных оценок координат наземного ИРИ и первый, второй и третий выходы запоминающего устройства 13 дисперсий единичных оценок координат наземного ИРИ подсоединены соответственно к первому, второму, третьему и четвертому входам вычислителя 14 координат наземного ИРИ, выход которого подсоединен к входу модуля 15 картографирования и индикации. Выход модуля 16 программного управления через управляющую шину подсоединен к управляющим входам радиоприемного блока 2, вычислителя азимута 3, блока 4 определения угловой ориентации бортовой пеленгаторной антенны, блока 5 определения координат фазового центра бортовой пеленгаторной антенны, блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения, запоминающего устройства 7 параметров единичных векторов азимутальных линий полуплоскостей положения и нормалей к полуплоскостям положения, блока 8 вычисления параметров линии положения, запоминающего устройства 9 параметров линий положения, вычислителя 10 единичной оценки координат наземного ИРИ, вычислителя 11 дисперсий единичной оценки координат наземного ИРИ, запоминающего устройства 12 единичных оценок координат наземного ИРИ, запоминающего устройства 13 дисперсий единичных оценок координат наземного ИРИ, вычислителя 14 координат наземного ИРИ и модуля 15 картографирования и индикации. Выход модуля 15 картографирования и индикации служит выходной шиной устройства определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата для формирования координат наземного источника радиоизлучения в нормальной земной системе координат.

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

Координаты наземного ИРИ определяются в нормальной земной прямоугольной декартовой системе координат [6. ГОСТ 20058-80 Динамика летательных аппаратов в атмосфере. Термины, определения и обозначения]; [7. ГОСТ 22268-76 Геодезия. Термины и определения]. Отображение координат наземных ИРИ осуществляется на топографических картах местности, представляющих собой отображение поверхности Земли в проекции Гаусса-Крюгера [8. ГОСТ Р 51794-2008 Глобальные навигационные спутниковые системы. Системы координат. Методы преобразования координат определяемых точек]; [9. Бубнов И.А., Кремп А.И., Фолимонов С.И. Военная топография. - М.: Военное издательство Министерства обороны СССР, 1958. - 376 с.]; [10. Сиваков А.М. Военная топография для связистов. - М.: Воениздат, 1974. - 172 с.]. При этом для отображения на топографических картах поверхности Земли используется плоская прямоугольная геодезическая система координат, представляющая собой плоскую левую прямоугольную декартовую систему координат Σ_пл={O, X, Y}, начало О которой фиксировано по отношению к Земле и находится на поверхности земного эллипсоида, оси абсцисс ОХ и ординат OY которой, называемые соответственно вертикальной и горизонтальной линиями координатной сетки плоской прямоугольной геодезической системы координат, находятся на поверхности земного эллипсоида, координаты x и y точек которой отображают по определенному математическому закону поверхность земного эллипсоида. Следует также отметить, что положительное направление оси абсцисс ОХ плоской прямоугольной геодезической системы координат, называемое северным направлением вертикальной линии координатной сетки плоской прямоугольной геодезической системы координат, совпадает с северным направлением истинного меридиана или отличается от него на угол сближения меридианов [9], [10]. Учитывая вышеупомянутые особенности отображения координат точек и обозначения осей в плоской прямоугольной геодезической системе координат для определения координат наземного источника радиоизлучения с борта летательного аппарата целесообразно использовать нормальную земную левую прямоугольную систему координат Σ_нз={O, X, Y, Z}, представляющую собой левую пространственную прямоугольную декартовую систему координат, начало О которой фиксировано по отношению к Земле, ось абсцисс ОХ которой, находящаяся в горизонтальной плоскости, совпадает с северным направлением истинного или магнитного меридиана или вертикальной линии координатной сетки плоской прямоугольной геодезической системы координат, ось аппликат OZ которой перпендикулярна горизонтальной плоскости и направлена вверх по вертикали, а ось ординат OY которой, находящаяся в горизонтальной плоскости XOY, дополняет систему до левой пространственной прямоугольной декартовой системы координат. Отметим также, что под местной вертикалью понимают прямую, совпадающую с направлением силы тяжести в районе рассматриваемой точки Земли, и, соответственно, плоскость, проходящую через местную вертикаль, называют вертикальной плоскостью, а вышеупомянутая горизонтальная плоскость XOY перпендикулярна местной вертикали [7].

Отсчет углов при радиопеленговании ИРИ и, в частности, наземных ИРИ, осуществляется в связанной системе координат бортовой пеленгаторной антенны (БПА) Σ_св={O', X', Y', Z'}, представляющей собой левую пространственную прямоугольную декартовую систему координат, начало О' которой помещено в фазовом центре БПА, осями которой являются продольная, поперечная и нормальная осевые линии БПА, являющиеся осями абсцисс О'Х', ординат O'Y' и аппликат O'Z' соответственно. При этом фазовый центр БПА представляет собой точку БПА, относительно которой алгебраическая сумма фазовых диаграмм, формируемых БПА для осуществления радиопеленгования ИРИ, постоянна или кусочно-постоянна с наличием скачков на 180° [11. ГОСТ 23282-91 Решетки антенные. Термины и определения]. Продольная осевая линия БПА, являющаяся осью абсцисс О'Х' связанной системы координат БПА Σ_св={O', X', Y', Z'}, расположена в азимутальной плоскости БПА, проходит через фазовый центр О' БПА и ориентирована в направлении, совпадающем с положительным направлением оси симметрии совокупности пересекающихся диаграмм направленности или фазовых диаграмм, формируемых БПА для осуществления радиопеленгования ИРИ, принятая за начало отсчета углов в азимутальной плоскости БПА. Поперечная осевая линия БПА, являющаяся осью ординат O'Y' связанной системы координат БПА Σ_св={O', X', Y', Z'}, также расположена в азимутальной плоскости БПА, проходит через фазовый центр О' БПА, перпендикулярна продольной осевой линии БПА и направлена к правой части БПА. То есть, азимутальная плоскостью БПА является плоскость X'O'Y', содержащая взаимно перпендикулярные продольную и поперечную осевые линии БПА, пересекающиеся в ее фазовом центре О' и соответствующие осям соответственно абсцисс О'Х' и ординат O'Y' левой плоской прямоугольной декартовой системы координат . И, наконец, нормальная осевая линия БПА, являющаяся осью аппликат O'Z' связанной системы координат БПА Σ_св={O', X', Y', Z'}, проходит через фазовый центр О' БПА, перпендикулярна к азимутальной плоскости БПА X'O'Y' и дополняет взаимно перпендикулярные продольную О'Х' и поперечную O'Y' осевые линии БПА до левой пространственной прямоугольной декартовой системы координат Σ_св={O', X', Y', Z'}. Отметим также, что плоскость, проходящая через нормальную осевую линию O'Z' БПА называется угломестной плоскостью БПА.

Пространственное положение БПА, размещенной на ЛА, в нормальной земной левой прямоугольной системе координат Σ_нз={O, X, Y, Z}, необходимое для определения координат наземных ИРИ, характеризуется: во-первых, координатами x_ц, y_ц и z_ц фазового центра О' БПА в нормальной земной левой прямоугольной системе координат Σ_нз={O, X, Y, Z}; во-вторых, углами курса ψ, тангажа μ и крена ϑ, определяющими угловую ориентацию продольной О'Х', поперечной O'Y' и нормальной O'Z' осевых линий БПА в нормальной земной подвижной системе координат Σ_нзп={O', X, Y, Z}, представляющей собой левую пространственную прямоугольную декартовую систему координат, начало О' которой помещено в фазовом центре О' БПА, о оси абсцисс О'Х, ординат O'Y и аппликат O'Z которой направлены также как и оси соответственно абсцисс ОХ, ординат OY и аппликат OZ нормальной земной левой пространственной прямоугольной декартовой системы координат Σ_нз={O, X, Y, Z} [6]; [8]; [12. Ишлинский А.Ю. Ориентация, гироскопы и инерциальная навигация. - М.: Наука, 1976. - 622 с.]. При этом угол курса ψ БПА представляет собой угол между проекцией на горизонтальную плоскость XO'Y нормальной земной подвижной системы координат Σ_нзп={O', X, Y, Z} положительного направления продольной осевой линии О'Х' БПА и линией О'Х, расположенной в горизонтальной плоскости XO'Y, проходящей через фазовый центр О' БПА, принятой за начало отсчета, совпадающей с северным направлением истинного или магнитного меридиана или вертикальной линии ОХ координатной сетки плоской прямоугольной геодезической системы координат Σ_пл={O, Х, Y}, измеряемый в градусах в пределах от 0° до 360° и отсчитываемый от вышеупомянутой линии О'Х по ходу часовой стрелки если смотреть в направлении вниз по вертикали, проходящей через фазовый центр О' БПА. Угол тангажа μ БПА представляет собой угол между продольной осевой линией О'Х' БПА и горизонтальной плоскостью XO'Y нормальной земной подвижной системы координат Σ_нзп={O', X, Y, Z}, проходящей через фазовый центр О' БПА, измеряемый в градусах в пределах от 0° до ±90°, отсчитываемый от горизонтальной плоскости XO'Y с положительным или отрицательным знаками в случае если положительное направление продольной осевой линии О'Х' БПА находится соответственно выше или ниже горизонтальной плоскости XO'Y. Угол крена ϑ БПА представляет собой угол между поперечной осевой линией O'Y' БПА и осью ординат O'Y_ψ смещенной нормальной земной подвижной системы координат Σ_нзпψ={O', X_ψ, Y_ψ, Z}, представляющую собой нормальную земную подвижную левую прямоугольную декартовую систему координат, оси абсцисс O'Х_ψ и ординат O'Y_ψ которой смещены относительно осей соответственно абсцисс О'Х и ординат O'Y нормальной земной подвижной системы координат Σ_нзп={O', X, Y, Z} в положение, при котором угол курса ψ БПА равен нулю, измеряемый в градусах в пределах от 0° до ±90°, отсчитываемый от оси ординат O'Y смещенной нормальной земной подвижной системы координат Σ_нзпψ={O', X_ψ, Y_ψ, Z} с положительным или отрицательным знаками в случае если положительное направление оси ординат O'Y_ψ смещенной нормальной земной подвижной системы координат Σ_нзпψ={О', X_ψ, Y_ψ, Z} совмещается с положительным направлением поперечной осевой линии O'Y' БПА поворотом вокруг продольной осевой линии О'Х' БПА соответственно по ходу или против хода часовой стрелки если смотреть в положительном направлении продольной осевой линии О'Х' БПА.

Единичный вектор S направления от фазового центра О' БПА на ИРИ в нормальной земной подвижной системе координат Σ_нзп={О', X, Y, Z} характеризуется двумя углами - углом в горизонтальной плоскости и углом в вертикальной плоскости, называемыми радиопеленгом θ и углом возвышения β соответственно [2]. При этом радиопеленг θ представляет собой угол между проекцией на горизонтальную плоскость XO'Y единичного вектора S направления от фазового центра О' БПА на ИРИ и линией, расположенной в горизонтальной плоскости XO'Y, проходящей через фазовый центр О' БПА и принятой за начало отсчета, совпадающей с северным направлением истинного или магнитного меридиана или вертикальной линией координатной сетки плоской прямоугольной геодезической системы координат, измеряемый в градусах в пределах от 0° до 360° и отсчитываемый от вышеупомянутой линии по ходу часовой стрелки если смотреть в направлении вниз по вертикали, проходящей через фазовый центр О' БПА. Угол возвышения β представляет собой угол между единичным вектором S направления от фазового центра О' БПА на ИРИ и проекцией вышеупомянутого направления на горизонтальную плоскость XO'Y, измеряемый в градусах в пределах от 0° до ±90°, отсчитываемый от вышеупомянутой проекции на горизонтальную плоскость XO'Y с положительным или отрицательным знаками в случае если ИРИ находится соответственно выше или ниже горизонтальной плоскости XO'Y, проходящей через фазовый центр О' БПА.

Кроме того, единичный вектор S направления от фазового центра О' БПА на ИРИ в ранее упомянутой Связанной системе координат БПА Σ_св={О', X', Y', Z'}, начало О' которой совпадает с фазовым центром О' БПА, характеризуется также двумя углами - углом в азимутальной плоскости БПА и углом в угломестной плоскости БПА, называемыми азимутом α и углом места ε соответственно. При этом азимут α представляет собой угол между проекцией единичного вектора S направления от фазового центра О' БПА на ИРИ на азимутальную плоскость БПА X'O'Y', содержащую продольную О'Х' и поперечную O'Y' осевые линии БПА, и положительным направлением продольной осевой линии О'Х' БПА, измеряемый в градусах в пределах от 0° до 360° и отсчитываемый от положительного направления продольной осевой линии О'Х' БПА по ходу часовой стрелки если смотреть в направлении навстречу положительному направлению нормальной осевой линии O'Z' БПА. Таким образом, параметром, характеризующим линию радиопеленга в азимутальной плоскости бортовой пеленгаторной антенны является азимут α ИРИ.

Угол места ε представляет собой угол между единичным вектором S направления от фазового центра О' БПА на ИРИ и проекцией вышеупомянутого направления на азимутальную плоскость X'O'Y' БПА, в которой размещены продольная О'Х' и поперечная O'Y' осевые линии БПА, измеряемый в градусах в пределах от 0° до ±90°, отсчитываемый от вышеупомянутой проекции на азимутальную плоскость X'O'Y' БПА с положительным или отрицательным знаками в случае если ИРИ находится соответственно выше или ниже азимутальной плоскости X'O'Y' БПА, проходящей через фазовый центр О' БПА.

В процессе полета летательного аппарата, на котором размещены БПА и устройство определения координат фазового центра БПА и угловой ориентации продольной, поперечной и нормальной осевых линий БПА, координаты x_ц, y_ц и z_ц фазового центра О' БПА, единичный вектор S направления от фазового центра О' БПА на наземный ИРИ и углы курса ψ, тангажа μ и крена ϑ БПА, определяющие угловую ориентацию продольной О'Х', поперечной O'Y' и нормальной O'Z' осевых линий БПА в нормальной земной подвижной системе координат Σ_нзп={O', X, Y, Z}, непрерывно изменяются. При этом для определения координат наземного ИРИ с борта летательного аппарата регистрацию полученных данных об азимуте ИРИ, координатах фазового центра БПА и углах курса, тангажа и крена БПА осуществляют периодически отсчетами, имеющими порядковые номера n, находящимися в пределах 1≤n≤n_max, где n_max - общее число отсчетов, для которого выполняется соотношение n_max≥2:

На фиг. 2 приведены соответствующие n-му отсчету вышеупомянутых данных (где n=1, 2, 3, … - порядковый номер отсчета) связанная система координат БПА , нормальная земная подвижная система координат и смещенная нормальная земная подвижная система координат , начала которых совпадают и размещены в фазовом центре БПА при n-м отсчете данных, где, кроме того, обозначены соответствующие n-му отсчету вышеупомянутых данных: z_n=0 - горизонтальная плоскость нормальной земной подвижной системы координат ; =0 - азимутальная плоскость связанной системы координат ; S_n - единичный вектор направления от фазового центра БПА на наземный ИРИ; и - проекции единичного вектора S_n на горизонтальную плоскость и азимутальную плоскость соответственно; α_n и ε_n - соответственно азимут и угол места ИРИ; θ_n и β_n - соответственно радиопеленг и угол возвышения ИРИ; x_n, y_n и z_n - прямоугольные декартовые координаты единичного вектора S_n в нормальной земной подвижной системе координат , являющиеся косинусами углов между единичным вектором S_n и положительными направлениями осей соответственно абсцисс , ординат и аппликат нормальной земной подвижной системы координат ; и - прямоугольные декартовые координаты единичного вектора S_n в связанной системе координат , являющиеся косинусами углов между единичным вектором S_n и положительными направлениями осей соответственно абсцисс , ординат и аппликат связанной системы координат ; ψ_n, μ_n и ϑ_n - углы соответственно курса, тангажа и крена БПА, определяющие угловую ориентацию продольной , поперечной и нормальной осевых линий БПА в нормальной земной подвижной системе координат .

Взаимосвязь между прямоугольными декартовыми координатами x_n, y_n и z_n единичного вектора S_n направления от фазового центра БПА на наземный ИРИ в нормальной земной подвижной системе координат с одной стороны и радиопеленгом θ_n и углом возвышения β_n наземного ИРИ с другой стороны определяется соотношениями

Аналогичная взаимосвязь между прямоугольными декартовыми координатами и единичного вектора S_n направления от фазового центра БПА на наземный ИРИ связанной системе координат с одной стороны и азимутом α_n и углом места ε_n наземного ИРИ с другой стороны определяется соотношениями

Взаимосвязь между прямоугольными декартовыми координатами , и единичного вектора S_n направления от фазового центра БПА на наземный ИРИ в связанной системе координат с одной стороны и прямоугольными декартовыми координатами x_n, y_n и z_n вышеупомянутого единичного вектора S_n в нормальной земной подвижной системе координат с другой стороны определяется соотношениями [13. Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся втузов. - М.: Издательство «Наука», 1980. - 976 с.]

где и - направляющие косинусы оси абсцисс , равные значениям косинусов углов между положительным направлением оси абсцисс (оси ординат OY'; оси аппликат OZ') связанной системы координат и положительными направлениями соответственно оси абсцисс , оси ординат и оси аппликат нормальной земной подвижной системы координат ;

и - направляющие косинусы оси ординат , равные значениям косинусов углов между положительным направлением оси ординат (оси аппликат OZ') связанной системы координат и положительными направлениями соответственно оси абсцисс , оси ординат и оси аппликат нормальной земной подвижной системы координат ;

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

Аналогичная взаимосвязь между прямоугольными декартовыми координатами x_n, y_n и z_n единичного вектора S_n направления от фазового центра БПА на наземный ИРИ в нормальной земной подвижной системе координат с одной стороны и прямоугольными декартовыми координатами и вышеупомянутого единичного вектора S_n в связанной системе координат с другой стороны определяется соотношениями [13]

Входящие в соотношения (3) и (4) направляющие косинусы осей координат и связаны со значениями углов курса ψ_n, тангажа μ_n и крена ϑ_n БПА соотношениями [6]; [12]:

С учетом соотношений (5) взаимосвязи между радиопеленгом θ_n и углом возвышения β_n наземного ИРИ с одной стороны и азимутом α_n и углом места ε_n наземного ИРИ с другой стороны определяются соотношениями

- функция отношения действительных параметров a и b, имеющая область значений в пределах от 0° до 360°, определяемая в градусах соотношением

С использованием БПА при n-м отсчете данных осуществляется измерение азимута α_n наземного ИРИ, который определяет угловую ориентацию проекции единичного вектора S_n направления от фазового центра БПА на наземный ИРИ на азимутальную плоскость (см. фиг. 2), являющейся линией радиопеленга в азимутальной плоскости БПА при n-м отсчете данных. Так как с использованием БПА измерения угла места ε_n не осуществляется, то при априорной неопределенности угла места ε_n наземного ИРИ при n-м отсчете азимута α_n наземного ИРИ осуществляется определение в азимутальной плоскости угловой ориентации n-й полуплоскости положения наземного ИРИ α_n=const, ортогональной азимутальной плоскости БПА, проходящей через полученную n-тую линию радиопеленга в азимутальной плоскости БПА и ограниченную нормальной осевой линией БПА (см. фиг. 2). При этом вышеупомянутая полуплоскость положения наземного ИРИ α_n=const проходит через точку М расположения наземного ИРИ. В условиях априорной неопределенности угла места ε_n наземного ИРИ для определения координат точки М расположения наземного ИРИ требуется не менее - двух пересекающихся полуплоскостей положения наземного ИРИ, ортогональных азимутальной плоскости БПА, проходящих через каждую полученную n-тую линию радиопеленга и ограниченных нормальной осевой линией БПА. В аналоге [4] и прототипе [5] для минимизации ошибок определения координат наземного ИРИ, связанных со случайными флуктуациями пространственного положения азимутальной плоскости БПА во время полета ЛА, осуществляется определение не менее двух линий положения наземного ИРИ как линий пересечения каждой вышеупомянутой полуплоскости положения наземного ИРИ с поверхностью Земли, и определения координат наземного ИРИ как точки пересечения полученных пар линий положения наземного ИРИ. При минимизации ошибок получения координат наземного ИРИ, обусловленных погрешностями определения линий положения как линий пересечения каждой полуплоскости положения наземного ИРИ с поверхностью Земли, обусловленных отклонением поверхности Земли от плоскости в зоне радиопеленгации, и снижении требуемых для реализации вычислительных затрат предлагаемым способом введены операции формирования нормалей к пересекающимся полуплоскостям положения наземного ИРИ, определения не менее одной линии положения наземного ИРИ как линии пересечения пары пересекающихся полуплоскостей положения наземного ИРИ, параметры которой определяются из условия ортогональности к нормалям пары пересекающихся полуплоскостей положения наземного ИРИ, и определения координат наземного ИРИ как точки пересечения линии положения наземного ИРИ с поверхностью Земли.

Единичные векторы Q_n n-х азимутальных линий и нормали N_n к соответствующим n-м полуплоскостям положения наземного ИРИ в соответствующих n-х связанных системах координат могут быть представлены в виде матриц-столбцов соответственно Q_n(св) и N_n(св) направляющих косинусов, имеющих размер 3×1, получаемых по результатам измерений азимутов α_n и определяемых соотношениями

где n=1, 2, … - порядковые номера отсчетов данных.

Вышеупомянутые единичные векторы Q_n n-х азимутальных линий и нормали N_n к соответствующим n-м полуплоскостям положения наземного ИРИ в n-й нормальной земной подвижной системе координат с учетом (5), (9) и (10) и результатов измерений углов курса ψ_n, тангажа μ_n и крена ϑ_n БПА могут быть представлены в виде матриц-столбцов соответственно Q_n(нзп) и N_n(нзп) направляющих косинусов, имеющих размер 3×1 и определяемых соотношениями

где Q_Xn, Q_Yn и Q_Zn - прямоугольные декартовые координаты единичного вектора Q_n n-й азимутальной линии в n-й нормальной земной подвижной системе координат , являющиеся косинусами углов между единичным вектором Q_n и положительными направлениями осей соответственно абсцисс , ординат и аппликат нормальной земной подвижной системы координат ;

N_Xn, N_Yn и N_Zn - прямоугольные декартовые координаты нормали N_n к n-й полуплоскости положения наземного ИРИ в n-й нормальной земной подвижной системе координат , являющиеся косинусами углов между нормалью N_n и положительными направлениями осей соответственно абсцисс , ординат и аппликат нормальной земной подвижной системы координат ;

- квадратная матрица вращения размера 3×3 при переходе от n-й связанной системы координат к n-й нормальной земной подвижной системе координат , определяемая с учетом (4) соотношением

Радиус-вектор r_n любой точки, принадлежащей n-й полуплоскости положения наземного ИРИ, в n-й нормальной земной подвижной системе координат , можно представить в виде соотношения

где r - радиус-вектор точки, характеризуемый прямоугольными декартовыми координатами x, y и z в нормальной земной системе координат ; r_цn - радиус-вектор фазового центра БПА, характеризуемый прямоугольными декартовыми координатами х_цn, y_цn и z_цn в нормальной земной системе координат при n-м отсчете данных.

Каждая m-я пара полуплоскостей положения наземного ИРИ формируется с использованием соответствующих m-й паре несовпадающих n₁-го и n₂-го отсчетов данных, которые из общего числа отсчетов данных 1≤n≤n_max выбираются в соответствии с соотношениями

Из соотношений (15) следует, что из общего числа отсчетов данных 1≤n≤n_max может быть сформирована совокупность пар полуплоскостей положения наземного ИРИ, соответствующих несовпадающим n₁-м и n₂-м отсчетам данных, максимальное число m_max которых определяется соотношением

При этом следует отметить, что из всей совокупности m_max пар полуплоскостей положения наземного ИРИ, соответствующих несовпадающим n₁-м и n₂-м отсчетам данных, из-за погрешностей формирования n₁-й и n₂-й полуплоскостей положения наземного ИРИ, обусловленных в первую очередь ошибками измерений азимутов и наземного ИРИ, соответствующих несовпадающим n₁-м и n₂-м отсчетам данных, не все пары полуплоскостей положения наземного ИРИ являются пересекающимися. То есть, из общего числа отсчетов данных 1≤n<n_max может быть сформирована совокупность пар пересекающихся полуплоскостей положения наземного ИРИ, соответствующих несовпадающим n₁-м и n₂-м отсчетам данных, не превышающих максимальное число m_max пар полуплоскостей положения наземного ИРИ, определяемое соотношением (16).

На фиг. 3 показано построение m-й линии положения наземного ИРИ как линии пересечения пересекающихся n₁-й и n₂-й полуплоскостей m-й пары пересекающихся полуплоскостей положения наземного ИРИ и , соответствующих несовпадающим n₁-му и n₂-му отсчетам данных, где кроме того обозначено: и - связанные системы координат БПА соответственно при n₁-м и n₂-м отсчетам данных; Σ_нз={O, X, Y, Z} - нормальная земная система координат, горизонтальная плоскость XOY которой касается поверхности геоида Земли в точке О начала координат; z=z_ум - горизонтальная плоскость, аппликата которой в нормальной земной системе координат Σ_нз={О, X, Y, Z} равна среднему значению высоты z_ум над уровнем моря земной поверхности в районе предполагаемого нахождения наземного ИРИ; - азимутальная плоскость связанной системы координат ; - азимутальная плоскость связанной системы координат ; и - соответственно азимут и угол места ИРИ при n₁-м отсчете данных; и - соответственно азимут и угол места ИРИ при n₂-м отсчете данных; - единичный вектор n₁-й азимутальной линии, являющейся линией пересечения n₁-й полуплоскости положения наземного ИРИ с азимутальной плоскостью связанной системы координат ; - единичный вектор n₂-й азимутальной линии, являющейся линией пересечения n₂-й полуплоскости положения наземного ИРИ с азимутальной плоскостью связанной системы координат ; и - нормали к n₁-й и n₂-й пересекающимся полуплоскостям положения наземного ИРИ соответственно и ; γ_m - угол между нормалями и пары пересекающихся полуплоскостей положения наземного ИРИ и , соответствующих n₁-му и n₂-му отсчетам данных; K_m - вектор, находящийся на линии положения наземного ИРИ (линии пересечения m-й пары пересекающихся полуплоскостей положения наземного ИРИ и ); P_m - точка пересечения линии положения наземного ИРИ с горизонтальной плоскостью z=z_ум; z=F(x, y) - уравнение, описывающее поверхность Земли при ее кусочно-постоянной аппроксимации; M_m - точка пересечения линии положения наземного ИРИ с поверхностью Земли; G_m - нормаль к поверхности Земли при ее кусочно-постоянной аппроксимации в окрестности точки M_m пересечения линии положения наземного ИРИ с поверхностью Земли; χ_m - угол между линией положения наземного ИРИ и нормалью G_m к поверхности Земли в окрестности точки M_m; - угол между нормалями и G_m; - угол между нормалями и G_ь; и - единичные векторы направлений от фазового центра БПА на наземный ИРИ при n₁-м и n₂-м отсчетам данных соответственно; - расстояние от фазового центра БПА до точки M_m пересечения линии положения наземного ИРИ с поверхностью Земли при n₁-м отсчете данных; - расстояние от фазового центра БПА до точки M_m пересечения линии положения наземного ИРИ с поверхностью Земли при n₂-м отсчете данных; O_mZ - ось аппликат нормальной земной подвижной системы координат Σ_нзпm={O_m, X, Y, Z}, начало О_m которой находится на m-й линии положения наземного ИРИ; В - местная вертикаль поверхности Земли в точке M_m пересечения линии положения наземного ИРИ с поверхностью Земли, ориентированная в положительном направлении оси аппликат O_mZ нормальной земной подвижной системы координат Z_нзпm={O_m, X, Y, Z}; χ_вm - угол между линией положения наземного ИРИ и местной вертикалью В поверхности Земли; М_и - точка истинного положения наземного ИРИ на поверхности Земли; η_m - радиус-вектор, начало и конец которого находятся в точке М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли и точке М_и истинного положения наземного ИРИ на поверхности Земли соответственно.

Параметры m-й линии положения наземного ИРИ (линии пересечения m-й пары пересекающихся полуплоскостей положения наземного ИРИ и ) в нормальной земной системе координат Σ_нз={O, X, Y, Z} определяются прямоугольными декартовыми координатами K_Xm, K_Ym и K_Zm вектора K_m в нормальной земной подвижной системе координат Σ_нзпm={O_m, X, Y, Z}, начало О_m которой находится на m-й линии положения наземного ИРИ (см. фиг. 3). Так как вектор K_m ортогонален одновременно двум нормалям и к n₁-й и n₂-й пересекающимся полуплоскостям положения наземного ИРИ соответственно и , то он определяется в виде векторного произведения нормалей и в соответствии с соотношением

Из соотношения (17) с учетом (12) и (13) следует, что прямоугольные декартовые координаты K_Xm, K_Ym и K_Zm вектора K_m в нормальной земной подвижной системе координат Z_нзпm={О_m, Х, Y, Z} могут быть представлены в виде матрицы-столбца K_m(нзп), имеющей размер 3×1, в соответствии с соотношением

где и - прямоугольные декартовые координаты нормали к n₁-й полуплоскости положения наземного ИРИ в n₁-й нормальной земной подвижной системе координат , являющиеся косинусами углов между нормалью и положительными направлениями осей соответственно абсцисс , ординат и аппликат нормальной земной подвижной системы координат , которые одновременно являются соответствующими прямоугольными декартовыми координатами нормали в нормальной земной подвижной системе координат Σ_нзпm={O_m, X, Y, Z} и нормальной земной системе координат Σ_нз={О, X, Y, Z};

и - прямоугольные декартовые координаты нормали к n₂-й полуплоскости положения наземного ИРИ в n₂-й нормальной земной подвижной системе координат , являющиеся косинусами углов между нормалью и положительными направлениями осей соответственно абсцисс , ординат и аппликат нормальной земной подвижной системы координат , которые одновременно являются соответствующими прямоугольные декартовыми координатами нормали в нормальной земной подвижной системе координат Σ_нзпm={O_m, X, Y, Z} и нормальной земной системе координат Σ_нз={О, X, Y, Z}.

Таким образом, в предложенном способе определения координат наземного источника радиоизлучения с борта летательного аппарата, в отличие от аналога [4] и прототипа [5], линии положения наземного источника радиоизлучения определяются в аналитическом виде в соответствии с соотношениями (17) и (18) как линии пересечения соответствующих пар пересекающихся полуплоскостей положения наземного источника радиоизлучения, параметры которых не зависят от параметров поверхности Земли, за счет чего повышается точность их определения при одновременном снижении вычислительных затрат. При этом определенные вышеупомянутым образом линии положения наземного источника радиоизлучения, в отличие от аналога [4] и прототипа [5], характеризуются ненулевыми значениями угла их пересечения с горизонтальной плоскостью и, соответственно, для определения координат наземного источника радиоизлучения как точки пересечения линии положения наземного источника радиоизлучения с поверхностью Земли можно использовать всего одну линию положения наземного источника радиоизлучения.

Система соотношений, которым должны удовлетворять прямоугольные декартовые координаты и точки М_m расположения на земной поверхности наземного ИРИ, получаемые по результатам формирования m-й линии положения наземного ИРИ как линии пересечения n₁-й и n₂-й пересекающихся полуплоскостей m-й пары пересекающихся полуплоскостей положения наземного ИРИ и и результатам определения точки М_m пересечения m-й линии положения с земной поверхностью (см. фиг. 3) при отсутствии ошибок пеленгования и ошибок определения координат фазового центра БПА и угловой ориентации осевых линий БПА, определяется с учетом (14) и обозначений вышеупомянутых параметров и в виде совокупности следующих соотношений

где - радиус-вектор точки, принадлежащей n₁-й полуплоскости положения наземного ИРИ, в n₁-й нормальной земной подвижной системе координат и находящейся на m-й линии положения наземного ИРИ; - радиус-вектор точки, принадлежащей n₂-й полуплоскости положения наземного ИРИ, в n₂-й нормальной земной подвижной системе координат и находящейся на m-й линии положения наземного ИРИ; - радиус-вектор точки М_m, находящейся на m-й линии положения наземного ИРИ (см. фиг. 3), который характеризуется прямоугольными декартовыми координатами и в нормальной земной системе координат Σ_нз={О, X, Y, Z}; - радиус-вектор фазового центра БПА, характеризуемый прямоугольными декартовыми координатами и в нормальной земной системе координат Σ_нз={О, X, Y, Z} при n₁-м отсчете данных; - радиус-вектор фазового центра БПА, характеризуемый прямоугольными декартовыми координатами и в нормальной земной системе координат Σ_нз={О, X, Y, Z} при n₂-м отсчете данных; - уравнение, описывающее поверхность Земли при ее кусочно-постоянной аппроксимации в окрестности точки М_m с прямоугольными декартовыми координатами , и .

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

где и - прямоугольные декартовые координаты единичного вектора n₁-й азимутальной линии наземного ИРИ в n₁-й нормальной земной подвижной системе координат , являющиеся косинусами углов между единичным вектором и положительными направлениями осей соответственно абсцисс , ординат и аппликат нормальной земной подвижной системы координат , которые одновременно являются соответствующими прямоугольными декартовыми координатами единичного вектора в нормальной земной подвижной системе координат Σ_нзпm={O_m, X, Y, Z} и нормальной земной системе координат Σ_нз={О, X, Y, Z};

и - прямоугольные декартовые координаты единичного вектора n₂-й азимутальной линии наземного ИРИ в n₂-й нормальной земной подвижной системе координат , являющиеся косинусами углов между единичным вектором и положительными направлениями осей соответственно абсцисс , ординат и аппликат нормальной земной подвижной системы координат , которые одновременно являются соответствующими прямоугольными декартовыми координатами единичного вектора в нормальной земной подвижной системе координат Σ_нзпm={O_m, X, Y, Z} и нормальной земной системе координат Z_нз={О, X, Y, Z}.

Следует отметить, что выполнение двух неравенств, определяемых третьим и четвертым соотношениями, входящими в систему соотношений (20), является необходимым условием пересечения n₁-й и n₂-й полуплоскостей положения наземного ИРИ, и, соответственно, необходимым условием формирования m-й линии положения наземного источника радиоизлучения как линии пересечения m-й пары пересекающихся n₁-й и n₂-й полуплоскостей положения наземного ИРИ. Поэтому для определения соответствия m-й линии положения наземного источника радиоизлучения линии пересечения m-й пары пересекающихся n₁-й и n₂-й полуплоскостей положения наземного ИРИ осуществляется следующая процедура проверок.

В нормальной земной системе координат Σ_нз={O, X, Y, Z} определяются прямоугольные декартовые координаты и точки P_m пересечения линии положения наземного ИРИ с горизонтальной плоскостью z=z_ум (см. фиг. 3), являющейся точкой начального приближения координат наземного ИРИ при использовании m-й линии положения наземного ИРИ, путем определения абсциссы и ординаты точки Р_m, которые удовлетворяют при аппликате этой точки Р_m, равной априорно известной аппликате z_ум, первым двум уравнениям системы соотношений (20) и, соответственно при , являются решением системы двух вышеупомянутых уравнений, которые описывается соотношениями

Решение системы уравнений (21) относительно абсциссы и ординаты точки Р_m определяется соотношениями

где и - параметры, зависящие от априорно известной аппликаты z_ум точки Р_m (см. фиг. 3), соответствующие n₁-й и n₂-й полуплоскостям m-й пары полуплоскостей положения наземного ИРИ соответственно, определяемые соотношениями

С учетом априорно известного значения аппликаты полученные в соответствии с соотношениями (22) и (23) прямоугольные декартовые координаты и точки Р_m пересечения линии положения наземного ИРИ с горизонтальной плоскостью z=z_ум, являющейся точкой начального приближения координат наземного ИРИ при использовании m-й линии положения наземного ИРИ, проверяются на непротиворечивость неравенствам, определяемым третьим и четвертым соотношениями, входящими в систему соотношений (20), путем проверки двух следующих неравенств

Если хотя бы одно неравенство пары неравенств (26) не выполняется, найденная точка начального приближения координат наземного ИРИ при использовании m-й линии положения не находится одновременно на n₁-й и n₂-й полуплоскостях положения наземного ИРИ, что соответствует случаю, когда из-за возможных погрешностей формирования n₁-й и n₂-й полуплоскостей положения наземного ИРИ, обусловленных в первую очередь ошибками измерений азимутов и , n₁-я и n₂-я полуплоскости положения наземного ИРИ не пересекаются. В этом случае m-я линия положения наземного источника радиоизлучения, соответствующая m-й паре непересекающихся n₁-й и n₂-й полуплоскостей положения наземного ИРИ не используется при определении координат наземного ИРИ.

Если оба неравенства пары неравенств (26) выполняются, найденная точка начального приближения координат наземного ИРИ при использовании m-й линии положения находится одновременно на n₁-й и n₂-й полуплоскостях положения наземного ИРИ, что соответствует случаю, когда n₁-я и n₂-я полуплоскости положения наземного ИРИ пересекаются. В этом случае m-я линия положения наземного источника радиоизлучения, соответствующая m-й паре пересекающихся n₁-й и n₂-й полуплоскостей положения наземного ИРИ, используется при определении координат наземного ИРИ. При этом с учетом (18) и с учетом представления текущих значений прямоугольных декартовых координат x_m, y_m и z_m m-й линии положения и прямоугольных декартовых координат и точки Р_m пересечения линии положения наземного ИРИ с горизонтальной плоскостью z=z_ум в виде матриц-столбцов, имеющих размер 3×1, параметрическое уравнение m-й линии положения наземного ИРИ в нормальной земной системе координат Σ_нз={О, X, Y, Z} определяется соотношением [14. Просветов Г.И. Линейная алгебра и аналитическая геометрия: задачи и решения. Учебно-практическое пособие. 2-е изд., доп. - М.: Издательство «Альфа-Пресс», 2009. - 208 с.]

где b_m - соответствующая m-й линии положения параметрическая переменная, представляющая собой действительное число.

При этом в заявляемом способе определения координат наземного источника радиоизлучения с борта летательного аппарата прямоугольные декартовые координаты и точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли определяют как прямоугольные декартовые координаты точки пересечения m-й линии положения наземного источника радиоизлучения, описываемой параметрическим уравнением (27), с поверхностью Земли, описываемой пятым соотношением системы соотношений (20).

Таким образом, в предложенном способе определения координат наземного источника радиоизлучения с борта летательного аппарата, в отличие от аналога [4] и прототипа [5], определение координат наземного источника радиоизлучения осуществляется путем определения координат точки пересечения линии положения наземного источника радиоизлучения, параметры которой заданы в аналитическом виде, с поверхностью Земли, что при заданной априорно известной кусочно-постоянной аппроксимации поверхности Земли позволяет повысить точность определения координат наземного источника радиоизлучения при одновременном снижении вычислительных затрат.

Для определения координат наземного ИРИ с борта летательного аппарата регистрацию полученных данных об азимуте ИРИ, координатах фазового центра БПА и углах курса, тангажа и крена БПА осуществляют периодически отсчетами. Определение азимута ИРИ, координат фазового центра БПА и углы курса, тангажа и крена БПА в ходе полета ЛА производятся многократно. Для нахождения координат наземного ИРИ достаточно определения точки пересечения одной (любой) m-й линии положения наземного ИРИ с поверхностью Земли, получаемой на основе измерения соответствующей m-й паре азимутов ИРИ и соответствующих этим азимутам m-м парам координат фазового центра БПА и углов курса, тангажа и крена БПА. В случае если общее число отсчетов данных n_max больше двух, то с учетом (16) из общего числа отсчетов данных 1≤n≤n_max может быть сформирована совокупность пар пересекающихся полуплоскостей положения наземного ИРИ и, соответственно, совокупность линий положения наземного ИРИ с поверхностью Земли, позволяющих определить совокупность точек пересечения линий положения наземного ИРИ с поверхностью Земли. Если максимальное число точек пересечения линий положения наземного ИРИ с поверхностью Земли превышает единицу то, как и в прототипе [5], для повышения точности определения координат наземного ИРИ используется весовая обработка полученных абсцисс и ординат точек М_m пересечения m-х линий положения наземного ИРИ с поверхностью Земли с получением средневзвешенных оценок абсциссы и ординаты точки М расположения наземного ИРИ в соответствии с соотношениями

где w_Xm и w_Ym - весовые коэффициенты m-й единичной оценки соответственно абсцисс и ординат , которые, которые для получения минимальных погрешностей выбирают обратно пропорциональными дисперсиям соответственно и составляющих линейных ошибок определения координат наземного ИРИ по двум взаимно перпендикулярным осям соответственно абсцисс ОХ и ординат OY нормальной земной системы координат Σ_нз={О, X, Y, Z}, соответствующим единичной оценке прямоугольных декартовых координат точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли [15. Мельников Ю.П., Попов С.В. Радиотехническая радиоразведка. Методы оценки эффективности местоопределения источников излучения. - М.: Радиотехника, 2008. - 432 с.];

- максимальное число независимых пар пересекающихся полуплоскостей положения наземного ИРИ, соответствующих независимым n₁-м и n₂-м отсчетам данных, которые для уменьшения погрешности измерений выбирают из общей совокупности несовпадающих пар пересекающихся полуплоскостей положения наземного ИРИ обеспечивающими минимальные значения дисперсий круговых (радиальных) ошибок m-х единичных оценок координат и точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли, зависящих от соответствующих m-х дисперсий и составляющих линейных ошибок определения координат наземного ИРИ в соответствии с соотношением [15]

где .

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

а при нечетном общем числе n_max отсчетов данных, не меньшем трех, определяется соотношением

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

С учетом соотношения (28) и пятого соотношения системы уравнений (20), средневзвешенная оценка аппликаты точки М_и истинного расположения наземного ИРИ определяется в соответствии с соотношением

Дисперсии и составляющих линейных ошибок определения координат наземного ИРИ, соответствующие единичной оценке прямоугольных декартовых координат точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли, зависят не только от погрешностей определения m-й пары радиопеленгов, включающих погрешности измерения соответствующих m-й пары азимутов наземного ИРИ, m-х пар углов курса, тангажа и крена БПА, m-й пары расстояний от БПА до наземного ИРИ, но и от параметров взаимной пространственной ориентации БПА при m-й паре отсчетов данных и параметров пространственной ориентации БПА при m-й паре отсчетов данных относительно земной поверхности.

Для определения дисперсий и m-й единичной оценки координат наземного ИРИ с использованием полученных прямоугольных декартовых координат и точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли, формируемой по результатам n₁-го и n₂-го отсчетов данных, вычисляют:

расстояния и от фазового центра БПА до точки М_m при n₁-м и n₂-м отсчетах данных соответственно (см. фиг. 3);

радиопеленги и наземного ИРИ при n₁-м и n₂-м отсчетах данных соответственно (см. фиг. 2);

углы возвышения и наземного ИРИ при n₁-м и n₂-м отсчетах данных соответственно (см. фиг. 2);

углы места и наземного ИРИ при n₁-м и n₂-м отсчетах данных соответственно (см. фиг. 3);

единичные векторы и направлений от фазового центра БПА на наземный ИРИ в нормальной земной подвижной системе координат , начало которой размещено в точке М_m, при n₁-м и n₂-м отсчетах данных соответственно (см. фиг. 3);

единичный вектор нормали G_m к поверхности Земли при ее кусочно-постоянной аппроксимации в районе точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли в нормальной земной подвижной системе координат , начало которой размещено в точке М_m (см. фиг. 3).

При этом расстояния и от фазового центра БПА до точки М_m при n₁-м и n₂-м отсчетах данных определяют в соответствии с соотношениями

Радиопеленги и наземного ИРИ при n₁-м и n₂-м отсчетах данных с учетом (8) определяют в соответствии с соотношениями

Углы возвышения и наземного ИРИ при n₁-м и n₂-м отсчетах данных с учетом (33) определяют в соответствии с соотношениями

Углы места и наземного ИРИ при n₁-м и n₂-м отсчетах данных с учетом (1), (3), (5), (7) (34) и (35), определяют в соответствии с соотношениями

где и - направляющие косинусы оси аппликат связанной системы координат при n₁-м отсчете данных (см. фиг. 3), определяемые с использованием полученных при n₁-м отсчете данных значений углов курса , тангажа и крена БПА в соответствии с соотношениями

и - направляющие косинусы оси аппликат связанной системы координат при n₂-м отсчете данных (см. фиг. 3), определяемые с использованием полученных при n₂-м отсчете данных значений углов курса , тангажа и крена БПА в соответствии с соотношениями

Единичные векторы и направлений от фазового центра БПА на наземный ИРИ в нормальной земной подвижной системе координат , начало которой размещено в точке М_m, при n₁-м и n₂-м отсчетах данных с учетом (1), (34) и (35) определяют в виде матрицы-столбца размером 3×1 в соответствии с соотношениями

где и - прямоугольные декартовые координаты единичного вектора в нормальной земной подвижной системе координат при n₁-м отсчете данных;

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

Единичный вектор нормали G_m к поверхности Земли при ее кусочно-постоянной аппроксимации в районе точки М_m пересечения линии m-й линии положения наземного ИРИ с поверхностью Земли в нормальной земной подвижной системе координат , начало которой размещено в точке М_m, определяется соотношением

где и - соответственно абсцисса, ордината и аппликата единичного вектора нормали G_m к поверхности Земли в нормальной земной подвижной системе координат ;

- транспонированная матрица-строка размера 1×3, являющаяся матрицей-столбцом размера 3×1;

и - частные производные функции F(x, y), описывающей поверхность Земли по осям соответственно абсцисс и ординат нормальной земной подвижной системы координат в точке М_m пересечения линии m-й линии положения наземного ИРИ с поверхностью Земли с прямоугольными декартовыми координатами и в нормальной земной системе координат Σ_нз={O, X, Y, Z}, которые при кусочно-постоянной аппроксимации поверхности Земли определяют в соответствии с соотношениями

Δl - шаги вдоль осей абсцисс и ординат двухмерной координатной сетки горизонтальной плоскости XOY нормальной земной системы координат Σ_нз={O, X, Y, Z} при кусочно-постоянной аппроксимации поверхности Земли, значения которых выбирают в соответствии с соотношением

где σ_р - априорно известное требуемое значение круговой (радиальной) средней квадратической ошибки определения координат наземного ИРИ.

С использованием полученных данных измерений и расчетов вышеупомянутых параметров определение дисперсий и m-й единичной оценки координат наземного ИРИ осуществляют следующим образом.

При отсчетах данных ошибки радиопеленгования являются статистически независимыми. Поэтому случайные ошибки и определения азимута наземного ИРИ соответственно при n₁-м и n₂-м отсчетах данных, используемых для получения m-й единичной оценки координат и точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли, являются статистически независимыми. Пара статистически независимых случайных ошибок и определения азимута наземного ИРИ приводит к случайной ошибке определения координат наземного ИРИ, которая может быть представлена в виде случайного радиус-вектора η_m, начало которого находится в точке М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли, а конец - в точке М_и априорно неизвестного истинного положения наземного ИРИ (см. фиг. 3). Считая, что точка М_и истинного положения наземного ИРИ и точка М_m m-й единичной оценки координат наземного ИРИ находятся в плоскости, касательной к поверхности Земли в точке М_m пересечения линии положения наземного ИРИ с поверхностью Земли, характеризуемой единичным вектором нормали G_m к поверхности Земли (см. фиг. 3), для случайного радиус-вектора η_m в нормальной земной подвижной системе координат выполняются следующие соотношения

где и - соответственно абсцисса, ордината и аппликата случайного радиус-вектора η_m в нормальной земной подвижной системе координат .

Систему уравнений (44) можно представить в матричном виде следующим образом

где V_m - квадратная матрица размера, имеющая размер 3×3, определяемая соотношением

Тогда из соотношения (45) с учетом (46) для матрицы-столбца случайного радиус-вектора η_m, имеющего размер 3×1, получаем следующее соотношение

где - матрица, обратная матрице V_m.

Корреляционную матрицу ошибок оценивания истинных координат наземного ИРИ, которая, в общем случае, является квадратной матрицей, имеющей размер 3×3, зависящей от статистических характеристик случайного радиус-вектора η_m, в нормальной земной подвижной системе координат можно представить в виде соотношения

где и - дисперсий составляющих линейных ошибок определения координат наземного ИРИ по осям соответственно абсцисс, ординат и аппликат нормальной земной подвижной системы координат ;

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

При этом с учетом (47) и (48) корреляционную матрицу можно определить в соответствии с соотношением

где - операция усреднения по вероятности;

- операция транспонирования матрицы.

Так как случайные ошибки и определения азимута наземного ИРИ являются статистически независимыми, то выполняются следующие соотношения

где и - дисперсии ошибок определения азимута наземного ИРИ соответственно при n₁-м и n₂-м отсчетах данных.

С учетом (46) и (50) соотношение (49) может быть представлено в виде

С учетом (46), (48) и (51) элементы и корреляционной матрицы определяются соотношениями

где D_m - определитель матрицы V_m, определяемый соотношением

Зависимость элементов и корреляционной матрицы , определяемых соотношениями (52)-(54), от углов и (см. фиг. 2 и фиг. 3) можно определить следующим образом.

Для векторов и , получаемых путем векторного перемножения соответственно пары единичных векторов и G_m и пары единичных векторов и G_m, с учетом (12), (18) и (41)-(43) выполняются следующие соотношения

где и - модули векторов и , векторного произведения единичных векторов и от и единичных векторов , и G_m соответственно;

- угол между единичными векторами нормалей и G_m, находящийся в пределах ;

- угол между единичными векторами нормалей и G_m, находящийся в пределах .

Из соотношений (59)-(62) получаем следующие соотношения

Кроме того, с учетом (39) и (40) взаимосвязи между векторами и с одной стороны и единичными векторами соответственно и с другой стороны определяются соотношениями

Зависимость определителя D_m матрицы V_m от углов γ_m и χ_m (см. фиг. 3) можно определить следующим образом.

Для вектора K_m в нормальной земной подвижной системе координат , с учетом (17) и (18) выполняются следующие соотношения

где и - модули вектора K_m(нзп), единичных векторов и , векторного произведения единичных векторов и действительного параметра sin γ_m соответственно.

С учетом (17), (18), (41) и (67) для скалярного произведения вектора K_m(нзп) и единичного вектора нормали G_m к поверхности Земли в нормальной земной подвижной системе координат , начало которой размещено в точке М_m, выполняются следующие соотношения

где и - модули транспонированного вектора-столбца и единичного вектора G_m соответственно.

Из соотношений (58), (68) и (69) следует, что определитель D_m матрицы V_m может быть представлен соотношением

Дисперсии и ошибок определения азимута наземного ИРИ с использованием БПА зависят как от отношения сигнал/шум в каналах БПА, собственных инструментальных погрешностей радиопеленгатора и погрешностей, обусловленных ошибками определения в пространстве угловой ориентации БПА (ошибками определения углов курса, тангажа и крена), но и от углов места и соответственно при n₁-м и n₂-м отсчетах данных (см. фиг. 3). Для обеспечения наименьших ошибок определения азимута в качестве БПА используют плоские пеленгаторные антенные решетки, характеризуемые максимально возможной апертурой в азимутальной плоскости, то есть в плоскости, для которой угол места ε равен нулю. В случае отклонения угла места наземного ИРИ от значения ε=0 проекция апертуры БПА на направление распространения радиосигнала наземного ИРИ уменьшается, что приводит к появлению угломестной зависимости дисперсии ошибок пеленгования в азимутальной плоскости . Поэтому при идентичных условиях определения азимутов и наземного ИРИ с использованием БПА соответственно при n₁-м и n₂-м отсчетах данных угломестные зависимости соответствующих дисперсий и ошибок определения азимута наземного ИРИ определяются соотношениями [2]; [16. Белавин О.В., Зерова М.В. Современные средства радионавигации. - М.: Советское радио, 1965, с. 48-53]:

где - априорно известное значение дисперсии ошибок радиопеленгования в азимутальной плоскости БПА, характеризующее точность определения азимута α ИРИ, расположенного в азимутальной плоскости БПА;

и - углы места наземного ИРИ соответственно при n₁-м и n₂-м отсчетах данных, определяемые в соответствии с соотношениями (36).

С учетом (42)-(54), (65), (70) и (71) диагональные элементы , и корреляционной матрицы можно представить в виде соотношений

где m=1, 2, …, .

В заявляемом способе определения координат наземного ИРИ средневзвешенная оценка аппликаты точки М расположения наземного ИРИ определяется в соответствии с соотношением (32). То есть ошибки определения аппликаты наземного ИРИ зависят от ошибок определения абсциссы и ординаты наземного ИРИ, и, соответственно, полное представление о точности определения координат наземного ИРИ может быть получено с использованием корреляционной матрицы , имеющей размер 2×2, которая является подматрицей корреляционной матрицы ошибок оценивания истинных координат наземного ИРИ, определяемой с учетом (48), соотношением

Дисперсия круговой (радиальной) ошибки m-й единичной оценки координат и точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли в горизонтальной плоскости XOY нормальной земной системы координат Σ_нз={О, X, Y, Z} по определению равна сумме диагональных элементов корреляционной матрицы и, с учетом (75), определяется в соответствии с соотношением

где - след (сумма диагональных элементов) матрицы .

С учетом (72), (73) и (76) дисперсия круговой (радиальной) ошибки и единичной оценки координат наземного ИРИ определяется соотношением

где m=1, 2, …,.

Согласно изобретению с учетом (15), (16), (30) и (31) из общего числа несовпадающих пар пересекающихся полуплоскостей положения наземного ИРИ, соответствующих несовпадающим n₁-м и n₂-м отсчетам данных, выбирают максимальное число независимых пар несовпадающих полуплоскостей положения наземного ИРИ, соответствующих независимым n₁-м и n₂-м отсчетам данных, для которых дисперсии круговых (радиальных) ошибок m-х единичных оценок координат, определяемых соотношением (77), имеют минимальные значения. При этом, в отличие от прототипа [5], значения дисперсий круговых (радиальных) ошибок m-х единичных оценок координат, определяемых соотношением (77), во-первых, зависят не от значений модулей разностей между парами радиопеленгов, а от значений углов γ_m между нормалями и к парам пересекающихся полуплоскостей положения наземного ИРИ (см. рис. 3), и, во-вторых, дополнительно зависят от значений углов χ_m между линиями положения наземного ИРИ и единичным вектором нормали G_m к поверхности Земли, и от значений пар углов места и , пар углов возвышения и и пар углов и между единичными векторами нормалей и к пересекающимся полуплоскостям положения наземного ИРИ соответственно и единичным вектором нормали G_m к поверхности Земли G_m, соответствующих n₁-м и n₂-м отсчетам данных, используемых при формировании m-х единичных оценок координат наземного ИРИ. Учет вышеупомянутых зависимостей позволяет из общего числа несовпадающих пар пересекающихся полуплоскостей положения наземного ИРИ выбрать совокупность независимых пар пересекающихся полуплоскостей положения наземного ИРИ, соответствующих независимым n₁-м и n₂-м отсчетам данных, для которых дисперсии круговых (радиальных) ошибок m-х единичных оценок координат наземного ИРИ имеют минимальные значения, что позволяет повысить точность определения координат наземного ИРИ при многократном радиопеленговании с разных точек траектории движения ЛА.

Согласно изобретению при определении средневзвешенных оценок абсциссы и ординаты точки М_и истинного расположения наземного ИРИ в соответствии с соотношением (28) значения весовых коэффициентов w_Xm и w_Ym m-й единичной оценки соответственно абсцисс и ординат с учетом (15), (16), (30), (31) (72) и (73) определяются в соответствии с соотношениями

где m=1, 2, …, - порядковые номера независимых пар пересекающихся полуплоскостей положения наземного ИРИ, соответствующих независимым n₁-м и n₂-м отсчетам данных.

При этом, согласно (78) и (79), в отличие от прототипа [5], значения весовых коэффициентов w_Xm и w_Ym m-й единичной оценки соответственно абсцисс и ординат наземного ИРИ, во-первых, зависят не от значений модулей разностей между парами радиопеленгов, а от значений углов γ_m между парами полуплоскостей положения наземного ИРИ (см. рис. 3), и, во-вторых, дополнительно зависят от значений углов χ_m между линиями положения наземного ИРИ и единичным вектором нормали G_m к поверхности Земли, и от значений пар углов места и , пар углов возвышения и и пар углов и между единичными векторами нормалей и к полуплоскостям положения наземного ИРИ и единичным вектором нормали G_m к поверхности Земли G_m, соответствующих независимым n₁-м и n₂-м отсчетам данных, используемых при формировании m-х единичных независимых оценок координат наземного ИРИ. Учет вышеупомянутых зависимостей позволяет при проведении весовой обработки полученных данных увеличить вклад отсчетов, при которых соответствующие независимые пары пересекающихся полуплоскостей положения наземного источника радиоизлучения наиболее близки к взаимно ортогональным полуплоскостям и обеспечивают углы пересечения линий положения наземного источника радиоизлучения с поверхностью Земли, наиболее близкие к нормали к поверхности Земли в зоне радиопеленгования, при которых нормали к независимым парам пересекающихся полуплоскостей положения наземного источника радиоизлучения расположены в плоскостях, наиболее близких к горизонтальной плоскости, и, наконец, при которых соответствующие независимым парам отсчетов модули углов возвышения наземного ИРИ характеризуются максимальными значениями, а модули углов места наземного ИРИ характеризуются минимальными значениями, что позволяет повысить точность определения координат при многократном радиопеленговании с разных точек траектории движения ЛА.

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

В заявляемом способе определения координат наземного источника радиоизлучения с борта летательного аппарата прямоугольные декартовые координаты и точки M_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли определяют как прямоугольные декартовые координаты точки пересечения m-й линии положения наземного источника радиоизлучения, описываемой параметрическим уравнением (27), с поверхностью Земли, описываемой пятым соотношением системы соотношений (20), с использованием следующей итерационной процедуры.

Для однозначного определения с использованием соотношения (27) точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли (см. фиг. 3), определяемой пятым соотношением системы соотношений (20), параметрическую переменную b_m, соответствующую m-й линии положения , выбирают в соответствии с соотношением

где sgn(K_Zm) - функция знака, определяемая соотношением [13]

При кусочно-постоянной аппроксимации поверхности Земли, определяемой пятым соотношением системы соотношений (20), решение системы уравнений (27) с учетом (80) относительно прямоугольных декартовых координат и точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли осуществляют с использованием итерационной процедуры определения прямоугольных декартовых координат x_m,J, y_m,J и z_m,J точки М_m с использованием полученных прямоугольных декартовых координат х_m,0, y_m,0 и z_m,0 точки начального приближения координат наземного ИРИ в соответствии с соотношениями

где x_m,J, y_m,J и z_m,J - начальные значения соответственно абсциссы, ординаты и аппликаты прямоугольных декартовых координат точки М_m на (j+1)-м шаге (итерации) итерационной процедуры;

x_m,j+1, y_m,j+1 и z_m,j+1 - конечные значения соответственно абсциссы, ординаты и аппликаты прямоугольных декартовых координат точки М_m на (j+1)-м шаге (итерации) итерационной процедуры;

j=0, 1, …, (J-1) - порядковые номера начальных значений прямоугольных декартовых координат точки М_m на (j+1)-х шагах (итерациях) итерационной процедуры;

- абсцисса точки Р_m начального приближения координат наземного ИРИ, определяемая в соответствии с соотношением (22);

- абсцисса точки Р_m начального приближения координат наземного ИРИ, определяемая в соответствии с соотношением (23);

- аппликата точки Р_m начального приближения координат наземного ИРИ, равная аппликате точки Р_m пересечения линии положения наземного ИРИ с горизонтальной плоскостью z=z_ум (см. фиг. 3);

J=1, 2, … - минимальное число шагов (итераций) итерационной процедуры, обеспечивающее выполнение пары неравенств, определяемых соотношениями

где δ_X и δ_Y - априорно известные максимально допустимые погрешности определения соответственно абсциссы x_m,J и ординаты y_m,J прямоугольных декартовых координат точки размещения ИРИ на земной поверхности, обусловленные итерационной процедурой определения прямоугольных декартовых координат x_m,J, y_m,J и z_m,J точки М_m при наличии погрешностей кусочно-постоянной аппроксимации поверхности Земли, которые выбирают на порядок меньшими априорно известного требуемого значения круговой (радиальной) средней квадратической ошибки σ_р определения координат наземного ИРИ.

По завершении итерационной процедуры, (82), характеризуемой минимальным числом J ее шагов (итераций), обеспечивающим выполнение пары неравенств (83), прямоугольные декартовые координаты и точки М_m пересечения m-й линии положения наземного ИРИ с поверхностью Земли определяются в виде следующих равенств

Структурная схема алгоритма, иллюстрирующего заявляемый способ определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата в части вычисления единичной оценки координат наземного источника радиоизлучения в соответствии с соотношениями (27), (80)-(84), приведена на фиг.5.

Таким образом, в заявляемом способе определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата в отличие от прототипа [5] точку пересечения линии положения наземного источника радиоизлучения с поверхностью Земли определяют с использованием итерационной процедуры поиска вышеупомянутой точки при кусочно-постоянной аппроксимации поверхности Земли с шагом сходимости итерационной процедуры, и, соответственно, с минимальным числом шагов итерационной процедуры, которые выбирают в зависимости от изрезанности рельефа поверхности Земли в зоне радиопеленгации, определяемого формой функции F(x_m,J, y_m,J), в зависимости от априорно известной заданной точности определения координат наземного ИРИ, характеризуемой круговой (радиальной) средней квадратической ошибкой σ_р определения координат наземного ИРИ, и в зависимости от угла χ_вm между m-й линией положения наземного ИРИ и местной вертикалью В поверхности Земли, совпадающей с осью аппликат нормальной земной подвижной системы координат , который характеризуется совокупностью отношений и прямоугольных декартовых координат K_Xm, K_Ym и K_Zm вектора K_m в нормальной земной подвижной системе координат Σ_нзпm=[O_m, X, Y, Z], начало О_m которой находится на m-й линии положения наземного ИРИ, что позволяет дополнительно повысить точность определения координат наземного ИРИ при одновременном снижении вычислительных затрат.

Следует отметить, что взаимосвязь между углом χ_вm и отношений и прямоугольных декартовых координат вектора K_m определяется в соответствии с соотношением

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

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

Под воздействием электромагнитного поля радиосигнала, излучаемого на частоте ƒ наземным источником радиоизлучения, характеризуемого при n-м отсчете данных при радиопеленговании с борта ЛА амплитудой Е_n и фазой ς_0n в фазовом центре БПА, а также - направлением единичного вектора S_n распространения электромагнитной волны, характеризуемым в связанной системе координат азимутом α_n и углом места ε_n (см. фиг. 2), на выходах антенных элементов 1.1, 1.2, …, 1.H бортовой пеленгаторной антенны 1, формируются радиосигналы , которые описываются соотношениями

где h=1, 2, …, H - порядковый номер антенных элементов БПА;

U_h,n=D_h(α_n, ƒ)cos ε_nE_n - амплитуда радиосигнала наземного ИРИ, принятого h-м антенным элементом БПА;

- фаза радиосигнала наземного ИРИ, принятого h-м антенным элементом БПА;

D_h(α_n, ƒ) и ϕ_h(α_n, ƒ) - соответственно амплитудная диаграмма направленности (по полю) и фазовая диаграмма h-го антенного элемента БПА в азимутальной плоскости БПА, соответствующей плоскости связанной системы координат (см фиг. 2);

t - время;

i= - мнимая единица.

Радиосигналы, принятые антенными элементами 1.1, 1.2, 1.Н бортовой пеленгаторной антенны 1, описываемые соотношениями (86), с соответствующих выходов бортовой пеленгаторной антенны 1 поступают на входы соответствующих каналов радиоприемного блока 2.

Запуск функционирования радиоприемного блока 2, блока 4 определения угловой ориентации бортовой пеленгаторной антенны и блока 5 определения координат фазового центра бортовой пеленгаторной антенны в процессе полета ЛА осуществляется многократно периодически для получения совокупности отсчетов данных, с их общим числом n_max, не меньшем двух. При этом запуск функционирования радиоприемного блока 2, блока 4 определения угловой ориентации бортовой пеленгаторной антенны и блока 5 определения координат фазового центра бортовой пеленгаторной антенны при каждом n-м отсчете данных (1≤n≤n_max) осуществляется одновременно по команде модуля 16 программного управления.

При этом по команде, соответствующей n-му отсчету данных, поступающей с выхода модуля 16 программного управления на управляющий вход радиоприемного блока 2, сигналы , поступившие на входы соответствующих Н каналов радиоприемного блока 2 (h=1, 2, …, H), подвергаются типовым для современных многоканальных радиоприемных блоков со спектральной обработкой сигналов преобразованиям: синхронной фильтрации в полосе частот ΔF, усилению, преобразованию на промежуточную частоту с идентичными общими комплексными коэффициентами передачи каналов радиоприемного блока 2, получаемыми либо с использованием идентичных каналов радиоприемного блока 2, либо с использованием результатов калибровки неидентичности общих комплексных коэффициентов передачи каналов радиоприемного блока 2, синхронному преобразованию сигналов промежуточной частоты в цифровые сигналы путем дискретизации по уровню и времени принимаемых сигналов и осуществления дискретного преобразования Фурье (в частности - быстрого преобразования Фурье) с получением действительной и мнимой составляющих сигналов, являющихся цифровыми эквивалентами принимаемых радиосигналов

где Δt - интервал равномерной дискретизации сигналов по времени;

- объем обрабатываемого массива временных отсчетов радиосигналов;

ΔT - требуемое (допустимое при фильтрации сигналов в полосе частот ΔF) время наблюдения радиосигналов;

ν_д=0, 1, …, (V_д-1) - порядковые номера временных отсчетов радиосигналов;

l_ƒ=0, 1, …, (V_д-1) - номер спектральных составляющих радиосигналов.

При этом после осуществления преобразования Фурье согласно вышеупомянутому выражению принятые радиосигналы представляются в виде V_д комплексных спектральных составляющих, частотное расстояние между которыми, то есть фактически частотное разрешение, составляет . В дальнейшем для упрощения обработки используют только комплексных спектральных составляющих, а остальные спектральных составляющих, соответствующих отрицательным частотам, полагают равными нулю. Таким образом, прием совокупности радиосигналов, отфильтрованных в каждом из каналов радиоприемного блока 2 в широкой полосе частот ΔF с последующим разделением по частоте методом преобразования Фурье эквивалентен одновременному приему сигналов в соответствующих элементарных частотных каналах с шириной полосы и общим количеством . Указанное обстоятельство обеспечивает возможность частотной селекции радиосигналов и получение высокой чувствительности радиоприемного блока 2 (за счет уменьшения действующего значения напряжения внутреннего шума, обратно пропорционального величине ) при одновременном увеличении быстродействия радиопеленгатора и повышении частотной избирательности радиоприемного блока 2.

По полученным действительным и мнимым спектральным компонентам сигналов в каждом канале радиоприемного блока 2 в зависимости от метода радиопеленгования (амплитудного или фазового), реализуемого при радиопеленговании с использованием бортовой пеленгаторной антенны 1, на частоте ƒ излучения наземного ИРИ определяется параметр радиосигнала (амплитуда или фаза), значение которого, согласно (85), зависит от азимута α_n наземного ИРИ. При этом, в частности, зависящие от азимута α_n амплитуда U_h,n(α_n) и фаза Ψ_h,n(α_n) радиосигнала определяются е соответствии с соотношениями

где h=1, 2, …, H;

n=1, 2, … - порядковые номера отсчетов данных.

Сигналы, соответствующие значениям измеренных параметров радиосигналов (амплитуд U_h,n(α_n) или фаз Ψ_h,n(α_n)), с Н выходов радиоприемного блока 2 поступают на соответствующие входы вычислителя азимута 3.

По команде, соответствующей n-му отсчету данных, поступающей с выхода модуля 16 программного управления на управляющий вход блока 4 определения угловой ориентации бортовой пеленгаторной антенны в блоке 4 определения угловой ориентации бортовой пеленгаторной антенны производится измерение углов курса ψ_n, тангажа μ_n и крена ϑ_n бортовой пеленгаторной антенны, определяющие угловую ориентацию продольной О'Х', поперечной O'Y' и нормальной O'Z' осевых линий бортовой пеленгаторной антенны 1 в нормальной земной подвижной системе координат Σ_нзп={О', X, Y, Z} при n-м отсчете данных (см. фиг. 2). Сигналы, соответствующие полученным значениям углов курса ψ_n, тангажа μ_n и крена ϑ_n бортовой пеленгаторной антенны с первого, второго и третьего выходов блока 4 определения угловой ориентации бортовой пеленгаторной антенны поступают соответственно на объединенные второй вход блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения и четвертый вход вычислителя 11 дисперсий единичной оценки координат наземного ИРИ, объединенные третий вход блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения и пятый вход вычислителя 11 дисперсий единичной оценки координат наземного ИРИ и объединенные четвертый вход блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения и шестой вход вычислителя 11 дисперсий единичной оценки координат наземного ИРИ.

По команде, соответствующей n-му отсчету данных, поступающей с выхода модуля 16 программного управления на управляющий вход блока 5 определения координат фазового центра бортовой пеленгаторной антенны в блоке 5 определения координат фазового центра бортовой пеленгаторной антенны производится измерение прямоугольных декартовых координат х_цn, y_цn и z_цn фазового центра бортовой пеленгаторной антенны 1 в нормальной земной системе координат Σ_нз={O, X, Y, Z} при n-м отсчете данных (см. фиг. 3). Сигналы, соответствующие полученным значениям координат х_цn, y_цn и z_цn фазового центра бортовой пеленгаторной антенны, с выхода блока 5 определения координат фазового центра бортовой пеленгаторной антенны в виде кодированных значений поступают на объединенные третий вход блока 8 вычисления параметров линии положения и седьмой вход вычислителя 11 дисперсий единичной оценки координат наземного ИРИ.

По команде, соответствующей n-му отсчету данных, поступающей с выхода модуля 16 программного управления на управляющий вход вычислителя азимута 3, в вычислителе азимута 3 путем сравнения параметров (амплитуд U_h,n(α_n) или фаз Ψ_h,n(α_n)) радиосигналов и с использованием априорно заданных параметров вычислений, соответствующих реализуемому методу радиопеленгования и структуре бортовой пеленгаторной антенны 1, в пределах от 0 до 360° определяется оценка азимута α_n наземного ИРИ, соответствующая n-му отсчету данных. Сигнал, соответствующий полученному значению азимута α_n наземного ИРИ, с выхода вычислителя азимута 3, поступает на объединенные первые входы блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения и вычислителя 11 дисперсий единичной оценки координат наземного ИРИ.

По команде, соответствующей n-му отсчету данных, поступающей с выхода модуля 16 программного управления на управляющий вход блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения в блоке 6 определения полуплоскости положения и формирования нормали к полуплоскости положения в соответствии с соотношениями (5) и (11)-(13) производится вычисление прямоугольных декартовых координат Q_Xn, Q_Yn и Q_Zn единичного вектора Q_n азимутальной линии n-й полуплоскости положения и прямоугольных декартовых координат N_Xn, N_Yn и N_Zn нормали N_n к n-й полуплоскости положения в нормальной земной подвижной системе координат при n-м отсчете данных (см. фиг. 3). Сигналы, соответствующие полученным значениям координат Q_Xn, Q_Yn и Q_Zn единичного вектора Q_n азимутальной линии n-й полуплоскости положения, с первого выхода блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения в виде кодированных значений поступают на первый вход запоминающего устройства 7 параметров единичных векторов азимутальных линий полуплоскостей положения и нормалей к полуплоскостям положения. Сигналы, соответствующие полученным значениям координат N_Xn, N_Yn и N_Zn нормали N_n к n-й полуплоскости положения, с второго выхода блока 6 определения полуплоскости положения и формирования нормали к полуплоскости положения в виде кодированных значений поступают на второй вход запоминающего устройства 7 параметров единичных векторов азимутальных линий полуплоскостей положения и нормалей к полуплоскостям положения.

По завершении совокупности n_max отсчетов данных, в результате которых в вышеупомянутых блоках устройства определения координат наземного ИРИ при радиопеленговании с борта ЛА определяются совокупности азимутов α_n, углов курса Ψ_n, тангажа μ_n и крена ϑ_n бортовой пеленгаторной антенны 1, координаты х_цп, y_цn и z_цn фазовых центров бортовой пеленгаторной антенны 1, координаты Q_Xn, Q_Yn и Q_Zn единичных векторов Q_n азимутальных линий n-х полуплоскостей положения и координаты N_Xn, N_Yn и N_Zn нормалей N_n к n-м полуплоскостям положения для n-х отсчетов данных, где 1≤n≤n_max, по команде, поступающей с выхода модуля 16 программного управления на управляющий вход запоминающего устройства 7 параметров единичных векторов азимутальных линий полуплоскостей положения и нормалей к полуплоскостям положения сигналы, соответствующие кодированным значениям координат Q_Xn, Q_Yn и Q_Zn единичного вектора Q_n азимутальной линии n-й полуплоскости положения и сигналы, соответствующие кодированным значениям координат N_Xn, N_Yn и N_Zn нормалей N_n к n-м полуплоскостям положения соответственно с первого и второго выходов запоминающего устройства 7 параметров единичных векторов азимутальных линий полуплоскостей положения и нормалей к полуплоскостям положения поступают соответственно на первый вход блока 8 вычисления параметров линии положения и объединенные второй вход блока 8 вычисления параметров линии положения и третий вход вычислителя 11 дисперсий единичной оценки координат наземного ИРИ.

По команде, поступающей с выхода модуля 16 программного управления на управляющий вход блока 8 вычисления параметров линии положения в блоке 8 вычисления параметров линии положения в соответствии с соотношениями (18) и (22)-(26) производится вычисление координат K_Хm, K_Ym и K_Zm m-х векторов K_m в нормальной земной подвижной системе координат Σ_нзпm={O_m, X, Y, Z}, находящихся на m-х линиях пересечения несовпадающих пар пересекающихся полуплоскостей положения наземного ИРИ, и координат и m-х точек Р_m начального приближения координат наземного ИРИ в нормальной земной системе координат Σ_нз={О, X, Y, Z} (см. фиг. 3), где 1≤m≤, причем представляет собой максимальное число пар пересекающихся полуплоскостей положения наземного ИРИ, соответствующих несовпадающим n₁-м и n₂-м отсчетам данных, не превышающих максимальное число m_max пар полуплоскостей положения наземного ИРИ, определяемое соотношением (16).

Сигналы, соответствующие полученным значениям координат K_Хm, K_Ym и K_Zm m-х векторов K_m, с первого выхода блока 8 вычисления параметров линии положения в виде кодированных значений поступают на первый вход запоминающего устройства 9 параметров линий положения.

Сигналы, соответствующие полученным значениям координат и m-х точек Р_m начального приближения координат наземного ИРИ с второго выхода блока 8 вычисления параметров линии положения в виде кодированных значений поступают на второй вход запоминающего устройства 9 параметров линий положения.

По команде, поступающей с выхода модуля 16 программного управления на управляющий вход запоминающего устройства 9 параметров линий положения сигналы, соответствующие кодированным значениям координат K_Хm, K_Ym и K_Zm m-х векторов K_m, и сигналы, соответствующие кодированным значениям координат и m-х точек Р_m начального приближения координат наземного ИРИ соответственно с первого и второго выходов запоминающего устройства 9 параметров линий положения поступают соответственно на первый и второй входы вычислителя 10 единичной оценки координат наземного ИРИ.

По команде, поступающей с выхода модуля 16 программного управления на управляющий вход вычислителя 10 единичной оценки координат наземного ИРИ в вычислителе 10 единичной оценки координат наземного ИРИ с использованием итерационной процедуры поиска точек пересечения линий положения наземного источника радиоизлучения с поверхностью Земли (см. фиг. 5) в соответствии с соотношениями (27), (80)-(84) производится вычисление прямоугольных декартовых координат и точек М_m пересечения m-х линий положения наземного ИРИ с поверхностью Земли в нормальной земной системе координат Σ_нз={О, X, Y, Z}, где 1≤m≤.

Сигналы, соответствующие полученным значениям координат и точек М_m пересечения m-х линий положения наземного ИРИ с поверхностью Земли, с выхода вычислителя 10 единичной оценки координат наземного ИРИ в виде кодированных значений поступают на объединенные второй вход вычислителя 11 дисперсий единичной оценки координат наземного ИРИ и вход запоминающего устройства 12 единичных оценок координат наземного ИРИ.

По команде, поступающей с выхода модуля 16 программного управления на управляющий вход вычислителя 11 дисперсий единичной оценки координат наземного ИРИ в вычислителе 11 дисперсий единичной оценки координат наземного ИРИ в соответствии с соотношениями (33)-(38), (42), (43), (52), (53), (58), (71) и (76) производится вычисление дисперсий и составляющих линейных ошибок m-х единичных оценок координат наземного ИРИ по осям соответственно абсцисс и ординат нормальной земной подвижной системы координат и дисперсий круговых (радиальных) ошибок m-х единичных оценок координат наземного ИРИ в горизонтальной плоскости нормальной земной подвижной системы координат , где 1≤m≤.

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

По команде, поступающей с выхода модуля 16 программного управления на управляющий вход запоминающего устройства 12 единичных оценок координат наземного ИРИ, сигналы, соответствующие кодированным значениям координат и точек М_m пересечения m-х линий положения наземного ИРИ с поверхностью Земли с выхода запоминающего устройства 12 единичных оценок координат наземного ИРИ поступают на первый вход вычислителя 14 координат наземного ИРИ.

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

По команде, поступающей с выхода модуля 16 программного управления на управляющий вход вычислителя 14 координат наземного ИРИ, в вычислителе 14 координат наземного ИРИ в соответствии с соотношениями (28)-(32) осуществляют выбор координат и точек М_m пересечения m-х линий положения наземного ИРИ, соответствующих независимым парам пересекающихся полуплоскостей положения и характеризуемых наименьшими значениями дисперсий круговых (радиальных) ошибок m-х единичных оценок координат наземного ИРИ, вычисление весовых коэффициентов w_Xm и w_Ym m-х единичных оценок соответственно абсцисс и ординат точек М_m пересечения m-х линий положения наземного ИРИ, где 1≤m≤, и вычисление с учетом полученных значений весовых коэффициентов w_Xm и w_Ym прямоугольных декартовых координат и точки М расположения наземного ИРИ в нормальной земной системе координат Е_нз={O, X, Y, Z}.

Сигнал, соответствующий полученным значениям координат и точки М расположения наземного ИРИ, с выхода вычислителя 14 координат наземного ИРИ в виде кодированных значений поступает на вход модуля 15 картографирования и индикации.

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

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

Бортовая пеленгаторная антенна 1 может быть реализована, например, в соответствии с описанием к патенту Российской Федерации №2027261 [17. Пеленгационная антенная система. Патент Российской Федерации №2027261, МПК G01Q 21/24] или в соответствии с описанием к патенту Российской Федерации №2310956 [18. Высоконаправленная кольцевая фазированная антенная решетка. Патент Российской Федерации №2310956, МПК G01Q 21/06, Бюл. №32, опубл. 20.11.2007]. При этом варианты реализации антенных элементов 1.1, 1.2, …, 1.Н бортовой пеленгаторной антенны 1 приведены, например, в работах [19. Жук М.С., Молочков Ю.Б. Проектирование антенно-фидерных устройств. - М. - Л.: Энергия, 1966, с. 108-158]; [20. Саидов А.С., Тагилаев А.Р., Алиев Н.М., Асланов Г.К. Проектирование фазовых автоматических радиопеленгаторов. - М.: Радио и связь, 1997, с. 62-65].

Радиоприемный блок 2 может быть реализован по схеме супергетеродинного радиоприемного устройства с преобразованием частоты, с использованием параллельно-последовательного способа обзора рабочего диапазона частот и аналого-цифрового преобразования радиосигналов с получением их спектральных характеристик на основе быстрого преобразования Фурье, что в виде описаний схемо-технических решений изложено, например, в работах [21. Рембовский А.М., Ашихмин А.В., Козьмин В.А. Радиомониторинг: задачи, методы, средства / Под редакцией А.М. Рембовского. 3-е изд. перераб. и доп. М.: Горячая линия - Телеком. 2012. с. 28-111;], [22. Побережский Е.С. Цифровые радиоприемные устройства. - М.: Радио и связь. - 184 с.], [23. Ашихмин А.В. Сергеев В.Б., Сергиенко А.Р. Радиоприемные тракты комплексов автоматизированного радиоконтроля: особенности, решения и перспективы. - Специальная техника, специальный выпуск, 2002, с. 57-64] и [24. Литюк В.И., Литюк Л.В. Методы цифровой многопроцессорной обработки ансамблей радиосигналов. - М.: СОЛОН-ПРЕСС, 2007. - 592 с.]. При этом для осуществления дискретизации сигналов в варианте реализации радиоприемного блока 2 может использоваться 14-ти разрядный аналого-цифровой преобразователь AD6644 фирмы Analog Devices и микросхемы DIGITAL DOWN CONVERTER (DDC) AD6224 фирмы Analog Devices. Для получения спектральных характеристик радиосигналов на основе быстрого преобразования Фурье в варианте реализации радиоприемного блока 2 может использоваться цифровой сигнальный процессор ADSP TS-202 фирмы Analog Devices, который имеет суперскалярную архитектуру, оптимизирован для задач цифровой обработки сигналов и выполняет преобразование Фурье радиосигналов, например, на 1024 точки за 20 мкс.

Вычислитель азимута 3 может быть реализован, например, а соответствии с алгоритмом, приведенным в описании к патенту Российской Федерации №2521959 [25. Амплитудный способ радиопеленгования и радиопеленгатор для его осуществления. Патент РФ на изобретение №2521959, МПК G01S 3/28, G01S 5/04, Бюл. №19, опубл. 10.07.2014] на базе сигнального процессора ADSP TS-201 фирмы Analog Devices.

Блок 4 определения угловой ориентации бортовой пеленгаторной антенны может быть реализован, например, в виде бортовой инерциальной навигационной системы, описание которой приведено в работах [26. Матвеев В.В., Распопов В.Я. Основы построения бесплатформенных инерциальных навигационных систем / Под общей ред. В.Я. Распопова. - СПб.: ГНЦ РФ ОАО «Концерн «ЦНИИ «Электроприбор», 2009. - 280 с.], [27. Бромберг. П.В. Теория инерциальных систем навигации. - М.: Наука, 1979. - 291 с.] и [28. Ориентация и навигация подвижных объектов: современные информационные технологии / Под ред. Б.С. Алешина, К.К. Веремеенко, А.И. Черноморского. - М.: ФИЗМАТ ЛИТ, 2006. - 424 с.], или с использованием устройства МРК-11 [29. Аппаратура МРК-11. Руководство по эксплуатации. УЭ 2.517.006 РЭ. Научно-исследовательский институт радиотехники КГТУ. - Красноярск, 2004].

Необходимо отметить, что в качестве блока 4 определения угловой ориентации бортовой пеленгаторной антенны может использоваться также блок определения угловой ориентации осевых линий ЛА, входящий в состав ЛА, на котором размещена БПА. При этом в качестве значений углов курса, тангажа и крена БПА можно использовать значения углов соответственно курса, тангажа и крена ЛА при условии попарного совпадения положительных направлений продольных, поперечных и нормальных осевых линий БПА и ЛА.

Блок 5 определения координат фазового центра бортовой пеленгаторной антенны может быть реализован, например, в виде спутникового навигатора, описание которого приведено в работах [28], [30. Помыкаев И.И., Селезнев В.Л., Дмитроченко Л.А. Навигационные приборы и системы. - М.: Машиностроение, 1983. - 456 с.] и [31. Спутниковые радионавигационные системы / Под ред. В.С. Шебшаевича..- М.: Радио и связь, 1993. - 272 с.], или с использованием вышеупомянутого устройства МРК-11 [29].

Необходимо отметить, что в качестве блока 5 определения координат фазового центра бортовой пеленгаторной антенны может использоваться также блок определения координат центра масс ЛА, входящий в состав ЛА, на котором размещена БПА. При этом в качестве значений координат x_цn, y_цn и z_цn фазового центра бортовой пеленгаторной антенны в абсолютном большинстве практически важных случаев, определяемых пренебрежимо малыми размерами ЛА относительно расстояний от ЛА до наземного ИРИ, можно использовать значения координат соответственно х_ЛАn, у_ЛАn и z_ЛАn центра масс ЛА, соответствующих n-м отсчетам данных. Вместе с тем, в случае необходимости, как и в аналоге [4], с использованием результатов измерений координат х_ЛАn, у_ЛАn и z_ЛАn центра масс ЛА можно определять вышеупомянутые значения координат х_цn, y_цn и z_цn фазового центра бортовой пеленгаторной антенны в соответствии с соотношениями

где х_цсв, х_цсв и х_цсв - связанные прямоугольные декартовые координаты фазового центра бортовой пеленгаторной антенны относительно центра масс ЛА.

Блок 6 определения полуплоскости положения и формирования нормали к полуплоскости положения, запоминающее устройство 7 параметров единичных векторов азимутальных линий полуплоскостей положения и нормалей к полуплоскостям положения, блок 8 вычисления параметров линии положения, запоминающее устройство 9 параметров линий положения, вычислитель 10 единичной оценки координат наземного ИРИ, вычислитель 11 дисперсий единичной оценки координат наземного ИРИ, запоминающее устройства 12 единичных оценок координат наземного ИРИ, запоминающее устройство 13 дисперсий единичных оценок координат наземного ИРИ, вычислитель 14 координат наземного ИРИ, модуль 15 картографирования и индикации и модуль 16 программного управления могут быть реализованы с использованием схемо-технических решений, изложенных, например, в работах [32. Потехин Д.С., Тарасов И.Е. Разработка систем цифровой обработки сигналов на базе ПЛИС. - М.: Горячая линия - Телеком, 2007.], [33. Сергиенко А.Б. Цифровая обработка сигналов. - СПб.: Питер, 2003. - 608 с.], [34. Шевкоплес Б.В. Микропроцессорные структуры. Инженерные решения: Справочник. - 2-е изд. перераб. и доп. - М.: Радио и связь, 1990. - 512 с.] и [35. Лебедев О.М. Микросхемы памяти и их применение. - М.: Радио и связь, 1990. - 160 с.], на базе процессоров Texas Instruments TMS 320 С 6416/6713 и программируемых логических интегральных схем фирмы Altera.

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

Результаты статистического моделирования показали, что при одинаковых условиях полета ЛА предложенный способ определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата по сравнению с прототипом [5] обеспечивает повышение точности определения координат наземного источника радиоизлучения не менее чем в 2 раза при одновременном снижении вычислительных затрат от 20% до 70% в зависимости от изрезанности рельефа поверхности Земли.

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

Изобретение относится к радиотехнике и может быть использовано для определения координат наземных источников радиоизлучения (ИРИ) при радиопеленговании с борта летательного аппарата (ЛА). Достигаемый технический результат - повышение точности определения координат наземных ИРИ и снижение вычислительных затрат при радиопеленговании с борта ЛА. Указанный результат достигается за счет того, что осуществляют прием радиосигналов бортовой пеленгаторной антенной (БПА), частотную селекцию радиосигналов, определение линий радиопеленгов в азимутальной плоскости БПА, регистрацию полученных данных периодически отсчетами, формирование не менее одной независимой пары пересекающихся полуплоскостей положения наземного ИРИ, ортогональных азимутальной плоскости БПА, проходящих через каждую полученную линию радиопеленга, выбор и весовую обработку пар независимых отсчетов данных, учитывающих зависимости дисперсий оценок координат наземного ИРИ от взаимного расположения в пространстве ЛА и наземного ИРИ. При этом дополнительно введены операции формирования нормалей к полуплоскостям положения наземного ИРИ, определения не менее одной линии положения наземного ИРИ как линии пересечения независимой пары пересекающихся полуплоскостей положения наземного ИРИ, параметры которой определяют из условия ортогональности к вышеупомянутым нормалям, и определения координат наземного ИРИ как точки пересечения линии положения наземного ИРИ с поверхностью Земли с использованием итерационной процедуры ее поиска. Кроме того, при выборе и весовой обработке пар независимых отсчетов данных дополнительно учтены зависимости дисперсий оценок координат наземного ИРИ от параметров угловой ориентации БПА и от углов пересечения линии положения и нормалей к полуплоскостям положения наземного ИРИ с поверхностью Земли. 1 з.п. ф-лы, 5 ил.

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

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

2. Способ определения координат наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата по п. 1, отличающийся тем, что точку пересечения линии положения наземного источника радиоизлучения с поверхностью Земли определяют с использованием итерационной процедуры ее поиска, причем шаг сходимости итерационной процедуры и минимальное число шагов итерационной процедуры выбирают зависящими от изрезанности рельефа поверхности Земли в зоне радиопеленгования, априорно известной заданной точности определения координат наземного источника радиоизлучения и угла между линией положения наземного источника радиоизлучения и местной вертикалью поверхности Земли.