Результат интеллектуальной деятельности: МУЛЬТИПЛИКАТИВНЫЙ РАЗНОСТНО-ОТНОСИТЕЛЬНЫЙ СПОСОБ ОПРЕДЕЛЕНИЯ КООРДИНАТ МЕСТОПОЛОЖЕНИЯ ИСТОЧНИКА ИМПУЛЬСНОГО РАДИОИЗЛУЧЕНИЯ

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

Известны способы определения координат ИРИ, в которых используются пассивные пеленгаторы в количестве не менее трех, центр тяжести области пересечения выявленных азимутов которых на фронт прихода волны принимается за оценку местоположения. Основными принципами работы таких пеленгаторов являются амплитудные, фазовые и интерферометрические [1, 2]. Широко применяемым является амплитудный способ пеленгования, при котором используется антенная система, имеющая диаграмму направленности с ярко выраженным максимумом главного лепестка и минимальными задним и боковыми лепестками. К таким антенным системам относятся, например, логопериодические или антенны, имеющие кардиоидную характеристику и др. При амплитудном способе механическим вращением добиваются положения антенны, при котором выходной сигнал имеет максимальную величину. Такое направление принимают за направление на ИРИ. К недостаткам большинства пеленгаторов следует отнести высокую степень сложности антенных систем, коммутационных устройств и наличие многоканальных радиоприемников, а также необходимость в быстродействующих системах обработки информации.

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

Известен способ [3], в котором для определения координат местоположения ИРИ используют N, не менее четырех, стационарных радиоконтрольных постов, расположенных не на одной прямой, один из которых принимают за базовый, соединяя с остальными N-1 постами линиями связи, на всех постах осуществляют квазисинхронное сканирование по заданным фиксированным частотам настройки, усредняют измеренные значения уровней сигналов на каждой из сканируемых частот, а затем на базовом посту для каждого из сочетаний (сочетаний из N по 4) на основании обратно пропорциональной зависимости отношений расстояний от поста до источника радиоизлучения и соответствующих им разностей уровней сигналов, выраженных в дБ, составляют три уравнения, каждое из которых описывает окружность равных отношений, по параметрам двух любых пар которых и определяют текущее среднее значение широты и долготы местоположения источника радиоизлучения. Недостатком этого способа является большое количество стационарных постов радиоконтроля.

Известны способы и устройства пеленгования [4, 5], которые могут быть использованы для целей определения координат.

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

К недостаткам способа следует отнести:

1) Большое количество антенн.

2) Способ не ориентирован на использование радиоконтрольных постов.

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

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

Данному пеленгатору свойственны следующие недостатки:

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

2) Большое количество специализированных пеленгационных (но не радиоконтрольных) постов.

3) Необоснованное и не раскрытое (хотя бы до функциональной схемы) применение системы единого времени на ЦП и хронизаторов на ПП, синхронизированных с системой единого времени.

4) Необходимость в наличии радиоканалов с большой пропускной способностью (до 625 Мбод) для передачи даже фрагментов сигналов с ПП1 и ПП2 на ЦП.

5) Для организации радиоканала необходимы радиопередающие устройства и получение разрешения на их работу в определенных условиях эксплуатации.

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

Способ, основанный на приеме сигналов ИРИ четырьмя антеннами, образующими три независимые измерительные базы, в разнесенных точках A, B, C, D таким образом, что объем фигуры, образованный из этих точек, больше нуля (V_{A, B, C, D}>0). Сигнал одновременно принимается на все антенны, измеряют три независимые разности времени Δt_AC, Δt_BC, Δt_DC приема сигнала парами антенн, образующих измерительные антенные базы (АС), (ВС) и (DC). По измеренным разностям времен вычисляют разности дальностей от ИРИ до пар точек (А, С), (В, С), (D, C), для k-й тройки антенн, расположенных в точках А, В, С при k=1, B, C, D при k=2, D, C, A при k=3, вычисляют с помощью измеренных разностей дальностей значения угла γ_k, характеризующие угловое положение плоскости положения ИРИ Ω_k, k=1, 2, 3 относительно соответствующей измерительной базы, и координаты точки F_k, принадлежащей k-й плоскости положения ИРИ, вычисляют искомые координаты ИРИ как координаты точки пересечения трех плоскостей положения ИРИ Ω_k, k=1, 2, 3, каждая из которых характеризуется координатами точек расположения k-й тройки антенн и вычисленными значениями угла γ_k и координатами точки F_k, отображают результаты вычисления координат ИРИ в заданном формате.

