Способ пеленгации источников радиоизлучения на одной частоте



Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте
Способ пеленгации источников радиоизлучения на одной частоте

 


Владельцы патента RU 2556699:

Федеральное государственное унитарное предприятие "Ростовский-на-Дону научно-исследовательский институт радиосвязи" (ФГУП "РНИИРС") (RU)

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

 

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

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

Известные способы пеленгации нескольких источников радиоизлучения на одной частоте опираются на статистические методы проверки гипотез (критерий отношения правдоподобий [1 - Радзиевский В.Т., Уфаев В.А. Алгоритмы обнаружения и пеленгования совокупности частотно неразделимых радиосигналов // Радиотехника, 2005, №9, с.56-69]), на метод максимума правдоподобия [2 - Дзвонковская А.Л., Дмитриенко А.Н., Кузьмин А.В. Эффективность измерения углов прихода сигнала радиопеленгатора на основе метода максимального правдоподобия // Радиотехника и электроника, 2001, т.46, №10, с.1242-1247], на сверхразрешающие методы [3 - Мюнье Ж., Делиль Ж.Ю. Пространственный анализ в пассивных локационных системах с помощью адаптивных методов // ТИИЭР, 1987, т.75, №11, с.21-37], [4 - Manikas A., Ratnarajah Т., Lee Jinsock Evaluation of Super-resolution Array-Techniques as Applied to Coherent Sources // International Journals of Electronic, Jan. 1997, Vol.82, No.1, pp.77-105] и др. Однако задача пеленгации ИРИ, как некорректная задача, не может быть решена надежно [5 - А.А. Грешилов, А.Л. Лебедев, П.А. Плохута Многосигнальная пеленгация источников радиоизлучения на одной частоте как некорректная задача // Успехи современной радиоэлектроники, 2008, №3, с.30-46]: ни статистическими методами, достоверность результата которых определяется точностью полученных оценок параметров сигналов; ни методом наименьших квадратов (МНК) в силу нелинейности и плохой обусловленности решаемой системы уравнений; ни сверхразрешающими методами, которые дают приемлемые результаты лишь при высоких соотношениях сигнал/шум и не обеспечивают разрешение ИРИ, имеющих близкие по значениям пеленги.

В частности, известен способ пеленгации источника сигнала [6 - Патент РФ №2192651, МПК G01S 3/00, G01S 3/14. Способ пеленгования источника сигнала, войсковая часть 11135. Тынянкин С.И., Апульцына И.В., Бурцев С.Ю. Публ. - 10.11.2002 г.], включающий прием пеленгуемого сигнала элементами двух линейных эквидистантных антенных решеток, расположенных взаимно-перпендикулярно, вычисление пространственного спектра Фурье пеленгуемого сигнала, принятого элементами первой линейной эквидистантной антенной решетки и комплексно-сопряженного пространственного спектра Фурье пеленгуемого сигнала, принятого элементами второй линейной эквидистантной антенной решетки, преобразование масштабов обоих вычисленных пространственных спектров пеленгуемого сигнала по логарифмическому закону, корреляционный анализ и измерение относительного сдвига преобразованных пространственных спектров пеленгуемого сигнала и оценку угловых координат.

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

Известен способ пеленгации радиосигналов [7 - Патент РФ №2144200, МПК G01S 3/14. Способ пеленгации радиосигналов и многоканальный пеленгатор. Ашихмин А.В., Виноградов А.Д., Кондращенко В.Н., Рембовский A.M.. Публ. - 10.01.2000 г.], который включает прием радиосигналов антенной решеткой, состоящей из N антенных элементов, выполненных идентичными в количестве не менее трех и расположенных в плоскости пеленгования, измерение в каждом частотном поддиапазоне комплексных амплитуд пар сигналов, характеризующих фазы каждого радиосигнала, принимаемого в соответствующем частотном поддиапазоне одним из антенных элементов пары, выбранным в качестве сигнального, относительно фазы радиосигнала, принимаемом в том же частотном поддиапазоне другим из антенных элементов пары, выбранным в качестве опорного для всех используемых пар антенных элементов, формирование двухмерных угловых спектров каждого принимаемого в соответствующем частотном поддиапазоне радиосигнала по измеренным комплексным амплитудам пар сигналов для различных пар антенных элементов антенной решетки соответственно взаимному расположению этих антенных элементов в плоскости пеленгования, по которым судят об азимутах и углах места принятых радиосигналов.

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

Рассмотрение решения задач радиопеленгации с учетом их некорректности, например, в [5], [8 - Некорректные задачи цифровой обработки информации и сигналов: монография / А.А. Грешилов. - 2-е изд., доп. - М.: Университетская книга: Логос, 2009.- 360 с], [9 - Патент РФ №2382379, МПК (2006.01) G01S 5/04. Способ многосигнальной пеленгации источников радиоизлучения на одной частоте, ГОУ ВПО "МГТУ им. Н.Э. Баумана". Грешилов А.А., Плохута П.А. Публ. - 20.02.2010 г.], [10 - Патент РФ №2380719, МПК (2006.01) G01S 5/04. Способ пеленгации источников радиоизлучения на одной частоте, ГОУ ВПО "МГТУ им. Н.Э. Баумана". Грешилов А.А. Публ. - 27.01.2010 г.], позволило повысить достоверность получаемых результатов пеленгования, сократить количество необходимых при реализации сложных математических операций, а также получить интервальные оценки результатов пеленгации.

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

1. Радиосигналы источников принимают посредством антенной системы (АС), состоящей из M слабонаправленных элементов (вибраторов), расположенных линейно на расстояниях друг от друга, кратных заданной величине d, как правило, не превышающей значения половины длины волны излучения ИРИ.

2. Получают комплексные амплитуды сигналов на выходах вибраторов ym, m = 1, M ¯ (вектор амплитудно-фазового распределения (АФР) - y = ( y 1 y 2 y M ) T ). При этом математическая модель сигнала для m-го элемента вектора АФР имеет вид

где uk - амплитуда сигнала k-го ИРИ;

φm - фаза сигнала k-го ИРИ на m-м вибраторе, зависящая от азимутального и угломестного пеленгов k-го ИРИ θk и βk соответственно;

nm - шум, имеющий место на m-м вибраторе, включающий в себя шум мирового фона и аппаратуры (с нулевым математическим ожиданием и заданным средним квадратичным отклонением (СКО): σm);

j - мнимая единица.

3. Формируют и решают систему линейных алгебраических уравнений (СЛАУ) вида

запоминают полученный вектор C = ( C 0 C 1 C K 1 ) T ,

где T - знак операции транспонирования.

4. Полагая элементы вектора C коэффициентами полинома

находят и запоминают его корни - экспоненциальные функции, показатели которых содержат информацию об угловом положении ИРИ:

где λ - длина волны сигналов ИРИ.

5. По формуле

находят азимутальные пеленги источников радиоизлучения (если угломестные пеленги равны нулю) или (если угломестные пеленги не равны нулю) по формуле

находят величины произведений косинусов азимутальных и угломестных пеленгов, а затем по тригонометрическим формулам, приведенным, например, в [5] находят значения азимутальных и угломестных пеленгов, запоминают значения пеленгов:

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

;

.

6. Находят диагональные элементы ковариационной матрицы DC. решения СЛАУ (2) - дисперсии элементов вектора коэффициентов полинома (3) при заданном СКО σy элементов матрицы системы (2)

на основе, например, упрощенной формулы [5], [8]:

7. Находят дисперсии пеленгов по формуле, приведенной в [5], [10] для K=2:

либо как скалярные функции случайных аргументов на основе дисперсий (8) и выражений (3)-(7).

8. На основе дисперсий, полученных в п.7 строят соответствующие доверительные интервалы.

9. Производят визуализацию результатов.

Указанный способ обладает следующими недостатками.

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

а. Применение выражений вида (5) обуславливает отсутствие единственности решения, поскольку математическая функция arccos{ϑ} единственному значению (отличному от ϑ=±1) ставит в соответствие два угла в интервале (0°, 360°). Поэтому требуется, например, предложенное в [5], явное задание рабочего сектора пеленгации по каждому из углов - его ширина не должна превышать 180° (что при применении слабонаправленных элементов АС ведет к необходимости введения дополнительных операций в рассматриваемый способ пеленгации).

б. Ограниченность оценки азимутальных углов при использовании формул вида (7) диапазоном от -90° до 90°.

в. Непосредственное применение выражения (5) для оценки азимутального положения ИРИ возможно только в частных случаях получения его действительного значения, когда выполняется условие | ξ k | = 1 , т.е. при некоторых частных наборах реализации вектора АФР. В противном случае

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

2. Алгоритмические и вычислительные сложности получения интервальных оценок, связанные:

- со сложностями оценки дисперсий азимутальных пеленгов ИРИ при их большом числе (основанные на выражениях вида (9), просто вычислимы только при возможности получения аналитических выражений для частных производных θ k C i , k = 1, K ¯ , i = 0, K 1 ¯ , что ограничено числом ИРИ, а для численного определения θ k C i , k = 1, K ¯ , i = 0, K 1 ¯ требуется многократное нахождение корней полинома (3));

- с получением дисперсий оценок азимутальных и угломестных пеленгов как дисперсий существенно нелинейных скалярных функций от комплекснозначных сигналов ξk, k = 1, K ¯ , являющихся в свою очередь функциями случайных величин - оценок решений СЛАУ вида (2).

3. Получаемые интервальные оценки пространственного углового положения ИРИ, основанные только на дисперсиях, не учитывают факта зависимости случайных величин оценок пространственного углового положения ИРИ (оценки азимутального и угломестного положения ИРИ). При формировании интервальной оценки в виде доверительной области (с заданной доверительной вероятностью) это приводит к тому, что она будет представлять собой прямоугольник, стороны которого определяются границами доверительных интервалов оценок азимутального и угломестного положения ИРИ полученных раздельно: прямоугольник, как форма доверительной области системы двух случайных величин, получается в том случае, когда они являются независимыми. Однако, поскольку каждая из случайных величин (оценки азимутального и угломестного положения ИРИ), входящих в рассматриваемую систему, является по сути функцией одного и того же векторного случайного аргумента (составленного из шумов, имеющих место на вибраторах), то в общем случае они являются зависимыми. В этом случае форма доверительной области представляет собой фигуру, вписанную в прямоугольник. Следовательно, интервальные оценки не адекватны реальным условиям. В практическом плане это приводит к увеличению области неопределенности местоопределения ИРИ. Так, например, если форма сечения пространственной фигуры, образованной точкой положения пеленгатора и областью неопределенности оценки азимутального и угломестного положения ИРИ, плоскостью его возможного положения (поверхность Земли) - эллипс, полуоси (a, b) которого расположены параллельно сторонам описывающего его прямоугольника, то площадь поиска ИРИ будет необоснованно увеличена на ( ( 2 a ) ( 2 b ) π a b 1 ) 100 % = 27,3 % .

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

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

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

- упрощение процесса получения интервальных оценок углового положения ИРИ;

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

при сохранении повышенного быстродействия (скорости) обработки сигналов.

Указанный технический результат достигают за счет:

1) определенной топологии слабонаправленных элементов АС;

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

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

Согласно изобретению, до приема сигналов из элементов антенной системы формируют L линейных антенных подсистем (L≥2), расположенных под разными углами относительно друг друга, каждая из которых содержит Ml, l = 1, L ¯ , элементов, выполняют вышеперечисленные операции с сигналами для каждой из антенных подсистем, получая оценки значений экспоненциальных функций для каждого k-го источника радиоизлучения k = 1, K ¯ (K - количество источников радиоизлучения), получают интервальную оценку углового положения каждого k-го источника радиоизлучения как двумерную область в пространстве азимутальных углов и углов места, являющуюся пересечением L областей, полученных для каждой l-й ( l = 1, L ¯ ) антенной подсистемы, при этом каждая из этих Z, областей формируется как отображение областей в комплексных пространствах значений экспоненциальных функций, накрывающих соответствующие им оценки, в пространство дискретных значений угловых координат по азимуту и углу места, при необходимости выделяют точечные оценки углового положения источников радиоизлучения как середины соответствующих интервальных оценок.

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

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

Радиосигналы источников принимают посредством антенной системы (АС), состоящей из L линейных подсистем AC - групп вибраторов (L≥2), каждая из которых содержит Ml, l = 1, L ¯ , элементов, расположенных на расстоянии кратном d между соседними элементами «своей» подсистемы. При этом некоторые элементы могут одновременно принадлежать нескольким группам. В качестве фазового центра, например, выбирается общий элемент, относительно него и направления, определяемого первой подсистемой АС, остальные подсистемы расположены под углами γ2, …, γL в плоскости определения азимутальных направлений θ = ( θ 1 θ 2 θ K ) T . Для определенности пусть таким элементом будет первый элемент каждой антенной подсистемы. - 1

Получают комплексные амплитуды сигналов на выходах вибраторов ym, m = 1, M ¯ (вектор АФР - y = ( y 1 y 2 y M ) T ). При этом, математическая модель сигнала для m-го элемента вектора АФР имеет вид, аналогичный (1). - 2

В соответствии с топологией АС из y = ( y 1 y 2 y M ) T формируют L векторов z 1 = ( z 11 z 12 z 1 M 1 ) T , …, z l 1 = ( z l 1 z l 2 z l M 1 ) T …, z L = ( z L 1 z L 2 z L M L ) T , при этом в связи с выбранным элементом - фазовым центром: z11=z21=…=zL1=y1. - 3

Формируют и решают L СЛАУ вида

,

запоминают полученные векторы коэффициентов полиномов Cl=(Cl0 Cl1 Cl2 … Cl,K-1)T, l = 1, L ¯ . - 4

Находят корни полиномов

- оценки η ^ l k экспоненциальных функций ηlk, k = 1, K ¯ , l = 1, L ¯ , содержащие информацию об угловом положении ИРИ, запоминают их. - 5

На основе интервального анализа с заданной точностью (и скоростью) на основе выражений

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

7. Производят визуализацию результатов. - 7

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

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

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

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

Полагаем, что излучение на одной частоте осуществляет K ИРИ. Требуется по сигналу комплексного амплитудно-фазового распределения y = ( y 1 y 2 y M ) T определить параметры ИРИ: θ = ( θ 1 θ 2 θ K ) T , β = ( β 1 β 2 β K ) T и u = ( u 1 u 2 u K ) T . В зависимости от пеленгационной обстановки и возможностей АС формируют L линейных антенных подсистем - групп вибраторов (L≥2), каждая из которых содержит Ml, l = 1, L ¯ элементов, расположенных на расстоянии кратном d между соседними элементами «своей» подсистемы. При этом некоторые элементы могут одновременно принадлежать нескольким группам. В качестве фазового центра, например, выбирается общий элемент, относительно него и направления, определяемого первой подсистемой АС, остальные подсистемы расположены под углами γ2, …, γL в плоскости определения азимутальных направлений θ = ( θ 1 θ 2 θ K ) T .

Таким образом, из элементов вектора y = ( y 1 y 2 y M ) T формируют L векторов z 1 = ( z 11 z 12 z 1 M 1 ) T , …, z l 1 = ( z l 1 z l 2 z l M 1 ) T …, z L = ( z L 1 z L 2 z L M L ) T , l = 1 L M l = M . Тогда выражения (2)-(5) принимают вид:

а СКО σm, m = 1, M ¯ , представляются в виде σlml, l = 1, L ¯ , m l = 1, M l ¯ .

В соответствии с прототипом решение (10), (11) позволяет найти в общем случае комплекснозначные оценки η ^ l k значений экспоненциальных функций ηlk, k = 1, K ¯ , l = 1, L ¯ . В отличие от прототипа, определение интервальных оценок искомых параметров k-го ИРИ осуществляется не на основе формирования точечной оценки по тригонометрическим выражениям, аргументами которых является только пара значений η ^ l k из L с последующим определением (с большими вычислительными затратами) ковариационной матрицы и формированием интервальных оценок, а путем интервального анализа результатов для всех L антенных подсистем на основе оценок η ^ l k , k = 1, K ¯ , l = 1, L ¯ . Отметим, что интервальные оценки формируются на основе данных, в которых впервые в порядке обработки сигналов от элементов АС, представлена информация об угловом положении каждого ИРИ по отдельности (т.е. каждый набор оценок η ^ l k , l = 1, L ¯ для фиксированного k содержит информацию об угловом положении только k-го ИРИ). Поэтому формирование интервальных оценок угловых положений ИРИ упрощается. Кроме того, в процесс их получения не вносятся погрешности и отсутствуют ограничения, связанные с применением вспомогательных операций над сигналами.

Согласно предлагаемому способу (θkk)инт ((θkk)инт - обозначение интервальной оценки углового положения k-го ИРИ) формируют на основе сравнения оценок η ^ l k со значениями ηlk, полученными по выражению (12). Рассмотрим следующий вариант построения схемы сравнения, основанный на введении системы событий и применении аппарата алгебры событий.

Пусть Ωlkkk) - событие, заключающееся в одновременном выполнении для некоторой точки (θkk) пространства угловых координат θ, β следующих условий:

