ПОМЕХОУСТОЙЧИВОСТЬ АЛГОРИТМОВ ПРИЁМА СИГНАЛОВ С МНОГОИМПУЛЬСНОЙ ПОЗИЦИОННО-ИМПУЛЬСНОЙ МОДУЛЯЦИЕЙ
В.И. Парфенов , ДЮ. Голованов 'Воронежский государственный университет, Воронеж, Россия
Аннотация
В работе предложен алгоритм оценки временных положений и амплитуд импульсов периодической последовательности по незначительному количеству отсчётов, число которых определяется только числом этих импульсов. Рассмотрена работа этого алгоритма в случае, когда фильтром, ограничивающим спектр исходного сигнала, является идеальный фильтр низких частот, а также фильтр Найквиста и определены условия перехода от одного фильтра к другому. Исследована эффективность предложенного алгоритма на основе анализа зависимости рассеяния оценки временных положений и амплитуд от отношения сигнал/шум и числа импульсов последовательности. Показано, что с этой точки зрения эффективность работы алгоритма снижается с ростом числа импульсов последовательности и уменьшением отношения сигнал/шум.
Отмечено, что предложенный алгоритм, в отличие от классического алгоритма максимального правдоподобия, не требует поиска максимума функции большого числа переменных, и в то же время характеристики оценок оказываются практически одинаковыми для обоих этих методов. Показано также, что точность оценки предложенного алгоритма может быть повышена за счёт незначительного увеличения количества отсчётов сигнала.
Полученные результаты могут быть использованы при практическом проектировании лазерных систем связи, в которых в качестве метода передачи сообщений используется многоимпульсная позиционно-импульсная модуляция.
Ключевые слова: многоимпульсная позиционно-импульсная модуляция, идеальный фильтр низких частот, фильтр Найквиста, алгоритм оценки, метод MUSIC, отношение сигнал/шум, рассеяние оценки, метод максимального правдоподобия, функция правдоподобия.
Цитирование: Парфенов, В.И. Помехоустойчивость алгоритмов приёма сигналов с многоимпульсной позиционно-импульсной модуляцией / В.И. Парфенов, Д.Ю. Голованов // Компьютерная оптика. - 2018. - Т. 42, № 1. - С. 167-174. - DOI: 10.18287/2412-6179-2018-421-167-174.
Введение
Лазерная техника открывает большие возможности для совершенствования систем связи, локации и радиоуправления. Эти возможности связаны, в первую очередь, с громадным коэффициентом усиления передающих оптических антенн, что позволяет получить большое отношение сигнал/шум в приемнике в широкой полосе частот при маломощном передатчике, а также с возможностью использования очень широкой полосы частот при передаче и приёме оптических сигналов. Мощность передатчика во многом определяется используемым методом модуляции. Одним из наиболее выгодных методов передачи сообщений в лазерных системах связи с прямым детектированием, позволяющих использовать передачу с минимальной мощностью, является цифровая позиционно-импульсная модуляция (ПИМ или PPM в английской аббревиатуре (PulsePosition Modulation)) [1]. Как показано в [2], для классических оптических каналов при ограничениях на пиковую и среднюю мощности, PPM может обеспечить близкую к оптимальной ёмкость канала. Подобный метод модуляции широко используется в оптоволоконных и спутниковых системах [3, 4], кроме того, он применяется и в сверхширокополосных радиокоммуникационных системах [5].
При классическом способе формирования сигналов с PPM символьный интервал [0; T] разбивается на
L подынтервалов (слотов) длительностью т0 каждый. Передаче l-го символа соответствует передача импульса с временным положением (l - 1) т0. Очевидно, максимальное количество возможных передаваемых сообщений в этом случае не превосходит L. Интересной альтернативой PPM является многоимпульсная позиционно-импульсная модуляция (Multipulse PulsePosition Modulation (MPPM)) [6], которая имеет преимущества по сравнению с классической PPM в уменьшенном отношении пиковой мощности к средней. При таком подходе передаётся уже не один, а K импульсов в течение тех же L слотов. В этом случае количество передаваемых сообщений существенно увеличивается по сравнению с классической PPM. Например, если количество слотов L = 10, то при классической PPM можно осуществить передачу 10 различных сообщений (отличающихся временными положениями передаваемых импульсов), в то время как при использовании MPPM количество разных сообщений будет определяться числом сочетаний Ск (например, при L = 10 и K = 2 это количество становится равным 45). Следует отметить, что MPPM превосходит PPM также по частотной эффективности [7].
В дальнейшем будут рассматриваться многоимпульсные PPM (MPPM) с точки зрения разработки новых методов их демодуляции. Существующие методы приёма (демодуляции) таких сигналов основы-
ваются на методах турбодекодирования с использованием алгоритма Витерби с мягким принятием решения [8], что бывает не просто реализовать на практике. Учтём, что информация о передаваемых сообщениях закодирована во временных положениях (позициях) передаваемых К импульсов. Следовательно, демодуляция (приём) может заключаться в оценке этих временных положений. Однако при цифровой обработке таких сигналов могут возникнуть сложности, связанные с необходимостью обработки огромного количества отсчётов и их перебора. Отметим, что аналогичные проблемы возникают везде, где используется поток коротких импульсов, например, при исследовании нейронной активности, при формировании и обработке биоизображений [9], при сверхширокополосной связи [10] и пр.
1. Алгоритм оценки временных положений и амплитуд периодической последовательности сверхкоротких импульсов с использованием идеального фильтра нижних частот
Постановка задачи
Предположим, что излучаемый сигнал представляет собой сумму К достаточно коротких импульсов, в идеале - дельта-функций. Представим такой сигнал в виде
s(t) = IX £ SCí - tk - nT),
(1)
где Ак и % - амплитудные множители (в дальнейшем просто амплитуды) и временные положения отдельных импульсов на периоде Т, предполагаемые априори неизвестными, а число таких импульсов К на периоде предполагаем известным. Теоретически спектр такого сигнала бесконечен, что при дискретизации потребует бесконечно малого интервала дискретизации (бесконечно большого количества отсчётов). Пройдя среду распространения с известной импульсной характеристикой, дельта-импульсы расширяются во времени так, что вместо (1) полезный сигнал будет иметь вид
s(t) = £Л £ F(t- tk - nT),
(2)
где - известная функция времени, описывающая форму импульсов, прошедших канал распространения и получающихся в результате вычисления свёртки входного сигнала с импульсной характеристикой канала. Подобный сигнал будет иметь вид также на выходе преобразователя лазерного (инфракрасного) излучения в электрические импульсы. Для такого преобразования обычно используется детектор сигнала (рш-фотодиод) [11].
Учитывая периодический характер, представим сумму по п в (2) в виде ряда Фурье
£ Г(г - пТ) = £ ехр(2п]т1 / Т), 1 Т
где /т = -1Г (г)ехр(-2п]т1/Т)& . Т 0
Тогда (2) перепишется как
где
s(t) = £ Ak £ fm exp(2njm(t - tk )/T) =
k=0 ш=-ж « K -1
= £ fm exp(2njmt/T)£ Ak exp(-2njmtk /T) =
ш=-ж k=0
= £ fmXm exp(2njmt/T),
K-1
Xn = £ At exp(-2njntk /T).
(3)
Учтём, что в оптической связи длительности импульсов, описываемых функцией F(t), чрезвычайно малы (пикосекунды и менее). Следовательно, при их цифровой обработке потребуется огромное количество отсчётов, что в настоящее время может привести к существенным аппаратным и финансовым затруднениям. Поэтому для уменьшения ширины спектра подадим такой сигнал на идеальный фильтр нижних частот (ФНЧ) с импульсной характеристикой вида h(t) = B sine (nBt), где sine (x) = sin (x)/x. Частота среза такого фильтра равна nB рад/с, в дальнейшем, не умаляя общности, будем считать величину nB чётной. В результате на выходе такого фильтра сигнал примет вид
s(t) = J s(v)h(t -v)dv =
от »
= X fmXm J exp(2njmv/T)h(t-v)dv =
= £ fmXmH
2nm
exp(2njmt/T).
Здесь Н(ш) - частотный коэффициент передачи рассматриваемого идеального ФНЧ с импульсной характеристикой Н(г).
С учётом того, что Н (2пт / Т) = 0 при |т| > ВТ/ 2 = М, имеем
s(t) = £ fmXm exp (2njmt/T) =
m=-M
M
= £ X 1m exp (2njmt/T),
(4)
Где X1 m fmXm.
Несложно заметить, что максимальная частота в спектре сигнала (4) равна пБ, что и следовало ожидать, т.к. широкополосный сигнал (2) был пропущен через идеальный ФНЧ с частотой среза пБ. Следовательно, его можно продискретизировать в соответствии с теоремой отсчётов с шагом At < п / пБ = 1/B. Причём количество необходимых отсчётов на периоде N становится весьма небольшим: N> T/At + 1 = 2M + 1 . Учтем также, что Xn - отсчёты суммы K комплексных экспонент (3). Как известно [12], для нахождения их амплитуд и фаз достаточно иметь N> 2K отсчётов. Однако M=BT/ 2. Следовательно,
Б>Л. (5)
k=0
k=0
k=0
m=-M
Таким образом, оказывается, что параметр В идеального ФНЧ определяется количеством импульсов на периоде К, а число необходимых отсчётов выходного сигнала фильтра
s n = £ X1m exp(2njmnAt / T), n = 0, N -1,
(6)
должно удовлетворять условию
N > 2K +1, (7)
и, следовательно, зависит лишь от количества импульсов на периоде, а не от их длительности.
Итак, если сигнал (2) предварительно подать на ФНЧ с указанными характеристиками, то дискретная обработка выходного сигнала фильтра потребует незначительного количества равноотстоящих отсчётов, зависящих не от спектральных характеристик входного сигнала, а от количества этих сигналов. Таким образом, с учётом вышеизложенного, задача определения временных положений tk и амплитуд Ak импульсов в (2) сводится к следующему. Выбираем идеальный ФНЧ с параметром (5) и подаём на него сигнал (2). Выходной сигнал этого фильтра равномерно дискретизируем на периоде T, получая N отсчётов (6), количество которых удовлетворяет условию (7). Формула (6) представляет собой, по сути, систему из N уравнений с 2М + 1 неизвестными Xn. Если отношение T/At выбрать целым, то эта система уравнений легко решается. Например, если T/At = N, то (6) - это обычное обратное дискретное преобразование Фурье. Следовательно, Xn может быть найдено как прямое преобразование Фурье от sn. Если T/At < N, то её решение может быть найдено через псевдообращение матрицы. Таким образом, находим N отсчётов Xn, которые в соответствии с (3) представляют собой сумму K комплексных экспонент с амплитудами Ak и временными положениями tk. Как известно [12], при выполнении условия (7) эти параметры (Ak и tk, n = 0, K -1 ) могут быть найдены любыми известными способами. В частности, для решения подобной задачи могут быть использованы методы Прони [13], Annihilating Filter [14], MUSIC [13] и пр.
Отметим, что рассмотренная выше задача является идеализированной, прежде всего, из-за отсутствия в наблюдаемых данных всегда присутствующих в реальности случайных искажений - шумов. Представим наблюдаемые данные в виде аддитивной смеси сигнала (2) и гауссовского белого шума n(t) с односторонней спектральной плотностью мощности N0: x(t) = s(t) + n(t). Тогда на выходе идеального ФНЧ будем иметь y(t) = s (t) + n (t), где s (t) описывается формулой (4), а n (t) - гауссовский центрированный случайный процесс
с корреляционной функцией К(т) = (N0 / 2) B sine (пВт). Тем не менее, и в этом случае для нахождения оценок
параметров Ak и tk, n = 0, K -1 , сигнала (2) можно воспользоваться аналогичным подходом, описанным выше, и применить перечисленные методы анализа [13, 14], а также метод Cadzow [15], специально разработанный для работы в шумах.
Результаты моделирования В среде Matlab было выполнено моделирование алгоритмов оценки временных положений tk и амплитудных множителей Ak сигнала (2), наблюдаемого на фоне аддитивного гауссовского шума. С этой целью формировались отсчёты процесса y(t) с шагом, обеспечивающим среднеквадратическую погрешность аппроксимации, не превышающую 10 %. По каждой сформированной реализации с использованием алгоритмов Annihilating Filter, Cadzow и MUSIC находились оценки tk и Ak (k = 1, К) при разных значениях К и отношениях сигнал/шум (ОСШ) z2, выраженного в децибелах (дБ). Далее определялись нормированные рассеяния оценок временных положений и амплитуд импульсов
1 К-1^ 42
=
KT k=0 1
£ (4 - h)
К-!/- \ 2 £ (Ak - Ak)
Va к £A;
k=0
для всех алгоритмов и для всех (не менее 10000) испытаний. Окончательные значения рассеяний Vt и VA получались путём усреднения по всем испытаниям (опытам). При моделировании использовались специальные функции, реализующие алгоритмы Annihilating Filter, Cadzow и MUSIC. Первые две из них были написаны самостоятельно, а для реализации третьего алгоритма использована функция rootmusic из библиотеки Matlab. На рис. 1 и 2 изображены зависимости этих нормированных рассеяний от ОСШ z2 (дБ). Кривые 1 построены с использованием метода Annihilating Filter, кривые 2 - с использованием метода Cadzow и кривые 3 - с использованием метода MUSIC. Сплошными линиями изображены соответствующие зависимости при K = 2, а штриховыми -при K = 7 (для рис. 1) и K = 5 (для рис. 2). ю->
ю-2
Ю-3 Ю-4 Ю-5 10~6
5 10 15 20 25 30 35 40
Рис. 1. Зависимость нормированного рассеяния оценки временных положений от отношения сигнал/шум
Анализ приведённых рисунков свидетельствует об ухудшении точности оценок рассматриваемых параметров с увеличением количества сигналов K. Кроме того, приведённые результаты свидетельствуют о том, что используемые методы оценки (Annihilating Filter, Cadzow, MUSIC) обладают приблизительно одинаковой точностью.
Vt
,1,2
3 Л. -W. 3
z2M
m=-M
5 10 15 20 25 30 35 40
Рис. 2. Зависимость нормированного рассеяния оценки амплитуд от отношения сигнал/шум
Действительно, на рис. 1 и рис. 2 кривые 1-3 отличаются друг от друга в пределах погрешности измерений. Поэтому в дальнейших исследованиях использовался только один из этих методов - алгоритм MUSIC. Следует, однако, отметить, что при моделировании временные положения импульсов выбирались хотя и случайным образом, но при некотором ограничении.
Это ограничение выглядит следующим образом: минимальное расстояние между двумя соседними импульсами должно быть не меньше, чем Г/10. Если уменьшить это минимальное расстояние между соседними импульсами, то точность измерения резко падает.
Приведённые результаты были получены в соответствии с (7) при наименьшем возможном количестве отсчётов N = 2K + 1. В работе было выполнено также моделирование аналогичных алгоритмов с увеличенным количеством отсчётов N. При этом было показано, что увеличение этого количества отсчётов приводит к повышению точности оценок соответствующих параметров. Например, при K = 5 увеличение количества отсчётов с минимально необходимых N = 2K+ 1 = 11 до 17 позволило уменьшить рассеяние оценок временных положений приблизительно в полтора раза при не очень малых значениях ОСШ.
2. Алгоритм оценки временных положений и амплитуд периодической последовательности сверхкоротких импульсов с использованием фильтра Найквиста
Постановка задачи
Одним из возможных недостатков изложенного подхода к оценке временных положений и амплитуд импульсов в (2) является использование идеального ФНЧ, который, как известно, не реализуем на практике. В связи с этим рассмотрим возможность замены этого фильтра на другие. В частности, используем для подобной замены фильтр, относящийся к классу фильтров Найквиста с амплитудно-частотной характеристикой (АЧХ) вида [16]
1, 1ш|<хВ(1 -Р),
H (ш) = <
(\
1 + cos
ш|-хВ(1 -Р) 2р8
(8)
п8(1 -р) <1ш|<хВ(1 + р),
0,|ш| > п8(1 + р).
Такой фильтр широко используется в качестве сглаживающего фильтра, предназначенного для компенсации искажений, вносимых каналом связи. Изменение параметра р этого фильтра приводит к изменению формы его АЧХ. Так, при р = 0 из (8) получаем рассмотренный ранее идеальный ФНЧ, а при р = 1 - это будет известный фильтр типа приподнятого косинуса.
Результаты моделирования Выполним моделирование рассмотренных алгоритмов оценки временных положений и амплитуд импульсов, заменив в них идеальный ФНЧ на фильтр Найквиста (8). При расчётах, очевидно, необходимо будет изменить характеристики шума на выходе этого фильтра, так как они будут отличаться от аналогичных характеристик на выходе идеального ФНЧ. Результаты моделирования при использовании фильтра Найквиста приведены на рис. 3-5.
10^
5 10 15 20 25 30 35 40
Рис. 3. Зависимость нормированного рассеяния оценки временных положений от отношения сигнал/шум 10°\
5 10 15 20 25 30 35 40
Рис. 4. Зависимость нормированного рассеяния оценки амплитуд от отношения сигнал/шум
ю-Я
10*
10-
10
-Р-
0,2 0,4 0,6 0,8 1,0
Рис. 5. Зависимость нормированного рассеяния оценки временных положений от параметра фильтра Найквиста в
Для исследования степени влияния параметров этого фильтра на характеристики оценок временных положений и амплитуд импульсов на этих же рисунках приведены соответствующие характеристики оценок и для идеального ФНЧ.
Так, на рис. 3 сплошными линиями изображены нормированные рассеяния оценок временных положений при использовании идеального ФНЧ, штриховыми линиями - нормированные рассеяния оценок временных положений при использовании фильтра Найквиста при в = 0,1, а штрихпунктирными линиями - аналогичные зависимости при р = 1. Кроме того, кривые 1 на этом рисунке построены при К = 2, а кривые 2 - при К = 5.
На рис. 4 приведены аналогичные зависимости, только для рассеяний оценок амплитуд импульсов.
Анализ приведённых зависимостей свидетельствует о следующем. Применение фильтра Найквиста, вместо идеального ФНЧ, вполне допустимо при выборе не очень больших значений параметра р. Так, независимо от количества импульсов и при любых (не очень малых) ОСШ, рассеяния оценок временных положений при использовании фильтра Найквиста незначительно отличаются от аналогичных характеристик при использовании идеального ФНЧ при р = 0,1. Увеличение параметра р приводит к ухудшению точности оценок временных положений импульсов, начиная примерно с г > 20 дБ. Что касается оценок амплитуд импульсов, то в этом случае выбор параметра фильтра Найквиста р оказывает несколько более существенное влияние на характеристики этих оценок. Замена идеального ФНЧ на фильтр Найквиста при р = 0,1 целесообразна при не очень малых значениях ОСШ (если К = 2 - при г2 > 20 дБ, в то же время, если К = 5, то уже при г2 > 30 дБ).
Рассмотрим более подробно влияние параметра фильтра Найквиста р на характеристики оценок временных положений импульсов (как наиболее информативных параметров по сравнению с амплитудами). Это влияние можно проследить из анализа рис. 5, на котором изображены нормированные рассеяния оценок временных положений импульсов при использовании как идеального ФНЧ (сплошные линии), так и фильтра Найквиста (штриховые линии). На этом рисунке кривые 1 построены при К = 2 и г2 > 20 дБ, кривые 2 - при К = 5 и г2 > 20 дБ, кривые 3 - при К = 2 и г2 > 40 дБ, кривые 4 - при К = 5 и г2 > 40 дБ.
Анализ приведённых зависимостей показывает, что замена идеального ФНЧ на фильтр Найквиста действительно целесообразна лишь при выборе малых значений параметра р. Например, при К = 2 параметр р < 0,3 ^ 0,5, а при К = 5 параметр р не должен превышать 0,1.
3. Оценка временных положений и амплитуд
импульсов последовательности классическим методом максимального правдоподобия
Рассмотрим теперь возможность решения поставленной задачи (оценка временных положений и амплитуд импульсов последовательности (2)) классиче-
ским методом максимального правдоподобия. Для этого в соответствии с (3) предположим, что наблюдению доступны N отсчётов аддитивной смеси полезного сигнала (6) и гауссовского центрированного случайного процесса с элементами корреляционной матрицы вида Kij = (N0 / 2) B sine (nBAt (i-j)), i, j = 1, N : Y = st + nt. Необходимо определить оценки параметров tk и Ak (k = 1, К). В соответствии с методом максимального правдоподобия ищется логарифм функции правдоподобия, который в нашем случае выглядит следующим образом:
N-1 N -1
Л[1, A] = £ £ (Y - m (t, A)) (Yj - mj (t, A)) Cjj , (9)
i=0 j=0
где щ (t, A) = s , Cij - элементы матрицы, обратной корреляционной матрице шума с элементами Kij.
Оценки амплитуд импульсов можно достаточно просто найти из решения системы уравнений
(ЭЛ / ЭД )|» = 0, l = 1, К. Действительно, вектор оценок амплитуд можно представить в виде A = uR-1, где u - вектор размера 1хК с элементами
N-1 N-1
и1 = ££ CjYG (jAt),
=0 j =0
R
Pg.1 PG.2... PQ.K -1
Pu... P1.K
Рк-1,0 РК-1,1 РК-1,2... a
К-1
матрица размера КхК с элементами
N -1 N -1
а = ЕЕ №) (
¡=0 ]=0
N-1 N -1
р^=ЕЕ №) ^ (т.
¡=0 ]=0
Значительно сложнее дело обстоит с оценкой временных положений импульсов 4 сигнала (2). Для этого найденные оценки амплитуд Ак необходимо подставить в (9), в результате получится функция Л, зависящая от К переменных к Необходимо осуществить поиск абсолютного максимума этой многомерной функции. Подобная задача может быть решена без существенных вычислительных затрат лишь для К = 1 и К = 2. Увеличение количества импульсов в наблюдаемых данных на единицу приводит к тому, что всякий раз размерность функции Л будет увеличиваться в N раз, где Nt - количество отсчётов по возможным временным положениям функции правдоподобия. В связи с вышеизложенным в работе была выполнена оценка по методу максимального правдоподобия параметров гк и
Ак (к = 0, К -1) только для случая К = 2. Как показали расчёты, точность этих оценок не существенно отличается от точности оценок анализируемых параметров, производимых предлагаемым в работе методом, что
подтверждает его высокую эффективность. Действительно, при K = 2 и z2 < 40 дБ рассеяние оценки временных положений при использовании метода максимального правдоподобия отличается от рассеяния оценки по предложенному методу не более, чем на 30 %. Отличие рассеяний оценок амплитуд для этих методов для тех же условий не превышает 50 %. Однако ещё раз подчеркнем, что практическая реализация алгоритма оценки временных положений и амплитуд импульсов по методу максимального правдоподобия уже при K > 3 весьма затруднительна (временные затраты возрастают настолько, что оценки не могут быть найдены в реальном режиме времени). Следовательно, несмотря на некоторый проигрыш в эффективности измерений, в реальных условиях целесообразно использовать предлагаемый в работе алгоритм вместо алгоритма максимального правдоподобия.
Отметим, что приведённые результаты справедливы в случае, если количество импульсов на периоде K априори известно. В то же время значительный интерес может представлять аналогичная задача при неизвестном количестве импульсов. Решение этой задачи может быть выполнено, например, с использованием теории Compressive Sensing [17-19], однако решение подобной задачи оставим на будущее.
Заключение
В работе показано, каким образом можно оценивать временные положения и амплитуды периодической последовательности импульсов по незначительному, зависящему лишь от числа этих импульсов количеству отсчётов. Определены условия, при которых вместо идеального ФНЧ можно использовать фильтр Найквиста. Показано, от чего зависит выбор параметров этого фильтра. Отмечено, что в отличие от метода максимального правдоподобия практическая реализация предложенного алгоритма значительно более простая, не требующая формирования и поиска максимума функции большого числа переменных. В то же время характеристики точности оценок временных положений и амплитуд импульсов при использовании этих двух методов отличаются незначительно. Показано, что, если точность оценки временных положений при использовании минимально необходимого числа отсчётов N = 2K + 1 недостаточна, можно несколько повысить эту точность путём незначительного увеличения количества этих отсчётов. Действительно, как показано ранее, увеличение количества отсчётов приблизительно в 1,5 раза (от 11 до 15-17) от минимально необходимого приводит к повышению точности измерения временных положений не менее чем в 1,5 раза (при не очень малых значениях ОСШ z2 > 20 дБ). Однако дальнейшее увеличение количества используемых отсчётов не приведет к существенному выигрышу в эффективности оценки временных положений. В то же время подобное увеличение количества отсчётов, как показывают исследования, практически не оказывает никакого влияния на эффективность оценки амплитуд. Таким образом, в
оптических системах с MPPM, в которых количество K передаваемых импульсов определяет не только энергетическую и частотную эффективности, но и скорость передачи информации, использование предлагаемого подхода позволяет осуществлять детектирование без существенных временных и аппаратных затрат с достаточно высокой точностью.
Литература
1. Радиосистемы передачи информации / И.М. Тепляков, Б.В. Рощин, А.И. Фомин, В.А. Вейцель; под ред. И.М. Теплякова. - М.: Радио и связь, 1982. - 264 с.
2. Гальярди, Р.М. Оптическая связь / Р.М. Гальярди, Ш. Карп; пер. с англ. - М.: Связь, 1978. - 424 с.
3. Wada, A. Enhancement of dynamic range of optical fiber sensor using fiber Bragg grating Fabry-Pérot interferometer with pulse-position modulation scheme: Compensation of source wavelength-sweep nonlinearity / A. Wada, S. Tanaka, N. Takahashi // IEEE Photonics Journal. - 2013. - Vol. 5, Issue 4. - P. 1-8. - DOI: 10.1109/JPHOT.2013.2276971.
4. Sahin, S. Improvement of required SNR for satellite-to-ground optical communication via aperture diversity by using PPM / S. §ahin // Proceedings of 10th International Symposium on Communication Systems, Networks and Digital Signal Processing (CSNDSP). - 2016. - P. 175179. - DOI: 10.1109/CSNDSP.2016.7573929.
5. Esman, D.J. Frequency-hopping pulse position modulation ultrawideband receiver / D.J. Esman, V. Ataie, B.P.-P. Kuo, N. Alic, S. Radic // Journal of Lightwave Technology. -2017. - Vol. 35, Issue 10. - P. 1894-1899. - DOI: 10.1109/JLT.2017.2672527.
6. Peppas, K.P. Performance of underwater optical wireless communication with multi-pulse pulse-position modulation receivers and spatial diversity / K.P. Peppas, A.C. Boucouvalas, Z. Ghassemloy // IET Optoelectronics. - 2017. - Vol. 11, Issue 5. - P. 180-185. - DOI: 10.1049/iet-opt.2016.0130.
7. Hassan, M. Error performance study of MPPM optical communication systems with finite extinction ratios / M. Hassan, S. Shapsough, T. Landolsi, A.F. Elrefaie // Proceedings of International Conference on Industrial Informatics and Computer Systems (CIICS). - 2016. - P. 1-4. -DOI: 10.1109/ICCSII.2016.7462445.
8. Zhang, P. Wagner-like decoding for noncoherent PPM based ultra-low-power communications / P. Zhang, F.M.J. Willems, L. Huang // Proceedings of IEEE 24th International Symposium on Personal Indoor and Mobile Radio Communications (PIMRC). - 2013. - P. 932-937. -DOI: 10.1109/PIMRC.2013.6666271.
9. Tur, R. Innovation rate sampling of pulse streams with application to ultrasound imaging / R. Tur, Y.C. Eldar, Z. Friedman // IEEE Transactions on Signal Processing. - 2011. - Vol. 59, Issue 4. - P. 1827-1842. - DOI: 10.1109/TSP.2011.2105480.
10. Katsu, K. Performance evaluation on correlation detection and energy detection for ultra wideband-impulse radio communication with multi-pulse position modulation scheme in implant body area networks / K. Katsu, D. Anzai, J. Wang // IET Communications. - 2013. - Vol. 7, Issue 13. - P. 1430-1436. - DOI: 10.1049/iet-com.2013.0015.
11. Зубаченко, В.Л. Разработка систем кодирования для оптоэлектронных атмосферных каналов (ОАК) и волоконно-оптических линий связи (ВОЛС) / В.Л. Зубаченко, В.П. Дмитриев, П.А. Бобович, И.В. Коршунов // Компьютерная оптика. - 2006. - № 30. - С. 80-91.
12. Оппенгейм, А. Цифровая обработка сигналов / А.В. Оп-пенгейм, Р.В. Шафер; под ред. А.С. Ненашева; пер. с англ. - М.: Техносфера, 2006. - 856 с.
13. Марпл.-мл., С.Л. Цифровой спектральный анализ и его приложения / С.Л. Марпл.-мл.; пер. с англ. - М.: Мир, 1990. - 584 с.
14. Ye, J.C. Compressive sampling using annihilating filter-based low-rank interpolation / J.C. Ye, J.M. Kim, K.H. Jin, K. Lee // IEEE Transaction on Information Theory. - 2016. - Vol. 63, Issue 2. - P. 777-801. - DOI: 10Л109/ГТГ.2016.2629078.
15. Condat, L. Cadzow denoising upgrated: A new projection method for the recovery of Dirac pulses from noisy linear measurements / L. Condat, A. Hirabayashi // Sampling Theory in Signal and Image Processing. - 2015. - Vol. 14, Issue 1. - P. 17-47.
16. Çafak, M. Digital communications / M. §afak. - Chichester: John Wiley and Sons Ltd., 2017. - 906 p. - ISBN: 9781-11909125-7.
17. Bellasi, D.E. Energy-efficiency analysis of analog and digital Compressive Sensing in wireless sensors / D.E. Bellasi,
L. Benini // IEEE Transactions on Circuits and Systems I: Regular Papers. - 2015. - Vol. 62, Issue 11. - P. 27182729. - DOI: 10.1109/TCSI.2015.2477579.
18. Donoho, D.L. Compressed sensing / D.L. Donoho // IEEE Transactions on Information Theory. - 2006. - Vol. 52, Issue 4. - P. 1289-1306. - DOI: 10.1109/TIT.2006.871582.
19. Парфенов, В.И. Принципы построения и анализ эффективности функционирования беспроводных сенсорных сетей на основе теории Compressive Sensing / В.И. Парфенов, Д.Ю. Голованов // Цифровая обработка сигналов. - 2016. - № 2. - С. 14-19.
20. Парфенов, В.И. Определение незанятых частотных поддиапазонов по сжатым измерениям / В.И. Парфенов, Д.Ю. Голованов // Инфокоммуникационные технологии. - 2015. -Т. 13, №3. - С. 305-312.
Сведения об авторах
Парфенов Владимир Иванович, 1961 года рождения, в 1983 окончил Воронежский государственный университет, в 1989 году защитил кандидатскую, а в 2002 - докторскую диссертацию. Доктор физико-математических наук, профессор, профессор кафедры радиофизики Воронежского государственного университета. Область научных интересов: применение методов статистической радиофизики к анализу систем передачи информации. E-mail: vip@phys.vsu.ru .
Голованов Дмитрий Юрьевич, 1988 года рождения, в 2012 году окончил Воронежский государственный университет. Аспирант кафедры радиофизики Воронежского государственного университета. Область научных интересов: разработка и исследование новых методов сжатия цифровой информации. E-mail: golovanov@amm.vsu.ru .
ГРНТИ: 49.03.05
Поступила в редакцию 6 ноября 2017 г. Окончательный вариант - 9 декабря 2017 г.
NOISE STABILITY OF SIGNAL RECEPTION ALGORITHMS WITH MULTIPULSE PULSE-POSITION MODULATION
V.I. Parfenov1, D. Y. Golovanov1
1Voronezh State University, Voronezh, Russia
Abstract
An algorithm for estimating time positions and amplitudes of a periodic pulse sequence from a small number of samples was proposed. The number of these samples was determined only by the number of pulses. The performance of this algorithm was considered on the assumption that the spectrum of the original signal is limited with an ideal low-pass filter or the Nyquist filter, and conditions for the conversion from one filter to the other were determined. The efficiency of the proposed algorithm was investigated through analyzing in which way the dispersion of estimates of time positions and amplitudes depends on the signal-to-noise ratio and on the number of pulses in the sequence. It was shown that, from this point of view, the efficiency of the algorithm decreases with increasing number of sequence pulses. Besides, the efficiency of the proposed algorithm decreases with decreasing signal-to-noise ratio.
It was found that, unlike the classical maximum likelihood algorithm, the proposed algorithm does not require a search for the maximum of a multivariable function, meanwhile characteristics of the estimates are practically the same for both these methods. Also, it was shown that the estimation accuracy of the proposed algorithm can be increased by an insignificant increase in the number of signal samples.
The results obtained may be used in the practical design of laser communication systems, in which the multipulse pulse-position modulation is used for message transmission.
Keywords: optical communications, modulation, optical systems, lasers, pulses, filters, estimation algorithms, method MUSIC, signal to noise ratio, dispersion of estimates, maximum likelihood method, likelihood function.
Citation: Parfenov VI, Golovanov DY. Noise stability of signal reception algorithms with multipulse pulse-position modulation. Computer Optics 2018; 42(1): 167-174. DOI: 10.18287/24126179-2018-42-1-167-174.
References
[1] Teplyakov IM, Roshcin BV, Fomin AI, Weitzel VA. Radio communication systems [In Russian]. Moscow: "Radio i svjaz" Publisher; 1982.
[2] Gagliardi RM, Karp S. Optical communications. New York: John Wiley & Sons, Inc.; 1974.
[3] Wada A, Tanaka S, Takahashi N. Enhancement of dynamic range of optical fiber sensor using fiber Bragg grating Fabry-Pérot interferometer with pulse-position modulation scheme: Compensation of source wavelength-sweep nonlinearity. IEEE Photonics Journal 2013; 5(4): 1-8. DOI: 10.1109/JPHOT. 2013.2276971.
[4] §ahin S. Improvement of required SNR for satellite-to-ground optical communication via aperture diversity by using PPM. Proceedings of 10th International Symposium on Communication Systems, Networks and Digital Signal Processing (CSNDSP) 2016; 175-179. DOI: 10.1109/CSNDSP.2016.7573929.
[5] Esman DJ, Ataie V, Kuo BP-P, Alic N, Radic S. Frequency-hopping pulse position modulation ultrawideband receiver. Journal of Lightwave Technology 2017; 35(10): 1894-1899. DOI: 10.1109/JLT.2017.2672527.
[6] Peppas KP, Boucouvalas AC, Ghassemloy Z. Performance of underwater optical wireless communication with multipulse pulse-position modulation receivers and spatial diversity. IET Optoelectronics 2017; 11(5): 180-185. DOI: 10.1049/iet-opt.2016.0130.
[7] Hassan M, Shapsough S, Landolsi T, Elrefaie AF. Error performance study of MPPM optical communication systems with finite extinction ratios. Proceedings of International Conference on Industrial Informatics and Computer Systems (CIICS) 2016; 1-4. DOI: 10.1109/ICCSII.2016.7462445.
[8] Zhang P, Willems FMJ, Huang L. Wagner-like decoding for noncoherent PPM based ultra-low-power communications. Proceedings of IEEE 24th International Symposium on Personal Indoor and Mobile Radio Communications (PIMRC) 2013; 932-937. DOI: 10.1109/PIMRC.2013.6666271.
[9] Tur R, Eldar YC, Friedman Z. Innovation rate sampling of pulse streams with application to ultrasound imaging. IEEE Transactions on Signal Processing 2011; 59(4): 1827-1842. DOI: 10.1109/TSP.2011.2105480.
[10] Katsu K, Anzai D, Wang J. Performance evaluation on correlation detection and energy detection for ultra wideband-impulse radio communication with multi-pulse position modulation scheme in implant body area networks. IET Communications 2013; 7(13): 1430-1436. DOI: 10.1049/iet-com.2013.0015.
[11] Zubachenco VL, Dmitriev VP, Bobovich PA, Korshunov IV. Development of coding systems for optoelectronic atmospheric channels (OAC) and fiber-optic communication lines (FOCL) [In Russian]. Computer Optics 2006; 30: 80-91.
[12] Oppenheim AV, Schafer RW. Digital signal processing. Englewood Cliffs, New Jersey: Prentice Hall, Inc.; 1999.
[13] Marple SL. Digital spectral analysis: With applications. Upper Saddle River, New Jersey: Prentice Hall, Inc.; 1987. ISBN: 0-132-14149-3.
[14] Ye JC, Kim JM, Jin KH, Lee K. Compressive Sampling using annihilating filter-based low-rank interpolation. IEEE Transaction on Information Theory 2016; 63(2): 777-801. DOI: 10.1109/TIT.2016.2629078.
[15] Condat L, Hirabayashi A. Cadzow denoising upgrated: A new projection method for the recovery of Dirac pulses from noisy linear measurements. Sampling Theory in Signal and Image Processing 2015; 14(1): 17-47.
[16] §afak M. Digital communications. Chichester: John Wiley and Sons Ltd.; 2017. ISBN: 978-1-11909125-7.
[17] Bellasi DE, Benini L. Energy-efficiency analysis of analog and digital Compressive Sensing in wireless sensors. IEEE Transactions on Circuits and Systems I: Regular Papers 2015; 62(11): 2718-2729. DOI: 10.1109/TCSI.2015.2477579.
[18] Donoho DL. Compressed Sensing. IEEE Transactions on Information Theory 2006; 52(4): 1289-1306. DOI: 10.1109/TIT.2006.871582.
[19] Parfenov VI, Golovanov DY. Principles of construction and effectiveness analysis of operation of wireless sensor networks based on compressive sensing theory [In Russian]. Digital signal processing 2016; 2: 14-19.
[20] Parfenov VI, Golovanov DY. Detection of vacant unoccupied frequency sub-bands by compressed measurements [In Russian]. Infocommunication Technologies 2015; 13(3): 305-312.
Author's information
Vladimir Ivanovich Parfenov (b. 1961) graduated from Voronezh State University in 1983. In 1989 he defended his candidate's dissertation and in 2002 - his doctoral dissertation. Doctor of Physical and Mathematical sciences, professor, professor of Radiophysics chair, Voronezh State University. Research interests are application statistical radio-physics methods to the analysis of information transmission systems. E-mail: vip@phys.vsu.ru .
Dmitry Yurievich Golovanov (b. 1988) graduated from Voronezh State University in 2012. Graduate student of Radiophysics chair, Voronezh State University. Research interests are development and investigation of new compression methods of digital information. E-mail: golovanov@amm.vsu.ru .
Received November 6, 2017. The final version - December 9, 2017.
Дизайн: Я.Е. Тахтаров. Оформление и вёрстка: М.А. Вахе, Е.В. Семиколенных, С.В. Смагин и Я.Е. Тахтаров. Лит. редактор и корректор Ю.Н. Литвинова. Консультант по оформлению англоязычного блока М.И. Котляр. E-mail: ko@smr. ru, http://www.computeroptics.smr.ru
Подписано в печать 26.02.2018 г. Усл. печ. л. 20,10. Заказ № 11/1. Тираж 218 экз. Печать офсетная. Формат 62х84 1/8. Цена: 550 рублей / Price of 550 rubles (6+)
Редакция: Институт систем обработки изображений РАН - филиал ФНИЦ «Кристаллография и фотоника» РАН, (443010, г. Самара, ул. Молодогвардейская, 151) Соучредители: Федеральное государственное автономное образовательное учреждение высшего образования «Самарский национальный исследовательский университет имени академика С.П. Королева» (443086, г. Самара, Московское шоссе, д.34), Федеральное государственное учреждение «Федеральный научно-исследовательский центр «Кристаллография и фотоника» Российской академии наук» (117342, г. Москва, ул. Бутлерова, д 17А) Отпечатано в типографии ООО «Предприятие «Новая техника» (443013 г. Самара, пр-кт. Карла Маркса, 24-76)