CC BY
41
Доклады БГУИР
2015 № 3 (89)
УДК 621.395
РАСЧЕТ ЗАЩИЩЕННОСТИ СИГНАЛА ТИПА OFDM ОТ ПОМЕХ НЕЛИНЕЙНОГО ПРОИСХОЖДЕНИЯ
В.И. КИРИЛЛОВ*, А.А. ПИЛЮШКО, Е.К. КАРПУК
*Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Военная академия Республики Беларусь Независимости, 220, Минск, 220057, Республика Беларусь
Поступила в редакцию 16 декабря 2014
Рассмотрена методика расчета защищенности сигнала типа OFDM от помех нелинейного происхождения, основанная на использовании универсального метода определения продуктов нелинейности на выходе трактов прохождения электрических сигналов со сложными (негладкими) мгновенными динамическими характеристиками и полигармоническом и/или модулированном входном воздействии.
Ключевые слова: OFDM-сигнал, защищенность от помех нелинейного происхождения, функциональный преобразователь электрических сигналов, метод определения продуктов нелинейности, мгновенная динамическая характеристика.
Введение
В настоящее время в системах телекоммуникаций широко используются цифровые методы обработки и передачи информации (например, OFDM, DMT и др.), что позволяет существенно улучшить такие технико-экономические показатели систем передачи, как пропускная способность и дальность связи [1-3]. При этом стремление к снижению стоимости цифровых магистралей заставляет зачастую между обслуживаемыми регенераторными пунктами (ОРП) включать подряд несколько необслуживаемых пунктов, работающих в режиме аналоговых усилительных пунктов. Последние, как впрочем и ОРП, могут являться источниками специфических нелинейных помех, которые приводят к возрастанию вероятности ошибок цифрового сигнала. Поэтому актуальность борьбы с нелинейными искажениями и помехами нелинейного происхождения (продуктами нелинейности) в системах телекоммуникаций до сих пор объективно сохраняется [3-7]. Анализ научно-технической литературы показывает, что существующие способы формирования многочастотного OFDM-сигнала, например в «Режиме 2К» [1] при известном выборе сетки несущих частот: fh = f0 + h / т, где f0 - нижняя частота выбранного диапазона; h - номер несущей (h = 0..K — 1; K - количество несущих); т - временной интервал передачи одного многоуровневого символа OFDM на каждой несущей; не позволяют избавиться от продуктов нелинейности [2-3].
Данное обстоятельство изначально (уже на этапе формирования сигнала) снижает потенциальные возможности OFDM-сигнала по его ретрансляции и передаче на максимальное расстояние. Для того, чтобы повысить защищенность OFDM-сигнала от продуктов нелинейности, возникающих в результате прохождения такого сигнала через нелинейный функциональный преобразователь (НФП) или группу НФП, входящих в состав тракта прохождения электрических сигналов (ТПС), прежде всего необходимо ответить на следующие вопросы: какое количество продуктов нелинейности присутствует в спектре сигнала, каковы особенности их накопления, амплитуда и частота, какова будет суммарная вероятность ошибки передаваемого сигнала в зависимости от вероятности ошибки на отдельной несущей.
Нахождение спектра OFDM-сигнала на выходе НФП
В общем случае OFDM-сигнал представляет собой сумму нескольких частотных составляющих (их количество может достигать нескольких тысяч), каждая из которых квадратурно модулирована (QAM) [1]. Исследование такого сложного сигнала весьма затруднительно, поэтому целесообразно вначале остановиться на рассмотрении сравнительно простого сигнала с двумя амплитудно-модулированными несущими вида UK (О = Uml [l + щ sin(Q1t)]sin(ffll0 + Um2 [l + m2 sin(Q2t)]sin(©2i), (1)
где Umi, Um2 - амплитуды несущих сигнала типа OFDM; щ, щ - коэффициенты амплитудной
модуляции; Qx, Q2 - частоты модулирующих сигналов; , ю2 - частоты несущих.
Для нахождения спектра сигнала (1) на выходе НФП целесообразно воспользоваться универсальным методом определения продуктов нелинейности, рассмотренным в [5-7]. Сущность метода заключается в аппроксимации, как правило, нелинейной мгновенной динамической характеристики (МДХ) НФП/Т ПС рядом Фурье с последующим разложением каждой гармоники этого ряда во временной ряд с помощью функций Бесселя. Использование универсального метода обусловлено тем, что он выгодно отличается от традиционных методов своей простотой и высокой точностью определения продуктов нелинейности (прежде всего комбинационных продуктов) в тех случаях, когда НФП обладает сложной («негладкой») МДХ и при полигармоническом и/или модулированном входном воздействии.
Пусть нелинейная МДХ НФП UBbIX г = ф(ивх г) представлена в виде ряда Фурье:
N
ивых = A + X A cos(k^UвХ) + Bk sin(кШвХ)], (2)
k=1
где N - количество точек (дискретных пар ив^х, =ф(ивх i) , где ивх i е[ивх m,nUвх max]),
полученных в результате эксперимента по снятию МДХ НФП; UBX ^,^^^ - минимальное и
максимальное значения входного сигнала;
1 n-1 2 Г 1
A0 = ^ 2 Uвых< ; Ak = Тг 2 Uвых¿ COS(k^UBKi )J;
0 AT 2 вых ' ' k N
N i=0 N i=0
2л 2л „2
*. = ^ = ^-^-; Bk = 2 2 [и вых i sin( kA, и^вх i)].
T U - U 1 k N
вх max вх min i=0
Если подставить исследуемый сигнал (1) в выражение (2) и учесть известные из теории функций Бесселя соотношения
ад ад
cos(a • sin x) = J0 (a) + 22 [J2n (a) cos(2nx)], sin(a • sin x) = 22 [J2n-1 (a) sin( [2n -1]x)],
n=1 n=1
где J (a) - табличная функция Бесселя порядка n от аргумента a; то можно получить математическую модель спектра этого сигнала на выходе НФП. Она позволяет рассчитать любые (в данном классе) нелинейные продукты на выходе анализируемого НФП. Например, амплитуды гармоник i-го (Um гГ) и j-го (U г) порядка - это коэффициенты при множителях
cos('raí) и sin(jrat) соответственно.
Из всех составляющих спектра сигнала (1) на выходе НФП наибольший интерес представляют комбинационные продукты вида +П2 и 2ю1±ю2, поскольку они в наибольшей мере воздействуют на полезные гармоники [3].
В общем случае амплитуда комбинационного продукта на частоте ±П2 (называемого часто перекрестным искажением) равна [5-7]:
Um k 21 = 4ÍtBk [ J 02( X1 >] - 4ÍtBk ^(Х1 m^1 (3)
k=1 2 k=1 2
m m
где с учетом обозначений (1) X1 = kkUm 1 ; S1 = J0(X2 J1(X1)J1(X2)J1(X2
Количество продуктов перекрестных искажений, воздействующих на отдельную несущую, равно М = К — 1, где К - количество несущих [5-6].
Амплитуда комбинационного продукта на частоте 2юх +ю2 (называемого часто интермодуляционным искажением) равна [5-6]:
N N
Umk3 = Bk[Jo2(Xi mi)Jo2(X2 m^)Jl(X2)J2(Xl)] + Bk[J02(X2 rm^-)Jl(X2)J2(Xl) x k=1 22 k=1 2
22(XI m-)]—2]TBk [ Jo2 (X 2) Jo2 (X 2 ^ Ji(X 2) Ji2( X1 m-)]+2]TBk [Jo2(Xi m 2 k=1 2 2 k=1 2
N - m
x J12 (X2 J2 (X1)] + 2^ Bk [ J (X2 )J1 (X2 m2)J2 (X1)J22 (X v)]. (4)
2 k=1 2 2 В общем случае количество продуктов интермодуляционных искажений, попадающих на каждую конкретную полезную гармонику, зависит от количества несущих в сигнале типа OFDM и конкретных значений их частот.
Влияние рассмотренных помех на качество передачи полезного сигнала зависит от соотношения их амплитуд. В частности, амплитуда первой гармоники полезного сигнала, например, на частоте roj, равна [5-6]:
N ™ „„ N
2,v mx.T 2^ m2
Um 1Г»1 = 2^Bk[Jo(X2)Jo2(X1 -^)Jo2(X2 m2)Jl(Xl)] — 2^Bk[J0(X2) x k=1 2 2 k=1
N
x Jo2(X -^(X —t)J 1 (X1)] + 2^ Bk [J o2( x 1 -^(хшхуЛх m-)] —
N
2^Bk[J1(X1)J12(X1 m)J12(X2 J2(X2)]. (5)
к=1 2 2
Коэффициенты гармоник К1Г и Кг, а также затухания нелинейности по гармоникам агГ и ауГ можно рассчитать по известным формулам [3]:
К,Г = ит адг / Ц, 1Г и аиТ =—201в К,,г = 201в и„ 1Г / Ц, ^, (6)
где ит1Г - амплитуда первой (основной) гармоники на выходе НФП/Т ПС.
Авторами был произведен расчет затухания нелинейности по перекрестным а и интермодуляционным искажениям а в зависимости от мгновенного значения входного сигнала х = ивх / ивх т, ивх т = 1В при различных исходных данных (соотношениях коэффициентов модуляции т1 /т2 и амплитуд составляющих входного сигнала ит1/ит2). Здесь ивхт - максимальная амплитуда переменной составляющей сигнала на входе скорректированного тракта. При этом зависимость нормированной МДХ для тракта в данном случае имеет вид у = х + Лу(х), где Ау = ах2 + Ъх3, а = 0,007, Ъ = 0,01, что соответствует затуханию нелинейности по 2-й и 3-й гармоникам равным 55 дБ и 73 дБ при выходном уровне испытательного сигнала на выходе тракта, равном нулю [5].
По результатам расчета сделаны следующие выводы:
1. Мощность интермодуляционных и перекрестных искажений существенно зависит от указанных соотношений: так, при итХ / ит2 е [0,1;0,9] и щ / т2 = 0,9 или щ / т2 = 0,1 мощность интермодуляционных искажений больше на 12-15 дБ, чем мощность перекрестных искажений. Если же щ /щ = 0,5, то при этих же условиях мощность перекрестных искажений больше на 2-3 дБ, чем мощность интермодуляционных искажений.
Защищенность от интермодуляционных и перекрестных искажений (А3) в любом из рассмотренных случаев составляет не менее 60 дБ, что соответствует требованиям по защищенности от помех нелинейного происхождения, предъявляемым, например, в цифровом
телевидении (А3 феб > 50 дБ) [8-9]. Однако следует учитывать, что при рассмотрении реального OFDM-сигнала с большим количеством несущих на каждую его полезную гармонику может воздействовать значительное количество таких помех, в результате чего А может стать меньше 50 дБ. Кроме того, после прохождения OFDM-сигнала через ряд последовательно соединенных НФП в его спектре будет появляться все больше и больше новых паразитных составляющих, как результат комбинирования не только полезных гармоник OFDM-сигнала, но и добавления помех нелинейного происхождения, образованных от предыдущих НФП. Очевидно, именно данный факт может быть основной проблемой при передаче OFDM-сигнала с ретрансляцией по высоким частотам.
Следует отметить, что порядок нахождения спектра сигнала (1) на выходе НФП и представленные выше по тексту формулы для расчета продуктов нелинейности применимы при исследовании НФП со сложной («негладкой») МДХ и при более сложных входных воздействиях. Они позволяют определить любые (в данном классе) продукты нелинейности, дать исчерпывающий ответ об их количестве и мощностях. Однако рассмотренный порядок нахождения спектра сигнала (1) на выходе НФП, если применить его к OFDM-сигналам с произвольным количеством и номиналами частот несущих, не позволяет получить выражение для расчета защищенности А3 OFDM-сигнала в общем виде. Это обстоятельство связано с тем, что в формулу для расчета А3 OFDM-сигнала необходимо подставлять значения мощности ряда комбинационных продуктов, которые попадают в конкретную полосу пропускания одной из составляющих OFDM-сигнала. При изменении номиналов частот несущих количество таких комбинационных продуктов тоже изменится. Например, пусть входной сигнал состоит из трех несущих с частотами f, f2, f3, равными соответственно 1, 2 и 3 МГц в первом случае и равными 1, 3 и 5 МГц во втором случае. Тогда, например, комбинационный продукт вида f + f2 попадет в полосу пропускания несущей f3 только в первом случае, поскольку f + f2 = f3 = 3 МГц. Во втором случае комбинационный продукт вида f + f2 не попадет в полосу пропускания несущей f, поскольку f + f2 = 4 МГц Ф f = 5 МГц.
В связи с этим целесообразно использовать общий подход, который применяется, в частности, для расчета нелинейных продуктов в аналоговых системах с частотным разделением каналов (ЧРК) [3]. Но поскольку технология OFDM подразумевает передачу цифрового сигнала, то необходимо предварительно найти связь между мощностью нелинейных продуктов и обусловленных ими ошибками в цифровом сигнале.
Расчет вероятности ошибки сигнала типа OFDM
Как известно [3], суммарная вероятность ошибки сигнала OFDM находится, как отношение общего количества ошибок на отдельных i-х несущих Nonn к скорости передачи исходного двоичного сигнала B^ :
K-1
Лкж.х=ЕNошi / B2, (7)
i=0
где K - количество несущих.
Тогда количество ошибок, возникающих на отдельной i-й несущей, можно найти как
Кшi = Рош.Д-, (8)
где p , B - вероятность ошибки и исходная битовая скорость передачи информации соответственно на отдельной i-й несущей (в первом приближении Bt = B^ / K).
Скорость передачи информации Bt известна, поскольку параметры OFDM-сигнала формируются программно непосредственно перед осуществлением его передачи по конкретной линии связи.
Вероятность ошибки рошг можно найти с использованием формул, применяемых для расчета защищенности А3 для /-го многоуровневого цифрового сигнала, передаваемого на i-й несущей с помощью mi -уровневого линейного кода [3]:
А = 10,65 + 11,42lglg(1 /Рош./) + 20lg((m, —1)/2). (9)
Отсюда
mZ m — 1
Рош./ = 10—10 , где Zt = А —10,65 — 20lg^)/11,42. (10)
Известно, что расчет суммарной защищенности полезного сигнала от основных источников помех сигнала осуществляется по следующей формуле [3]:
АзЕ = —10lg(Рпом.Е /Pc) = —10lg((Рпом.с + Рпом.лп + Рпом.нп.)/Pc) , (11)
где Pc - мощность полезного сигнала типа OFDM, прошедшего через НФП; Рпом.2 - суммарная мощность помехи, состоящая из суммы мощностей собственных (тепловых) шумов Рпом с, помех линейных переходов Рпомлп. и помех нелинейного происхождения Рпом.нп (прежде всего интермодуляционных и перекрестных искажений) соответственно.
Соотношение этих мощностей в общем случае может быть любым. Для определения конкретных значений необходимо проведение многочисленных измерений в течение длительного времени и статистическая обработка результатов. В настоящее время в современной научно-технической литературе нет данных о допустимом (или типовом) соотношении мощностей помех при передаче OFDM-сигналов. Здесь для простоты можно принять соотношение (по аналогии, как это принято для многопарных симметричных кабелей) Рпомс : Рпомлп.: РПом.нп. = 1:1:1. Учитывая данное обстоятельство, формулу (11) можно представить в следующем виде:
Аз^нп. = — 10М3РпоМШ1. /Pc) = — 10М3(РпоМНИ + Рпомин. + Рпомпер.)/Pc) , (12)
где PmMm , PII0MИН - мощности продуктов нелинейности 2-го и 3-го порядка (нелинейных и
интермодуляционных искажений соответственно); P0№nep - мощность перекрестных искажений.
Мощность перекрестных искажений находится с помощью универсального метода определения продуктов нелинейности [5-7]. Однако, как уже говорилось выше по тексту, применение этого метода для нахождения количества и мощностей продуктов нелинейности 2-го и 3-го порядков на i-ой несущей при использовании OFDM-сигналов с большим количеством несущих весьма затруднительно. Данное обстоятельство связано с тем, что выбор частот несущих зависит от полосы пропускания предоставленного канала связи, а также от условий распространения каждой i-ой несущей в этой полосе. В любом случае перед началом работы происходит анализ предоставленного канала на соответствие его предъявляемым требованиям, на основании чего, как правило, программно осуществляется окончательный выбор количества несущих OFDM-сигнала и закон ее модуляции (по сути, скорость передачи сигнала на i-й несущей). При изменении параметров канала передачи расчетные выражения, полученные для OFDM-сигнала первоначально, изменяются, поскольку изменится количество и номиналы частот несущих.
Для расчета мощности продуктов нелинейности 2-го и 3-го порядков на i-ой несущей целесообразно воспользоваться формулами для нахождения спектральной плотности мощности многоканальных сигналов, известными из теории многоканальных систем передачи с частотным разделением каналов [3]:
Pj.2(8л) = 4PZ210—01а- Д^^*), (13)
Д F
Pn3 (8*) = — 24P| 10—01^° Y3(5k), (14)
где Р^ - суммарная мощность основных гармоник входного сигнала на выходе НФП; агтс а3Г0 - затухания нелинейности по 2-й и 3-й гармоникам при выходном уровне
испытательного сигнала на выходе НФП, равном нулю (или при амплитуде выходного напряжения, равной эталонному напряжению, развивающему на эталонном сопротивлении эталонную мощность 1 мВт); А^ - полоса пропускания г-ой несущей; А/ = /в — /н - полоса пропускания ОРБЫ-сигнала с граничными частотами /в и /н ; 5к = (/ — /н)/ А/ - средняя нормированная частота к-ой несущей в спектре частот ОРБЫ-сигнала (/н < /к < /в). Очевидно, 5к е [0;1,0].
Функция У^ (5) определяется одним из приведенных ниже выражений в зависимости от величины нормированной частоты 5 = 5к и коэффициента перекрытия частотного диапазона
Р = /в / /н [3]:
Р-2 —5,0 <5<Р — 2
Y2(5) Н
р-1 ' р-1'
( 1 >
8 —
Р-1/ Р-1
1 ^8< 1, (15)
2Р-5 5 1
2(р-1) 2' р-1 р-1
Функция 73(8) определяется выражением [3]:
3
73(8) = -[1 + 28(1 -8)] ,0<8< 1. (16)
8
Следует отметить, что формулы (13)-(16) применимы в случаях, когда исследуемый НФП обладает МДХ со «слабой» нелинейностью (МДХ можно аппроксимировать полиномом не старше 3-й степени). Если исследуемый НФП обладает сложной («негладкой») МДХ, то вначале ее необходимо линеаризировать, т.е. привести к виду со «слабой» нелинейностью. Для этого целесообразно воспользоваться методикой, приведенной в [7], после чего дальнейший расчет вероятности ошибки сигнала OFDM проводится по методике, рассмотренной выше по тексту.
Таким образом, для расчета вероятности ошибки OFDM-сигнала необходимо:
1) знать исходные параметры OFDM-сигнала (амплитуды, частоты, закон модуляции и количество несущих);
2) с помощью универсального метода определения продуктов нелинейности по формулам (3)-(6) определить амплитуды несущих OFDM-сигнала, прошедшего через НФП, и мощности перекрестных искажений;
3) найти мощности продуктов нелинейности 2-го, 3-го порядка, оказывающих влияние на каждую несущую, подставив поочередно в (13), (14) вместо AF конкретные значения полосы частот, занимаемой данной несущей;
4) на основании известных мощностей продуктов нелинейности найти поочередно соответствующие защищенности несущих от этих помех, вероятности ошибок и количество ошибок на отдельных несущих по формулам (8)-(12);
5) используя формулу (7), найти суммарную вероятность ошибки OFDM-сигнала.
Заключение
Расчет защищенности сигнала типа OFDM от помех нелинейного происхождения представляет собой весьма сложную задачу, которая требует проведения большого количества математических вычислений. В первую очередь необходимо найти спектр исследуемого сигнала, прошедшего через НФП. Для решения этой задачи лучше всего подходит универсальный метод определения продуктов нелинейности, который позволяет с высокой точностью определять продукты нелинейности на выходе НФП/ТПС с МДХ сложной
(«негладкой») формы и полигармоническом и/или модулированном входном воздействии. Кроме того метод является сравнительно простым по вычислительной процедуре (опирается на известное, стандартное программное обеспечение по расчету коэффициентов Фурье и функций Бесселя) и не требует проведения многочисленных и трудоемких натурных экспериментов по определению продуктов нелинейности в случае изменения параметров входного сигнала (уровней, соотношения частот и мощностей составляющих входного сигнала при полигармоническом входном воздействии, параметров модуляции при модулированном входном воздействии).
Результаты, полученные с помощью универсального метода, могут быть использованы для осуществления расчетов суммарной вероятности ошибки сигналов типа OFDM в соответствии с предложенной здесь методикой.
В свою очередь, расчет суммарной вероятности ошибки сигналов OFDM позволит на этапе проектирования техники связи ответить на вопрос, какова будет помехозащищенность исследуемого НФП/ТПС при воздействии на него OFDM сигнала с определенными параметрами и какой она должна быть, чтобы обеспечить его ретрансляцию по высоким частотам.
Таким образом, комплексное применение универсального метода определения продуктов нелинейности и формул для расчета защищенности многочастотных сигналов типа OFDM (7)—(16) позволяет произвести анализ таких сигналов и НФП на предмет защищенности их от помех нелинейного происхождения. По результатам анализа становится возможным выдать конкретные практические рекомендации по структурно-параметрической оптимизации НФП и сигналов типа OFDM.
CALCULATION OF IMMUNITY TYPE OFDM SIGNAL FROM INTERFERENCE
NON-LINEAR ORIGIN
V.I. KIRILLOV, A.A. PILIUSHKO, E.K. KARPUK
Abstract
The method of calculation of the security OFDM signal from interference nonlinear origin is described in this article. Such method based on the use of a generic method definition products nonlinearity output paths passing electrical signals with complex (non-smooth) instant dynamic characteristics and polyharmonic and / or modulated input action.
Список литературы
1. ГОСТ Р 55694-2013. Телевидение вещательное цифровое. Наземное цифровое телевизионное вещание.
2. Шахнович И.В. Современные технологии беспроводной связи. М., 2006.
3. Кириллов В.И. Многоканальные системы передачи: учебник для вузов. М., 2003.
4. Фалько А.И., Омуралиева С.С. // Вест. СибГУТИ. 2011. № 1. С. 3-16.
5. Пилюшко А.А. Анализ и структурно-параметрический синтез волоконно-оптических модемов для ведомственных сетей связи: Дис. ... канд. тех. наук. Минск, 2007.
6. Карпук Е.К. Методика расчета помех нелинейного происхождения на выходе функциональных преобразователей сигналов с негладкими мгновенными динамическими характеристиками при сложных входных воздействиях: Дис. ... магистра тех. наук. Минск, 2012.
7. Кириллов В.И., Пилюшко А.А., Карпук Е.К. // Электросвязь. 2014. № 12. С. 5-10.
8. ГОСТ Р 52722-2007. Каналы передачи цифровых телевизионных сигналов аппаратно-студийного комплекса и передвижной телевизионной станции цифрового вещательного телевидения.
9. ГОСТ Р 52595-2006. Линии соединительные цифровые для передачи телевизионных программ.