,

где ε l k Re , ε l k Im ( l = 1, L ¯ , k = 1, K ¯ ) - параметры, характеризующие размеры ε - окрестностей оценок η ^ l k в области действительных и мнимых их частей и следовательно, интервальных оценок пеленгов в пространстве угловых координат θ, β.

Значения ε l k Re , ε l k Im определяются на этапе настройки реализации предлагаемого способа и в общем случае зависят от топологии АС и отношения сигнал/шум. В качестве вариантов обоснования их значений можно использовать методы статистического моделирования или конфлюэнтного анализа для получения СКО Re [ η ^ l k ] и Im [ η ^ l k ] на основе значений σlml, l = 1, L ¯ , m l = 1, M l ¯ . Кроме того обоснование выбора их значений можно выполнить на основе построения доверительных областей с заданными доверительными вероятностями в соответствии с подходами многомерного статистического анализа, изложенными в [11 - Химмельблау Д. Анализ процессов статистическими методами. - М.: Мир, 1973. - С.208], [12 - Чанышева А.Ф. Нахождение интервальной оценки комплексного уравнения регрессии // Экономическое прогнозирование: модели и методы: материалы Международной научно-практической конференции. 5-6 апреля 2009 г.: в 2 ч., под ред. Давниса В.В. Воронеж.: Изд-во ВГУ, 2009].

Тогда событие, заключающееся в одновременном появлении всех событий Ωlkkk), k = 1, K ¯ , l = 1, L ¯ представимо в виде:

При этом каждому событию Ωlkkk) соответствует точка в пространстве угловых координат θ, β, принадлежащая интервальной оценке

где (θkk)l - множество значений азимутальных углов и углов места, характеризующих положение k-го ИРИ, полученных на основе оценок η ^ l k l-й антенной подсистемы, а из (13), (14) следует:

Практическую реализацию приближенного поиска (θkk)инт можно представить в виде процедуры формирования сети точек (θkk) в пространстве угловых координат θ, β, проверки их принадлежности допустимому множеству по выражениям:

с последующим определением (θkk)инт как оболочки множества допустимых точек.

Таким образом, интервальная оценка (θkk)инт углового положения k-го ИРИ есть двумерная область в пространстве азимутальных углов и углов места, являющаяся пересечением L областей, полученных для каждой l-й ( l = 1, L ¯ ) антенной подсистемы. При этом каждая из этих L областей формируется как отображение областей в комплексных пространствах значений ηlk, накрывающих соответствующие им оценки η ^ l k , в пространство дискретных значений угловых координат θ, β.

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

Критерий определения количества сигналов, принимаемых АС, относительно способа-прототипа модифицируется следующим образом: их количество равно обобщенному (по одному из возможных правил) числу «больших» собственных чисел матриц A C l T A C l ,

