Оптимальный алгоритм приема М-ичных частотно-манипулированных сигналов при использовании модели несущей в пространстве состояний

Размещено на http://www.allbest.ru/

ФГБОУ ВПО «Российский государственный университет туризма и сервиса»

Оптимальный алгоритм приема МЧМН-сигналов при использовании модели несущей в пространстве состояний

Вадим Михайлович Советов, д.т.н., с.н.с., проф., каф. «Информационные системы и технологии»

Дарья Анатольевна Рыженкова, аспирант

Москва

Аннотация

Синтезирована схема оптимального приема М-ичных частотно-манипулированных (МЧМн) сигналов при представлении несущих сигналов динамической моделью в пространстве состояний и использовании вектора начального состояния в качестве информационного вектора; показано, что оптимальный приемник МЧМн-сигналов представляет собой М каналов, в которых для вычисления апостериорной плотности вероятности М гипотез применяются фильтры Калмана, и данный алгоритм может использоваться также для приема сигналов с несинусоидальной несущей, при представлении их динамической моделью; получена формула расчета вероятности ошибки и приведены результаты этого расчета и моделирования.

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

The authors synthesized a scheme of an optimal reception of M frequency-manipulated signals in submitting carrier signals of dynamic model in the space of states and with the using of the initial state vector as an informational vector. It is shown that the optimal receiver of M frequency-manipulated signals represents M channels in which for calculating the posterior of probability of the M hypotheses, the Kalman filters are applied. Except that, the algorithm can also be used for receiving signals from the non-sinusoidal carrier, during the presentation of their dynamic model. The authors obtained a formula for calculating the probability of error and presented the results of this calculation and simulation.

Keywords: dynamic model, the Kalman filter, the optimal reception, the state space, frequency shift keying.

При М-ичной частотной манипуляции (МЧМн) для передачи информации битовый поток разбивается на пакеты из n символов и каждый пакет передается на отдельной частоте из M возможных, при этом бит. Все M возможных пакетов и, соответственно, частот обозначим как M сообщения: , .

При МЧМн каждый из М сигналов разной частоты можно записать как [1]

, для , ,

где Т - длительность пакета из n символов.

Если начальная фаза одинакова для всех i, то сигналы когерентны. Для ортогональности когерентных сигналов минимальное частотное разделение между двумя любыми соседними когерентными сигналами должно быть , а соседними некогерентными . В общем случае когерентные сигналы могут отстоять друг от друга на , а некогерентные на . Обычно выбирается одинаковое частотное разделение между двумя соседними частотами [1].

Цель статьи - провести анализ помехоустойчивости приема МЧМн-сигналов при использовании динамической модели несущих сигналов в пространстве состояний (ПС).

Известно, что гармоническое колебание можно представить динамической моделью в ПС в виде уравнения смены состояний и уравнения наблюдения состояний [2]

,

,

где - матрица перехода состояний (МПС), - круговая частота; - интервал дискретизации; - вектор состояния (ВС); - вектор наблюдения; - матрица выхода уравнения наблюдения.

В работе [3] была рассмотрена задача оптимального приема одного из М сигналов, представленного моделью в ПС в общем виде в дискретном времени. В данном же случае начальный двухмерный вектор xi(0), представляет собой случайный вектор сообщения, передаваемого на одной из М частот, и имеет детерминированные, но неизвестные на приемной стороне значения элементов. Тогда имеем случай выбора одной среди М гипотез , ,…, в любой момент времени k пока длится эксперимент:

,

, .

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

- математическое ожидание;

- ковариационная матрица.

Решение о выборе той или иной гипотезы принимается по k последовательным наблюдениям:

.

В результате решения задачи оптимального приема М-ичных сигналов при представлении несущей моделью в ПС (1), (2) показано, что оптимальный приемник вместо коррелятора или согласованного фильтра содержит фильтр Калмана (ФК) [3]:

, (3)

, (4)

, (5)

, (6)

где - оценка вектора состояния в i-м канале приема; - матрица ковариации ошибки оценки в i-м канале приема; - весовая матрица.

Данный фильтр стоит в каждом канале приема, и для каждой гипотезы вычисляются в отдельном канале приема по алгоритму Калмана (3) - (6) оценка и матрица ковариации ошибки оценки . На рис. 1 показаны результаты моделирования работы алгоритма (3) - (6) для двух каналов приема при совпадении (рис. 1, а) и несовпадении (рис. 1, б) частоты несущей с частотой канала. Как видно, в отличие от коррелятора, в котором амплитуда сигнала постоянно возрастает, ФК выдает оценку несущей с амплитудой близкой к амплитуде несущего сигнала, при этом отношение сигнал/шум (ОСШ) возрастает за счет уменьшения дисперсии оценки, а не роста амплитуды.

Рис. 1. Сигнал на выходе ФК: а - при совпадении частоты канала с частотой несущей; б - при несовпадении частоты канна с частотой несущей (кривая 1 - несущая частота; кривая 2 - оценка несущей; кривая 3 - разностная частота)

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

имеет максимальное значение [3].

