СПОСОБ ПЕЛЕНГАЦИИ УЗКОПОЛОСНЫХ РАДИОСИГНАЛОВ КВ ДИАПАЗОНА
LOCATION FINDING METHOD OF NARROW-BAND RADIO SIGNALS OF SHORT-WAVE RANGE

Изобретение относится к радиотехнике и может быть использовано для повышения точности пеленгации радиосигналов в коротковолновом (КВ) диапазоне частот. Характерная особенность КВ диапазона - это многолучевость, т.е. возможность распространения радиоволн по нескольким несовпадающим траекториям.
Известен способ широкополосного обнаружения и пеленгации [1], при котором из выходных сигналов каждого элемента антенной решетки выделяются цифровые сигналы, характеризующие спектры принятых сигналов, и для каждой выбранной частоты в полосе приема, используя фазу сигналов, производится прямое вычисление пространственного ряда Фурье, дискретно описывающего угловой спектр мощности на выбранной частоте. После восстановления углового спектра на всех частотах определяется пеленг любого источника, излучающего сигналы на любой из частот в пределах текущей полосы приема. Этот способ из максимально возможной амплитудно-фазовой информации использует только фазу сигнала, обладает низким быстродействием при определении азимутального пеленга и не позволяет разрешать многолучевые сигналы, особенно при их близком угловом положении.
Известен также способ [2] определения пеленга и частоты источников излучения, основанный на двумерном Фурье-преобразовании (по времени и пространству) сигналов, принятых кольцевой решеткой из N элементов на каждой временной частоте. При фиксированных значениях наведения по углу места каждую сформированную в результате преобразования пространственную частотную составляющую q раз умножают на Фурье-образ функции наведения, зависящей от частоты, угла места, параметров антенной решетки и азимутального циклического сдвига 2πq/M. Далее производят q пространственных обратных Фурье-преобразований полученной совокупности скорректированных пространственных частотных составляющих, определяют наличие сигнала и получают оценки азимута и угла места его источника на каждой частоте по квадрату модуля двумерного комплексного углового спектра. К недостаткам данного способа можно отнести то, что при малом количестве излучателей в антенной решетке точность определения пеленга оказывается не достаточно высокой. Также он не позволяет разделить многолучевые сигналы.
Наиболее близким к заявленному изобретению является способ пеленгования источников радиоизлучений в условиях многолучевости, описанный в [3]. Поставленная задача достигается тем, что указанный способ включает прием сигналов на N-элементную антенную решетку (AP), формирование сигналов сканирования в заданном угловом секторе, расчет весовых коэффициентов и последующее перемножение весовых коэффициентов с соответствующими сигналами сканирования, запоминание полученных сигналов для всех дискретных значений угла сканирования. Затем производится выбор из них величины угла, соответствующего сигналу с максимальной амплитудой, сигналы повторно принимаются, рассчитываются их весовые коэффициенты, после чего их перемножают с соответствующими сигналами сканирования для каждого дискретного значения угла сканирования во время передачи информационного пакета и защитного интервала соответственно. Далее вычитают результирующий сигнал, рассчитанный во время защитного интервала из результирующего сигнала, рассчитанного в интервале передачи информационного пакета, запоминают разностные сигналы для всех дискретных значений угла сканирования, из которых выбирают значение угла, соответствующее сигналу с максимальной амплитудой. Однако изложенный метод требует знания защитного интервала для каждого из сигналов, большого объема вычислений и не обеспечивает требуемой точности определения направлений прихода, особенно для близкорасположенных сигналов, а также определения амплитуд этих сигналов. Кроме того, точность пеленгации сигналов в вышеперечисленных способах существенно ухудшается при малых отношениях мощности пеленгуемых сигналов к мощности шумов и помех.
Целью изобретения являются обеспечение высокой точности пеленгования сигналов с полностью неизвестными структурой и параметрами и определения их комплексных амплитуд при приходе этих сигналов по многолучевым радиотрассам с существенно различными амплитудами и при малом отношении сигнал/(помеха + шум).
x
x
k
Поставленная задача достигается тем, что в способе пеленгации узкополосных радиосигналов КВ диапазона, заключающемся в том, что узкополосные радиосигналы КВ диапазона, распространяющиеся по многолучевым радиотрассам, каждый из которых лежит в соответствующем k-ом (k=1,…,K) частотном поддиапазоне заданного интервала К частотных поддиапазонов, принимают каждым излучателем многоэлементной антенны и вычисляют корреляционную матрицу принятых сигналов, согласно изобретению в качестве многоэлементной антенны используют линейную антенную решетку, состоящую из N независимых приемных излучателей, осуществляют прием узкополосных радиосигналов, падающих в угломестной плоскости, в течение интервала наблюдения, принятые сигналы разделяют с помощью K узкополосных фильтров, установленных на входе каждого независимого приемного излучателя и перекрывающих заданный интервал K частотных поддиапазонов, корреляционную матрицу  вычисляют по принятым сигналам в каждом k-ом (k=1,…,K) частотном поддиапазоне, раскладывают каждую из k-ых корреляционных матриц  по собственным сигнальным  и шумовым  (j=1,…,(N-1)) векторам, определяют собственные диаграммы направленности антенны  и  в каждом k-ом частотном поддиапазоне как диаграммы направленности в режиме возбуждения независимых приемных излучателей соответствующими комплексно-сопряженными собственными векторами:  и  где , k - номер частотного поддиапазона, λ - средняя длина волны сигнала в k-ом частотном поддиапазоне, n - номер независимого приемного излучателя, d - расстояние между независимыми приемными излучателями, а направления прихода  и комплексные амплитуды  принимаемых сигналов в k-ом частотном поддиапазоне определяют путем решения систем уравнений:
, j=1,…,(N-1),
m=1,…, M, l=1,…,M.
где .
l
1
2
pl
l
Поставленная задача также достигается тем, что из найденных решений систем уравнений выбирают p решений, для которых направления прихода сигналов  совпадают, где k=k, k,…,k; l=1,…,L, с последующим объединением p частотных поддиапазонов в полосу частот одного l-го широкополосного сигнала.
1
2
pl
Поставленная задача также достигается тем, что спектр широкополосного сигнала формируют из спектральных составляющих , где k=k, k,…,k, l=1,…,L, объединенных частотных поддиапазонов.
Поставленная задача также достигается тем, что для каждого из найденных направлений прихода  сигналы от каждого излучателя суммируют с весовыми коэффициентами, пропорциональными комплексно-сопряженному сигнальному собственному вектору  для этого направления, умноженному на обратную корреляционную матрицу помех .
y
x
Поставленная задача также достигается тем, что дополнительно осуществляют прием узкополосных радиосигналов, падающих с произвольного направления  в угломестной и  азимутальной плоскостях, в течение интервала наблюдения с помощью дополнительной линейной антенной решетки, состоящей из N независимых приемных излучателей, расположенной перпендикулярно оси основной линейной антенной решетки из N независимых приемных излучателей, определяют обобщенные направления прихода  и  сигналов из решения систем уравнений:
x
y
, j=1,…,(N-1) и , j=1,…,(N-1),
где  и  - собственные шумовые диаграммы направленности для основной и дополнительной линейных антенных решеток соответственно, после чего восстанавливают истинные направления прихода сигналов  и  по соотношениям:
Изобретение поясняется чертежами.
c
ш
c
ш
c
ш
c
ш
На фиг.1 схематично представлена антенная решетка из N излучателей для одномерной пеленгации; на фиг.2 а, б приведены собственные амплитудные диаграммы направленности (ДН) антенны, состоящей из двух излучателей, при условии одного падающего на антенну сигнала при соотношении мощности сигнала к мощности шума P/Р=10 дБ и 0 дБ; на фиг.3а, б показаны собственные амплитудные ДН антенны, состоящей из четырех излучателей, при условии одного падающего на антенну сигнала при соотношении мощности сигнала к мощности шума P/Р=10 дБ и 0 дБ; на фиг.4а, б приведены собственные амплитудные ДН антенны, состоящей из двух излучателей, при условии двух падающих на антенну сигналов при соотношении суммарной мощности сигналов к мощности шума P/Р=10 дБ и 0 дБ; на фиг.5а, б приведены собственные амплитудные ДН антенны, состоящей из четырех излучателей, при условии двух падающих на антенну сигналов при соотношении суммарной мощности сигналов к мощности шума P/Р=10 дБ и 0 дБ; на фиг.6 приведена схема антенной части КВ пеленгатора для реализации способа; на фиг.7 схематично представлена антенная решетка для двумерной пеленгации; на фиг.8 и 9 приведены таблицы с данными об ошибках определения параметров приходящих сигналов.
1
k
K
Антенная часть КВ пеленгатора (фиг.6) содержит линейную антенную решетку 1, состоящую из N (1,…n,…N) независимых приемных излучателей; N блоков 2, каждый из которых состоит их K узкополосных фильтров, настроенных соответственно на частоты (ω…ω…ω), при этом цифрой 3 обозначена контрольная плоскость, в которой снимаются комплексные амплитуды сигналов  блок 4 управления весовыми коэффициентами; блок 5 сумматоров.
k
Предлагаемый способ основан на разбиении заданного частотного диапазона с помощью узкополосных аналоговых фильтров, которые имеют большой динамический диапазон, на K более узкополосных частотных поддиапазонов. Разбиение на K частотных поддиапазонов производится таким образом, что в каждом из k-ых частотных поддиапазонов (k=1,…K) в данный момент времени существует только один падающий на многоэлементную антенну сигнал с некоторой средней гармонической частотой ω. Этот сигнал из-за многолучевости распространения может приходить на антенную решетку (AP) КВ пеленгатора одновременно с одного или большего числа направлений с неизвестными комплексными амплитудами  для каждого m-го направления. В процессе пеленгации должны быть определены направления  и амплитуды  приходящих сигналов в каждом k-ом частотном поддиапазоне.
k
Рассмотрим способ пеленгации сигнала на примере одномерной пеленгации в фиксированной плоскости (φ=0) для антенны из N независимых приемных излучателей, расположенных на расстоянии d друг от друга вдоль оси ох (фиг.1). При падении плоской электромагнитной волны с комплексной амплитудой  и с гармонической частотой ω с направления  на входе n-го излучателя возникает напряжение
k
k
где  - набег фазы между антенными элементами, определяемый направлением на источник радиосигнала,  - начальная фаза сигнала частоты ω с m-го направления, а  - комплексная амплитуда напряжения на входе n-го излучателя на частоте ω.
ш
Как правило, наряду с гармоническими сигналами на излучатели падает шумовой сигнал с мощностью Р в полосе приема. Тогда усредненную за период наблюдения T комплексную амплитуду сигнала + шума на входе излучателей можно записать в виде
По комплексным амплитудам и  определяется корреляционная матрица напряжений на входах излучателей , которая в предположении некоррелированности принимаемых шумов в соседних излучателях как с полезным сигналом, так и между собой имеет вид:
где  - корреляционная матрица сигналов в отсутствии шумов,  обозначает вектор-столбец,  обозначает вектор-строку); [E] - единичная матрица.
ш
Матрица  при Р>0 является неособенной. По корреляционной матрице определяется N отличных от нуля собственных значений и N соответствующих им собственных векторов. Собственные значения сигнально-шумовой корреляционной матрицы имеют вид:
где  - суммарная мощность сигнала, принимаемая всеми излучателями AP.
2
N
Структура собственного вектора, соответствующего собственному числу  и называемого сигнальным , при любом количестве падающих волн и количестве излучателей совпадает с сигнальным вектором входных напряжений на клеммах излучателей. Остальные N-1 собственных векторов, соответствующих λ,…,λ, называются шумовыми и обозначаются  (j=1,…,(N-1)). С учетом указанных обозначений разложение сигнально-шумовой корреляционной матрицы по собственным векторам имеет вид:
k
с
ш
Определим собственную ДН Fq(θ), соответствующую вектору  как ДН в режиме возбуждения излучателей Ap комплексно-сопряженным вектором . Тот факт, что возбуждение излучателей является комплексно-сопряженным вектору , при условии одного падающего сигнала обеспечивает формирование собственной сигнальной ДН AP с максимумом и собственных шумовых ДН с глубоким провалом (в идеале - нулевым значением) в направлении прихода сигнала (фиг.2) для двух излучателей и (фиг.3) - для четырех излучателей на частоте  при соотношении мощности сигнала к мощности шума Р/Р=10 дБ и 0 дБ). Поэтому для одного сигнала (M=1) направление его прихода определяется по совпадающим нулям собственных шумовых ДН.
Собственные сигнальная и шумовые ДН имеют вид:
k
где , где λ - длина волны сигнала в k-ом частотном поддиапазоне.
m
По сигнальным и шумовым собственным ДН определяется уровень полезного сигнала , принимаемого излучателями AP. Для M плоских волн, падающих с направлений θ на AP из излучателей, уровень полезного сигнала, принимаемого по j-й шумовой ДН, определяется выражением:
c
ш
m
Как показывает анализ, в случае двух и большего количества сигналов, приходящих с различных направлений, направления минимумов шумовых ДН уже не совпадают с направлениями прихода сигналов (фиг.4 для двух излучателей и фиг.5 для четырех излучателей на частоте  при соотношении мощности сигналов к мощности шума P/Р=10 дБ и 0 дБ), как для случая одного сигнала. Поэтому по минимумам собственных шумовых ДН нельзя определить направления двух и более сигналов. Для определения направлений θ нескольких (M>1) сигналов решается система функциональных уравнений:
которое следует из того факта, что уровень полезного сигнала, принимаемого по собственным шумовым ДН, равен нулю.
Для определения абсолютных значений амплитуд  (m=1,…,M) приходящих сигналов при найденных из (9)  используется М дополнительных уравнений:
или, учитывая, что структура собственного сигнального вектора совпадает с сигнальным вектором, система уравнений (10) принимает вид:
c
ш
c
ш
Кривые для сигнальной и шумовых ДН, приведенные на фиг.2, 3, 4, 5, очень слабо зависят от соотношения P/Р, что указывает на возможность применения изложенного ниже способа высокоточной пеленгации и при малых соотношениях P/Р.
После определения из систем уравнений (9), (11) направления прихода и амплитуды m-й плоской волны осуществляется операция максимизации отношения С/(П+Ш) для этого направления  с помощью весовых коэффициентов , вводимых устройствами управления амплитудой и фазой сигнала на выходе каждого n-го излучателя (фиг.6). С этой целью вектор весовых коэффициентов  выбирается из условия
k
где  - корреляционная матрица помех, в которой учитываются как шумовые составляющие, так и сигнальные составляющие со всех направлений за исключением направления ;  - собственный сигнальный вектор для m-го направления прихода сигнала с частотой ω:
После операции (12) сигнал отправляется на дальнейшую обработку.
Для реализации описанного выше способа определения направлений прихода сигналов и их амплитуд используется КВ пеленгатор (фиг.6).
Способ заключается в следующем. Сигнал, приходящий с произвольного m-го направления, попадает на антенную решетку, состоящую из N независимых приемных излучателей, на выходе каждого из которых стоит система узкополосных фильтров. С выхода каждого k-го фильтра n-го излучателя снимается гармонический сигнал .
После узкополосной фильтрации по соотношению (3) вычисляется корреляционная матрица комплексных амплитуд сигналов, , пришедших в каждом k-ом частотном поддиапазоне. Направления прихода сигналов находятся из условия равенства нулю уровня полезного сигнала, принимаемого по шумовым собственным диаграммам направленности, определяемым по собственным векторам корреляционной матрицы, соотношение (9). Для максимизации отношения сигнал/(помеха + шум) для направлений прихода сигналов, найденных на предыдущем этапе, производится суммирование сигналов с каждого n-го излучателя для каждого k-го поддиапазона с весовыми коэффициентами , выбираемыми в соответствии с законом (12).
Описанный способ определения направлений прихода и амплитуд гармонических сигналов можно использовать не только для гармонических сигналов, но и для широкополосных сигналов или сигналов со скачкообразно изменяющейся частотой. В этом случае алгоритм должен быть дополнен следующими операциями: после того как определены все возможные направления прихода сигналов во всех частотных поддиапазонах, необходимо объединить в один частотный интервал те из частотных поддиапазонов, в которых направления прихода сигналов одинаковы. Что касается амплитудного спектра широкополосного сигнала, то он может быть найден как совокупность спектральных составляющих с комплексными амплитудами , найденными из (11) в каждом из k-ых частотных поддиапазонов, входящих в его частотный интервал.
Для двумерной пеленгации используются две ортогонально расположенные линейные антенные решетки (фиг.7). При этом для нахождения угломестных  и азимутальных  углов прихода сигналов должна быть выполнена процедура определения обобщенных координат  и  в соответствии с алгоритмом нахождения направлений прихода сигналов на линейную антенную решетку из N излучателей для случая одномерной пеленгации.
Отличие двумерной пеленгации от одномерной состоит в том, что вначале определяются обобщенные координаты  и , а затем в соответствии с алгоритмом для одномерной пеленгации определяются истинные углы прихода сигналов и их амплитуды.
k
x
y
При падении плоской электромагнитной волны с комплексной амплитудой  и с гармонической частотой ω с направления  на входах n-го и n-го излучателей возникают комплексные амплитуды напряжений:
x
y
x
y
где  N и N - количество излучателей вдоль осей ox и oy соответственно, d и d - расстояние между излучателями в первой и второй линейных антенных решетках соответственно.
Соотношения (14) и (15) можно рассматривать независимо друг от друга и таким образом определить  и  в соответствии с алгоритмом определения углов, справедливым для одномерной пеленгации.
Обобщенные координаты связаны с искомыми параметрами сигналов соотношениями
Используя (16) и (17), по найденным обобщенным координатам восстанавливаются истинные направления прихода сигналов  и  в пространстве.
Таким образом, использование дополнительной линейной антенной решетки, ортогональной первой, позволяет снизить размерность сигнально-шумовой матрицы  для двумерного пеленга до размерности двух независимых одномерных матриц соответствующих линейных AP. Данный способ существенно сокращает объем вычислений при определении углов прихода  и .
01
02
1
2
1
2
01
02
01
02
1
2
1
2
c01
c02
01
02
В качестве наглядного примера приведены результаты моделирования для одномерной пеленгации одного и двух гармонических узкополосных сигналов с частотой  которые показывают, что ошибка определения угловых координат источников радиосигналов и их амплитуд зависит от точности решения систем уравнений (9), (11) и может быть сколь угодно малой. В частности, при задании шага решения уравнений (9), (11) в интервале 0.1°÷1° максимальная ошибка определения угловых координат сопоставима с величиной шага, а относительная точность определения амплитуд составляет несколько процентов, см. табл.1-6 (фиг.8, 9). В этих таблицах введены следующие обозначения: θ и θ - истинные направления прихода сигналов, θ и θ - направления, определяемые из решения системы (9), Δθ и Δθ - ошибки определения направлений прихода сигналов, A и A - истинные значения амплитуд сигналов, принимаемых AP с направлений θ и θ, A и A - значения амплитуд сигналов, полученные в результате решения системы (11), ΔA и ΔA - величины относительной ошибки определения амплитуды первого и второго сигналов соответственно, P и P - мощности сигналов, принимаемых AP с направлений θ и θ.
c
ш
c
ш
Результаты моделирования в табл.1-6 (фиг.8, 9) получены для разных отношений P/Р. Как следует из этих результатов, точность пеленгации многолучевых сигналов сохраняется очень высокой даже при P/Р≥-10 дБ.
Источники информации
1. Патент US 4626859, кл. G01S 5/04.
2. Патент РФ 2150122, кл. G01S 3/14.
3. Патент РФ 2141675, кл. G01S 3/00.