.

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

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

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

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

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

Таблица 1
Последовательность реализации способа-прототипа Последовательность реализации предлагаемого способа
0. Предварительный этап: при пеленгации ИРИ, расположенных в секторе по азимуту большем от 0° до 180°, а также при отличии угломестного положения ИРИ от 0° - выделение в составе линейной АС отдельного элемента - вибратора, смещенного относительно линии расположения остальных. 0. Предварительный этап: формируют из M слабонаправленных элементов (вибраторов) L (L≥2) линейных антенных подсистем с минимум одним общим элементом.
1. Принимают радиосигналы источников излучения посредством АС, состоящей из М слабонаправленных элементов (вибраторов), расположенных линейно на расстояниях друг от друга, кратных заданной величине d, как правило, не превышающей значения половины длины волны излучения ИРИ. 1. Принимают радиосигналы источников излучения посредством АС, состоящей из M слабонаправленных элементов (вибраторов), расположенных линейно на расстояниях друг от друга, кратных заданной величине d, как правило, не превышающей значения половины длины волны излучения ИРИ.
2. Получают комплексные амплитуды сигналов на выходах вибраторов ym,
m = 1, M ¯
(вектор АФР -
y = ( y 1 y 2 y M ) T ) .
2. Получают комплексные амплитуды сигналов на выходах вибраторов ym,
m = 1, M ¯
(вектор АФР -
y = ( y 1 y 2 y M ) T ) .
3. Формируют и решают систему линейных алгебраических
уравнений вида (2), запоминают полученное решение - вектор
C = ( C 0 C 1 C K 1 ) T .
3. Формируют и решают L систем линейных алгебраических
уравнений вида (10), запоминают полученные решения - Clk,
l = 1, L ¯ , k = 1, K ¯ .
4. Полагая элементы вектора C коэффициентами полинома (3) находят и запоминают его корни ξk,
k = 1, K ¯ .
4. Полагая элементы вектора C коэффициентами полиномов (11) находят и запоминают их корни η ^ l k ,
l = 1, L ¯ , k = 1, K ¯ .
5. Находят точечные оценки азимутальных и угломестных пеленгов по выражениям (5)-(7), содержащие функции, ограничивающие применение способа по рабочему сектору углового положения ИРИ и совокупности полученных реальных измерений. 5. Находят интервальные оценки (θkk)инт углового положения каждого k-го ИРИ в виде двумерных областей в пространстве азимутальных углов и углов места, являющихся пересечением L областей, полученных для каждой l-й
( l = 1, L ¯ )
антенной подсистемы. При этом каждая из этих L областей формируется как отображение областей в комплексных пространствах значений ηlk, накрывающих соответствующие им оценки η ^ l k , в пространство дискретных значений угловых координат θ, β.
6. Находят дисперсии элементов вектора коэффициентов полинома (3) при заданном СКО σy элементов матрицы системы (2).
7. Находят дисперсии пеленгов по формуле (9) или как скалярные функции случайных аргументов на основе дисперсий п.6 и выражений (3)-(7), являющиеся источником алгоритмической и вычислительной сложности.
8. На основе полученных в п.5 - п.7
результатов формируют интервальные оценки в виде соответствующих доверительных интервалов, не учитывающих зависимость случайных величин - точечных оценок азимутальных и угломестных пеленгов.
6. При необходимости выделяют
точечные оценки углового положения источников радиоизлучения как середины соответствующих интервальных оценок.
9. Производят визуализацию результатов. 7. Производят визуализацию результатов.

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