Этот способ ближе к заявляемому, но также обладает рядом существенных недостатков:

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

2) Необходимость сведения сигналов ИРИ с разнесенных на оптимальное расстояние до 0,6-0,7 R ЭМД антенн согласно [2] в одну точку, что реализовывать экономически нецелесообразно.

3) Для измерения разности времени приема сигнала непосредственно с антенн используются двухвходовые измерители.

4) Сложность технической реализации, обусловленная большим количеством различных вычислителей.

5) Неопределенность в построении поверхности положения в виде плоскости, перпендикулярной плоскости расположения антенн, так как антенны в точках A, B, C, D не располагаются в одной плоскости, о чем свидетельствует условие V_{A, B, C, D}>0 в формуле изобретении.

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

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

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

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

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

- измеряют две разности времен Δt_AC и Δt_BC приема сигнала ИРИ антеннами;

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

- вычисляют значения разностей дальностей от ИРИ до пар антенн с использованием выражений Δr_AC=CΔt_AC, Δr_BC=CΔt_BC, Δr_AB=Δr_AC-Δr_BC, где С - скорость распространения электромагнитной волны;

- вычисляют координаты по полученной формуле.

В соответствии с [7] в состав устройства, реализующего способ, входит:

- три антенны;

- два измерителя разности времен;

- два измерителя плотности потока мощности;

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

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

Прототипу свойственны следующие недостатки:

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

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

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

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

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

приводит к появлению ошибки сингулярности (когда знаменатель может быть близок к нулю).

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

Эта цель достигается с помощью признаков, указанных в формуле изобретения:

Мультипликативный разностно-относительный способ определения координат местоположения источника импульсного радиоизлучения, основанный на измерении величин запаздывания сигналов радиоизлучений на назначенных частотах в нескольких точках пространства радиоприемными устройствами, отличающийся тем, что для измерения моментов прихода сигналов используют предварительно откалиброванный мобильный пост радиоконтроля в качестве базового, который перемещают по М≥1 точкам нелинейной траектории, и два ведомых одинаковых стационарных поста радиоконтроля, соединенные линиями связи с базовым, измеренные значения моментов прихода сигналов с которых передают на базовый пост, где вычисляют разность времени прихода сигналов, формируют определитель Кэли-Менгера размерностью 5×5, составляют на его основе полное уравнение четвертой степени относительно одного из неизвестных расстояний от источника до поста, по нему вычисляют расстояния до других постов и отношения расстояний, равнозначные отношениям величин времени распространения сигналов до постов, составляют и мультипликативных функций, представляющих сочетания, взятые по два и по три, из вычисленных парных сочетаний (М+2) разностей отношений расстояний, рассчитанных от точек измерения до местоположения искомого ИРИ по заданным его координатам, и вычисленных парных сочетаний (М+2) отношений соответствующих измеренных величин запаздывания прихода сигналов, расстояния от точек измерения до источника получают при этом для заданных из известного диапазона значений широт и долгот местоположения искомого ИРИ, изменяют дихотомически или методом наискорейшего спуска последовательно каждый из параметров местоположения искомого источника радиоизлучения при неизменном значении другого и находят точки экстремумов и точки перегиба мультипликативных функций с усреднением и последовательной фиксацией каждого искомого параметра местоположения источника в этих точках, как окончательного. Перед проведением измерений и вычислений координат местоположения искомого ИРИ выполняют калибровку измерителей на постах.

