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

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

Способ поиска шумоподобных сигналов с минимальной частотной манипуляцией, заключающийся в том, что осуществляют разделение входного шумоподобного сигнала на квадратурные составляющие путем его перемножения с опорными гармоническими сигналами несущей частоты, сдвинутыми относительно друг друга на фазовый угол /2, вычисляют в каждом квадратурном канале на каждом цикле поиска функцию взаимной корреляции входного и опорного сигналов, выделяют модуль функции взаимной корреляции, отличающийся тем, что на каждом цикле поиска осуществляют накопление отсчетов квадратурных составляющих входного шумоподобного сигнала на тактовых интервалах, равных половине длительности элемента кодовой последовательности, запоминают на время, равное n периодов повторения входного шумоподобного сигнала, результаты накопления квадратурных составляющих входного шумоподобного сигнала, формируют 4N отсчетов квадратурных компонент входного шумоподобного сигнала z_1k и z_2k , формируют 4N отсчетов квадратурных опорных видеосигналов I _k=cos (t_k) и Q_k=sin (t_k) при фиксированном времени задержки элемента кодовой последовательности t_k, производят знаковую двухступенчатую аппроксимацию формы квадратурных опорных видеосигналов I_k и Q_k, причем отсчеты квадратурных опорных видеосигналов при знаковой двухступенчатой аппроксимации равны I_k=(±a₁, ±a₂) и Q _k=(±a₂, ±a₁), где a ₁ и a₂ - весовые коэффициенты, на каждом цикле поиска в каждом квадратурном канале осуществляют перемножение отсчетов квадратурных компонент входного шумоподобного сигнала z_1k и z_2k с отсчетами аппроксимированных квадратурных опорных видеосигналов I_k=(±a ₁, ±a₂) и Q_k=(±a₂ , ±a₁), попарно объединяют результаты перемножения квадратурных компонент z_1k и z_2k с отсчетами аппроксимированных квадратурных опорных видеосигналов I_k =(±a₁, ±a₂) и Q_k(±a ₂, ±a₁), вычисляют 2N отсчетов квадратурных составляющих функции взаимной корреляции путем интегрирования объединенных результатов перемножения при фиксированной задержке элемента кодовой последовательности t_k, выделяют N значений модуля функции взаимной корреляции, определяют значение задержки элемента кодовой последовательности, соответствующее максимальному значению модуля функции взаимной корреляции, используют найденное значение задержки кодовой последовательности для установки генератора кода в состояние синхронизма с принятым шумоподобным сигналом с минимальной частотной манипуляцией.

Описание изобретения к патенту

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

Известен способ параллельного поиска шумоподобных сигналов по времени запаздывания [Радиотехнические системы / Ю.П.Гришин, В.П.Ипатов, Ю.М.Казаринов и др.; под ред. Ю.М.Казаринова. - М.: высш. шк., 1990, с.64 (рис.3.16), с.99 (рис.4.6)], заключающийся в перемножении принятого сигнала с N парами опорных шумоподобных сигналов, являющихся квадратурными копиями принятого сигнала для N дискретных значений времени запаздывания, интегрировании результатов перемножения в 2N квадратурных каналах на интервале, равном длительности шумоподобного сигнала, выделении N значений модуля взаимной корреляционной функции (ВКФ) и выборе в качестве оценки времени запаздывания значения задержки опорного сигнала в канале с максимальным значением ВКФ.

Способ поиска обеспечивает потенциально достижимую помехоустойчивость, минимально возможное время поиска, однако трудно реализуем при числе каналов N>>1.

Известен способ быстрого поиска шумоподобного сигнала (ШПС) [патент RU № 2206180, МКИ H04L 7/10, опубл. 10.06.2003], заключающийся в последовательной оценке символов псевдослучайной последовательности (ПСП), позволяющий по любому неискаженному сегменту ПСП длиной в m символов синтезировать в приемном устройстве сигнал с требуемой задержкой. В регистр опорного генератора ШПС корреляционного накопителя записывают m оценок принятых двоичных символов псевдослучайной последовательности и затем проверяют принадлежность данной m-значной комбинации символов искомому шумоподобному сигналу. Если принятая m-значная комбинация символов не принадлежит шумоподобному сигналу (ложная тревога), то регистр обнуляется и вновь заполняется очередными оценками символов псевдослучайной последовательности искомого шумоподобного сигнала для последующей проверки с помощью корреляционного накопителя.

Способ быстрого поиска обеспечивает сокращение времени поиска ШПС, но только при достаточно большом отношении сигнал/шум для элемента ШПС.