Порядок выполнения интервального анализа рассмотрим на примере пеленгации двух ИРИ, работающих на частоте 1000 МГц. Характеристики углового положения ИРИ: θ1=35°, β1=20°, θ2=195°, β2=40°; амплитудные характеристики ИРИ: u1=10 мВ, u2=8 мВ.

Помеха имеет математическое ожидание, равное нулю, и СКО σm=0,1 мВ. Пеленгацию будем осуществлять посредством АС, состоящей из двух линейных подсистем по 7 вибраторов в каждой, отстоящих друг от друга на расстоянии, равном 0,15 м. Угол между подсистемами в азимутальной плоскости: γ2=90°.

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

С целью сравнения результатов зафиксируем реализацию шумового воздействия. В таблице 2 представлены значения шумов в сигналах z 1 = ( z 11 z 12 z 17 ) T , z 2 = ( z 21 z 22 z 27 ) T :

Таблица 2
Номер элемента в подсистеме на элементы первой подсистемы на элементы второй подсистемы
1 -0,040598 -0,040598
2 -0,104730 -0,153490
3 0,153570 0,022137
4 0,043443 -0,137450
5 -0,191710 -0,083929
6 0,046994 -0,020864
7 0,127440 0,075591

Для представленной реализации шумов значения оценок экспоненциальных функций принимают вид: η ^ 11 = 0,74939 + j 0,66402 , η ^ 12 = 0,68497 j 0,73201 , η ^ 21 = 0,1298 j 0,99151 , η ^ 22 = 0,8134 + j 0,58394 .

