Системы управления,связи и безопасности №1. 2020
Systems of Control, Communication and Security ISSN 2410-9916
УДК 621.396
Обнаружение слабых сигналов методом взаимной корреляции с компенсацией фазовых нестабильностей при радиоконтроле частотного ресурса спутниковых систем связи
Кулакова В. И.
Постановка задачи: космические аппараты (КА) широко используются для ретрансляции сигналов Земных станций (ЗС) в целях их передачи на большие расстояния. Обнаружение всех ЗС, использующих спутниковые системы связи, необходимо для контроля частотного ресурса КА связи и его защиты от несанкционированного доступа. Для обнаружения сигналов ЗС используется метод взаимной корреляции. Сложность заключается в том, что сигналы несанкционированных ЗС зачастую наблюдаются станцией мониторинга при значительном отрицательном отношении сигнал/шум (ОСШ) (менее минус 20 дБ). В связи с этим для получения пика огибающей взаимной корреляционной функции (ВКФ) требуется длительное время когерентного накопления обнаруживаемого сигнала (десятки секунд). При этом важным фактором становятся фазовые нестабильности, вносимые в сигнал приемо-передающими трактами в составе задействованных КА и наземного радиоприемного устройства. Целью работы является исследование возможности обнаружения слабых сигналов ЗС в реальных условиях применения. При этом для компенсации фазовых нестабильностей в каналах используется сигнал реперной ЗС. Результат: Приведены формулы расчета ОСШ на выходе коррелятора в зависимости от ОСШ в каждом канале, полосы приема и длительности интервала интегрирования. Для обнаружения пика огибающей ВКФ на фоне выходных шумов необходимо обеспечить выходное ОСШ не менее 10...12 дБ. Это означает, что если сигнал ЗС принимается при ОСШ около -30 дБ в каждом канале, то потребуется 16 с обработки в полосе частот 1 МГц. При вычислении взаимной корреляции на таких длительных интервалах необходимо учитывать присутствие в каналах частотных отстроек и фазовых нестабильностей. При этом фазовая нестабильность на входе коррелятора является похожей для всех сигналов ЗС, которые ретранслируется данным набором КА, в то время как сдвиг по частоте для каждой ЗС будет отличным из-за эффекта Доплера. На основании этого предлагается использовать эквалайзер для компенсации фазовых нестабильностей, который настраивается по известному сигналу реперной ЗС. После применения эквалайзера, алгоритм обнаружения всех ЗС сводится к расчету ВКФ сигналов, принятых в соседних каналах, для разных значений отстроек по частоте и сравнению ее огибающей с заданным порогом. Приведены результаты эксперимента, в ходе которого были приняты и зарегистрированы сигналы с двух КА. В данных присутствовал сигнал реперной ЗС, у которого расширение спектра происходило за счет использования псевдослучайной последовательности. Длительность записи сигнала составляла 19 с. Оцененные фазовые нестабильности по амплитуде достигали 300°, что разрушало корреляцию сигналов в разных каналах. Показано, что при использовании эквалайзера, корректном выборе шага поиска максимума огибающей ВКФ по времени и частоте, экспериментальное выходное ОСШ с высокой точностью согласуются с теоретическими расчетами. Помимо сигнала репер-ной ЗС было обнаружено еще три ЗС, которые отличались от сигнала реперной ЗС отстройкой по частоте.
Ключевые слова: спутниковые системы связи, обнаружение сигналов, метод взаимной корреляции, отношение сигнал/шум.
Библиографическая ссылка на статью:
Кулакова В. И. Обнаружение слабых сигналов методом взаимной корреляции с компенсацией фазовых нестабильностей при радиоконтроле частотного ресурса спутниковых систем связи // Системы управления, связи и безопасности. 2020. № 1. С. 33-48. DOI: 10.24411/2410-9916-202010102.
Reference for citation:
Kulakova V. I. Detection of Weak Signals by Cross-Correlation with Phase Instabilities Compensation for Spectrum Monitoring in Satellite Communication Systems. Systems of Control, Communication and Security, 2020, no. 1, pp. 33-48. DOI: 10.24411/2410-9916-2020-10102 (in Russian).
DOI: 10.24411/2410-9916-2020-10102
Systems of Control, Communication and Security
ISSN 2410-9916
Введение
Спутниковые системы связи широко используются для ретрансляции сигналов земных станций (ЗС) в целях их передачи на большие расстояния [1]. Информация о всех передатчиках, использующих спутниковые системы связи, необходима для контроля частотного ресурса телекоммуникационного космического аппарата (КА) и его защиты от несанкционированного доступа [2].
Одна из распространенных схем контроля частотного ресурса КА представлена на рис. 1 [2-6]. Принцип ее работы основан на том, что сигналы нескольких ЗС одновременно принимаются двумя (или более) КА, которые затем ретранслируют их в другое место на поверхности Земли, где они регистрируется станцией мониторинга. Если местоположение ЗС и классы излучаемых ею сигналов известны, то ее называют реперной станцией (РС) [3, 4]. На станции мониторинга вычисляется корреляция сигналов, принятых с двух КА, для различных временных и частотных сдвигов, до тех пор, пока не появится пик в огибающей взаимной корреляционной функции (ВКФ). Временные и частотные сдвиги, при которых получен данный пик огибающей ВКФ, могут быть использованы для поиска местоположения соответствующей им ЗС.
КА№1 КА№2
Рис. 1. Принцип работы схемы контроля спутниковой системы связи
На практике возникает сложность с использованием метода ВКФ, поскольку коррелируемые сигналы являются слабыми, т.е. наблюдаемыми при большом отрицательном значении отношения сигнал/шум (ОСШ) по мощности, которое может достигать значений -35...-20 дБ [2, 3]. В связи с этим для получения пика огибающей ВКФ требуется длительное время когерентного накопления обнаруживаемого сигнала [7]. При этом важным фактором становятся фазовые нестабильности, вносимые в сигнал приемо-передающими трактами в составе задействованных КА и наземного радиоприемного устройства [3, 8]. Для компенсации фазовых нестабильностей в каналах приема может быть использован сигнал РС [3, 5, 6]. Подход основан на том, что фазовая нестабильность на входе коррелятора является похожей для всех сигналов ЗС, которые ретранслируется данным набором КА, в то время как сдвиг по частоте
DOI: 10.24411/2410-9916-2020-10102
Системы управления,связи и безопасности №1. 2020
Systems of Control, Communication and Security ISSN 2410-9916
для каждой ЗС отличается из-за эффекта Доплера. В работах [5, 6] предполагается, что данная компенсация с использованием РС выполняется, но не предоставляется описание алгоритма компенсации или результатов экспериментов. В работе [3] предложен алгоритм компенсации фазовых нестабильностей и продемонстрирована его работа для повышения выходного ОСШ и улучшения точности оценки частотного сдвига. Данный алгоритм был также применен в [9]. Согласно [3] после получения пика огибающей ВКФ для сигналов, принятых от РС, вычисляются временной сдвиг и примерное значение частотного сдвига сигналов. Это позволяет рассчитать подынтегральную функцию в выражении для ВКФ для оцененных сдвигов и, таким образом, определить сигнал рассогласования, который содержит информацию о фазовых нестабильностях в каналах и остаточном частотном сдвиге. Корректирующий сигнал представляет собой сигнал рассогласования после применения к нему фильтрации. При вычислении ВКФ для неизвестных сигналов подынтегральная функция умножается на комплексно-сопряженный корректирующий сигнал. Недостатком данной схемы компенсации является условие получения пика огибающей ВКФ для сигнала РС, что не всегда выполнимо, если сигнал РС сам является слабым. В противном случае необходима априорная информация о сдвигах сигналов РС в каналах.
Компенсация фазовых нестабильностей возможна и без применения РС [10]. В этом случае интервал наблюдения разбивается на несколько коротких сегментов, внутри которых оценивается изменение частотного сдвига. Изменение частотного сдвига моделируется линейной функцией, параметры которой подбираются из условия максимизации пика ВКФ. После компенсации изменения сдвига частоты от сегмента к сегменту, результаты корреляций по каждому сегменту сращиваются, что позволяет увеличить выходное ОСШ. Данный подход предполагает, что на коротком сегменте, на котором справедлива квадратичная модель фазовой нестабильности, возможно получить пик огибающей ВКФ с достаточным выходным ОСШ, чтобы оценить динамику изменения частотного сдвига. Для слабых сигналов это условие зачастую не выполняется.
Алгоритм компенсации фазовых нестабильностей, предлагаемый в данной работе, во многом похож на [3], но использует корреляцию принятых сигналов РС с известным излученным сигналом, что позволяет работать даже со слабыми сигналами РС. Это связано с тем, что в данном случае один из коррелируемых сигналов является чистым (без шумов), что позволяет получить пик огибающей ВКФ на коротких интервалах интегрирования. Отличительной особенностью предлагаемой работы является также сопоставление практических результатов с теоретическими расчетами с подробным описанием всех проводимых вычислений.
Работа в дальнейшем представлена следующим образом. В первом разделе приводятся формулы для расчета выходного ОСШ. Во втором разделе рассматривается вопрос компенсации при расчете ВКФ фазовых нестабильностей в принятых сигналах. Затем представлены результаты эксперимента.
DOI: 10.24411/2410-9916-2020-10102
Systems of Control, Communication and Security
ISSN 2410-9916
1. Отношение сигнал шум на выходе коррелятора
Пусть спектральная плотность шумов в пределах полосы пропускания по частоте радиоприемного устройства сохраняет постоянное значение. Пусть есть две выборки сигнала ЗС х(и) и у(и), принятые в двух независимых каналах приема (полученные с двух КА), содержащие по N комплексных отсчетов на частоте дискретизации
х(и) = Sl(n) + д\(п), (1.1)
у(и) = ^(и) + ^(и) = а-^и+Ад) + ^(и), п = 0,1,...N-1, (1.2)
где s1(n) - комплексная огибающая сигнала ЗС, принятого в первом канале приема, ^2(и) - комплексная огибающая сигнала ЗС, принятого во втором канале приема с коэффициентом затухания а, ^(и), д2(и) - аддитивные гауссовские центрированные шумы в первом и втором каналах соответственно, Ад - задержка сигнала s1(n) относительно s2(n), выраженная в отсчетах сигнала. Отметим, что сигнал s1(n) считается некоррелированным с шумами д1(п) и д2(и).
Тогда взаимная корреляция сигналов х(и) и у(и), вычисленная для максимальной возможной задержки сигнала ЗС Отс1х [11]:
Kxy (m) =
Гу (n m) ht;, . Г"-т-ч , (3)
Z х(п + т) ■ y(«)*, т = 0,1,...,Dmax, 1 1 " n = 0 (2)
N-\m|-1 (2)
NM Z Х(П) ■ У(П - ^ m = Dmax ~Dmax + U."1,
1 1 n = 0
где ()* - обозначает комплексное сопряжение. Введем функцию, описывающую произведение сигналов х и y с учетом их сдвига по времени на m отсчетов:
x(n + т) ■ y(n)*, m = 0,1,...,Dn
[x(n) ■ y(n - m)\ m = -Dmax, -Dmax + -1.
Используя (3), выражение (2) можно записать в виде:
N- \ m \-1
Kxy(m) = Nm\ Z Гху (П m)- (4)
1 1 n = 0
Используя (1), раскроем выражение (4) для Kxy(m):
N-\m\-1 N-\m\-1
Kxy(m) = ЩМ\ Z rs1s2 (П m) + Nqmf Z rsiq2 (П m) + 1 1 n = 0 1 1 n = 0
N-\m\-1 N-\m\-1
+ ЩМ\ Z rqis2 (П m) + Nm Z rqiq2 m) = Kss(m) + vsiq2(m) +
1 1 n = 0 1 1 n = 0
+ Vqis2(m) + Vqiq2(m) = Kss(m) + vs(m), (5)
N-\m\-1
где Kss(m)= Z rs\si(n, m) - ВКФ для сигнала ЗС, а все остальные состав-
\ \ n = 0
ляющие vs1q2, vq1s2, vq\q2 являются выходным шумом vs, полученным за счет усреднения произведения шумовых отсчетов, а также произведения шума с сигналом. Выходной шум vs стремится к нулю при увеличении времени интегрирования T = NFs [7].
Примем Dmax << N. ОСШ на выходе коррелятора определятся выражением [12]:
Rout = |Kss(Di,2)| / STD{vs} = Sout / ov, (6)
N-l ml-1
DOI: 10.24411/2410-9916-2020-10102
Systems of Control, Communication and Security
ISSN 2410-9916
где STD{x} - операция вычисления среднеквадратического отклонения (СКО), Оу - СКО выходного шума = Vslq2 + Vqls2 + Vqlq2 , Sout = \К^1,2)\ - уровень выходного сигнала. При этом Sout = Ой = , где Osl, as2 - СКО сигналов ¿1 и s2 соответственно, Р^, Р^ - их средние мощности в полосе частот приема В. Здесь и далее частотные полосы определяются как двусторонние (в диапазоне от -да до +да), т.е. в данном случае В = Ец.
СКО выходного шума уе уменьшается при усреднении пропорционально 4ЁТ [7]:
= Л/а2я1аJ2 + а\а2Л + а>J2 l4ßT , (7)
где oq1, oq2 - СКО шумов q1 и q2 на входе коррелятора. Таким образом, с учетом (6) и (7), ОСШ на выходе коррелятора:
Rout = Osi ■ OsiJBT / д/a+ a2qla2s2 + a2qla2q2 . (8)
При условии, что полезный сигнал принимается под шумами, т.е.
Os << Oq:
Rout ~ Os1 ■ Os2 VBT / [Oq1 ■ Oq2]= Rini Rin2' 4BT, (9)
где Rin1 = os1/oq1, Rin2 = os2/oq2 - ОСШ для сигналов в первом и втором канале соответственно, посчитанное в полосе частот B. Выражение (9) удобно рассматривать в дБ:
R'out = R'in1 + R'in2 + 10-log(BT) = R'in1 + R'in2 + Radd (дБ), (10.1)
R'in1= 20^log(Rin1), R'in2 = 20^log(Rin2), Radd = 10-log(BT) (дБ), (10.2) где Radd - выигрыш в дБ за счет когерентного накопления сигнала ЗС. Например, Radd увеличивается на 3 дБ при удвоении времени интегрирования в корреляторе.
Необходимое ОСШ на выходе коррелятора выбирается исходя из заданных характеристик правильного обнаружения и ложной тревоги [12]. Как правило, алгоритм обнаружения сигнала сводится к поиску такого значения огибающей ВКФ Z = \Kxy\, которое больше некоторого заданного порога: Z > Zo. Согласно работам [7, 12] для обнаружения пика огибающей ВКВ на фоне выходных шумов необходимо иметь R'out не менее 10...12 дБ, а для точной оценки задержки и частоты около 20 дБ [2]. Например, если Rin1 = Rin2 = -30 дБ, то согласно (10) потребуется 16 с обработки в полосе частот 1 МГц, чтобы получить R'out = 12 дБ (4 ед.).
Из выражения (8) следует, что увеличение Rout можно достигнуть также за счет увеличения полосы частот B, в которой ведется обработка сигналов. Однако при увеличении B растет СКО шумов на входе коррелятора пропорционально 4B. Если при этом не увеличится СКО сигналов, также пропорционально 4B (например, если полоса частот сигнала Bs меньше нового значения B), то Rout будет ниже, чем до увеличения B. Это означает, что при условии Fs > Bs, необходимо использовать фильтр до коррелятора, полоса пропускания которого Bf согласована с полосой сигнала, т.е. Bf = Bs. В противном случае проигрыш по выходному ОСШ будет равен ■sjFJB.
DOI: 10.24411/2410-9916-2020-10102
Systems of Control, Communication and Security
ISSN 2410-9916
За оценку выходного ОСШ Rout по измеренной ВКФ, что потребуется для сопоставления теоретических и экспериментальных расчетов, примем:
Rout = Sout / av ~ Zmax / (7v , (11)
где Zmax - максимум огибающей измеренной ВКФ, который принимается за оценку уровня выходного сигнала Soui, а av - оценка СКО шумовых отсчетов ВКФ. Важно отметить, что оценка Sout = Zmax будет смещенной, поскольку плотность вероятности случайной величины Z(D^2) = \Kxy(D12)\ = = \KssD12) + vs(A,2)\ будет определяться законом Райса [12], а ее среднее значение будет больше значения Sout. При увеличении значения Rout распределение Z(DU) будет стремиться к нормальному закону распределения со средним значением равным Sout. Для интересующего диапазона Rout > 3 ед. и целей работы смещением среднего значения Z(D1,2) относительно значения Sout можно пренебречь.
На точность оценки Rout будет также влиять выбор Fs, поскольку частота дискретизации определяет точность нахождения максимума огибающей ВКФ. Так при условии, что Fs > Bf > Bs, интервал корреляции сигнала на входе коррелятора тс = 1/Bs. Выбор Fs = Bs приводит к ухудшению характеристик обнаружения, которое в худшем случае составит двукратное уменьшение мощности выходного сигнала. Это означает, что Rout упадет в 0,7 раза. Поэтому на практике для точного нахождения максимума огибающей ВКФ часто полагают Fs = (2...4)Bs.
2. Коррекция частотной отстройки и фазовых искажений при расчете ВКФ
Выражение для достижимого выходного ОСШ (8) справедливо для модели (1 ), когда мгновенные фазы сигналов в каналах совпадают (без учета начальных фаз), т.е. в каналах отсутствует частотная отстройка и фазовые нестабильности. В реальных условиях данное допущение не выполняется, и модель принятого в k-м канале сигнала необходимо дополнить [3, 10, 13]:
^ (n) =ak • s0 (n - D) • e [2%'¥k 'n'Ts +SPk (nWk], (12.1)
Afk = A/пер k + /доп k , (12.2)
5фк(n )= 5p ПрД (n)+ 5p Прм k(n)+ 5p KA k(n)+ 5p атм k(n)+ SP доп k(n X (12.3) где k - номер канала, S0(n) - комплексная огибающая сигнала, переданного ЗС, Dk - задержка сигнала ЗС в канале, A/k - сдвиг частоты в канале, A/пер k - суммарная отстройка по частоте, вносимая в сигнал посредством передатчика, приемника и задействованного КА, А/доп k - сдвиг частоты за счет допплеровского эффекта, вызванного движением КА относительно ЗС и приемника, 5фk - фазовые искажения в канале, 5фпрд, 5фпрм k, 5фКА k - фазовые искажения, вносимые в сигнал посредством передатчика, приемника и задействованного КА соответственно, 5фатм k - фазовые искажения, вносимые в сигнал из-за атмосферных явлений в канале, 5фдоп k - фазовые искажения, вызванные изменением А/доп k за время наблюдения, фk - неизвестная и неинформативная начальная фаза.
DOI: 10.24411/2410-9916-2020-10102
Системы управления,связи и безопасности №1. 2020
Systems of Control, Communication and Security ISSN 2410-9916
Наличие частотных отстроек А/1, А/2 и фазовых искажений (нестабильно-стей) 5ф1, 5ф2 ограничивает время корреляции сигналов s1(n) и s2(n). Частотная отстройка учитывается путем одновременного поиска сдвига сигналов как по задержке, так и по частоте на интервале времени T [8]. Поиск по частоте выполняется путем перебора по конечному числу значений /Ае [-/а max, /а max] с заданным шагом 5/
N- | m | -1
кхуЫа) = Nlm Z 'ху (птУ e""2,
n = 0
/А = - /А max5 7 A max + 5/,..., /а max-
Максимум огибающей ВКФ достигается при частоте /А = А/1 - А/2. Максимальное допустимое значение шага поиска по частоте составляет 5/ = 2/(3 T). В худшем случае оно приведет к двукратному уменьшению мощности выходного сигнала. На практике для точного нахождения максимума огибающей ВКФ |Kxy(m/A)| часто полагают 5/ = 1/(3T)...1/(6T).
Исходя из (12), искажения фазы, которые будут присутствовать в сигнале rs1s2(n, Ад) при корреляции сигналов x(n) и у(п) по формуле (4):
5фА(п) = 5ф1(п+Ад) - 5ф2(п) ~ 5ф1(п) - 5ф2(п) = 5фпрм i (п) + + 5фКА 1 (п) + 5фатм 1 (п) + 5фдоп 1 (п) - 5фпрм 2 (п) - 5фКА 2 (п) -- 5фатм 2 (п) - 5фдоп 2 (п), (13)
где упрощение 5ф1(п ± Ад) ~ 5ф1(п) принято в виду того, что нестабильности фазы являются медленными процессами и при их анализе можно пренебречь влиянием задержки Ад, которая обычно составляет миллисекунды. Из (13) следует, что составляющие фазы 5фпрм 1, 5фпрм 2, 5фКА 1, 5фКА 1 будут одинаковыми для всех сигналов ЗС, которые ретранслируется данным набором КА. Это заключение также в целом справедливо по отношению к искажениям фазы 5фатм 1, 5фатм 2, вызванным атмосферными явлениями [3]. Таким образом, зная искажения фазы для одного сигнала ЗС можно построить эквалайзер для компенсации этих составляющих 5фА в сигналах других ЗС. Отметим, что при этом сдвиг по частоте /А для каждой ЗС будет отличным. Некомпенсированными останутся искажения 5фдоп 1, 5фдоп 2, которые начинают играть роль на интервалах наблюдения длительностью в несколько сотен секунд [10]. В данном случае, после применения эквалайзера к сигналу ЗС можно использовать подход, предложенный в [10], когда при вычислении ВКФ оценивается изменение /А на интервале наблюдения в виде линейной функции от времени.
Искажения в каналах оцениваются с помощью известного переданного РС сигнала sR, который принимается одновременно с сигналом интересующей ЗС:
х(п) = ^(п) + sr (n) + 41(п), (14.1)
у(п) = ^(п) + s2r (n) + 42(п), (14.2)
где sR, sR связаны с sR через (12).
Получить сигнал s0 на станции мониторинга можно несколькими способами. Например, в определенные моменты времени РС отправляет заданную битовую последовательность, преобразованную в модулированный сигнал из-
DOI: 10.24411/2410-9916-2020-10102
Системы управления,связи и безопасности №1. 2020
Systems of Control, Communication and Security ISSN 2410-9916
вестным образом. Тогда на станции мониторинга возможно смоделировать сигнал sR с точностью до частотных нестабильностей передатчика, которые роли в дальнейшей обработке не играют, так как будут скомпенсированы при вычислении разности фаз принятых с двух КА сигналов. Возможен также вариант, когда сигнал с выхода модулятора РС до его переноса на несущую частоту вместе с метками единого координированного времени пересылается на станцию мониторинга.
Поскольку сигнал sR является чистым (без шумов), то пик огибающей ВКФ принятых сигналов x, y и сигнала s0 можно получить с высоким выходным ОСШ (выше 20 дБ) на коротком интервале интегрирования (несколько секунд) и, таким образом, оценить сдвиги сигнала sR по времени DDR и частоте А/R для каждого канала. С учетом (12) и (14) сигнал рассогласования для РС в к-м канале после компенсации сдвигов DR и А/0:
eR(n)= [ sR(n)+sk (n)+q (n )]•sR (n - DR) • e~"2К f n'Ts =
= aR •e'•[2К */0nTs+8pR(n)+фк 1 + ^ (n), где 8/R = А/0 - А/0 - остаточный сдвиг частоты, 8фR - фазовые искажения в сигнале РС в к-м канале, vk - выходной шум. Для получения оценки 8pR сигнал eR подвергается низкочастотной фильтрации, после чего извлекается его фаза, в фазе компенсируется линейный набег, вызванный наличием 8/0.
При расчете ВКФ для обнаружения интересующей ЗС выполняется компенсация фазовых искажений, т.е. применяется эквалайзер:
N-\m\-1 Г „„ „ 1
F / гч 1 У т. / ч 0-1 [2 К •fа • n • Ts +2 К • fA •n • Ts + 8ф А (n )J
Kxy(m,f) = N-rmr ^ rxy (n тУe , (15)
1 1 n = 0
где 8фА (n) = 8pR (n )-8pR (n) и fR = А/°.
3. Результаты эксперимента
В ходе эксперимента были приняты и зарегистрированы сигналы x(n) и
y(n) с двух КА. В данных присутствовал сигнал РС sR (n) с шириной полосы частот 1,2 МГц, у которого расширение спектра происходило за счет использования псевдослучайной последовательности. Обработка велась на частоте дискретизации Fs = 3 МГц при использовании фильтра низких частот с полосой Bf = 1,2 МГц. Длительность записи сигнала составляла 19 с, что позволяет согласно (10) достигнуть в данной полосе выигрыша по ОСШ Radd = 73,6 дБ.
Прежде всего было измерено входное ОСШ в каждом канале с помощью корреляции принятого сигнала с сигналом sR с единичным СКО (gs0 = 1) на интервале длительностью T = 3 с. При расчете корреляции выполнялся поиск отстройки принятого сигнала по частоте с шагом f = 1/(4T). На рис. 2 и 3 представлены полученные огибающие ВКФ в двух каналах приема. Поскольку шум на входе по уровню значительно превосходит полезный сигнал, то СКО вход-
DOI: 10.24411/2410-9916-2020-10102
Systems of Control, Communication and Security
ISSN 2410-9916
ного шума можно принять равным СКО принятого сигнала, т.е. оч1 ~ ах, ~ ау. В то время как за СКО полезного сигнала при условии большого выходного ОСШ (больше 20 дБ) можно принять значение максимума огибающей ВКФ. В результате было полечено Я'м = -33,9 дБ, Я'п2 = -23,3 дБ, при этом согласно (10) ожидаемое выходное ОСШ при накоплении сигнала в течение 19 с составит Я'апг = 16,4 дБ (6,6 ед).
Рис. 2. Огибающая ВКФ принятого в Рис. 3. Огибающая ВКФ принятого
r r
первом канале сигнала с s0 во втором канале сигнала с s0
На рис. 4 представлены оцененные по РС реализации фазовых искажений в каждом канале 8ф R (t), 8ф R (t), а также разностные фазовые искажения 8ф R (t).
400
300
S
ö 200
0
и
Л
В 100
«о
С5 H
8 0 к
а
1 100
о
а
е
200
-300
Ô9f(0
10 15 20
Время I. с
Рис. 4. Реализации фазовых искажений 8ф R (t), 8ф R (t), 8ф R (t)
Видно, что фазовые искажения 8ф Д ) достигает по амплитуде 300°, что нарушает условия когерентности сигналов в соседних каналах при вычислении ВКФ. На рис. 5 и 6 представлены полученные огибающие ВКФ для Т = 19 с и
DOI: 10.24411/2410-9916-2020-10102
Systems of Control, Communication and Security
ISSN 2410-9916
5/ = 1/(47) для случаев без использования эквалайзера и с применением эквалайзера (15). В первом случае сигнал РС не наблюдается, так как уровень выходного сигнала находится на уровне выходных шумов. После применения эквалайзера наблюдается сигнал РС с ОСШ Rout = 6,5 ед. (16,2 дБ), что соответствует теоретически достижимому значению Rout = 6,6 ед.
Рис. 5. Огибающая ВКФ без эквалайзера (T = 19 с)
Рис. 6. Огибающая ВКФ с эквалайзером (T = 19 с)
Оценки изменения уровня выходного сигнала §ош (Т) и СКО выходного шума (гу (Т) по мере увеличения интервала интегрирования Т с 4 до 19 с при условии применения эквалайзера представлены на рис. 7 и 8. Поскольку при расчете ВКФ выполняется нормирование выходного значения к длительности
интервала интегрирования, значение 8ои{ должно оставаться неизменным, в то время как уровень выходных шумов должен уменьшатся со скоростью пропорциональной
4Т , что и видно на представленных графиках. График (Т) показан на рис. 9. Здесь же пунктирной линией обозначено ожидаемое выходное ОСШ Яоиг (Т), рассчитанное по формуле (9). Видно, что результаты эксперимента с высокой точностью согласуются с теоретическими расчетами.
Помимо сигнала РС в зарегистрированных сигналах х(п) и у(п) также присутствуют сигналы других ЗС, которые отличаются от сигнала РС отстройкой по частоте /д. Для обнаружения этих передатчиков выполнялся расчет ВКФ после применения эквалайзера (15) с учетом сдвига по частоте =А/* , соответствующего РС. Значения сдвига по частоте / перебирались в диапазоне -1 до 1 Гц с шагом 5/ = 1/(3Т). Для каждой ВКФ искался максимум ее огибающей и сравнивался с выбранным порогом 20. Приведенный к оценке СКО шумов максимум огибающей ВКФ 2 = 2Пах / с для набора частот/ представлен на рис. 10. Здесь же обозначен нормированный порог обнаружения И = 20/ ссу.
DOI: 10.24411/2410-9916-2020-10102
Systems of Control, Communication and Security
ISSN 2410-9916
Время интегрирования Г, с
Рис. 7. Оценка уровня выходного сигнала Sout (T)
4M)
350 300 250 200 15(1
10 12 14 16 18 20
Время интегрирования Г, с
Рис. 8. Оценка уровня выходного шума ст 0 (Т)
[ интегрирования
Рис. 9. Оценка выходного ОСШ Rout (T) и его ожидаемое значение Rout (T)
У
1 \
У
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Отстройка по частоте Гц
Рис. 10. Приведенный к оценке СКО шумов макодмум огибающей ВКФ
S 461S
46 s/ 46 IS
5 А 50 ЗС №2 1
5150/ \5150 \ Л\ S7QS
ЗС №3
/ PC w8795 г - г; х чД4618
/ //
ЗС №1
Отстройка по частоте /д, Гц
Рис. 11. Обнаруженные ЗС
DOI: 10.24411/2410-9916-2020-10102
Systems of Control, Communication and Security
ISSN 2410-9916
Группа обнаруженных сигналов в укрупненном масштабе показана на рис. 11. В дополнение обозначены измеренные в отсчетах задержки для каждого значения 2, превысившего заданный порог И = 4. Передатчик, обнаруженный при / = 0, соответствует РС. Также присутствуют еще три ЗС, обозначенные на рис. 11. Огибающие ВКФ для каждой обнаруженной ЗС для максимального уровня выходного сигнала представлены на рис. 12-14.
1000
N ©
g еоо
pq «
S 600
УО 400
ЗС №2
Zo
200
Время г. с Время г, с
Рис. 12. Огибающая ВКФ для первой Рис. 13. Огибающая ВКФ для второй обнаруженной ЗС обнаруженной ЗС
Рис. 14. Огибающая ВКФ для третей обнаруженной ЗС
На рис. 15 показаны графики изменения выходного ОСШ (Т) для трех обнаруженных ЗС. Видно, что выходное ОСШ увеличивается по мере накопления сигнала для всех трех ЗС, что говорит о выполнении условий когерентного их накопления после применения эквалайзера.
DOI: 10.24411/2410-9916-2020-10102
Systems of Control, Communication and Security
ISSN 2410-9916
Выводы
В работе рассмотрена задача обнаружения слабых сигналов ЗС в спутниковых системах связи в целях контроля частотного ресурса КА. Для обеспечения необходимых для нахождения пика огибающей ВКФ длительностей интегрирования (десятков секунд) выполняется компенсация фазовых нестабильно-стей в принятых сигналах с помощью эквалайзера, настроенного по известному сигналу РС. После применения эквалайзера, алгоритм обнаружения всех передатчиков сводится к расчету ВКФ сигналов, принятых в соседних каналах, для разных значений отстроек по частоте и сравнению ее огибающей с заданным порогом. Приведены результаты эксперимента, которые подтвердили возможность длительного когерентного накопления сигнала ЗС при использовании предложенного эквалайзера. Отличительной особенностью предлагаемого эквалайзера является коррелирование принятых сигналов РС с известным излученным сигналом, что позволяет работать со слабыми сигналами РС.
Литература
1. Roddy D. Satellite Communications. - McGraw-Hill Telecommunications, 2001. - 569 p.
2. Haworth D. P., Smith N. G., Bardelli R., Clement T. Interference localization for Eutelsat satellites - the first European transmitter location system // International Journal of satellite communications. 1997. Т. 15. С. 155-183.
3. Haworth D. P. Locating the source of an unknown signal. U.S. Patent 6018312, Jan. 25. 2000.
4. Xue Y., Li X., Xu L., Ren Y. Research on position differential method of dual-satellites TDOA and FDOA in passive location system // IEEE International Frequency Control Symposium. - Baltimore, USA, 2012. - С. 1-5.
5. Агиевич С. Н., Бердников А. В., Буев С. Г., Ватутин В. М., Глуздов А. Н., Климов С. А., Полтавец Ю. И., Пономарев А. А., Севидов В. В., Смирнов М. А., Топорков И. С., Ширшов М. В. Способ определения
DOI: 10.24411/2410-9916-2020-10102
Системы управления,связи и безопасности №1. 2020
Systems of Control, Communication and Security ISSN 2410-9916
местоположения абонентского терминала с помощью не менее двух спутников -ретрансляторов на низкой околоземной орбите // Патент на изобретение RU 2684740 C1, опубл. 15.04.2019, бюл. № 11.
6. Богомолов А. В., Кельян А. Х., Севидов В. В., Чемаров А. О., Эконом В. П. Способ определения местоположения земной станции спутниковой связи по ретранслированному сигналу // Патент на изобретение RU 2663193 C1, опубл. 02.08.2018, бюл. № 22.
7. Stein S. Algorithms for ambiguity function processing // IEEE Trans. on acoustics, speech and signal process. 1981. T. ASSP-29. № 3. С. 588-599.
8. ГЛОНАСС. Принципы построения и функционирования / под. ред. А. И. Перова. - М.: Радиотехника, 2010. - 800 с.
9. Chan M. H. Application of a dual satellite geolocation system on locating sweeping interference // World Academy of Science, Engineering and Technology. 2012. Т. 6. С. 939-944.
10. Wu R., Zhang Y., Huang Y., Xiong J., Deng Z. A Novel Long-Time Accumulation Method for Double-Satellite TDOA/FDOA Interference Localization. // Radio Science. 2018. № 53. С. 129-142.
11. Proakis J. G., Manolakis D. G. Digital Signal Processing: Principles, Algorithms and Applications. - Upper Saddle River, New Jersey: Pearson Prentice Hall, Inc., 2007. - 1084 c.
12. Перов А. И. Статистическая теория радиотехнических систем. - М.: Радиотехника, 2003. - 400 с.
13. Волков Р. В., Саяпин В. Н., Севидов В. В. Модель измерения временной задержки и частотного сдвига радиосигнала, принятого от спутника-ретранслятора при определении местоположения земной станции // T-Comm: Телекоммуникации и транспорт. 2016. Т. 10. № 9. С. 14-18.
References
1. Roddy D. Satellite Communications. McGraw-Hill Telecommunications, 2001. 569 p.
2. Haworth D. P., Smith N. G., Bardelli R., Clement T. Interference localization for Eutelsat satellites - the first European transmitter location system. International Journal of satellite communications, 1997, vol. 15, pp. 155-183.
3. Haworth D. P. Locating the source of an unknown signal. U.S. Patent 6018312, Jan. 25, 2000.
4. Xue Y., Li X., Xu L., Ren Y. Research on position differential method of dual-satellites TDOA and FDOA in passive location system. IEEE International Frequency Control Symposium. Baltimore, USA, 2012, pp. 1-5.
5. Agievich S. N., Berdnikov A. V., Buev S. G., Vatutin V. M., Gluzdov A. N., Klimov S. A., Poltavets Yu. I., Ponomarev А. А., Sevidov V. V., Smirnov М. А., Toporkov I. S., Shirshov М. V. Method of determining location of subscriber terminal by means of at least two satellite converters on low oriental orbit. Patent Russia, no. RU 2684740 C1. Publish. 15.04.2019, bul. no.11 (in Russian).
6. Bogomolov А. V., Kelyan А. K., Sevidov V. V., Chemarov А. О., Ekonom V. P. Method of determining location of satellite earth station according to a
DOI: 10.24411/2410-9916-2020-10102
Системы управления,связи и безопасности №1. 2020
Systems of Control, Communication and Security ISSN 2410-9916
repeated signal. Patent Russia, no. RU 2663193 C1. Publish. 02.08.2018, bul. no. 22 (in Russian).
7. Stein S. Algorithms for ambiguity function processing. IEEE Trans. on acoustics, speech and signal process, 1981, vol. ASSP-29, no. 3, pp. 588-599.
8. GLONASS. Principy postroeniya i funkcionirovaniya [GLONASS. Principles of construction and operation]. By edit. A. I. Perov. Moscow, Radiotehnika Publ., 2010. 800 p. (in Russian).
9. Chan M. H. Application of a dual satellite geolocation system on locating sweeping interference. World Academy of Science, Engineering and Technology, 2012, vol. 6, pp. 939-944.
10. Wu R., Zhang Y., Huang Y., Xiong J., Deng Z. A Novel Long-Time Accumulation Method for Double-Satellite TDOA/FDOA Interference Localization. Radio Science, 2018, no. 53, pp. 129-142.
11. Proakis J. G., Manolakis D. G. Digital Signal Processing: Principles, Algorithms and Applications. Upper Saddle River, New Jersey, Pearson Prentice Hall, Inc., 2007. 1084 p.
12. Perov A. I. Statisticheskaya teoriya radiotekhnicheskih sistem [Statistical theory of radio engineering systems]. Moscow, Radiotehnika Publ, 2003. 400 p. (in Russian).
13. Volkov R. V., Saypin V. N., Sevidov V. V. Model measuring the time delay and frequency shift of the radio signal received from the satellite repeater in locating ground terminal. T-Comm: telecommunications and transport, 2016, vol. 10, no. 9, pp. 14-18 (in Russian).
Статья поступила 12 ноября 2019 г.
Информация об авторе
Кулакова Вероника Игоревна - кандидат технических наук. Начальник отдела - главный конструктор ООО "Специальный технологический центр". Область научных интересов: радиопеленгация, определение местоположения источников радиоизлучений, геолокация, навигация. E-mail: nika_kulakova@mail.ru
Адрес: 195220, Россия, г. Санкт-Петербург, Гжатская ул., дом 21, лит. Б, офис 53.
Detection of Weak Signals by Cross-Correlation with Phase Instabilities Compensation for Spectrum Monitoring in Satellite Communication Systems
V. I. Kulakova
Purpose. Satellites are widely used to retranslate signals from Earth stations over long distances. Detection of all Earth stations that use satellite communication systems is necessary to control the frequency resource of the satellite link and to protect it from unauthorized access. In this case a cross-correlation method is used for signal detection. The difficulty lies in the fact that the signals of unauthorized transmitters are often observed by a monitoring station at the significant negative signal-to-noise ratio (SNR) (less than
DOI: 10.24411/2410-9916-2020-10102
Системы управления,связи и безопасности №1. 2020
Systems of Control, Communication and Security ISSN 2410-9916
minus 20 dB). In this regard, a long time of coherent accumulation of the detected signal (tens of seconds) is required to obtain the peak in the cross-correlation function (CCF) envelope. The purpose of the present paper is to investigate the possibility of detecting of weak signals of Earth stations in real conditions of application. Paper uses the reference Earth station for correcting the phase distortions introduced into the signals in each channel. Results. Paper presents the equations for calculating the SNR at the output of the cross-correlator, depending on the SNR in each channel, the receiver frequency band and the duration of the integration interval. In order to detect the peak of the envelope of the CCF in the presence of output noise, the output SNR should be greater than 10...12 dB. This means that if an Earth station signal is received with an SNR about -30 dB in each channel, it should be integrated during 16 s in a 1 MHz frequency band. When calculating the cross-correlation for such long intervals, it is necessary to take into account the presence of frequency shifts and phase distortions in different channels. It is assumed that the phase instability is the same for all signals that are retransmitted by this set of satellites, while the frequency shift for each signal is different due to the Doppler effect. Based on this, the equalizer is used to compensate for the phase distortions, which is adjusted by the known signal from the reference Earth station. Paper presents the experimental results when signals from two satellites were received and registered. In the data there was a signal of the reference Earth station, whose spectrum was expanded with the use of a pseudo-random sequence. The duration of the recorded signal was 19 s. The phase distortions reached 300° in amplitude, which destroyed the correlation of signals in different channels. It is shown that with the use of equalizer, adequate steps for finding the maximum of the envelope of the cross-correlation function in time and frequency, the experimental output SNR with high accuracy corresponds to theoretical calculations. In addition to the signal of the reference station, three more transmitters were found which differ from the signal of the reference station by frequency shifts.
Key words: satellite communication system, signal detection, cross-correlation, signal to noise ratio.
Information about Author
Kulakova Veronika Igorevna - Ph.D. of Engineering Sciences. Head of Department - Chief Designer. «Special Technological Center» ltd. Field of research: the radio direction finding, determining the location of radio emission sources, geoloca-tion, navigation. E-mail: nika_kulakova@mail.ru
Address: Russia, 195220, Saint-Petersburg, Gzhatskaya, 21/B.
DOI: 10.24411/2410-9916-2020-10102