Известен способ циклического поиска шумоподобных сигналов [Е.П.Петров, Д.Е.Прозоров, И.Е.Петров, А.В.Смирнов. Быстрый поиск шумоподобных сигналов / Успехи современной радиоэлектроники, 2008, № 8, с.47-48], заключающийся в разделении входного сигнала на квадратурные составляющие путем перемножения входного сигнала с опорными гармоническими сигналами несущей частоты, сдвинутыми друг относительно друга на фазовый угол /2. Формируют в каждом квадратурном канале на каждом цикле поиска функции взаимной корреляции принятого и опорного шумоподобных сигналов. Выделяют на каждом цикле модуль функции взаимной корреляции и производят сравнение значения модуля функции взаимной корреляции с порогом обнаружения. По результатам сравнения принимают решение о наличии или отсутствии наличия шумоподобного сигнала.

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

Наиболее близким техническим решением к заявляемому изобретению является способ поиска шумоподобных сигналов с минимальной частотной манипуляцией [патент RU № 2353064, МКИ H04L 27/14, опубл. 20.04.2009], заключающийся в том, что на каждом цикле поиска осуществляют разделение входного шумоподобного сигнала на квадратурные составляющие, дискретизацию, оцифровку и интегрирование квадратурных составляющих на тактовых интервалах, в m раз меньших длительности элемента кодовой последовательности. Разделение входного шумоподобного сигнала осуществляют путем его перемножения с опорными гармоническими сигналами несущей частоты, сдвинутыми друг относительно друга на фазовый угол /2. Результаты интегрирования запоминают на время, равное длительности элемента кодовой последовательности, в течение которого их перемножают с отсчетами соответствующих опорных шумоподобных сигналов, формируемых с частотой, в М раз большей тактовой частоты кодовой последовательности входного сигнала. При этом за время, равное периоду повторения кодовой последовательности, формируют М значений модуля функции взаимной корреляции, которые затем накапливают на интервале, равном фиксированному числу периодов повторения кодовой последовательности. Решение о значении времени запаздывания входного сигнала принимают путем выбора максимального из М значений накопленных модулей функции взаимной корреляции, запоминания максимального значения, его адреса и номера цикла и повторения процедуры поиска фиксированное число раз со сдвигом на М элементов последовательности при переходе на каждый последующий цикл поиска.

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

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

Поставленная задача решается тем, что в способе поиска шумоподобных сигналов с минимальной частотной манипуляцией, заключающемся в том, что осуществляется разделение входного сигнала на квадратурные составляющие путем перемножения входного сигнала с опорными гармоническими сигналами несущей частоты, сдвинутыми друг относительно друга на фазовый угол /2, формирование в каждом квадратурном канале на каждом цикле поиска значения функции взаимной корреляции входного и опорного шумоподобных сигналов, выделение модуля функции взаимной корреляции, согласно изобретению на каждом цикле поиска осуществляют накопление отсчетов квадратурных составляющих входного шумоподобного сигнала на тактовых интервалах, равных половине длительности элемента кодовой последовательности, запоминают на время, равное n периодов повторения входного шумоподобного сигнала, результаты накопления квадратурных составляющих входного шумоподобного сигнала, формируют 4N отсчетов квадратурных компонент входного шумоподобного сигнала z_1k и z_2k, формируют 4N отсчетов квадратурных опорных видеосигналов I_k=cos (t_k) и Q_k=sin (t_k) при фиксированном времени задержки элемента кодовой последовательности t_k, производят знаковую двухступенчатую аппроксимацию формы квадратурных опорных видеосигналов I_k и Q_k, причем отсчеты квадратурных опорных видеосигналов при знаковой двухступенчатой аппроксимации равны I_k=(±а₁, ±а₂) и Q _k=(±a₂, ±a₁), где а ₁ и а₂ - весовые коэффициенты, на каждом цикле поиска в каждом квадратурном канале осуществляют перемножение отсчетов квадратурных компонент входного шумоподобного сигнала z_1k и z_2k с отсчетами аппроксимированных квадратурных опорных видеосигналов I_k=(±а ₁, ±а₂) и Q_k=(±a₂ , ±a₁), попарно объединяют результаты перемножения квадратурных компонент z_1k и z_2k с отсчетами аппроксимированных квадратурных опорных видеосигналов I_k =(±а₁, ±а₂) и Q_k=(±a ₂, ±a₁), вычисляют 2N отсчетов квадратурных составляющих функции взаимной корреляции путем интегрирования объединенных результатов перемножения при фиксированной задержке элемента кодовой последовательности t_k, выделяют N значений модуля функции взаимной корреляции, определяют значение задержки элемента кодовой последовательности, соответствующее максимальному значению модуля функции взаимной корреляции, используют найденное значение задержки кодовой последовательности для установки генератора кода в состояние синхронизма с принятым шумоподобным сигналом с минимальной частотной манипуляцией.

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