Реализация интервального анализа в пространстве угловых координат θ ∈ [0°, 360°], β ∈ [0°, 90°] с шагом в один градус для εlk=ε ( l = 1, L ¯ , k = 1, K ¯ ) на основе выражений (14), (15) позволила получить следующие результаты, представленные в таблицах 3-7.

Результаты интервальной оценки при ε=0,05

Таблица 3
Номер положительного результата при сравнении Значение азимутального угла, град. Значение угла места, град.
1 34 19
2 34 20
3 35 17
4 35 18
5 35 19
6 35 20
7 35 21
8 35 22
9 35 23
10 36 20
11 36 21
12 194 39
13 194 40
14 194 41
15 195 39
16 195 40
17 195 41
18 196 40
19 196 41

Результаты интервальной оценки при ε=0,04

Таблица 4
Номер положительного результата при сравнении Значение азимутального угла, град. Значение угла места, град.
1 35 18
2 35 19
3 35 20
4 35 21
5 35 22
6 195 39
7 195 40
8 195 41

Результаты интервальной оценки при ε=0,03

Таблица 5
Номер положительного результата при сравнении Значение азимутального угла, град. Значение угла места, град.
1 35 19
2 35 20
3 35 21
4 35 22
5 195 40

Результаты интервальной оценки при ε=0,02

Таблица 6
Номер положительного результата при сравнении Значение азимутального угла, град. Значение угла места, град.
1 35 19
2 35 20
3 35 21
4 195 40

Результаты интервальной оценки при ε=0,01

Таблица 7
Номер положительного результата при сравнении Значение азимутального угла, град. Значение угла места, град.
1 35 20
2 195 40

