Математические модели сигналов систем БШД, использующих технологию MIMO Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

Математические модели сигналов систем БШД использующих технологию MIMO

Технология MIMO (Multiple Input Multiple Output - так называемая технология множественного входа и множе-ственного выхода) нашла широкое практическое применение в сетях беспроводного широкополосного доступа стандарта IEEE 802.11 n (Wi-Fi), IEEE 802.16 (WiMAX), LTE и недавно вышедшего стандарта 802.22 (White Space). Главное достоинство технологии - возможность увеличения пропускной способности, ёмкости и надёжности сети пропорционально увеличению количества передающих и приёмных антенн. При этом технология позволяет отчасти, либо полностью решить такие проблемы, присущие всем беспроводным сетям, как интерференция, потери при распространении радиоволн, быстрые и медленные замирания сигнала, затенение, экранирование. Технология благодаря своим достоинствам имеет большое будущее, но должна постоянно совершенствоваться, поэтому вопросам её изучения, в том числе моделирования передаваемых и принимаемых сигналов, необходимо уделять должное внимание. В вышеназванных современных беспроводных системах технология Ключевые слова: математическая MIMO используется вместе с технологией ортогонального кодирования OFDM (Orthogona

модель, беспроводной Frequency Division Multiplexing), которая предполагает как можно более эффективное

широкополосный доступ, Wi-Fi, использование спектра за счёт параллельной передачи данных по множеству близко

WiMAX, LTE, MIMO, OFDM расположенных друг к другу узкополосных частотных подканалов.

Передатчик

Прнсмиик

Бухарев ИА,

ведущий инженер ФГУП ЦНИИС, bukharev@zniis.ru

Математические модели принимаемых сигналов и помех систем БШД, использующих технологию MIMO

Разработка математических моделей (ММ) принимаемого полезного сигнала и помех - это первый этап решения задачи синтеза алгоритмов приёма и обработки OFDM-сигналов в абонентском терминале (АТ), работающем с использованием технологии MIMO, методами марковской теории оптимального оценивания случайных процессов [I]. Разрабатываемые математические модели должны в полной мере соответствовать параметрам реальных сигналов и помех, наблюдаемых на входе абонентского терминала, работающего в сетях Wi-Fi/WiMAX/LTE/White Space с использованием технологий MIMO и OFDM.

Согласно рисунку 1, MIMO-сигнал на приёмной стороне в общем виде записывается следующим образом: Х(1)=Н S(t)+Z(t) (1),

где Х(1) — матрица принятого сообщения; Н - канальная матрица размера М x N (М - число приёмных антенн, N - число передающих антенн); S(t) - матрица передаваемых OFDM-сигналов; Z(t) - матрица собственных шумов приёмных элементов антенны.

В формуле (1) матрица передаваемых OFDM-сигналов S(t) имеет размер /Vxl, зависящий только от числа передающих антенн:

«,(/)

