Разностно-дальномерный способ определения координат источника радиоизлучения
DIFFERENCE-RANGEFINDER METHOD FOR DETERMINATION OF RADIATION SOURCE COORDINATES

Изобретение относится к радиотехнике и может быть использовано в многопозиционных радиотехнических системах для определения координат (ОК) заданных источников радиоизлучения (ИРИ) с кодовым или временным разделением каналов.
i,k
i,k
Известен способ измерения взаимной задержки MSK сигналов пакетных радиосетей в разностно-дальномерной системе местоопределения (см. Пат, РФ №2623094, МПК G01S 5/06, опубл. 22.06.2017, бюл. №18). Способ предполагает прием периферийными пунктами приема (ППП) разностно-дальномерной системы местоопределения (РДСМО) сигналов пакетных радиосетей, измерение времени прихода сигналов относительно единой шкалы времени, передачу по линии связи измеренных значений на центральный пункт приема и обработки (ЦППО), где вычисляют взаимные задержки сигналов r, причем модуль взаимно корреляционной функции (ВКФ) R(τ) вычисляют с помощью парциальных ВКФ, чем достигается повышение точности измерения τ.
Однако аналогу присущи недостатки, ограничивающие его применение. Основными из них являются недостаточная точность ОК ИРИ, сложность реализации ППП. Способ предусматривает использование на всех ППП системы единого времени при измерении задержки сигналов. Кроме того, отсутствие привязки измеренных значений задержки сигнала на ППП к заданному источнику радиоизлучения при оценке их координат в условиях, когда работа нескольких ИРИ осуществляется на одной частоте (ИРИ с временным или кодовым разделением) приводит к ошибкам оценки координат ИРИ.
i,k
Известны Пат. РФ №№2258242, 2309420, 2521084 и др., эффективность применения которых резко падает при оценке координат ИРИ с временным или кодовым разделением каналов. В зависимости от длительности оценки задержек сигналов τ возникают ошибки ОК различной природы.
Известен разностно-дальномерный способ определения координат источника радиоизлучения (см. Пат. №2539968, МПК G01S 3/46, опубл. 27.01.2015, бюл. №3).
i,k
i,k
Способ-аналог заключается в приеме сигналов ИРИ разнесенными в пространстве периферийными пунктами приема, связанными с ЦППО командными линиями связи и линиями аналоговой ретрансляции сигналов, причем по командным линиям связи с ЦППО на ППП передают команды настройки на частоту сигнала ИРИ, а по линиям аналоговой ретрансляции принятые в ППП сигналы ИРИ передают на ЦППО, где измеряют разности времени приема сигналов в ППП и ЦППО τ, i и k - номера ЦППО и ППП соответственно, при этом задержки τ определяют как аргумент максимизации модуля взаимнокорреляционной функции, а на их основе вычисляют координаты ИРИ.
Аналог обеспечивает уменьшение количества вычислений при оценке временных задержек, принимаемых на ППП сигналов ИРИ в реализуемой ЦППО процедуре. Кроме того, выполнение ППП максимально упрощено, что способствует широкому применению способа-прототипа в РДСМО. Каждый ППП представляет собой совокупность устройств, выделяющих радиосигналы ИРИ на фоне помех, а также устройств, организующих аналоговые линии ретрансляции.
k
Однако аналогу присущ недостаток, ограничивающий его применение: недостаточная точность определения координат ИРИ, а в ряде случаев он теряет работоспособность. Причиной ошибок ОК прежде всего является многолучевость распространения радиоволн. В результате последние приходят на измеритель с разных направлений и с разной задержкой. Ошибочное определение координат возникает и в тех случаях, когда на частоте искомого ИРИ работают несколько корреспондентов (ИРИ) с временным или кодовым разделением каналов, например, сети Wi-Fi, GSM 2G, UMTS и др. Кроме того, в способе-аналоге сделано допущение о априорной известности с заданной точностью координат ЦППО (х,у)цппо и ППП (х,у), k=0, 1, 2, …, K. Однако в условиях развертывания РДСМО на стадионах, спортивно-концертных залах и др. это условие во многих случаях не выполняется. Сигналы глобальной навигационной спутниковой системы (ГНСС) GPS недоступны. В результате аналог теряет свою работоспособность.
Известен способ определения местоположения объекта навигации (см. Патент РФ №2676862, МПК G01C 21/20. опубл. 11.01.2019, бюл. №2), реализующий локальную систему определения координат. Способ обеспечивает определение координат объектов в заданном ограниченном районе в условиях отсутствия доступности сигналов ГНСС. Однако аналог не обеспечивает измерение координат ИРИ, а его использование для определения местоположения ЦППО и ППП экономически нецелесообразно в условиях развертывания на стадионах, помещениях спортивно-концертных залов.
Наиболее близким по технической сущности к заявляемому техническому решению является разностно-дальномерный способ определения координат источника радиоизлучения (см. Патент РФ №2719770, МПК G01S 3/46, опубл. 23.04.2020, бюл. №13).
k
k
i
j
i
j
i
j
k
i,j,k
k
i
j
i,j,k
k
k
k
k
k
k
k
i,j,k
k
i,j
k
i
j
а
a
a
-1
-1
Способ-прототип заключается в приеме сигналов ИРИ пространственно разнесенными ППП с известными координатами, связанными с ЦППО командными линиями связи и линиями аналоговой ретрансляции сигнала, причем командными линиями связи с ЦППО на ППП передают команды настройки на частоту сигнала ИРИ, а по линиям аналоговой ретрансляции принятые в ППП сигналы ИРИ передают на ЦППО, где измеряют разность времени приема этих сигналов в ППП и ЦППО τ, где k ∈ 1, 2, …, K - номера ППП и измерительных баз ЦППО - ППП, при этом на подготовительном этапе на ЦППО в прямоугольном виде задают контролируемый район (КР), границы которого по осям 0х и 0у обозначают через А и В соответственно, КР делят на элементарные участки со сторонами и b по осям 0х и 0у соответственно, площадь которых выбирают исходя из заданной точности Δd определения координат ИРИ, определяют координаты центров элементарных участков (Χ, Y), на основе которых формируют матрицу координат N размерности А' × В', где Χ = i, i=1, 2, …, Α/=А', Y=bj, j=1, 2, …, B/b=В', для каждого элемента (X, Y) матрицы координат N всех K измерительных баз определяют эталонные значения разности времени прихода сигнала на ЦППО и ППП τ, формируют K эталонных матриц M, k=1, 2, …, K, элементами каждой из которых является соответствующее координатам (X, Y) эталонное значение разности задержки прихода сигналов τ для k-й измерительной базы ЦППО - ППП, а в процессе работы для учета разности времени приема сигнала вычисляют взаимнокорреляционную функцию (ВКФ) сигналов R[τ], принятых на k-м ППП, k=1, 2, …, K, и на ЦППО, для значений временного сдвига в интервале τ ∈ (-cd, cd), где d - расстояние между ЦППО и k-ы ППП, с заданным шагом Δτ, соответствующим частоте дискретизации сигнала, на основе K эталонных матриц M формируют корреляционные матрицы Ф путем замены элементов задержки τ на соответствующие им измеренные значения ВКФ R[τ], значения матриц Ф, k=1, 2, …, K, суммируют по всем K измерительным базам  а за наиболее вероятное расположение ИРИ принимают координаты точки (Х, Y), соответствующей максимальному значению элемента матрицы
Прототип обеспечивает повышение точности определения координат заданного источника радиоизлучения в условиях работы на одной частоте нескольких ИРИ с временным или кодовым разделением каналов.
Однако способу-прототипу присущ недостаток, ограничивающий его применение. Прототип и известные аналоги сохраняют свою работоспособность при наличии достоверных данных о координатах всех элементов РДСМО. Местоположение ЦППО и ППП может быть определено с использованием сигналов ГНСС или в ручном режиме. В условиях недоступности сигналов ГНСС лишь второй подход сохраняет работоспособность.
Ручной режим предполагает определение координат элементов РДСМО путем измерения удаления и угла отклонения ЦППО и ППП от объектов с известными координатами. Данному подходу присущи существенные временные издержки и зависимость качества выполняемых операций от степени подготовленности операторов. Метод находит применение при стационарном размещении РДСМО на объекте. Для оперативного развертывания и уточнения позиций ППП при необходимости их перемещения в ходе выполняемых измерений этот подход теряет свою работоспособность.
Целью заявляемого технического решения является разработка разностно-дальномерного способа (РДС) определения координат ИРИ, обеспечивающего повышение оперативности при развертывании измерителя и его перемещениях в условиях отсутствия доступности сигналов ГНСС.
k
k
i
j
i
j
i
j
k
i,j,k
k
i
j
i,j,k
k
k
k
k
k
k
k
i,j,k
k
i,j
k
i
j
1
1
1
1
ЦППО,k
k
k
l,k
1
1
k
1
k
1
1
k
1
k
ЦППО
1
1
а
a
a
-1
-1
Поставленная цель достигается тем, что в известном РДС определения координат ИРИ, основанном на приеме сигналов пространственно разнесенными периферийными пунктами приема, связанными с центральным пунктом приема и обработки командными линиями связи и линиями аналоговой ретрансляции сигнала, причем командными линиями связи с ЦППО на ППП передают команды настройки на частоту сигнала ИРИ, а по линиям аналоговой ретрансляции принятые в ППП сигналы ИРИ передают на ЦППО, где измеряют разность времени приема этих сигналов в ППП и ЦППО τ, k ∈ 1, 2, …, K - номера ППП и измерительных баз ЦППО - ППП, при этом на подготовительном этапе на ЦППО в прямоугольном виде задают контролируемый район (КР), границы которого по осям 0х и 0у обозначают через А и В соответственно, КР делят на элементарные участки со сторонами  и b по осям 0х и 0у соответственно, площадь которых выбирают исходя из заданной точности Δd определения координат ИРИ, определяют координаты центров элементарных участков (Χ, Y), на основе которых формируют матрицу координат N размерности А' × В', где X=i, i=1, 2, …, Α/=А', Y=bj, j=1, 2, …, B/b=B', для каждого элемента (X, Y) матрицы координат N всех K измерительных баз определяют эталонные значения разности времени прихода сигнала на ЦППО и ППП τ, формируют K эталонных матриц M, k=1, 2, …, K, элементами каждой из которых является соответствующее координатам (X, Y) эталонное значение разности задержки прихода сигналов τ для k-й измерительной базы ЦППО - ППП, а в процессе работы для учета разности времени приема сигнала вычисляют взаимно-корреляционную функцию сигналов R[τ], принятых на k-м ППП, k=1, 2, …, K, и на ЦППО, для значений временного сдвига в интервале τ ∈ (-cd, cd), где d - расстояние между ЦППО и k-м ППП, с заданным шагом Δτ, соответствующим частоте дискретизации сигнала, на основе К эталонных матриц М формируют корреляционные матрицы Ф путем замены элементов задержки τ на соответствующие им измеренные значения ВКФ R[τ], значения матриц Ф, k=1, 2, …, K, суммируют по всем K измерительным базам  а за наиболее вероятное расположение ИРИ принимают координаты точки (X, Y), соответствующей максимальному значению элемента матрицы  на подготовительном этапе в условиях отсутствия доступности сигналов ГНСС определяют координаты местоположения лишь ЦППО и опорного ППП в ручном режиме посредством измерения удаления и угла отклонения ЦППО и ППП от точки или объекта с известными координатами, а в качестве ППП используют один из K ППП, в состав ЦППО и ППП дополнительно вводят по передатчику, которые предназначены для излучения контрольных сигналов, определяют расстояние d от ЦППО до всех K ППП путем излучения ЦППО контрольных сигналов в направлении каждого ППП и приеме от них ретранслированного сигнала с последующим определением их задержки Т,  где с - скорость света, t - внутренняя задержка сигнала в k-м ППП, на ЦППО определяют расстояние d между опорным ППП и остальными K - 1 ППП путем измерения задержки  на трассе ППП - ППП, k=2, 3, …, K, на основе излучения передатчиком ППП контрольного сигнала в направлении ППП, приеме на ППП от него ретранслированного сигнала с последующей передачей его на ЦППО  на основе теоремы косинусов определяют углы в формируемых треугольниках ЦППО - ППП - ППП на основе расстояний между ЦППО, опорным ΠΠΠ и искомым ППП, определяют координаты k-то ППП (х, у) на основе полученных расстояний и углов с использованием координат ЦППО (х, у) и опорного ΠΠΠ (x, y).
Благодаря новой совокупности существенных признаков за счет реализации с минимальными издержками локальной навигационной системы в заявляемом способе обеспечивается повышение оперативности при развертывании РДСМО и ее перемещениях в условиях отсутствия доступности сигналов ГНСС.
Заявляемый способ поясняется чертежами, на которых:
на фиг. 1 - показан обобщенный алгоритм РДСМО на подготовительном этапе работы в соответствии с заявляемым изобретением;
на фиг. 2 - приведен обобщенный алгоритм РДСМО в процессе измерения координат ИРИ в соответствии с заявляемым изобретением;
k
на фиг. 3 - приведен алгоритм поиска максимального значения взаимнокор-реляционной функции R[τ];
k
1
на фиг. 4 - иллюстрируется рисунок, поясняющий порядок определения угла ω ППП - ЦППО - ППП в соответствующем треугольнике;
на фиг. 5 - порядок задания контролируемого района:
а) задание границ контролируемого района;
б) формирование элементарных участков;
в) определение координат центров элементарных участков;
на фиг. 6 - порядок формирования матрицы координат N;
на фиг. 7 - порядок формирования первой из K:
k
а) эталонной матрицы M;
k
б) корреляционной матрицы Ф;
на фиг. 8 - приведена обобщенная структурная схема устройства, реализующая данный способ;
на фиг. 9 - приведено измеренное значение ВКФ для одной измерительной
базы;
на фиг. 10 - представлена визуализация суммы проекций ВКФ трех измерительных баз;
1
на фиг. 11 - приведена зависимость среднеквадратического отклонения погрешности (СКО) при измерении координат ППП от погрешности в определении ме-таположения ЦППО и ППП;
на фиг. 12 - показана зависимость СКО погрешности определения координат ППП от погрешности измерения момента прихода сигнала;
1
2
3
на фиг. 13 - вариант пространственного размещения ППП (А-ЦППО, И - ППП, С - ППП, D - ППП);
на фиг. 14 - приведена зависимость СКО погрешности измерения координат ИРИ от погрешности в определении местоположения ППП.
Реализация заявленного способа поясняется следующим образом. Подготовительный этап (см. фиг. 1) предполагает задание границ контролируемого района, выполнение анализа по оптимизации размещения элементов РДСМО (ЦППО и K ППП) для достижения максимальной точности определения координат контролируемых ИРИ.
Известно, что точностные характеристики РДС в основном зависят от: точности измерения временных параметров на ЦППО;
геометрии размещения элементов РДСМО;
точности определения координат элементов РДСМО (ЦППО и K ППП).
Реализация последнего из названных требований связана с наибольшими трудностями в условиях размещения КР в помещениях спортивно-концертных залов, стадионах и др., в условиях отсутствия доступности сигналов ГНСС. Известные способы определения координат (см. Патенты РФ №№2704793, 2719770, 2717231, 2740640) предполагают наличие достоверной информации о местоположении ЦППО и всех ППП.
1
В названных условиях используют ручной режим определения координат элементов РДСМО. Последний предполагает измерение удаления и угла отклонения ЦППО и опорного ППП от точки или объекта с известными координатами (угол или окно задания и т.д.). Данный процесс трудоемкий, предполагает определенный уровень подготовленности оператора и главное - существенные временные затраты. Последние составляют 10-15 минут на один элемент РДСМО при выполнении работ подготовленным специалистом. В связи с этим оперативное развертывание оборудования РДСМО или перенос его элементов в процессе работы как правило невозможен.
Для описания пространственного размещения всех элементов РДСМО (получения координат ППП) необходимо знание:
ЦППО
1
1
координат размещения двух ее элементов ЦППО (х, у) и опорного ППП (х, y);
1
взаимных расстояний от ЦППО до всех ППП и от опорного ППП до всех остальных ППП;
1
углов в треугольниках между ЦППО, опорным ППП и искомым ППП;
k
k
внутренних задержек сигналов ППП t.
1
В предлагаемом способе на подготовительном этапе (на фиг. 1) в ручном или ином доступном режиме определяют координаты двух точек, предназначенных для развертывания ЦППО и опорного ППП.
ЦППО
1
1
Определяют тип локальной системы координат (СК), в которой будет выполняться определение координат ИРИ. Выбор СК зависит от конкретных условий ее применения. Для названных выше ограничений предпочтительно использование метровых координат с их началом в левой нижней зоне контролируемого района. В качестве последнего может быть использован план строения или местности. Зададим позиции ЦППО (х, у) и опорного ППП в этой системе координат. Далее для упрощения расчетов определим новую локальную СК с центром в точке развертывания ЦППО и определим позицию опорного ППП
В качестве опорного ППП может быть назначен любой из их совокупности K, k=1, 2, …, K.
k
Для каждой полученной измерительной базы, где источником контрольного сигнала является ЦППО, определяют значение временного сдвига Т (см. фиг. 3)
а
k
где τ - аналитически определенный отсчет с максимальным значением ВКФ R,
D
D
D
- количество отсчетов ВКФ, T=1/F - период дискретизации, F - частота дискретизации.
k
Цифровая обработка сигналов приводит к снижению точностных характеристик. Причина состоит в том, что нет объективных оснований для того, чтобы время задержки между парой сигналов было кратно периоду дискретизации или сколь либо близко к кратному значению. В этом случае положение во времени наибольшего по абсолютному значению отсчета дискретной ВКФ не является точной оценкой времени задержки Т.
k
2
Для повышения точности оценки времени задержки используют метод, основанный ца интерполяции абсолютных значений ВКФ R в окрестностях ее максимума двух соседних отсчетов. Для этой цели подходит полином второго порядка A(τ)
a
k
2
1
1
3
3
2
Для нахождения τ определяют абсолютное значение максимального по модулю отсчета ВКФ Е(х, y) и определяют две соседние ближайшие точки, располагающиеся слева (x, y) и справа (x, y). Далее определяют коэффициенты параболы и находят ее максимум:
2
k
1
2
D
3
2
D
2
k
2
1
k
1
3
k
3
х=maxE[τ], х=х - T, х=х+T, у=Е[х], у=Е[х], y=Е[х],
k
После измерения временного сдвига для каждой базы определяют удаление ППП от ЦППО
где с - скорость света.
1
k
1
1
k
k
1
1,k
На следующем этапе источником контрольного сигнала является опорный ППП. С его помощью (аналогично сигналам ЦППО) определяют значения временного сдвига  (см. фиг. 1). Для этого принятый от ППП ретранслированный сигнал ППП далее поступает на ЦППО для измерения временного сдвига  на трассе ППП - ППП. После нахождения  для каждой измерительной базы определяют удаление ППП, k=1, 2, …, K, от опорного ППП d
k
Зная удаление ППП от ЦППО становится возможным определение его координат
а
k
k
где  - угол между вектором ЦППО - ППП и осью абсцисс,  - позиция ППП в системе координат с центром в точке расположения ЦППО.
Х
1
а
В связи с тем, что координаты опорного ППП известны (x, y), угол  определяют как сумму углов (см. фиг. 4)
а
=ϕ+ω,
1
k
1
где ϕ - угол между вектором ЦППО – ППП и осью абсцисс; ω - угол ППП - ЦППО – ППП в соответствующем треугольнике.
Угол ϕ находят в соответствии со следующим выражением
где
k
Для каждого ППП угол ω определяют на основе теоремы косинусов следующим образом
1
В связи с тем, что существует два возможных треугольника, симметричных относительно стороны ЦППО – ΠΠΠ. Оператору необходимо определить расположение ППП - по часовой стрелке или против. Если ППП располагают по часовой стрелке, то
После определения координат всех элементов РДСМО на ЦППО приступают к делению контролируемого района на элементарные участки по аналогии со способом-прототипом. Для этого границы КР по осям 0х и 0у обозначают соответственно через А и В (см. фиг. 5 а, б и в).
k
k
Далее определяют позицию ППП (х, у) в системе координат КР, заданного на подготовительном этапе. Для этого выполняем следующее
а
a
a
а
i
j
i
j
i
j
k
i,j,k
Контролируемый район делят на элементарные участки со сторонами  и b по осям 0х и 0у соответственно, площадь которых выбирают исходя из заданной точности Δd определения координат ИРИ. На практике, как правило, =b. Определяют координаты центров элементарных участков (X, Y), на основе которых формируют матрицу координат N размерности A'× В', где Χ=i, i=1, 2, …, А/=Α', Y=bj, j=1, 2, …, B/b=В' (см. фиг. 6). Для каждого элемента (Х, Y) матрицы координат N всех K измерительных баз определяют эталонные разности времени прихода сигнала на ЦППО и ППП τ в соответствии с выражением
0
0
k
k
где с - скорость света, [Х, Y) - координаты ЦППО, (Х, Y) - координаты k-го ППП, k - 1, 2, 3 - номера ППП.
k
i
j
i,j,k
k
н
Далее формируют K эталонных матриц M, k=1, 2, …, K, элементами каждой из которых является соответствующее координатам (X, Y) эталонное значение разности задержке прихода сигнала τ для k-й измерительной базы ЦППО - ППП (см. фиг. 7). Далее задают МАС-адрес пользователя  длительность интервала анализа Δt и его начало t', частоту сигнала ƒ.
Обобщенная структурная схема устройства, реализующего данный способ, приведена на фиг. 8. Здесь:
1.1-1.3 - периферийные пункты приема (ППП);
2 - центральный пункт приема и обработки (ЦППО);
3 - центральный пост обработки (ЦПО);
4.1-4.3 - тракты приема радиосигналов ППП;
5.1-5.3 - тракты ретрансляции радиосигналов ППП;
6 - первая входная шина ЦППО;
7 - радиопередатчик командной линии связи;
8.0-8.3 - тракты приема радиосигналов ЦППО;
9.0-9.3 - аналого-цифровые преобразователи;
10 - вторая входная шина ЦППО;
11 - блок управления;
12 - тракт анализа;
13.0-13.3 - блоки памяти;
14 - тактовый генератор;
15 - вычислитель;
k
16 - блок формирования корреляционных матриц Ф;
k
17 - блок формирования эталонных матриц М;
18 - блок формирования матрицы координат N;
19 - сумматор;
20 - блок принятия решения;
21 - третья входная шина ЦППО;
22 - четвертая входная шина ЦППО;
23 - пятая входная шина ЦППО;
24 - выходная шина ЦППО;
25.0-25.1 - передатчики контрольных каналов.
Таким образом устройства содержит три ППП (k ∈ 1, 2 и 3) и один ЦППО (k ∈ 0).
k
k
Каждый ППП сигналов ИРИ с временным или кодовым разделением в общем виде представляет совокупность последовательно соединенных трактов приема сигналов ИРИ названных классов и трактов, реализующих аналоговую ретрансляцию излучений. ЦППО содержит радиопередатчик, предназначенный для реализации командной линии связи, совокупность из K+1 трактов приема радиосигналов (из них K трактов, предназначенных для приема ретранслированных ППП сигналов) и центральный пункт обработки в составе: блок управления, тракт анализа, K+1 аналогово-цифровых преобразователей (АЦП), K+1 блоков памяти, вычислитель, блок формирования матриц эталонных значений M, блок формирования матрицы координат N, блок формирования корреляционных матриц Ф, сумматор и блок принятия решения. При этом все ППП настраиваются на заданную частоту сигнала, значение которой поступает по командным каналам связи от ЦППО.
Сигналы ИРИ, принятые на ЦППО и ППП имеют вид
соответственно.
Последние с выходов приемных трактов ППП поступают на входы соответствующих трактов аналоговой ретрансляции и далее излучают в эфир.
На ЦППО эти излучения принимают с помощью соответствующих K приемных трактов и запоминаются на заданное время Δt. Одновременно с ними осуществляют прием сигналов этого ИРИ и на собственный K+1-вый приемный тракт в течение того же интервала времени Δt, в течение которого осуществляют их анализ.
0
0
a
Для этого принятые в K+1-м приемном тракте ЦППО сигналы с помощью АЦП и демодулятора (снимают DBPSK модуляцию), превращая аналоговый сигнал x(t) в битовую последовательность (t) в тракте анализа. На его основе осуществляют поиск преамбулы и заголовка PLCP. В преамбуле в поле SFD содержится значение 0xF3A0, которое служит индикатором начала фрейма. В заголовке PLCP в поле Signal содержится информация о способе цифровой модуляции самих данных (МАС-фрейма). Значение 0х0А соответствуют DBPSK модуляции, a 0x14-DQPSK модуляции. Кроме того, в заголовке PLCP содержится информация о длине данных в МАС-фрейме.
Известны МАС-фреймы различных типов: CTS, PS-Poll, CF-End, CF-End+CF-Ack, Block AckReq. Помимо перечисленных, существуют и другие типы МАС-фреймов, однако в них отсутствует информация о МАС-адресах передающего устройства. Поэтому их использование затруднено для местоопределения ИРИ.
У всех фреймов в их начале есть поле Frame Control, где содержится информация о типе данного фрейма - это поля Туре и Subtype. В поле ТА у фреймов содержится МАС-адрес передающего устройства. Для получения названных данных необходимо предварительно выполнить цифровую демодуляцию МАС-фрейма в соответствии с содержимом поля Signal заголовке PLCP.
1
k
k
В случае совпадения МАС-адреса Ρ обнаруженного фрейма с заданным адресом Pi устройства, участок сигнала с границами от начала преамбулы PLCP t*до конца МАС-фрейма t используют для нахождения взаимнокорреляционной функции принятого сигнала R[τ]. При этом свертка принятых ЦППО и ППП сигналов R может осуществляться как в аналоговом, так и в дискретном виде. В случае использования дискретных сверток (предлагаемый вариант устройства) запомненные сигналы от ППП дополнительно оцифровывают.
k
k
Дополнительная синхронизация при измерении τ не требуется в связи с тем, что обработку всех K+1 сигналов осуществляют в одном месте. Нахождение задержек τ с последующим определением координаты ИРИ выполняют в соответствии со способом-прототипом.
k
k
k
0
-1
-1
Вычисление ВКФ сигналов R[τ], принятой k-й измерительной базой, k=1, 2, …, K, осуществляют для значений временного сдвига τ ∈ (-cd; cd], где  - расстояние между ЦППО и k-м ППП с заданным шагом τ по формуле:
0
k
где τ - временной сдвиг между сигналами (задержка), u(t),  - сигнал, принятый на ЦППО и сопряженный с ним сигнал, u(t),  - сигнал, принятый на k-м ППП и сопряженный с ним сигнал. На фиг. 9 приведено значение ВКФ для одной измерительной базы.
k
k
k
i,j,k
i
j
i,j,k
i
j)
k
i,j
k
После этого полученные значения каждой рассчитанной ВКФ R[τ] проецируют на элементы соответствующей k-й эталонной матрицы M (элементы зоны поиска). В матрице M каждый элемент τ соответствует определенной точке (X, Y). По ранее полученным значениям τ для каждой точки (Х, Y контролируемого района рассчитывают значения R[τ] в соответствии с (16) и, следовательно, соответствующие этой точке. Названная операция соответствует формированию соответствующей (из набора в K) корреляционной матрицы Ф.
k
k
i,j
k
i
j
Формирование K корреляционных матриц Ф осуществляют путем проецирования полученных значений ВКФ R[τ] в соответствующих местах эталонной матрицы M, а следовательно, в соответствующих точках (Х, Y).
k
k
i
j
На следующем этапе на зону поиска (матрицы Ф) значения R[τ] суммируют по всем K измерительным базам:  За наиболее вероятное местоположение заданного ИРИ принимают координаты точки (X, Y), соответствующее максимальному значению элемента матрицы
i,j
i
j
k
Такой подход к ОК заданного ИРИ обеспечивает повышение точности выполняемых измерений. Известно, что во взаимнокорреляционных функциях R[τ] присутствуют побочные максимумы от отраженных сигналов, иногда с максимальным значением. Последние имеют местоположение на временной оси не соответствующее положению ИРИ. В известных способах и устройствах это приводит к ошибкам измерений координат. В предложенном способе значение ВКФ R[τ] принимает максимальное значение (глобальный максимум) благодаря тому, что в точке (X, Y) оно равно сумме величин корреляций, соответствующих местоположению истинным локальным максимумам R[τ] всех трех измерительных баз. В свою очередь локальные максимумы, полученные от переотраженных сигналов всех трех измерительных баз, будут распределены случайным образом. Их суммирование не дает глобального максимума ВКФ. На фиг. 10 представлена визуализация суммы проекций ВКФ на контролируемый район. Наблюдается пересечение трех линий, образованных истинными пиками ВКФ. Побочные пики не внесли ошибки при определении местоположения ИРИ.
Таким образом имеет место следующая последовательность действий на подготовительном этапе:
1. Задание контролируемого района в прямоугольном виде, в котором будет осуществляться поиск и местоопределение заданного ИРИ. Границы по 0х и 0у обозначают через А и В соответственно (см. фиг. 5а).
2. Выбор позиций для ЦППО и K ППП, обеспечивающих максимальную точность определения координат ИРИ.
ЦППО
1
1
3. Определение координат ЦППО (х, у) и ППП [х, у) одним из доступных методов.
ЦППО,k
k
4. Определение расстояния d от ЦППО до всех K ППП путем излучения ЦППО контрольных сигналов в направлении каждого из ППП и приема от них ретранслированного сигнала с последующим определением их задержки Т (см. выражение 7).
1
1,k
5. Определение расстояния между опорными ППП и остальными K - 1 ППП d в соответствии с (8).
1
k
1
6. На основе теоремы косинусов определение углов в формируемых треугольниках ЦППО - ППП - ППП на основе расстояний между ЦППО, опорным ППП и k-м искомым ППП (выражения 10, 11 и 12).
k
k
ЦППО,k
k
7. Определение координат k-го ППП (х, у) на основе полученных данных о расстоянии d и угле между вектором ЦППО - ППП и осью абсцисс.
а
8. КР делят на элементарные участки со сторонами  и b по осям 0х и 0у соответственно, площадь которых выбирают исходя из заданной точности Δd определения координат ИРИ (см. фиг. 5б).
i
j
9. Определение координат центров элементарных участков (X, Y) (см. фиг. 5в).
i
j
i
i
a
а
10. Формирование матрицы координат N размерности А' × В', на основе координат (X, Y), где X=i, i=1, 2, …, А/=A', Y=bj, j=1, 2, …, B/b=B' (см. фиг. 6).
i
j
k
i,j,k
11. Определение для каждого элемента [X, Y) матрицы N и всех K измерительных баз эталонных значений разности времени прихода сигнала на ЦППО и ППП τ.
k
i
j
i,j,k
k
12. Формирование K эталонных матриц M, k=1, 2, …, K, элементами каждой из которых является соответствующее координатам (X, Y) эталонное значение разности задержки прихода сигнала τ для k-й измерительной базы ЦППО - ППП (см. фиг. 7).
В процессе работы по аналогии с прототипом осуществляют:
1. Накопление в ЦППО сигналов за некоторый промежуток времени Δt от всех ППП и ЦППО.
2. Демодуляцию накопленного сигнала, принятого собственно ЦППО и модулированного посредством DBPSK модуляции. В результате получают массив битов.
2
3. Поиск в массиве битов SFD кода - значения 0xF3A0 (1111001110100000).
4. Чтение заголовка PLCP- следующие 48 бит после SFD кода.
2
2
5. Определение способа модуляции МАС-фрейма: в поле Signal заголовка PLCP: значение 0х0А (00001010) соответствуют DBPSK модуляции, 0x14 (00010100) - DQPSK - модуляции.
6. Определение длины МАС-фрейма - поле Length заголовка PLCP.
7. Демодуляцию МАС-фрейма.
8. Определение МАС-адреса устройства.
9. Если МАС-адрес устройства соответствует МАС-адресу искомого устройства, то выделяют часть сигнала, соответствующего началу преамбулы PLCP и окончанию МАС-фрейма, в противном случае начиная с конца текущего фрейма, осуществляется поиск следующего фрейма (пункт 3).
k
k
k
k
-1
-1
10. Вычисление взаимнокорреляционной функции сигналов R[τ], принятых на k-м ППП, k=1, 2, …, K, и на ЦППО для значений временного сдвига в интервале τ ∈ (-cd, cd), где d - расстояние между ЦППО и k-м ППП, с заданным шагом Δτ.
k
i,j,k
k
k
i,j
11. Формирование К корреляционных матриц Ф путем замены элементов задержки τ матриц M на соответствующие им измеренные значения ВКФ R[τ].
k
12. Суммирование матриц Ф по всем K измерительным базам (см. фиг. 10)
i
j
13. Принятие за наиболее вероятное местоположение заданного ИРИ точки с координатами (X,Y), соответствующей максимальному значению элемента  матрицы
Выполнена оценка эффективности предлагаемого способа определения координат ИРИ. Определим суммарные временные затраты на развертывание измерителя в соответствии с предлагаемым изобретением и прототипа в составе ЦППО и трех ППП. Практические испытания показали, что для ручного режима определения координат ЦППО и каждого ППП требуется примерно по 15 минут.На определение координат второго и третьего ППП в предлагаемом способе дополнительно требуется примерно 15+15=30 миллисекунд. Временем на анализ КР (деление, определение координат центров элементарных участков) можно пренебречь. Оно составляет десятки миллисекунд. В результате имеем: прототип - 60 минут, предлагаемый способ - 30 минут.Тогда выигрыш по скорости развертывания (оперативность) η для названных условий составляет
Следует отметить, что при увеличении количества ППП значение ηη линейно возрастает.
Оценка точностных характеристик предлагаемого способа определения координат ИРИ выполнена на основе математического моделирования в среде Visual Studio 2019 с использованием языка программирования С#. Для отображение результатов моделирования использована среда MATLAB. Моделирование выполнено на компьютере со следующими характеристиками:
процессор - Intel Core™i7-7500UCHU@ 2.70 GHz;
оперативная память - 32 ГБ;
видеоадаптер - Intel HD Graphics 620.
Погрешность оценивания определилась как второй начальный момент распределения модуля вектора ошибок
2
2
где δx, δy - среднее значение квадратов соответствующих отклонений от истинных координат в произвольно выбранной системе координат.
Точность определения координат ИРИ находится в зависимости от величины погрешности в измерении координат ЦППО и ППП, взаимного местоположения ЦППО и ППП, их пространственного положения относительно КР (анализ последнего выполнен в рамках Пат. РФ №2717231, опубл. 19.03.2020). В свою очередь точность определения координат ППП в предлагаемом способе находится в зависимости от погрешностей при измерении временных параметров на основе ВКФ.
1
1
1
1
В выполненном анализе погрешность определения координат ЦППО и ППП задана в пределах от 0 до 50 нс. Погрешность определения моментов прихода сигнала на ЦППО и ΠΠΠ задана в пределах от 0 до 100 нс. Зависимость среднеквадратической ошибки (СКО) оценки координат от погрешности координат ЦППО и опорного ΠΠΠ получены для четырех значений погрешности определения моментов прихода сигнала: 0 не, 25 нс, 50 нс и 100 нс. Для каждой погрешности определения момента прихода сигнала и погрешности в координатах ЦППО и ППП выполнено 10 тысяч испытаний. Соответствующие зависимости приведены на фиг. 11 и 12.
Выполнен анализ зависимости точности оценки координат ИРИ от погрешности определения местоположения ППП, результаты которого приведены на фиг. 14. Ошибки в определении координат ППП заданы в пределах от 0 до 20 м. Для каждой погрешности было проведено 10000 испытаний. Исходная расстановка ППП при моделировании приведена на фиг. 13.
1
Из фиг. 11-14 следует, что СКО оценки координат в 5 метров можно обеспечить при погрешности определения координат ЦППО и опорного ППП менее 5 метров и погрешности измерения временных параметров сигнала менее 25 нс.
1
При погрешностях определения координат ЦППО и ППП 1 метр и определения момента прихода сигнала 25 нс. (реализация этих требований подтверждена практическими испытаниями) значение СКО оценки координат ИРИ составляет 2-3 метра.
1
Результаты моделирования и практические испытания позволяют сделать вывод о том, что предлагаемый способ определения координат ИРИ при погрешностях определения координат ЦППО и ΠΠΠ менее 22 м и измерении временных параметров сигнала с погрешностью менее 25 нc дает погрешность в определении координат ППП менее 18 м. Данная величина меньше погрешности, получаемой при использовании сигналов ГНСС GPS (см. Электронный ресурс https://ru.wikipedia.org/wiki/GPS). Повышение точности определения координат ППП в предлагаемом способе влечет за собой повышение точности определения местоположения ИРИ.
Таким образом, предлагаемый способ с минимальными издержками (два передатчика контрольных сигналов) одновременно позволяет определять местоположение объектов навигации в виде ЦППО и K ППП (реализуют локальную систему навигации) и определение координат заданных ИРИ в условиях отсутствия доступности сигналов ГНСС. При этом получаемый выигрыш η по оперативности (скорости развертывания и перемещениям элементов измерителя) зависит от количества ППП и составляет η≥2.