Оценки, полученные на основе способа-прототипа:

θ1=35,03°, β1=18,29°, θ2=15,00°, β2=139,91°.

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

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

Отметим, что при отсутствии шумов значения оценок экспоненциальных функций принимают вид: η11=-0,7507+j0,66064, η12=-0,68559-j0,72799, η21=-0,12333-j0,99237, η22=0,81195+j0,58372, т.е. абсолютные погрешности определения их действительных и мнимых частей, соответственно равны (таблица 8).

Таблица 8
l, k | Re [ η l k ] Re [ η ^ l k ] | | Im [ η l k ] Im [ η ^ l k ] |
1, 1 0,00131 0,00338
1, 2 0,00062 0,00402
2, 1 0,00035 0,00086
2, 2 0,00145 0,00022

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

При этом оценки, полученные на основе способа-прототипа, существенно не улучшаются: θ1=35,00°, β1=13,09°, θ2=15,00°, β2=148,41°. Результаты по предлагаемому способу для этого случая представлены в таблицах 9-13.

Результаты интервальной оценки при ε=0,05

Таблица 9
Номер положительного результата при сравнении Значение азимутального угла, град. Значение угла места, град.
1 34 19
2 34 20
3 35 17
4 35 18
5 35 19
6 35 20
7 35 21
8 35 22
9 35 23
10 36 20
11 36 21
12 194 39
13 194 40
14 194 41
15 195 39
16 195 40
17 195 41
18 196 40
19 196 41

Результаты интервальной оценки при ε=0,04

Таблица 10
Номер положительного результата при сравнении Значение азимутального угла, град. Значение угла места, град.
1 35 18
2 35 19
3 35 20
4 35 21
5 35 22
6 195 39
7 195 40
8 195 41

Результаты интервальной оценки при ε=0,03

Таблица 11
Номер положительного результата при сравнении Значение азимутального угла, град. Значение угла места, град.
1 35 18
2 35 19
3 35 20
4 35 21
5 195 40

Результаты интервальной оценки при ε=0,02

Таблица 12
Номер положительного
результата при сравнении
Значение азимутального угла, град. Значение угла места, град.
1 35 19
2 35 20
3 35 21
4 195 40

Результаты интервальной оценки при ε=0,01

Таблица 13
Номер положительного результата при сравнении Значение азимутального угла, град. Значение угла места, град.
1 35 20
2 195 40

Таким образом, предлагаемый способ обладает устойчивостью к действию рассматриваемого класса помех. Кроме того, приведенные результаты показывают, что способ-прототип даже при отсутствии помех дает существенную погрешность в определении углов места и азимутальных углов, больших 180°.

Способ пеленгации источников радиоизлучения на одной частоте, при котором принимают многолучевой сигнал посредством многоэлементной антенной системы, синхронно преобразовывают ансамбль принятых сигналов, зависящих от времени и номера элемента антенной системы, в цифровые сигналы, преобразовывают цифровые сигналы в сигнал амплитудно-фазового распределения, описывающий распределение амплитуд и фаз на элементах антенной системы, формируют из амплитуд комплексных сигналов амплитудно-фазового распределения системы линейных алгебраических уравнений, первая строка которой представляет собой в левой части уравнения сумму произведений коэффициентов полинома экспоненциальных функций и амплитуд сигналов с линейно расположенных антенных элементов, взятых последовательно, начиная с амплитуды, соответствующей антенному элементу, выбранного в качестве фазового центра, число элементов суммирования равно количеству источников радиоизлучения, а правая часть представляет невыбранное еще значение амплитуды сигнала со следующего последовательно расположенного антенного элемента, при этом каждое последующее уравнение в системе сформированных линейных алгебраических уравнений является сдвигом коэффициентов уравнения вправо на одну позицию по отношению к предыдущему уравнению, определяют решения этой системы уравнений и определяют корни полинома - экспоненциальные функции, аргументами каждой из которых являются азимутальные и угломестные пеленги только одного из источников радиоизлучения, отличающийся тем, что до приема сигналов из элементов антенной системы формируют L линейных антенных подсистем, L≥2, расположенных под разными углами относительно друг друга, каждая из которых содержит Ml, l = 1, L ¯ , элементов, выполняют вышеперечисленные операции с сигналами для каждой из антенных подсистем, получая оценки η l k значений экспоненциальных функций ηlk=ехр(j(2π/λ)dcos(θkl)cosβk), k = 1, K ¯ , l = 1, L ¯ , где K - количество источников радиоизлучения, j - мнимая единица, λ - длина волны сигналов источников радиоизлучения, d - расстояние между соседними элементами антенных подсистем, θk, βk - азимутальный и угломестный пеленги k-го источника радиоизлучения, γl - угловое положение в азимутальной плоскости линейных антенных подсистем относительно первой из них, l = 1, L ¯ , γl=0, получают интервальную оценку (θk, βk)инт углового положения каждого k-го источника радиоизлучения как двумерную область в пространстве азимутальных углов и углов места, являющуюся пересечением L областей, полученных для каждой l-й, l = 1, L ¯ , антенной подсистемы, при этом каждая из этих L областей формируется как отображение областей в комплексных пространствах значений ηlk, накрывающих соответствующие им оценки η l k , в пространство дискретных значений угловых координат по азимуту и углу места, при необходимости выделяют точечные оценки углового положения источников радиоизлучения как середины соответствующих интервальных оценок (θk, βk)инт.



 