Устройство для реализации способа поиска шумоподобных сигналов с минимальной частотной манипуляцией содержит блок поэлементной обработки 1, соединенный с блоком 6 декодирования. Блок поэлементной обработки 1 включает в себя входные перемножители 2₁ и 2₂, сигнальные входы которых объединены, а опорные входы подключены к соответствующему выходу опорного генератора 5. К выходам первого и второго перемножителей 2₁ и 2₂ подключены сигнальные входы соответственно первого и второго накапливающих сумматоров 3₁ и 3₂ , выходы которых соединены с информационными входами первого и второго оперативных запоминающих устройств 4₁ и 4₂ соответственно. Тактируемые входы оперативных запоминающих устройств 4₁ и 4₂ соединены с соответствующими выходами тактового генератора 12 (на фиг.1 не показано). Выходы запоминающих устройств 4₁ и 4₂ соединены соответственно с сигнальными входами третьего и четвертого накапливающих сумматоров 3₃ и 3₄, опорные входы которых подключены к соответствующим информационным выходам блока хранения данных 16.

Блок 6 декодирования (см. фиг.2) состоит из седьмого и восьмого, девятого и десятого перемножителей 2 ₇ и 2₈, 2₉ и 2₁₀, выходы которых попарно объединены через первый и второй вычитатели 17 ₁ и 17₂ соответственно, а также из одиннадцатого и двенадцатого, тринадцатого и четырнадцатого перемножителей 2₁₁ и 2₁₂, 2₁₃ и 2₁₄ , выходы которых попарно соединены через третий и четвертый сумматоры 7₃ и 7₄ соответственно. Объединенные сигнальные входы перемножителей 2₇, 2₉, 2₁₁ , 2₁₃ и 2₈, 2₁₀, 2₁₂ , 2₁₄ подключены к выходу третьего и четвертого накапливающего сумматора 3₃ и 3₄ соответственно. Опорные входы перемножителей 2₇ и 2₁₂, 2₉ и 2₁₄, 2₈ и 2₁₁, 2₁₀ и 2₁₃ попарно объединены и подключены к соответствующим выходам блока 15 формирования кодов. Выходы первого и второго вычитателей 17₁ и 17₂ соединены соответственно с входами пятого и шестого накапливающих сумматоров 3₅ и 3₆, а выходы третьего и четвертого сумматоров 7 ₃ и 7₄ подключены к входам соответственно седьмого и восьмого накапливающих сумматоров 3₇ и 3₈ .

Выходы пятого, шестого, седьмого и восьмого накапливающих сумматоров 3₅, 3₆, 3 ₇ и 3₈ подключены к сигнальным входам соответственно третьего, четвертого, пятого и шестого перемножителей 2₃ , 2₄, 2₅ и 2₆. Выходы третьего и четвертого, пятого и шестого перемножителей 2₃ и 2₄, 2₅ и 2₆ объединены соответственно через первый и второй сумматоры 7₁ и 7₂ , к выходам которых подключены соответствующие входы вычислителя 8 модуля. К выходу вычислителя 8 модуля подключены последовательно соединенные решающий блок 9 и блок 10 управления, выход которого подключен к сигнальному входу элемента 11 управляемой задержки. Тактируемый вход элемента 11 управляемой задержки подключен к соответствующему выходу тактового генератора 12, а выход элемента 11 управляемой задержки - к входу генератора 13 кодов. Выход генератора 13 кодов, являющийся выходом устройства поиска (см. фиг.1), соединен с входом синтезатора 14 отсчетов, первый и второй выходы которого подключены к соответствующим сигнальным входам блока 15 формирования кодов. Попарно объединенные опорные входы третьего и пятого, четвертого и шестого перемножителей 2 ₃ и 2₅, 2₄ и 2₆, соединенные с соответствующими сигнальными входами блока 15 формирования кодов, подключены к соответствующим информационным выходам блока хранения данных 16. Опорный вход блока хранения данных 16 соединен с соответствующим выходом тактового генератора 12 (на фиг.1. не показано).

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

На вход устройства поиска (фиг.1) поступают с шагом дискретизации T_д отсчеты принимаемого периодического шумоподобного сигнала (ШПС) с минимальной частотной манипуляцией (индекс k не указан):

