Способ обнаружения линейно-частотно-модулированного сигнала с неизвестными параметрами

 

Изобретение относится к радиолокации и может быть использовано в радиолокационной технике для обнаружения маневрирующей цели в импульсно-доплеровских радиолокационных станциях. Способ основан на определении максимума модуля корреляционной суммы выборок сигнала, отраженного от цели, и опорного сигнала в пространстве параметров: частота сигнала и ее производная. Достигаемый технический результат - уменьшение вычислительных затрат. Выборки сигналов разбивают на S сегментов одинаковой длины, причем величина S должна быть не менее восьми и кратна четырем, вычисляются отсчеты дискретного преобразования Фурье сегментов, корреляционные суммы вычисляются в узлах сетки гексагонального типа с периодами 0,8/Т и 6/T² по частоте и производной частоты соответственно, где Т - длительность выборки сигнала. Корреляционные суммы вычисляются путем суммирования отсчетов дискретного преобразования Фурье сегментов, умноженных на комплексные коэффициенты, при этом от каждого сегмента в сумме учитывается один или два отсчета, комплексные коэффициенты в корреляционной сумме для узлов сетки, соответствующих первым 5S/4 номерам частотных дискретов, заранее вычисляются и сохраняются в памяти, остальные комплексные коэффициенты вычисляются путем использования соответствующих значений коэффициентов с первыми 5S/4 номерами. 4 ил.

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

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

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