Похожие патенты:

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

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

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

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

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

Способ предназначен для мониторинга радиоэлектронной обстановки при многолучевом распространении радиоволн, воздействии преднамеренных и непреднамеренных помех, отражениях сигнала от различных объектов и слоев атмосферы.

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

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

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

Изобретение относится к измерительной технике, в частности к пеленгаторам. .

Изобретение относится к радиотехнике и может быть использовано для определения местоположения и скорости априорно неизвестного источника радиоизлучения (ИРИ). Достигаемый технический результат - определение за один этап обработки одновременно координат и скорости ИРИ. Способ основан на вычислении количества N элементарных зон привязки возможного расположения ИРИ, определении координат местоположения центров элементарных зон привязки, присвоении каждой элементарной зоне привязки порядкового номера n=1, 2, …, N, задании полосы частот ΔF, в которой ведется прием сигналов, разбиении заданной полосы частот ΔF на P поддиапазонов шириной Δf, присвоении каждому поддиапазону порядкового номера p=1, 2, … P, определении для R взаимосвязанных периферийных и центрального пеленгаторных пунктов (ПП) с известным их местоположением, каждый из которых включает M антенных элементов, значений эталонных первичных пространственно-информационных параметров (ППИП), на выходах всех антенных элементов, которые рассчитывают для средних частот всех частотных поддиапазонов, приеме сигналов ИРИ в заданной полосе частот ΔF всеми ПП, измерении ППИП для каждого антенного элемента всех ПП и передаче их с периферийных ПП на центральный ПП, при этом перед определением эталонных ППИП задают диапазоны возможных значений составляющих скорости ИРИ, разбивают заданные диапазоны на G каналов каждый, присваивают каждому каналу порядковый номер c=1, 2 … G, d=1, 2 … G, задают интервал обработки, определяемый шириной полосы частотного поддиапазона, и время накопления сигнала, для каждого канала составляющих скорости, каждой элементарной зоны привязки, каждого антенного элемента всех ПП определяют значения эталонных ППИП, измеряют ППИП принятых сигналов исходя из ожидаемого положения ИРИ в каждом частотном поддиапазоне, на каждом антенном элементе всех ПП, далее для каждой элементарной зоны привязки, каждого антенного элемента всех ПП, каждого частотного поддиапазона и для каждого канала составляющих скорости определяют произведение измеренных и эталонных ППИП, полученные произведения суммируют по всем антенным элементам, находят абсолютное значение суммы и результат суммируют по частотным каналам, и по положению максимума результирующей суммы, определяемого по совокупности значений суммы в дискретных точках n, c, d, определяют координаты ИРИ и его скорость. 7 ил.

Изобретение относится к области локационной техники и может быть использовано в системах поиска объектов. Достигаемый технический результат - повышение точности определения направления на импульсные излучатели. Указанный результат достигается за счет того, что устройство содержит две разнесенные антенны с полями зрения 180°, два амплитудных селектора, дешифратор, блок определения малого временного интервала, постоянное запоминающее устройство, при этом выходы первой и второй антенн с углами поля зрения 180° соответственно соединены через первый и второй приемники, через первый и второй амплитудные селекторы с первым и вторым входами блока определения малого временного интервала, имеющего группу выходов, соединенную с группой входов постоянного запоминающего устройства и с группой входов дешифратора, первый и второй выход которого соединен с первым и вторым входом постоянного запоминающего устройства, имеющего группу выходов, соединенную с группой входов вычислителя и со второй группой входов блока вторичной обработки. 2 ил.
Наверх