Калибровку измерителя момента времени прихода сигналов на посты выполняют, используя эталонные РЭС с известными параметрами сигналов и координатами местоположения. Каждая эталонная РЭС должна находиться в зоне ЭМД всех трех постов. Количество эталонных РЭС и распределение в зоне ЭМД постов должно быть достаточным для обеспечения заданной точности калибровки как по расстоянию, так и азимуту относительно постов. В основе определения координат лежит концептуальный отказ от прменения каких-либо сложных линий положения ИРИ, например парабол, гипербол, окружностей Аполлония Пергского, овалов Кассини, гиперболоидов вращения и других, и применение универсального численного способа последовательного определения параметров местоположения ИРИ. При этом вместо физических величин (расстояний, величин запаздываний прихода сигналов радиоизлучений в точки их приема) предлагается использовать мультипликативные разности отношений этих физических величин, при этом использовать критерий минимума разностей отношений расстояний от местоположения ИРИ до постов радиоконтроля и соответствующих им отношений величин времени распространения сигналов. Координаты могут вычисляться по методу дихотомии, например методу поразрядного уравновешивания. Для его использования априори должны быть известны диапазоны D значений искомых величин. Эти диапазоны обычно известны, исходя из параметров общей зоны электромагнитной доступности используемых трех постов радиоконтроля. В соответствии с алгоритмом поразрядного уравновешивания первоначально задают среднее из диапазона D значение определяемой величины (например, широты) при фиксированной, но лежащей в известных диапазонах значений, долготы. Вычисляют расстояния от i-го местоположения ИРИ до каждого j-го поста, . Затем вычисляют парные отношения этих расстояний . Эти отношения позволяют исключить зависимость вычисления координат местоположения от мощности ИРИ. Полученные отношения сравнивают путем вычитания с измеренными отношениями величин запаздывания прихода сигналов или с соответствующими им расстояниями от ИРИ до постов , где r_a, r_b, r_c - фактические расстояния до постов, указанные на фиг. 5.

Например, для постов А и В эту разность определяют как f_1abi,ab=(n_abi-n_ab). Для В и С - как f_1bci,bc=(n_bci-n_bc) и т.д.

Если разность отношений меньше нуля, то к первоначальному значению широты добавляют 1/4 часть диапазона. В противном случае из первоначального значения широты вычитают 1/4 часть диапазона ее значения. Затем опять производят вычисление расстояний до постов и оценку результатов сравнения, как описано выше. При этом добавляют (или вычитают) уже 1/8 часть диапазона, затем 1/16 часть и т.д. Такие итерации продолжают до тех пор, пока результат сравнения не окажется по модулю меньше заранее заданного значения погрешности дискретизации каждого параметра местоположения , где m - количество итераций.

После этого фиксируют полученное значение параметра. Затем аналогично вычисляют значение долготы при найденной широте, Минимум любой из разностей свидетельствует о нахождении местоположения ИРИ в точке с выбранными координатами. Но так как координаты ИРИ находятся на перпендикуляре к линии баз, то каждая из отдельных разностей f_1abi,ab, f_1bci,bc, f_1cai,ca будет иметь минимальное значение при нахождении ИРИ по обе стороны от баз. Возникает неоднозначность в определении местоположения ИРИ. Неоднозначность снимают путем нахождения экстремумов или точек перегиба мультипликативных функций разностей отношений. Так как в способе рассматриваются две физических величины (см. выше), то берут разности отношений расстояний от источника радиоизлучений до точек их приема и соответствующих величин времени распространения сигналов. Так как измерители величин моментов прихода сигналов расположены согласно способу на всех постах, то будет три таких разности: f_1abi,ab=(n_abi-n_ab), f_1bci,bc=(n_bci-n_bc), f_1cai,ca=(n_cai-n_ca).

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

Примеры мультипликативных функций, включающих две разности отношений: f_2abc= Общее число таких функций равно .

Пример мультипликативных функций, включающих три разности отношений функций: f_3abc=f_1abi,abf_1bci,bcf_1aci,ac. Общее число таких функций равно .

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

По отметкам времени прихода сигналов ИРИ на синхронизированные сигналами эталонных РЭС измерители моментов времени прихода [8] всех трех постов на ведущем посту вычисляются разности времени запаздывания сигналов искомого РЭС . Выражая время запаздывания прихода сигналов T_c, T_b, T_a через соответствующие проходимые сигналами расстояния, получим , или в другом виде: , где С - скорость распространения электромагнитной волны.

Из полученных соотношений выразим расстояния: r_a и r_b через r_c (можно делать и другие замены): r_a=r_c+Δr_ca, r_b=r_c-Δr_bc.

Для вычисления неизвестного расстояния r_c составляют и раскрывают определитель Кэли-Менгера размерностью 5×5, на единицу большей, чем количество вершин описываемого им объема четырехвершинной симплекс-фигуры, представленной на фиг. 5. Поскольку объем проекции этой фигуры на плоскость равен нулю, то определитель Кэли-Менгера в соответствии с [9] представляется в виде:

В этом определителе в соответствии с поясняющей фиг. 5 через а, b и с обозначены базы (расстояния между постами). Раскрывая определитель, получим:

Определитель (1) с учетом введенных соотношений расстояний r_a и r_b будет представлять полное уравнение 4-й степени относительно неизвестного расстояния r_c, имеющее вид:

где:

Уравнение (2) относительно r_c решают численным методом. После этого находят и расстояния r_a и r_b. А затем определяют отношения этих расстояний, составляют и решают мультипликативные функции разностей отношений. В результате решения последних и усреднения получают искомые координаты местоположения ИРИ.

Техническая реализация способа, адекватная [10], приведена на фиг. 4, где показаны три одинаковых радиоконтрольных поста - РКП А, РКП Б и РКПС, содержащие:

1. Всенаправленные антенны 1, 6, 11;

2. Сканирующие радиоприемники (РП) 2, 7, 12;

3. Измерители величин запаздывания сигналов (ИВЗ) 3, 8, 13;

4. Компьютеры 4, 9, 14;

5. Устройства связи 5,10,15.

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

1) Калибруют измеритель величины запаздывания прихода на посты сигналов (ИВЗ), используя массив эталонных РЭС с известными параметрами сигналов и координатами местоположения. Каждая эталонная РЭС должна находиться в зоне ЭМД всех трех постов. Количество эталонных РЭС и распределение в зоне ЭМД постов должно быть достаточным для обеспечения заданной точности калибровки как по расстоянию, так и азимуту относительно постов.

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

3) Сведения, полученные в пп. 1 и 2, пересылают по каналу связи устройства связи из ведомых компьютеров на ведущий.

4) Составляют определитель Кэли-Менгера размерностью 5×5.

5) Раскрывают определитель Кэли-Менгера, получая при этом полное уравнение четвертой степени относительно одного неизвестного.

6) Численным методом определяют по полному уравнению четвертой степени расстояния, а затем величины и отношение величин запаздывания прихода сигналов от ИРИ на основании измеренных измерителями [8] моментов прихода сигналов.

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

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

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

Из таблицы видно, что способ дает увеличение статистики по сравнению с прототипом в 165 раз.

В предлагаемом способе устранены недостатки прототипа:

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

2) Определение координат местоположения ИРИ выполняется радиоконтрольными постами, а не средствами, далекими от их практической реализации.

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

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

Литература

1. Корнеев И.В., Ленцман В.Л. и др. Теория и практика государственного регулирования использования радиочастот и РЭС гражданского применения.

Сборник материалов курсов повышения квалификации специалистов радиочастотных центров федеральных округов. Книга 2. - СПб.: СПбГУТ. 2003.

2. Липатников В.А., Соломатин А.И., Терентьев А.В. Радиопеленгация. Теория и практика. Спб.: ВАС, 2006. - 356 с.

3. Способ определения координат местоположения источников радиоизлучения. Заявка №2009138071, опубл. 20.04.2011 г. Б.И. №11. Авторы: Логинов Ю.И., Екимов О.Б., Рудаков Р.Н.

4. Разностно-дальномерный способ пеленгования источника радиоизлучения. Патент РФ №2325666 С2. Авторы: Сайбель А.Г., Сидоров П.А.

5. Разнесенный разностно-дальномерный пеленгатор. Патент РФ №2382378, С1. Авторы: Ивасенко А.В., Сайбель А.Г., Хохлов П.Ю.

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

7. Дальномерно-разностно-дальномерный способ определения координат источника радиоизлучения и реализующее его устройство. Патент РФ №2363010, С2, опубл. 27.10.2007 г. Авторы: Сайбель А.Г., Вайгель К.И.

8. Определение координат местоположения источников излучения при радиоконтроле. Труды 9-го международного симпозиума ЭМС. - 2011. Авторы: Логинов Ю.И., Екимов О.Б., Антипин Б.М.

9. Владимиров Ю.С. Пространство - время: явные и скрытые размерности. Изд. 2-е, перераб. и доп. - М.: Книжный дом «ЛИБРОКОМ», 2010, 208 с.

10. Дальномерно-разностно-дальномерный способ определения координат местоположения источника радиоизлученияи и реализующее его устройство. Заявка №2011134103/07, опубл. 28.02.2013 г. Авторы: Логинов Ю.И., Екимов О.Б., Антипин Б.М., Гриценко А.А., Павлов В.Н., Портнаго Л.Б.