где ₀ - средняя частота;

- начальная фаза (амплитуда полагается равной единице);

(t) - функция, определяющая закон угловой модуляции;

d(t) - двоичный сигнал, соответствующий кодовой псевдослучайной последовательности (ПСП) d₀, d₁, , d_N-1;

rect(t) - прямоугольная функция (импульс единичной амплитуды и длительности Т);

N - длина ПСП, определяющая период T_п=NT повторения ШПС,

М - число отсчетов на периоде.

Входные перемножители 2₁ и 2₂ осуществляют перемножение отсчетов сигнала (1) с отсчетами опорных сигналов несущей частоты cos ₀t_ij и sin ₀t_ij, вырабатываемыми опорным генератором 5. На выходах перемножителей 2₁ и 2₂ образуются отсчеты видеочастотных составляющих входного сигнала соответственно:

а также составляющие удвоенной частоты 2 ₀, которые отфильтровываются последующим трактом обработки.

Отсчеты видеочастотных составляющих (2) входного сигнала x_jk и y_jk поступают соответственно на сигнальные входы накапливающих сумматоров 3 ₁ и 3₂, на выходах которых формируются величины соответственно:

Суммирование по i в (3) ведется на интервалах интегрирования, равных половине длительности элемента ПСП. При этом число отсчетов на каждом k-м интервале интегрирования равно m=Т/(2Т_д). Накопленные величины поступают в запоминающие устройства (ЗУ) 4₁ и 4₂, где запоминаются на время наблюдения, равное n периоду T_п . По управляющему сигналу запоминающие устройства 4₁ и 4₂ соответственно отправляют отсчеты X_jk и Y_jk на сигнальные входы третьего и четвертого накапливающих сумматоров 3₃ и 3₄, на опорные входы которых поступают величины u_1j и u_2j соответственно.

В накапливающих сумматорах 3₃ и 3 ₄ осуществляется когерентное межпериодное накопление соответствующих квадратурных компонент (3) по следующему алгоритму:

Результаты накопления (4), полученные за n периодов повторения ШПС, используются для последовательного вычисления N дискретных значений модуля взаимной корреляционной функции (ВКФ) в режиме постобработки на интервале времени, равном:

nT_п<t<nT_п+T_доп ,

где Т_доп - дополнительное время, затраченное на постобработку.

На каждом интервале поиска в синтезаторе 14 отсчетов формируют отсчеты квадратурных опорных видеочастотных сигналов I_k=cos (t_k) и Q_k=sin (t_k) при фиксированном значении задержки ПСП, равном где k=1, 2, , 2N. Отсчеты квадратурных опорных видеочастотных сигналов I_k=cos (t_k) и Q_k=sin (t_k) поступают на соответствующие сигнальные входы блока 15 формирования кодов, на соответствующие опорные входы которого поступают значения весовых коэффициентов a ₁=cos( /8) и a₂=sin( /8), хранящиеся в блоке хранения данных 16. В результате знаковой двухступенчатой аппроксимации формы квадратурных опорных видеочастотных сигналов I_k=cos (t_k) и Q_k=sin (t_k), осуществляемой в блоке 15 формирования кодов, отсчеты квадратурных опорных видеочастотных сигналов принимают следующие значения: I_k=(±a₁, ±a ₂) и Q_k=(±a₂, ±а₁ ). На каждом интервале поиска на основе полученных в результате аппроксимации отсчетов квадратурных опорных видеочастотных сигналов I_k=(±а₁, ±a₂) и Q _k(±a₂, ±a₁) в блоке 15 формирования кодов получают элементы кодов b_1k и b _2k, c_1k и c_2k в соответствии со следующим алгоритмом:

где sign(x) - знаковая функция.

В соответствии с вышеприведенным алгоритмом (5) элементы кодов b_1k, b_2k, c_1k, c_2k принимают значения 0, 1 или -1.

Квадратурные составляющие входного сигнала z_1k и z_2k поступают соответственно (см. фиг.2) на объединенные сигнальные входы перемножителей 2₇, 2₉, 2₁₁ , 2₁₃ и 2₈, 2₁₀, 2₁₂ , 2₁₄ блока 6 декодирования. На попарно объединенные опорные входы перемножителей 2₇ и 2₁₂, 2₈ и 2₁₁, 2₉ и 2₁₄ , 2₁₀ и 2₁₃ подаются соответственно элементы кодов b_1k, c_1k, b_2k, c_2k , формируемые блоком 15 формирования кодов. Сигналы с выходов перемножителей 2₇ и 2₈, 2₉ и 2₁₀, 2₁₁ и 2₁₂, 2₁₃ и 2₁₄ попарно объединяются соответственно в вычитателях 17₁ и 17₂ и в сумматорах 7₃ и 7₄, образуя соответствующие квадратурные составляющие z'_1k, z''_1k, z'_2k , z''_2k.