Рассмотрим симметричный случай, когда М сигналов имеют одинаковую энергию и равную априорную вероятность. В этом случае можно рассчитать вероятность ошибки при справедливости любой гипотезы, например .

Ошибка произойдет, если , , :

Поскольку гипотезы независимы и имеют одинаковые плотности вероятности, когда , вероятность ошибки на частотный символ вычисляется в соответствии с (8) как

Для нахождения необходимо вычислить достаточную статистику по формуле (7) при :

.

Для нахождения необходимо вычислить по формуле (7) при :

.

Рассмотрим составляющую

, ,

входящую как в (10), так и в (11). На рис. 2 показаны результаты расчета для случая четырех гипотез при равной энергии частотных символов. Как видно, при использовании частот с одинаковой амплитудой энергии ВС примерно равны и при равных матрицах ковариации ошибки оценки можно принять, что для всех , и в выражениях (10), (11) опустить эти члены, так как они приводят к одинаковому смещению . Тогда можно записать выражения как

,

.

Рис. 2. Составляющие Sj(k) решающего правила на выходе четырех каналов приема МЧМн-сигнала при равной энергии частотных символов

Теперь приемник должен сравнить между собой все (12), (13) и выбрать максимальную, при этом для гипотезы он будет использовать решающее правило:

,

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

Рассмотрим левую часть неравенства (14) и представим оценку в виде

.

Подставив (15) в (12) получим:

.

Раскроем скобки и представим выражение в следующем виде:

.

Математическое ожидание (16)

.

По определению, дисперсия

.

Раскрыв скобки, получим:

.

По определению и с учетом симметрии матрицы ковариации , окончательно получим:

.

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

.

Рассмотрим далее правую часть неравенства (14)

.

Значение в точках ортогональности сигналов равно нулю. Более того, при увеличении k весовая матрица Калмана стремится к нулю и стремится к нулю (см. рис. 1, б). Поэтому математическое ожидание будет

.

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

.

Вероятность ошибки на частотный символ в соответствии с (9), (17) и (18) будет. частотный сигнал манипулированный информационный

.

Путем замены переменных выражение (20) можно представить в виде

Для сокращения записи используем функцию

,

Откуда

.

Верхнюю границу вероятности ошибки можно записать как

.

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

,

где - энергия, приходящая на символ; - энергия, приходящая на один бит (один элемент символа).

Вероятность символьной ошибки для равновероятных М-ичных ортогональных сигналов с некогерентным детектированием определяется следующим выражением [1]:

,

где - биномиальный коэффициент.

Для бинарного случая формула сокращается до такой:

.

На рис. 3 показаны результаты расчета вероятности ошибки на бит по формулам (21) - (25) при некогерентном (кривые 1) и когерентном приеме (кривые 2) в зависимости от отношения в течение времени, т.е. от количества накапливаемых отсчетов. Как видно, вероятность ошибки при использовании ФК колеблется около значений кривой при когерентном приеме для заданного значения . Колебания объясняются тем, что в выражении (22) ОСШ определяется функционалом в зависимости от числа отсчетов, в котором представляет сумму двух квадратов отсчетов синусоиды, а в пределе представляет диагональную матрицу, члены которой стремятся к обратной величине дисперсии оценки.

Рис. 3. Вероятность ошибки на символ для М=2 (верхние кривые) и М=16 (нижние кривые): кривая 1 - некогерентный корреляционный прием; кривая 2 - когерентный корреляционный прием; кривая 3 - некогерентный прием с ФК

В связи с этим с течением времени, т.е. с увеличением числа накапливаемых отсчетов (см. рис. 2), наблюдаются колебания ОСШ и, соответственно, колебания вероятности ошибки. Следует отметить также, что случаи, когда вероятность ошибки при использовании ФК становится меньше вероятности ошибки при когерентном приеме, не означают, что превышена потенциальная помехоустойчивость. На рис. 3 кривые 3 показывают процесс накопления энергии со временем. При одном и том же времени накопления в случае пропорционального (как это имеет место при прямоугольном импульсе) и непропорционального (как это имеет место при накоплении энергии синусоиды) увеличения энергии времени накопления отношения энергий будут различны. При одинаковых энергиях помехоустойчивость некогерентного приема с использованием ФК будет хуже, чем помехоустойчивость при использовании когерентного корреляционного приема, однако, как это видно из графиков, лучше, чем при использовании некогерентного корреляционного приема. Последнее объясняется тем, что ФК накапливает энергию отрезка синусоиды независимо от ее начальной фазы (см. рис.1, а), поэтому нет необходимости в использовании двух квадратурных каналов, которые удваивают дисперсию шума.

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

Литература

1. Xiong Fuqin, Digital Modulation Techniques, Second Edition. ARTECH HOUSE, INC. 2006.

2. Советов В.М. Оптимальный прием частотно манипулированных сигналов на основе динамической модели несущей частоты // Электротехнические и информационные комплексы и системы. 2010. Т.6. №3. С. 14 - 19.

3. Советов В.М., Медведев П.А. Оптимальный алгоритм приема М сигналов при использовании модели несущей в пространстве состояний // Электротехнические и информационные комплексы и системы. 2010. Т.6. №4. С. 7 - 12.

Размещено на Allbest.ru