S(D =

f2(0

sN(t)

Рис. 1. Модель MIMO-канала

Через каждую антенну передаётся свой модулированный OFDM-сигнал. OFDM-сигнал формируется путём преобразования последовательного потока бинарных символов S[p] в Р параллельных потоков, каждый из которых согласуется (в том числе, и комплексно) с исходным потоком S[p], используя различные виды модуляции, с помощью инверсного мультиплексирования (демультиплексирования) [2].

Независимо от способа модуляции (QPSK, QAM или OQAM) модулированное колебание представляет собой сумму синфазной компоненты (косинусоиды) и квадратурной компоненты (синусоиды) [2]. Амплитуда синфазной компоненты равна вещественной части Re{X'/=X* нормированного комплексного модуляционного символа X'. Амплитуда квадратурной компоненты равна мнимой части 1т/Х'}=Х** этого символа [2]. В общем случае модуляция разных потоков может быть различной, поэтому некоторые потоки могут нести битовую скорость больше, чем у других. Модуляционный символ X' представляет собой совокупность модулированных ортогональных несущих частот. Его нормировка как раз и вводится для того, чтобы средние мощности колебаний с разными способами модуляции были одинаковыми. Так, для модуляции QPSK нормированный комплексный модуляционный символ Х'=Х/(2)' а для модуляции 16QAM -Х'=Х/(10)‘2 [3].

Сигнал несущей с номером к и частотой у*, модулированной символом Х\, может быть представлен в виде вещественной части произведения Х\ и комплексного колебания (комплексной экспоненты) с частотой /* [2]:

2 Art—

1

Хк(!) = ?л{Х\е‘,жи') = Re^V* (2)

где j = V-T» ft = \/Тп -к-я гармоника основной частоты

(величины, обратной длительности полезной части символа и равной расстоянию между частотами соседних несущих).

Сигнал OFDM, записанный на интервале одного символа, представляет собой сумму всех несущих колебаний, модулированных своими нормированными модуляционными символами [2]:

*(о=X хк (о=Z Re{* v е 1Г*Т‘} • (3>

*=0 *=0

где К - общее количество поднесущих.

Учитывая, что цифровая система передачи данных является системой с дискретным временем I = пТ (Г — интервал дискретизации; п - номер отсчёта) [4], получаем из выражения (3) формулу:

X(nT) = X„ = Re'£{X\e

ZJXknlr

1

(4)

Полагая, что T</T=N, получаем обратное дискретное преобразование Фурье для вычисления модулированных сигналов и переводом их цифровых отсчётов во временную область:

Jf.-RegiJirV^)- (5)

*=0

К-I 2 irk—

Вычисление мнимой части выражения \{Х\ е ' }

*=о

наряду с вычислением его вещественной части по формуле (5) используется для переноса спектра сигнала OFDM в заданный диапазон частот и в общем виде даёт следующий общий результат:

X(t) = ^\X\elJ’f"} = a + jb-*=о

Данные составляющие подвергаются цифроаналоговому преобразованию (ЦАП). Полученные аналоговые сигналы используются для модуляции соответственно синусоиды и косинусоиды (получаемой сдвигом синусоиды на 90°) несущей частоты/. [2].

Умножим вещественную часть а на синфазное колебание cos(2irfct) с частотой /., а мнимую часть Ь — на квадратурное колебание sin(2nfct) (сдвинутое по фазе на 90° по отношению к синфазному). Тогда суммирование полученных произведений даёт сигнал OFDM S(l), спектр которого смещён на частоту fc по отношению к спектру сигнала X(t). Данный сигнал после суммирования передаётся в линию связи [3]:

S(0=a cos(2xft+%)+jb sin( 2 кfjt+%) =

, (о)

А-I » V /

= X № (r~'s~ rXcos[2^/, +f)t+%\+js\n[2n(fk +l)r + <p:,})} где Sk(t-i, -г) _ номинальная амплитуда принимаемого полезного сигнала, определяемая закономерностью распространения радиоволн в свободном пространстве; lg - защитный интервал, добавляемый в начало каждого передаваемого символа, Г - время задержки передаваемого сигнала, (р0 - случайная начальная фаза сигнала.

Используя формулу Эйлера [4], а также с учётом того, что частота дискретизации СОк = 2/г([к + /с), выражение (6) можно привести к следующему виду:

5(0 = £ 5* (/-*,-г)-(7) *=0

Передача ОРЭМ-символа начинается не в момент времени 1=О, а в момент /я. До этого, в интервале [0.../^], передаётся периодическое повторение символа. Для упрощения представления сигнала 8(1) значение защитного интервала 1К можно не учитывать. Тогда выражение (1.7) принимает следующий вид:

(8)

Вследствие разности промежуточных частот приёмника и передатчика возникает случайный сдвиг сигнала по частоте:

Да)к = ц ~(ог- ко)\ + с!(о, (9)

где к(0\ - сдвиг частоты на целое количество гармоник (0)\ - расстояние между поднесущими в спектре сигнала), с!(О - сдвиг сигнала по частоте на значение, которое не превосходит СОк. В этом случае выражение (8) записывается следующим образом:

5(0=Х5*('-Г)'е'<л‘ц'+Л>- (10)

*=0

Время задержки радиосигнала на трассе от передатчика до абонентского терминала может быть представлено в виде:

г^+Д^ + Д/,, (П)

где г, = — - задержка сигнала, обусловленная дально-

1 с

стыо трассы с/ (с - скорость распространения радиоволн); 1тр - тропосферная задержка сигнала; Д/(/ -

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

Канальная матрица Н сигнала из выражения (1) состоит из комплексных коэффициентов И\м> которые не являются постоянными, и записывается в следующем виде:

Н =

ч К \ h]\

К\ л,, | К

J’m А*л І h\t\

Для беспроводных систем характерны амплитудные и фазовые искажения, обусловленные большим количеством факторов, в том числе эффектом Доплера и уходом частоты опорного генератора [5]. Такие искажения приводят к возникновению явления «замираний» сигнала. Значения коэффициентов Иш случайно изменяются во времени именно вследствие их зависимости от случайных амплитуды и фазы, доплеровский сдвиг которых обусловлен параметрами передающего устройства, многолучевым распространением сигнала, расположением и перемещением абонентского терминала, характером рельефа местности, плотностью застройки, наличием атмосферных осадков, листвы на деревьях, затуханием

сигнала в туннелях и другими помехами искусственного и естественного происхождения [1].

При перемножении в выражении (1) матрицы Н и матрицы 3(1) [4] получим сигнал в общем виде, наблюдаемый на входе приёмного устройства без учёта шумовой составляющей:

бег \'и'тиш\и

Д,) = //•£(,) =

ч ь* ! К,'

К К \ А:л

К\ Кг \ А

5,(0

52(1)_

МО

А, ,5, (/) + Л1252(/)+...+А,лЛ'(0 , л, (/) + /?„ (/) + ... + /)2Л л-л (I)

(/)+Л*\ (0

XV, (о ХМЛО

Хл»в*до

(12)

(15)

где С ' - параметр релеевской плотности вероятности.

Применительно к плотности вероятности математическое ожидание /И, для огибающей £*(/) равно [4]:

т

«Л

Дисперсия <Т7, для огибающей ^к(1) равна [4]:

2 / ■) о\=и--)о-

Амплитуда импульса Д принимаемого сигнала представляет собой дискретный случайный процесс, области значений и определения которого представляют собой дискретные множества [5]. В процессе синтеза алгоритма приёма и обработки сигнала такой процесс аппроксимируется простой цепью Маркова на несколько положений [1].

Матрица собственных шумов приёмных элементов антенны 2(1) имеет размер Мх 1, зависящий только от числа приёмных антенн:

*,(0

(16)

2(1) =

_5з(_0

*>,(')

где N - число передающих антенн, М - число приёмных антенн.

По сути, значения коэффициентов Лдда представляют собой импульсную характеристику каждого передаваемого сигнала. Общее выражение для суммарной импульсной характеристики канала одной антенны записывается следующим образом [3]:

И(т,0 = '£/Зк(0е~т,)Я(т-тк(0), (13)

*=о

где /Зк, вк - случайные амплитуда и фаза импульса, Тк

- момент появления на приёмном конце импульса при распространении по одному из лучей, / - текущее время (так как параметры канала могут меняться со временем).

Для упрощения рассуждений импульсную характеристику канала, которая создаётся при распространении сигнала, можно обозначить И(1). Тогда принятый сигнал одной антенны с учётом выражения (10) и без учёта шума имеет вид:

5(/) = £ Д (1)е^{пд(т-тк (0)5* (/ - т)ен^'*п) =

*=° (14)

= X А (05* (/ - т)6( т - тк (0У<Да|'+лн9*

*=<>

Формула справедлива в том случае, если импульсная характеристика короче защитного интервала, то есть Тк < /, тогда символы не перекрываются, а также при

стационарности импульсной характеристики в течение символа, то есть /=со/и/ [6].

Амплитуда принимаемых ОРПМ-сигналов подчиняется закону распределения Релея, при этом одномерная плотность вероятности амплигуды 5* (/) имеет вид [4]:

Такие шумы возникают из-за случайного движения электронов вследствие тепловой активности в приёмнике и называются тепловыми [7]. Эти помехи имеют Гауссовское распределение и равную мощность для всех частотных компонент, то есть спектральная плотность

N

мощности шума постоянна для всех частот и равна —2.

2

(где NI) - интенсивность шума), при этом __2. называют

2

двухсторонней шумовой спектральной плотностью [7]. Таким образом, в общем случае можно считать, что на вход приёмника воздействует стационарный аддитивный белый гауссовский шум (АБГШ) с нулевым средним значением. Статистические характеристики такого шума считаются известными [4]. Математическое ожидание (МО):

М{2(/)} = 0

Автокорреляционная функция (АКФ):

М{2(!1 )2(/2)} = (У18(1, -/,)=^(/, -12),

где N0 — интенсивность АБГШ, 6 — символ дельтафункции.

Наблюдаемый на входе приёмного устройства сигнал в общем виде с учётом (1), (12) и (16) представляется следующим образом:

ЛГ(/) = //5(0 + 2(/) =

Г И

/=1

+ =

=й(0

X *,«*<(')

(17)

С учётом (14) выражение (17) переписывается следующим образом:

ДГ(;)=Я5(/)+г(/)=

X £ А (')5(, О -г)<?( г - г* + г, (/)

1-1 А-О 1*1 1-0

(18)

Вектор наблюдения Е(/) для абонентского терминала М1МО-системы имеет вид:

Е(0 = ...£л«(')]\ i = \,M, (19)

где - наблюдение от /-го антенного элемента;

М - общее число антенн.

Математическая модель информационных

и сопутствующих процессов непрерывного

вектора состояния

Из рассмотрения математической модели полезных OFDM-сигналов, наблюдаемых на входе многоантенного приёмника (18), следует, что, наряду с дискретными параметрами 6, (/) и pu{t), они зависят от непрерывных случайных параметров: случайного ’запаздывания сигнала Г(/) и случайной фазы сигнала <p(t) ■ Непрерывные параметры содержат информацию об условиях распространения радиоволн и стабильности опорной частоты передатчика и используются для обеспечения синхронизации устройства приёма и обработки сигналов.

Случайная фаза сигнала обуславливается естественной нестабильностью опорной частоты передатчика, фазовыми флуктуациями, вносимыми средой распространения радиоволн, и флуктуациями движения абонентского терминала [1]. Эти же факторы могут влиять и на амплитуду принимаемого сигнала. При этом, как правило, мгновенные значения амплитуды сигнала, вследствие быстрых замираний, распределены по закону Релея (частота замираний 0.3...0.5 Гц) с параметром С 4, а начальная фаза имеет равномерное распределение на отрезке [- Л,л].

Случайные процессы <p(t) и A(Oki могут быть описаны системой стохастических дифференциальных уравнений [5]:

т=дйх/)--«/(/)W1 /2(г«(о, <ао=(р ....

С (20)

AoKt)=-Y^aW+yjlYjjri'Jt), Aa{tu) = \co где nv{t) и ИД(Ы(/) - АБГШ с нулевым математическим ожиданием и единичной интенсивностью; A(0(t)

- уход частоты сигнала за счёт нестабильности опорной частоты передатчика и приёмника; уы - коэффициент, характеризующий ширину спектра шума опорного генератора; —d(t) - доплеровское приращение частоты. с

Типовые значения параметров составляют [8]:

(?1Ц = (5...125) Гц; ущ = (КГ4... КГ2) Л!/;

= (2 • 10Л..2 • КГ4).

Применительно к решаемой задаче синтеза алгоритма приёма и обработки сигналов OFDM-MIMO непрерывный вектор состояния может быть представлен в виде:

*Т(/) = [ФГ(0,Д£2Г(/)Ь (21>

где Ф(/) = [(px(t),(p2(t).<Pg(t)Y - вектор случайных

фаз принимаемых OFDM-сигналов;

AQ(0 = [A<q(/)» Ac&(f),...,Аеок(/)У - вектор уходов частоты принимаемых OFDM-сигналов.

Динамика компонент непрерывного вектора состояния X(t) описывается системой дифференциальных уравнений [5]:

Фк«) = Acoy(t)-^dK(t)+yj\/2^ /Ч(/), ф,(Г0) = фхЛ (0 = -У,*Щ(')+рУА"<*('), Асоу(/„) = Аа)кМ

Щ)=х0,

где Пф (/) и П(0 (/) - независимые формирующие белого гауссовского шума с нулевым математическим ожиданием и известными интенсивностями; Дй g(t) -уход частоты сигнала за счёт нестабильности частоты опорного генератора передатчика; у - коэффициент,

характеризующий ширину спектра шума опорного генератора; ——d (t) - доплеровское приращение частоты;

С *

Х0 - начальные значения компонент вектора состояния

X(t).

Для оценки вектора состояния динамической системы целесообразно использовать рекурсивный фильтр Калмана, имеющий способ гарантированно достигать результат за конечное число действий, используя ряд неполных и зашумлённых измерений [9]. Фильтр Калмана предназначен для рекурсивного дооценивания вектора состояния априорно известной динамической системы, то есть для расчёта текущего состояния системы необходимо знать текущее измерение, а также предыдущее состояние самого фильтра. Таким образом, фильтр Калмана, как и множество других рекурсивных фильтров, реализован во временном представлении, а не в частотном.

При использовании фильтра Калмана для получения оценок вектора состояния процесса по серии зашумлённых измерений необходимо представить модель данного процесса в соответствии со структурой фильтра, то есть непрерывный вектор состояния X(t) в общем виде может быть описан векторно-матричным стохастическим дифференциальным уравнением [10]:

X0)=F(t)X«) + C(t)U(t) + G(l)N(t), X(i0)= Х0, (22) где N(t) — вектор формирующих белого гауссовского шума с известными статистическими характеристиками; F(t) - матрица состояния системы; C(t) — матрица управления, которая прикладывается к вектору управляющих воздействий U(t); G(t) - матрица возмущения;

U(t) = [X,,Wt„,Wt.......Wx,,Wy,,Wd] - вектор управ-

ляющих воздействий (вектор управления).

Применительно к решаемой задаче синтеза полный вектор состояния является дискретно-непрерывным и в общем случае имеет вид [11 ][ 12]:

[0l(t),eL(t),XT(t)h (23)

где 9{t) = [#! ^ , (0]7 - вектор дискретных ин-

формационных параметров; /?(/) = [Д ДД/)]7 -вектор дискретных энергетических параметров.

Рассмотренные математические модели принимаемых ОРОМ-сигналов хорошо отражают реальные условия функционирования систем БШД, использующих технологию М1МО. Они являются в требуемой степени общими и достаточно детально учитывают различные факторы, влияющие на работу системы: фазовую манипуляцию несущих колебаний; фазовые флуктуации ОР1)М-сигналов; наличие дискретных информационных параметров; наличие внутрисистемных и широкополосных помех.

Параметрическая математическая модель

информационных и сопутствующих процессов

Реализация оптимальных алгоритмов приёма и обработки ОРОМ-сигналов М1МО-систем связана с решением стохастического дифференциального уравнения Стратоновича, что сопровождается значительными трудностями. Поэтому при синтезе оптимальных алгоритмов стремятся упростить математическую модель вектора состояния. Для малых отрезков времени это удаётся сделать путём замены вектора состояния (21) и (22) на его параметрическую модель [11].

Непрерывный вектор состояния А'(0 (21) является многокомпонентным марковским диффузионным процессом. При переходе к параметрической модели действительную кривую, определяемую уравнением (22), заменяют ломаной (в простейшем случае ступенчатой) линией с шагом, равным сравнительно малому отрезку времени Т = - !к [13]. Для этого на каждом полу-

интервале времени [ЛЛ +7’)’ где к = 0,1,2,..., в уравнении (22) пренебрегают вектором формирующих белого гауссовского шума ЛД/) как малым (что соответствует равенству нулю коэффициента диффузии процесса Х(0 па малых отрезках времени длительностью Т).

Отсюда следует, что процесс X(() становится марковским квазислучайным и на полуинтервалах времени

- Т ), где к = 0,1,2,..., может быть описан следующим выражением [13]:

Х(1) = РХ(0 + Си(1), Х«к) = Хк,1е [/*,/* - Г), к = 0,1,2,... (24)

Для параметрической модели вектора состояния Х(1) решение уравнения (24) имеет вид [ 13][ 14]:

Г

ХЦ) = Ф(1-1к)Хк+\ф(1-т)Си(т)с1т< /е [/*,/* - Г), '1

Аг=0,1,2 (25)

где Ф(/ — ) - фундаментальная матрица, которая в

стационарном случае является решением дифференциального уравнения вида [14]:

</Ф(г~,к)_ { ), ф(0) = /, (26)

dt

где / - единичная матрица.

При дальнейшей аппроксимации математической модели вектора состояния Х(() за счёт упрощения коэффициента сноса РХ(/) фундаментальная матрица Ф(/,/*) на сравнительно малых отрезках времени мо-

жет быть заменена матрицей A(t,tk), которая может быть получена путём поэлементного разложения матрицы Ф(!,1к) в ряды Тейлора и удерживания только линейных членов [15]. При этом в виде неступенчатых квазислучайных процессов сохраняются лишь те компоненты, для которых недопустимо с требуемой точностью пренебречь детерминированным изменением (сносом) на полуинтервалах \tk,tk - Т). Остальные компоненты Х(1) могут быть аппроксимированы ступенчатыми квазислучайными процессами. Таким образом, параметрическая модель вектора состояния X(t) в конечном итоге имеет вид [ 13][ 14][ 15]: t

X(t) = A(t,tk)Xk + J A(t,tk )C(T)U (T)dTi (27)

h

tz[tk,tk- T), к = 0,1,2... X(t0) = X0,

A(tK,tk) = I, Xk=X(tk).

Сформировав параметрическую модель вектора состояния Х(/), необходимо проверить обоснованность таких допущений для M1MO-OFDM систем, т.е. оценить, насколько повлияет пренебрежение вышеуказанными параметрами на точность оценивания синтезируемых компонент вектора состояния X(I).

Определённую оценку можно получить, рассчитав приращения компонент вектора математического ожидания и элементов матрицы вторых центральных моментов исходного процесса X(t) (21) на характерных

отрезках времени Г = tk+l - tk на основе следующих соотношений [13][14][15]:

Мх(к + 1) = Ф(Г)Л/Г(£), Мх(0) = МхО; (28)

D(k + \) = Ф(Т)0(к)Ф' \Т)+ Г(Т)Г\Т), D(0) = 0, (29)

где Мх(к) = М{Хк}. Д*)=Л/|[^-А/д(*)][^-Л/л.(*)]т};

А/ДО) и D(0) - вектор математического ожидания и

матрица вторых центральных моментов вектора состояния в начальный момент времени; Ф(Т) - фундаментальная матрица, Г(Т) - переходная матрица возмущений; Г(Т)Гт(Т)= D,(T) - матрица приращений вторых

центральных моментов, которую можно определить следующим образом:

^^ = га|(0 + Д(0/7Т + 5. D,(0) = 0, f е [0,7"), (30)

dt

где В = GQC1, G - матрица возмущений, Q — матрица интенсивностей ЛД/).

Полученные значения приращений компонент вектора математического ожидания Мх(к +1) и элементов

матрицы дисперсий D(£ + l) необходимо сравнить с типовыми значениями соответствующих погрешностей оценивания компонент вектора состояния X(t)-

Расчёт приращений математического ожидания и дисперсий для компонент вектора состояния Х(1) производился при следующих значениях параметров [16]: Т= 20 мс, fuec = 2,4 ГГц (на данной частоте функционируют системы Wi-Fi, WiMAX и LTE); д«=1,8 10'5.

На рисунке 2 показано приращение математического ожидания по фазе ДМ[(р\ радиосигнала. Приращение

математического ожидания этой компоненты достаточно велико, поэтому пренебречь детерминированным изменением компоненты ц. нельзя, и, следовательно, соответствующие элементы матрицы А(Т) должны быть сохранены.

Рис. 2. Приращение математического ожидания по фазе

Расчёт приращения априорной дисперсии по фазе вектора состояний Х(1) производился на ПЭВМ в программном пакете МАТЬАВ при помощи численного решения дифференциального уравнения (29) с использованием метода Рунге-Кутты четвёртого порядка [10]. График, представленный на рисунке 3, отражает характер изменений априорной дисперсии по фазе. Анализируя рисунок 3, можно увидеть, что приращение дисперсии на отрезке Т достаточно мало (по сравнению с погрешностями оценивания) и поэтому использование параметрической модели вида (24) и (27) вполне допустимо для решения задачи синтеза применительно к ОРЭМ-М1МО системам.

Ьо'[<А, рад1

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

Сформулированные модели наблюдений совместно с предложенными моделями информационных и сопутствующих процессов позволят корректно осуществить постановку задачи синтеза алгоритмов комплексной первичной обработки информации в абонентском терминале систем БШД с использованием технологии М1МО.

Математические модели информационных и сопутствующих процессов дискретной компоненты вектора состояния

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

собой вектор [pTkj (/), 0ki (/)]. Одна компонента 6(t)

такого вектора дискретных параметров является информационным параметром, предназначенным для передачи речи, служебной информации и т.п., а другая /?(/) -энергетическим, характеризующим, например, амплитуду принимаемого радиосигнала на входе абонентского терминала [17].

Основным информационным параметром, предназначенным для передачи речевых сообщений, данных, служебной информации и т.п. в системах OFDM-MIMO и содержащимся в полезных сигналах, является дискретный параметр 6(1).

Обычно изменения 6(t) во времени удаётся описать однородной цепью Маркова на два (в более общем случае на R) положения с известными статистическими характеристиками [17]:

= ('*)]' Iй-v~ 1, А/,-матрица вероятностей

перехода в(1),

где ) = p'Mtk) = ву !6(tk_x) = } - вероятность

перехода информационного параметра на к-ом шаге;

Ра =[^(/0). ^('о). •••• РМ)' ~ вектор вероятностей начального состояния 6(1).

Дискретный случайный процесс /?(/), характеризующий амплитуду принимаемого полезного сигнала, описывается простой цепью Маркова на L положений и может принимать одно из значений A., j = l,L. Возможные моменты перехода параметра J3(t) из одного состояния в другое являются дискретными, хотя в общем случае случайными, так как зависят от случайного запаздывания принимаемого сигнала [ 11 ][ 15]:

Тс- Г,(/), к = 0,1,2,... (1.31)

где т (t) = г(t,X(t)) - случайное запаздывание сигнала; Т = const — длительность тактового интервала.

Значения дискретного случайного процесса /?(/) на соседних тактовых интервалах образуют простую цепь Маркова с матрицами одношаговых вероятностей перехода [15]| 18]:

W = [/>(/*)]. /,у = и;

и вектором вероятностей начального состояния:

/>,.=[vaw--v'„>i ’ у=1>^•

На всех полуинтервалах времени [tll,tk+]), где к = 0,1,2 ..., дискретные параметры /3(t) и 6(t) остаются постоянными и описываются априорными уравнениями вида [ 13][ 15]:

£^ = о; = о, (32)

dt dt

где /€ [/*,/*+,), * = 0,1,2...

Применение представленных математических моделей полезных сигналов, помех, информационных и сопутствующих параметров OFDM-MIMO систем позволяет достаточно корректно сформулировать постановку задачи синтеза оптимальных алгоритмов оценивания. Располагая априорными сведениями о дискретных параметрах /?(/) и 6(/), сведениями о непрерывном векторе состояния X(t) и наблюдениями Е(1), требуется определить оценку непрерывного вектра состояния Х’(1), оптимальную по критерию минимума апостериорных дисперсий, и оценки информационного и энергетического параметров в'(!) и /?*(/), оптимальные по

критерию максимума апостериорной вероятности.

Литература

1. Левин Б.Р. Теоретические основы статистической радиотехники. М.: Радио и связь, 1989. - 656 с.

2. Лебедев В. Модуляция OFDM в радиосвязи. Радиолюбитель. Август, 2008. С.51-55; Сентябрь, 2008. С.36-39)

3. Гаранин М.В., Журавлев В.И., Кунегин С.В. Системы и сети передачи информации. Учебное пособие для вузов. — М.: Радио и связь, 2001.-336 с.

4. Корн Г. и Корн Т. Справочник по математике для научных работников и инженеров. - М., 1973. - 832 с.)

5. Варакин Л.Е. Теория систем сигналов. - М.: Советское радио, 1978. - 304 с.

6. Иванов А. Алгоритмы синхронизации в OFDM-системах. - 6 с.

7. Сетевые проекты “Применение современных методов модуляции, кодирования и разделения в WiMAX” -http://www.networkspro.ru/pages_64/index.htm I.

8. Зверев В.А., Стромков А.А. Выделение сигналов из

помех численными методами. - Нижний Новгород: ИПФ РАН» 2001. 188 с.

9. Kalman R. Е. A New Approach to Linear Filtering and Prediction Problems. - I960. -12 c.

10. Фильтр Калмана для “чайников”: матрица наблюдаемости и фундаментальная матрица. - 4 с.

И. Тяжев А.И. Оптимизация цифровых детекторов в приёмниках по минимуму вычислительных затрат. — Самара, 1994.-256 с.

12. Некинелов Н. Калмановская фильтрация. — http://www.basegroup.ru/library/cleaning/kalmanfilter.

13. Павлов Б.А. Синхронный приём. - М.: Энергия, 1977.-80 с.

14. Тихонов В.И., Харисов В.Н. Статистический анализ и синтез радиотехнических устройств и систем. - М.: Радио и связь, 2004. - 608 с.

15. Трифонов А.П., Шинаков Ю.С. Совместное различение сигналов и оценка их параметров на фоне помех. — М.: Радио и связь, 1986. - 264 с.

16. Волков Л.Н., Немнровский М.С., Шинаков Ю.С. Системы цифровой радиосвязи: базовые методы и характеристики. - М.: Эко-Трендз, 2005. - 392 с.

17. Френке Л. Теория сигналов. Перевод с английского под редакцией Д.Е. Вакмана. - М.: Советское радио, 1974. - 344 с.

18. Липкнн И.А. Статистическая радиотехника. Теория информации и кодирования.-М.: Вузовская книга,2002.-216с.

MATHEMATICAL MODELS OF SIGNAL OF BWA SYSTEMS BASED ON MIMO Bukharev I.A., FSUE ZNIIS, lead engineer, bukharev@zniis.ru

Annotation: MIMO technology (Multiple Input Multiple Output) has found wide practical applications in wireless broadband network of standard IEEE 802.11n (Wi-Fi), IEEE 802.16 (WiMAX), LTE and the recently released standard 802.22 (White Space). The main advantage of technology - increasing throughput, capacity and reliability of the network in proportion to the number of transmitting and receiving antennas. The technology allows you to partially or completely solve problems inherent in all wireless networks, such as interference with propagation loss, slow and fast fading, shadowing, shielding. The technology has a great future, but should be constantly improved, so the question of its study, including modeling of transmitted and received signals, must be given due attention. In modern wireless systems MIMO technology is used in conjunction with the technology of orthogonal coding OFDM.

Keywords: mathematical model, broadband wireless access, Wi-Fi, WiMAX, LTE, MIMO, OFDM.