С выходов вычитателей 17₁, 17₂ и сумматоров 7₃, 7 ₄ блока 6 декодирования квадратурные компоненты z' _1k, z''_1k, z'_2k, z'' _2k поступают соответственно на пятый, шестой, седьмой и восьмой накапливающие сумматоры 3₅, 3₆, 3₇, 3₈, формирующие величины:

Результаты накопления (6) поступают на сигнальные входы соответствующих перемножителей 2₃ , 2₄, 2₅, 2₆. На опорные входы соответствующих перемножителей 2₃ и 2₅, 2₄ и 2₆ соответственно подаются весовые коэффициенты а₁ и а₂, хранящиеся в блоке хранения данных 16. «Взвешенные» результаты перемножения объединяются в сумматорах 7₁ и 7₂ соответственно, образуя соответствующие квадратурные составляющие ВКФ входного и опорного сигналов:

Совокупность выражений (6) и (7) эквивалентна вычислению квадратурных составляющих ВКФ входного и опорного сигналов при знаковой двухступенчатой аппроксимации квадратурных опорных видеочастотных сигналов:

где I_k=(±a₁ , ±a₂) и Q_k=(±a₂, ±a₁) - отсчеты квадратурных опорных видеочастотных сигналов при знаковой двухступенчатой аппроксимации формы квадратурных опорных видеочастотных сигналов;

z_1k и z_2k - отсчеты квадратурных составляющих входного шумоподобного сигнала.

На каждом интервале поиска при фиксированной задержке ПСП _l=lT, где l=0, 1, , (N-1) - номер квадратурного канала, вычислитель 8 модуля формирует значение модуля ВКФ входного и опорного сигналов:

Решающий блок 9 определяет максимальное значение модуля ВКФ входного и опорного сигналов:

которое используется в блоке 10 управления, при формировании кода задержки опорной ПСП

где µ - адрес канала с Z_µ =max.

Код задержки опорной ПСП подается на управляющий вход элемента 11 управляемой задержки, на опорный вход которого поступают с частотой f_т тактовые импульсы, вырабатываемые тактовым генератором 12. Указанный код определяет оценку задержки принятого ШПС и используется для установки генератора 13 кода в состояние синхронизма с принятым ШПС с точностью не хуже ±T/2 (при условии, что аномальные ошибки отсутствуют).

Качественные показатели описанного способа поиска ШПС характеризуются вероятностью Р_ош аномальных ошибок, превышающих значение Т/2 (по абсолютной величине), а также временем поиска t_поиск . При длине кодовой ПСП N>>1 задачу поиска ШПС по времени запаздывания можно свести к задаче распознавания N ортогональных сигналов, применительно к которой вероятность ошибки можно оценить как [5]:

где Ф(х) - интеграл вероятности, q - отношение сигнал/шум на выходе "синхронного" канала (при относительной задержке принятого и опорного сигналов =0).

Проигрыш в отношении сигнал/шум из-за "равновесовой" поэлементной обработки по сравнению с оптимальной корреляционной обработкой составляет около 0,2 дБ [В.Н.Бондаренко / Оптимальный алгоритм поиска шумоподобного сигнала с минимальной частотной манипуляцией. - М., «Радиотехника и электроника», 2008, т.53, № 2. С.222-229], т.е. предлагаемый способ поиска обеспечивает помехоустойчивость, близкую к потенциально достижимой.

Время поиска для предлагаемого способа при тактовой частоте процессора f_т=1 ГТц, числе периодов накопления n=25, длительности ШПС T_п=40 мс превышает время наблюдения на величину дополнительного времени T_доп N²Т_слож, затраченного на постобработку результатов. При поиске в режиме реального времени с помощью одноканального устройства, реализующего циклический способ поиска, время поиска для тех же условий составляет t_поиск=nT _пN 4,5 мин.

Таким образом, предлагаемый способ поиска ШПС с минимальной частотной манипуляцией позволяет значительно сократить время поиска по сравнению с прототипом (более чем в 270 раз) при пренебрежимо малых потерях в помехоустойчивости (менее 0.2 дБ). В этом заключается технико-экономический эффект по сравнению с известными способами поиска шумоподобных сигналов.