Наиболее близким к предлагаемому изобретению является способ, заключающийся в том, что осуществляют аналого-цифровое преобразование (АЦП) сигнала, вычисляют корреляционные суммы выборки сигнала и опорных комплексных линейно-частотно-модулированных сигналов, фазы действительной и мнимой частей которых сдвинуты на девяносто градусов, в узлах сетки, покрывающей область возможных значений частоты сигнала и ее производной, определяют максимум модуля корреляционной суммы по узлам сетки, сравнивают значение максимума модуля корреляционной суммы с порогом, при этом корреляционные суммы вычисляют в узлах прямоугольной сетки (f_i, f'_j), где i = 0...N-1; j = 0...N²-1; f_i=i/T - частота опорного сигнала в узлах с индексом i; производная частоты опорного сигнала в узлах с индексом j; T - длительность выборки сигнала; N - размер выборки сигнала, путем умножения сигнальных отсчетов на фазовые множители соответствующие значению производной частоты f'_j, где время n-го отсчета; n = 0...N-1; t - период дискретизации выборки сигнала,
и применения быстрого преобразования Фурье (БПФ), вычисляющего значения корреляционной суммы во всех узлах с одинаковым индексом j [2].

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

Цель изобретения - снижение вычислительных ресурсов при обнаружении ЛЧМ-сигнала с широким диапазоном изменения его параметров: частоты и ее производной.

Предлагаемый способ заключается в следующем.

Осуществляют АЦП сигнала, разбивают выборку сигнала на S сегментов одинаковой длины L, причем величина S должна быть не менее восьми, кратна четырем и удовлетворять условию
f'_maxT² < S², (1)
где f'_max - максимальное абсолютное значение производной частоты сигнала,
вычисляют отсчеты дискретного преобразования Фурье (ДПФ) сегментов по алгоритму БПФ, вычисляют корреляционные суммы выборки сигнала и опорных комплексных ЛЧМ-сигналов, фазы действительной и мнимой частей которых сдвинуты на девяносто градусов, в узлах сетки гексагонального типа (f_ij, f'_ij), i = 0. . . N₁ - 1, j = 0...N₂ - 1, с периодами 0.8/T и 6/T² по частоте опорного сигнала и производной частоты соответственно, где f_ij, f'_ij - значения частоты и ее производной в узле с индексом (ij); N₁- количество отсчетов частоты; N₂ - количество отсчетов производной частоты, по формуле


s = 0...S-1, k = 0...L-1,


Ф_s = Ф(t_s),

(2)
f_s = f+ft_s,




где Z(f,f') - значение корреляционной суммы;
f - частота опорного сигнала;
f' - производная частоты опорного сигнала;
F_sk - k-й отсчет дискретного преобразования Фурье s-го сегмента;
g_s, g'_s - комплексные коэффициенты;
k_s, k'_s - номера отсчетов дискретного преобразования Фурье s-го сегмента, которые учитываются в корреляционной сумме;
x_n - значения отсчетов выборки сигнала, n=0...N-1;
Ф(t) - зависимость фазы опорного сигнала от времени,
предварительно вычисляют и сохраняют в памяти номера отсчетов k_s, k'_s и значения комплексных коэффициентов g_s, g'_s для узлов с индексами i = 0... 5S/4-1, j = 0...N₂-1, остальные номера отсчетов и комплексные коэффициенты вычисляют по формулам




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

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

1. Выборку сигнала разбивают на S сегментов одинаковой длины. Величина S должна быть не менее восьми и кратна четырем, а также удовлетворять условию f'_maxT² < S². На основе алгоритма БПФ вычисляют значения ДПФ сегментов.

2. Корреляционные суммы вычисляют в узлах сетки гексагонального типа (f_ij, f'_ij), i = 0...N₁ - 1, j = 0...N₂ - 1, с периодами d = 0.8/T, d' = 6/T² по частоте и производной частоты соответственно, путем суммирования отсчетов ДПФ сегментов, умноженных на комплексные коэффициенты, при этом от каждого сегмента в сумме учитывается один или два отсчета ДПФ сегмента.

3. Комплексные коэффициенты в корреляционной сумме для узлов с номерами i = 0...5S/4-1, j = 0...N₂-1 заранее вычисляют и сохраняют в памяти. Комплексные коэффициенты с номерами i 5S/4 вычисляют путем использования соответствующих значений коэффициентов с номерами i = 0...5S/4-1.

На фиг. 1 приведены пояснения к правилу выбора отсчетов ДПФ сегментов, учитываемых в корреляционной сумме; на фиг. 2 - линии уровня функции рассогласования ЛЧМ-сигнала в пространстве частоты и ее производной; на фиг. 3 - зависимости коэффициента выигрыша по вычислительным затратам от размеров области изменения параметров; на фиг. 4 - схема устройства обнаружения ЛЧМ-сигналов.

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

В ИД РЛС цели облучаются сигналами, представляющими из себя пачки радиоимпульсов, излучаемые с определенной частотой повторения. При отражении от воздушной цели несущая частота радиоимпульсов сдвигается на доплеровскую частоту. Частота повторения радиоимпульсов выбирается из условия превышения возможных величин доплеровских частот. В аналоговой части приемника (в каждом канале дальности) сигнал подвергается фильтрации в области спектра возможных доплеровских частот так, что на выходе импульсная структура разрушается, и сигнал можно рассматривать как протяженный одиночный радиоимпульс с длительностью, равной длительности всей пачки. При отражении от маневрирующей цели, имеющей радиальное ускорение, радиоимпульсы преобразуются в ЛЧМ-сигнал [3]. Таким образом, задача обнаружения маневрирующей цели сводится к задаче обнаружения ЛЧМ-сигналов на фоне шумов. В качестве исходных данных используется следующая выборка сигнальных отсчетов с выхода АЦП
x_n = s_n+_n, n = 0...N-1, (3)
где отсчеты сигнала;
_n - отсчеты шума, статистически независимые и распределенные по нормальному закону с нулевым средним и дисперсией ²;
U, f, f' и Ф₀ - неизвестные параметры: амплитуда, частота, производная частоты и фаза сигнала соответственно;
N - количество отсчетов;
моменты времени отсчетов;
T = tN - - длительность выборки.

Временной интервал расположен симметрично относительно нуля. Значения параметров f и f'находятся в известных интервалах [f_a, f_b] и [f'_a, f'_b] соответственно.

Обнаружение сигнала производится методом отношения правдоподобия (ОП) по правилу [4]: при Z > v принимается решение о наличии цели, иначе - о ее отсутствии. Величина Z определяется следующим выражением

где корреляционная сумма (комплексная), зависящая от вектора параметров (f, f');
u(t_n|f,f) = exp(-jФ(t_n)) - выборка опорного комплексного сигнала;
фаза опорного сигнала;
v - порог, значение которого выбирается из уровня вероятности ложных тревог.

Для вычисления корреляционных сумм выборка сигнала x_n, n = 0...N-1, разбивается на S сегментов одинаковой длины L. В s-ый сегмент входят отсчеты x_sL+1, 1 = 0...L-1, s = 0...S-1. На основе алгоритма БПФ вычисляются значения ДПФ сегментов

Используя обратное ДПФ, можно выразить отсчеты сигнала в сегментах через величины отсчетов ДПФ сегментов F_sk

Подставим выражение (6) в (4), получим

где Ф_s - значение фазы опорного сигнала в середине сегмента (в момент времени t_s, (соотношения (2))).

Величина f'_max удовлетворяет соотношению (1), которое означает, что девиация частоты сигнала на интервале одного сегмента не превышает шага частоты ДПФ сегмента, равного 1/(tL). B силу этого соотношения величину фазы опорного сигнала для s-того сегмента приближенно можно вычислять по формуле

С учетом соотношения (8) выражение (7) принимает вид

где коэффициенты g_sk определяются по соотношениям (2).

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

Если |k_s-f_stL| < 0,3, то в сумме по k в выражении (9) учитывается один член с номером k_s. В противном случае из двух номеров k_s 1 выбирается ближайший к величине f_stL номер, который обозначается как k'_s. Тогда в сумме по k в выражении (9) учитываются слагаемые с номерами k_s и k'_s. Пояснения к данному правилу показаны на фиг. 1, где приведена зависимость мгновенной частоты сигнала от времени (наклонная прямая линия). Ось времени разделена на сегменты, а ось частот - на частотные дискрета ДПФ сегментов так, что в середине прямоугольников оказываются отсчеты ДПФ (k/tL, k = 0... L-l). Затененные прямоугольники соответствуют номерам k, которые учитываются в сумме.

Показатели качества вычисления корреляционной суммы по формуле (2) будут следующими. Если в выражение для корреляционной суммы (4) подставить выборку чистого сигнала (без шума) с параметрами (f_c, f'_c), то, рассматривая зависимость модуля суммы от параметров, получим функцию рассогласования ЛЧМ-сигнала (f_c, f'_c). При больших значениях N нормированную функцию рассогласования можно представить в виде

Из выражения (11) следует, что функция (f_c, f'_c) зависит от разности векторов, Ha фиг. 2 представлены линии уровня функции рассогласования в одном квадранте (поскольку в остальных квадрантах картина симметрична). Откуда следует, что в точке (0,0) расположен глобальный максимум, а в точках (0.91/T, 5.2/T²) - максимальные боковые лепестки, уровень которых равен 0.53.

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

где числитель вычисляется при отсутствии в выборке шума, а знаменатель - при отсутствии в выборке сигнала;
(.) - среднее квадратическое отклонение (СКО), усредненное по реализациям шума.

Обозначим через q' величину отношения С/Ш, вычисленную по формуле (12), где Z(f_c, f'_c) определяется по выражению (2). Тогда потери q в отношении С/Ш равны
q = 20 log₁₀(q/q), дБ. (13)
В выражении для величины q' знаменатель вычисляется по формуле

Расчеты показателей проводились для следующего набора параметров N = 1024 и S = 8, 16; N = 768 и S = 12. Установлено, что величина L_бл не превышает 0.6, а потери в отношении С/Ш q - 1 дБ. Причем с ростом числа сегментов S характеристики улучшаются.

При обнаружении сигналов максимум модуля корреляционной суммы (выражение (4)) определяют в узлах сетки, покрывающей область параметров [f_a, f_b] [f'_a, f'_b] . Экономичное размещение узлов сетки на плоскости (f, f') следующее. Линии уровня тела неопределенности (функции рассогласования) на уровнях 0.6...1 являются симметричными овалами с осями симметрии, направленными вдоль координат f и f' (см. фиг. 2). Надлежащим масштабированием осей можно добиться равенства их длин так, что овалы будут по форме близки к окружностям. Следовательно, экономичными будут являться сетки гексагонального типа. Присвоим узлу номер, состоящий из двух индексов (ij), i = 0...N₁ - 1, j = 0...N₂ - 1. Параметры с номером (i, j) определяются по выражениям
f_ij = f_a + d(i+0.5), i = 0...N₁ - 1, j = 0,2,4,...,
f_ij = f_a+ di, i = 0...N₁ - 1, j = 1,3,5,...,
f'_ij = f'_a + d'(j + 0.5), i = 0..N₁- 1, j = 0...N₂ - 1.

N₁ = (f_b - f_a)/d + 0.5, N₂ = (f'_b - f'_a)/d'. (15)
Величины d и d' определяются из условия пересечения тел неопределенности, размещенных в трех соседних узлах с номерами (i, j), (i+1, j) (i, j + 1) по уровню 0.63. Расчеты показали, что d = 0.8/T, d' = 6/T².

Размеры области параметров, которая приходится на один узел, равны d d'. Это позволяет оценить количество узлов N_n, покрывающих всю область [f_a, f_b] [f'_a, f'_b],
N_n = V₁V₂/4.8, V₁ = (f_b - f_a)T, V₂ = (f'_b - f'_b)T². (16)
Рассмотрим узлы с номерами (i, j) и (i + 5k, j), где k = 1,2,3,.... Частоты, соответствующие этим узлам, согласно выражению (15) отличаются на величину 4k/T. Поскольку количество сегментов S выбрано кратным четырем, S = 8, 12, 16, ... , то при k = S/4 частоты будут отличаться на величину S/T, равную 1/tL - частотному дискрету ДПФ сегмента. Можно показать, что номера отсчетов k_s и k'_s, соответствующие этим узлам, при одинаковых s отличаются на единицу, а коэффициенты g_s и g'_s отличаются на постоянный фазовый множитель, который можно опустить. Такая закономерность вытекает непосредственно из фиг. 1, где переходу из узла (i, j) в узел (i + 5S/4, j) соответствует сдвиг графика мгновенной частоты сигнала вверх на величину частотного дискрета.

Таким образом, если заранее вычислены и сохранены значения номеров отсчетов k_s, k'_s и коэффициентов g_s, g'_s для узлов с номерами i=0...5S/4-1, j=0...N₂-1, то эти величины для остальных узлов вычислять не требуется.

Сравнительная оценка вычислительных затрат предлагаемого способа и способа прототипа. Затраты вычислительных ресурсов измеряются в элементарных операциях (ЭО) типа умножения или сложения над действительными числами.

Так, для получения одного значения корреляционной суммы (без учета предварительных вычислений ДПФ сегментов), если заранее вычислены и сохранены значения номеров отсчетов k_s, k'_s и коэффициентов g_s, g'_s для узлов с номерами i=0...5S/4-1, j=0...N₁-1, требуется 16S ЭО.

На вычисление БПФ выборки размера N требуется 2Nlog₂(N) действительных сложений и 3Nlog₂(N) умножений [5]. В сумме это составляет 5Nlog₂(N) ЭО. Следовательно, если выборка размера N разделена на S одинаковых сегментов, то вычисление БПФ сегментов требует 5Nlog₂(N/S) ЭО.

Суммарные вычислительные затраты С_об1, необходимые для обнаружения сигнала в области [f_a, f_b] [f'_a, f'_b], будут определяться затратами на вычисление ДПФ сегментов и корреляционных сумм, т.е.

С_об1 = 5Nlog₂(N/S) + 3.3SV₁V₂. (17)
Определим вычислительные затраты способа прототипа, при котором корреляционные суммы вычисляют в узлах прямоугольной сетки с периодами 1/T и 2/T² по частоте и производной частоты соответственно. Количество узлов N_n, накрывающих область параметров [f_a, f_b] [f'_a, f'_b] для этого способа равно
N_n = V₁V₂/2, (18)
что в 2.4 раза больше величины, соответствующей предлагаемому способу (см. выражения (16)). Обозначим слоем сетки узлов их совокупность с одинаковым значением f'. БПФ позволяет вычислять сразу все значения в слое узлов. Кроме собственно БПФ необходимо производить умножения значений выборки сигнала на комплексные экспоненты (фазовые множители). Таким образом
С_об2 = V₂N(5log₂(N)+18)/2. (19)
Коэффициент выигрыша по вычислительным затратам К_вз, равный отношению соответствующих величин С_об1 и С_об2, будет определяться следующим образом

На фиг. 3 представлены зависимости К_вз от величин V₁ и V₂ при N = 1024, S = 8. Например, при V₁ = 200, V₂ = 30 величина К_вз равна 5.4.

Техническая реализация способа возможна на основе устройства, представленного на фиг. 4. Устройство состоит из блока АЦП 1, S блоков регистров сдвига 2, S блоков БПФ 3, S блоков регистров 4, S коммутаторов 5, постоянных запоминающих устройств (ПЗУ) 6 и 8, 2S умножителей комплексных чисел 7, сумматора комплексных чисел 9, вычислителя модуля комплексного числа 10, регистра 11, счетчика 12, вычислителя максимума из двух чисел 13, ключа 14, регистра 15, порогового устройства 16.

Устройство работает следующим образом. На вход блока АЦП 1 поступает аналоговый сигнал. Цифровые отсчеты сигнала последовательно накапливаются в S блоках регистров сдвига 2, соответствующих S сегментам выборки. После накопления выборки отсчеты сигнала поступают на входы S блоков БПФ 3, которые параллельно вычисляют комплексные отсчеты ДПФ сегментов и засылают их в S блоков регистров. Далее производится N_n циклов вычисления модулей корреляционной суммы по количеству узлов сетки, покрывающей область значений параметров. В каждом таком цикле коммутаторы 5 подключаются к тем регистрам 4, которые содержат отсчеты ДПФ сегментов, используемые при вычислении значения корреляционной суммы текущего узла. Номера этих регистров хранятся в ПЗУ 6. С выходов коммутаторов 5 отсчеты ДПФ сегментов поступают на первые входы умножителей 7, на вторые входы которых из ПЗУ 8 поступают комплексные коэффициенты. Отсчеты ДПФ сегментов, умноженные на комплексные коэффициенты, суммируются в блоке 9, на выходе которого формируется комплексное значение корреляционной суммы текущего узла. В блоке 10 вычисляется модуль этой величины. В регистре 11 содержится значение текущего максимума модуля корреляционной суммы (по пройденным узлам), в первом цикле в регистре 11 содержится нулевое значение. В счетчике 12 содержится номер текущего цикла. На первый вход вычислителя максимума из двух чисел 13 поступает значение модуля корреляционной суммы текущего узла, на второй вход - значение текущего максимума модуля корреляционной суммы (по пройденным узлам). На первом выходе вычислителя 13 формируется единица, если значение на первом его входе превышает значение на втором входе, в противном случае формируется ноль. На втором выходе вычислителя 13 формируется значение текущего максимума модуля корреляционной суммы. С первого выхода вычислителя 13 код поступает на управляющий вход ключа 14. На первый вход ключа 14 подается номер цикла. Таким образом, в регистр 15 записывается номер узла, которому соответствует текущее максимальное значение модуля корреляционной суммы. После завершения всех Nn циклов на втором выходе устройства 13 присутствует максимум модуля корреляционной суммы по всем узлам, значение которого поступает на первый вход порогового устройства 16. На второй вход порогового устройства 16 поступает значение порога. На выходе порогового устройства 16 формируется единичное значение в случае обнаружения сигнала и нулевое значение - в противном случае. В случае обнаружения сигнала в регистре 15 содержится номер узла, соответствующего сигналу.

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

ИСТОЧНИКИ ИНФОРМАЦИИ
1. Обработка сигналов в многоканальных РЛС./Под. ред. А.П. Лукошкина. М. : Радио и связь, 1983, стр. 307, рис. 12.18.

2. Кузьменков В.Ю., Логинов В.М. Способы и устройства совместного измерения радиальной скорости и радиального ускорения.// Радиотехника и электроника, 1997, Т.42, N 12, С. 1465...1475.

3. Многофункциональные радиоэлектронные комплексы истребителей./ Под. ред. Г.С.Кондратенкова. - М.: Военное издательство, 1994.

4. Ширман Я.Д., Манжос В.Н. Теория и техника обработки радиолокационной информации на фоне помех. - М.: Радио и связь, 1981.

5. Блейхут Р. Быстрые алгоритмы цифровой обработки сигналов. - М.: Мир, 1989.

Формула изобретения

Способ обнаружения линейно-частотно-модулированного сигнала, заключающийся в том, что осуществляют аналого-цифровое преобразование сигнала, вычисляют корреляционные суммы выборки сигнала и опорных комплексных линейно-частотно-модулированных сигналов, фазы действительной и мнимой частей которых сдвинуты на девяносто градусов, в узлах сетки, покрывающей область возможных значений параметров: частоты сигнала и ее производной, определяют максимум модуля корреляционной суммы по узлам сетки, сравнивают значение максимума модуля корреляционной суммы с порогом, отличающийся тем, что при вычислении корреляционных сумм разбивают выборку сигнала на S сегментов одинаковой длины L, причем величина S должна быть не менее восьми, кратна четырем и удовлетворять условию
f'_maxT² < S²,
где f'_max - максимальное абсолютное значение производной частоты сигнала;
T - длительность выборки сигнала,
вычисляют отсчеты дискретного преобразования Фурье сегментов на основе алгоритма быстрого преобразования Фурье, вычисляют корреляционные суммы в узлах сетки гексагонального типа (f_ij, f'_ij), i = N₁ - 1, j=0 ... N₂ - 1, с периодами 0,8/T и 6/T² по частоте опорного сигнала и производной частоты соответственно, где f_ij, f'_ij - значения частоты и ее производной в узле с индексом (ij), N₁ - количество отсчетов частоты, N₂ - количество отсчетов производной частоты, по формуле


s=0...S-1, K = 0...L-1,


Ф_s=Ф(t_s),


f_s= f+ft_s,




где Z (f, f') - значение корреляционной суммы;
f - частота опорного сигнала;
f' - производная частоты опорного сигнала;
F_sk - k-й отсчет дискретного преобразования Фурье s-го сегмента;
g_s, g'_s - комплексные коэффициенты;
k_s, k'_s - номера отсчетов дискретного преобразования Фурье s-го сегмента, которые учитываются в корреляционной сумме;
x_n - значения отсчетов выборки сигнала, n = 0 ... N-1;
Ф(t) - зависимость фазы опорного сигнала от времени;
t - период дискретизации выборки сигнала,
предварительно вычисляют и сохраняют в памяти номера отсчетов k_s, k'_s и значения комплексных коэффициентов g_s, g'_s - для узлов с индексами i = 0 ... 5S/4-1, j= 0. . .N₂-1, остальные номера отсчетов и комплексные коэффициенты вычисляют по формулам




где после вертикальной черты указаны индексы узлов.

РИСУНКИ

Рисунок 1, Рисунок 2, Рисунок 3, Рисунок 4