The method of calculating the error probability per bit in the communication channels with block freeze-ups for the receiver with linear addition of soft solutions of the incoherent demodulator

Full Text

Известно [1-3], что при наличии в каналах связи глубоких замираний (например, релеевского типа) для достижения вероятности ошибки на бит 5 10-=bP требуется обеспечить на входе оптимального некогерентного приемника сигналов с двукратной частотной манипуляцией отношение энергии сигнала Eb, приходящейся на бит информации, к спектральной плотности мощности шума N0 (отношение «сигнал/шум», далее - СШ) порядка 50 дБ. Это указывает на необходимость увеличения СШ на входе приемника на 37 дБ по сравнению со случаем отсутствия замираний (когда достаточно обеспечить СШ на уровне 13 дБ). Кодирование с исправлением ошибок обеспечивает существенное повышение помехоустойчивости в каналах связи с замираниями [3; 4]. Определяющим при применении кодирования в каналах связи с замираниями является время когерентности канала связи. Когда время когерентности превышает длительность передачи символа кода, замирания становятся время-селективными (блоковыми, пакетированными). Для борьбы с такими замираниями в сочетании с кодированием используется перемежение. Идея перемежения сводится к тому, что кодовые символы передаются в канал связи не в том порядке, в котором они поступают с выхода кодера, а после перестановки, после которой соседние символы одной кодовой комбинации разделяются символами других кодовых комбинаций [21]. Расстояние между такими символами зависит от глубины перемежения m. Как правило, глубина перемежения выбирается таким образом, чтобы исключить группирование ошибок в пакеты [11]. ТЕХНОЛОГИИ ТЕЛЕКОММУНИКАЦИЙ УДК 681.3 МЕТОД РАСЧЕТА ВЕРОЯТНОСТИ ОШИБКИ НА БИТ В КАНАЛАХ СВЯЗИ С БЛОЧНЫМИ ЗАМИРАНИЯМИ ДЛЯ ПРИЕМНИКА С ЛИНЕЙНЫМ СЛОЖЕНИЕМ МЯГКИХ РЕШЕНИЙ НЕКОГЕРЕНТНОГО ДЕМОДУЛЯТОРА Шевченко В.А.1, Пашинцев В.П.2 1 Департамент информационных систем Министерства обороны РФ, Москва, РФ 2 Северо-Кавказский федеральный университет, Ставрополь, РФ E-mail: shevv67@mail.ru Определены выражения для верхней аддитивной границы вероятности ошибки на бит в «некогерентных» каналах связи с кодированием и псевдослучайным перемежением при использовании метрики с линейным сложением мягких решений демодулятора в условиях воздействия блочных замираний, вызывающих группирование ошибок. Выражения конкретизированы для замираний, подчиненных распределениям Райса, Накагами-m и Накагами-q. С использованием полученных выражений на примере двоичных расширенных кодов Боуза-Чоудхури-Хоквингхема показано, что при заданном ограничении на задержку передачи информации коды длиной 64, 128 обеспечивают допустимую вероятность ошибки 10-5 при отношении «сигнал/шум» на входе приемника 16-18 дБ и энергетический выигрыш 3-4 дБ по сравнению с кодом с восьмикратным повторением. Ключевые слова: вероятность ошибки, пакет ошибок, замирания, некогерентный канал связи, метрика, кодирование Согласно приближенным расчетам, приведенным в [3], допустимая вероятность ошибки на бит информации 5 10-=bP в канале с релеевскими замираниями при использовании кодирования с перемежением обеспечивается при отношении СШ в пределах 13…21 дБ. Это указывает на возможность получения энергетических выигрышей 37…29 дБ по сравнению со случаем применения одиночного (неразнесенного) приема замирающих сигналов (когда требуется обеспечить СШ на уровне 50 дБ). На практике, когда имеются ограничения на задержку передачи информации tD, следует учитывать, что увеличение глубины перемежения с целью обеспечить независимость ошибок предполагает уменьшение длины кода и его корректирующей способности. В условиях априорной неопределенности относительно времени когерентности канала исключение группирования ошибок может привести к неоправданному укорочению длины кода. Возможным выходом из сложившейся ситуации является использование блочных (n, k) кодов большой длины и с высокой корректирующей способностью в сочетании с перемежением, допускающим группирование ошибок. Для подтверждения реализуемости указанного метода повышения помехоустойчивости канала связи требуется знание зависимости вероятности ошибки на бит ( ) 0,,,,bb P E N n k m F =ψ от СШ на входе приемника (Eb / N0), параметров помехоустойчивых кодов (n, k), глубины перемежения m и параметров замираний F для каналов с время-селективными (блочными) замираниями затруднительно. В качестве параметров замираний понимается вероятностное распределение, которому они подчиняются, а также время когерентности канала связи, нормированное относительно длительности передачи символа кода. Для описания замираний в канале связи ограничимся рассмотрением распределений Накагами (-m, -q) и Райса. Как показали исследования [10], Накагами-m распределением наиболее хорошо описывается плотность распределения амплитуды сигнала в трансионосферном канале связи. При этом данное распределение достаточно близко аппроксимируется распределением Райса [11]. Распределение Накагами-q (или Хойта) обычно наблюдается в каналах спутниковой связи, подвергающихся сильной ионосферной сцинтилляции [11]. Получить зависимость ( ) 0,,,,bb P E N n k m F =ψ в виде точного аналитического выражения затруднительно. Известные [4; 6-8] методы оценки помехоустойчивости приема сигналов в каналах связи с кодированием и время-селективными замираниями ограничены допущением, что в пределах пакета ошибок амплитуда и начальная фаза сигналов не изменяются. Для данной модели известен метод частично-когерентного приема, предполагающий «когерентное» сложение отсчетов на выходах каналов приема синфазной и квадратурной составляющих сигналов с последующим их «некогерентным» объединением. Порядок комбинирования когерентных решений в некогерентные при замираниях по законам Райса (c параметром глубины замираний K) и Накагами-m определен соответственно в [6; 7]. В [4; 8] проанализированы характеристики приемника с «субоптимальным» квадратичным комбинированием когерентных «мягких решений» синфазного и квадратурного каналов приема в условиях блочных замираний. Однако частично-когерентный прием неприменим для каналов связи с псевдослучайной перестройкой рабочей частоты (ППРЧ). В данных каналах связи время, в течение которого начальная фаза остается неизменной, определяется значением, обратным скорости перестройки рабочей частоты, и может быть существенно меньше времени когерентности канала связи. Последнее будет определять только время, в течение которого остается постоянной амплитуда сигнала, при условии, что замирания не являются частотно-селективными в пределах полосы частот сигнала с ППРЧ. Назовем такие блочные замирания некогерентными. Известен общий метод расчета асимптотически точных верхних границ вероятности ошибки Pb с использованием границы Чернова для каналов связи с кодированием и псевдослучайным перемежением при наличии группирования ошибок, природа возникновения которых не конкретизирована [5]. Целью статьи является разработка метода расчета верхней границы вероятности ошибки на бит в «некогерентных» каналах связи при использовании метрики с линейным сложением мягких решений демодулятора в условиях воздействия блочных замираний. Модель канала связи Пусть в канале связи для передачи информации используется блочный (n, k) код, имеющий скорость / = r k n и минимальное расстояние dmin. Дистанционные свойства кода характеризуются коэффициентами /, = dd A dC n где Сd - число кодовых комбинаций весом d. Кодовая комбинация состоит из символов кода xi, где 1, , = in которые могут принимать одно из целых чисел от 0 до 1, -M где М - объем алфавита символов кода, и образует последовательность 12 ,( ,= xxx ),. … n x Символы кода передаются соответствующими ортогональными сигналами из ансамбля 0{ ( ) , St 1() , St …, -1( )}. M St Объем ансамбля используемых сигналов равен объему алфавита кодовых символов М. На вход приемника поступает смесь полезного сигнала и теплового шума. Кроме того, полезный сигнал подвергается замираниям. Энергия сигнала, приходящаяся на бит информации, на входе демодулятора определяется величиной Eb. Тепловой шум () nt представляет собой аддитивный белый гауссовский шум с односторонней спектральной плотностью N0 и является независимым от амплитуды замирания. Замирания представляют собой мультипликативную помеху, в результате воздействия которой мгновенная мощность принятого сигнала модулируется амплитудой замирания aj [7]. Амплитуда сигнала aj является случайной величиной со среднеквадратичным значением ( ) 2 .Ω= j Ea При расчете энергетического бюджета радиолинии величину Ω следует учесть в качестве дополнительных потерь (например, на поглощение энергии сигнала в среде распространения радиоволн). Флуктуационной составляющей амплитуды сигнала является нормированная интенсивность этого сигнала 2 . = Ω jj ca (1) Отношение энергии сигнала с замираниями к спектральной плотности шума имеет среднее значение 2 2 0 log 2 γΩ = sn b E rM N (2) и мгновенное 2 2. ξ = γ j j sn c (3) Время когерентности канала связи τF определяет количество кодовых символов B, в течение передачи которых амплитуда (интенсивность) сигнала остается постоянной. При псевдослучайном перемежении с глубиной , ≤ mB как показано в [5], возникают сдвоенные пакеты (блоки) либо одинаковой длины /,= b B m либо с длинами -1 b и 1. +b На интервале передачи одной кодовой комбинации возможно возникновение 2 =   P n b таких блоков. Кроме того, в случае, когда длина кода не делится на величину 2b без остатка, существует блок длины -2 . n bP В демодуляторе принятая канальная последовательность обрабатывается параллельно в детекторах сигналов 0() , St 1() , St …, -1 ). (M St На выходе детекторов по окончанию приема i-го символа кода формируются мягкие решения, нормированные значения которых относительно дисперсии теплового шума образуют вектор 0, , (= ii YY 1, , iY ..., -1, . ) MiY Примем, что точная и своевременная оценка начальной фазы передаваемых сигналов затруднительна и возможен только «некогерентный» прием. Тогда в случае воздействия замираний и теплового шума величина Yk,i, где 0, 1, = - kM имеет следующую плотность распределения [9]: ( ) , , ; = ki Y k i j i p y z x ( ) ( ) 2 , 2 0, , 1exp 22 , , 1exp , . 22   +γ   -×   = ×γ =    -≠   k i j sn k i j sn i ki i yc I y c kx y kx (4) Декодирование производится по максимуму правдоподобия на основе решений демодулятора с использованием метрики m(Yi, xi; cj). Под воздействием замираний и теплового шума возможна трансформация переданной кодовой последовательности x в последовательность ( ) 12 ˆ ˆ ˆ ˆ , ,..., ,= n x x xx отличающуюся от переданной d символами кода. Без потери общности примем, что 0 =ix и ˆ 1 =ix для 1, . = in Верхняя граница вероятности ошибки в канале связи с группированием ошибок Верхняя аддитивная граница вероятности ошибки на бит имеет вид [19] ( ) min ˆ, , = ≤ ∑ n b dd dd P A P Exx (5) где ( ) ˆ,d PE xx - вероятность ошибочного приема кодовой комбинации ˆ x с расстоянием Хэмминга d вместо переданной x. Для определения ( ) ˆ,d PE xx можно использовать производящие функции моментов разности метрик для принятых ошибочно и переданных символов в блоках. Пусть d - i ошибочных символов содержится в P блоках длиной 2b, а остальные i - в блоке длиной -2 n bP (в 1 +P блоке). Производящую функцию моментов разности метрик для i ошибочных символов в 1 +P блоке обозначим как ( ) 1 , , 2 , - D s i n bP производящую функцию разности метрик для d ошибочных символов при их распределении по P блокам длиной 2b - как ( ) , ,2 . P D s d b Производящие функции моментов разности метрик ошибок для блоков и пакетов назовем весовыми функциями. В этом случае вероятность ошибочного приема кодовой комбинации в (5) определяется согласно выражению [5] ( ) ( ) ( ) ( ) ( ) min , 2 ˆ, 22 max 0, 2 1 1 , , 2 , ,2 , - - - = - ≤× × - - ∑  d n bP d i i d bP n bP d i d bPn P PE CC C D s i n bP D s d i b xx (6) где ! !( )! = -d n C n d n d Весовая функция P блоков длиной 2b в (6) описывается выражением ( ) ( ) ( ) ( )( ) ( )( ) 1 2 1 2 2 ! ! , ,2 2! , ; , ,..., , ! -+ =  - = × = × -∑      P F dl l d b bF d F D s d b bF A d l g g g Fl g (7) где ( ) ( ) 1 2 1, ; , ,..., -+=    dlA d l g g g g - однородные полиномы Белла, l - количество блоков с «ненулевым» весом, максимальное число которых ( ) min , .= F d P Формула для расчета v-го элемента последовательности  g в выражении (7) имеет вид ( )( ) ( ) 1 , ,2 2 ! 2 !ν = -ν  g D s v b b b (8) Здесь весовая функция блока длиной 2b определяется как ( ) ( ) ( ( )( ( ) ( ) ( ) ( ) ( ) ( ))) 11 2 min , 1 , 1 1 max 0, 11 ,11 1, ,2 , ,2 1 , ; 1 , ; 1 , ; , ; , - - -+ = = - - = ∆ + + - × × - - + + +- ∑  v b vb i v i b b b b i k v b i v i b b b b D s v b D s v b C p C C D s i b D s v i b p C C D s i b D s v i b (9) где ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 ,1 , 1 1 ,1 ,1 1 1 ,1 1 1 , ,2 0, , , ; , , ,1; , ; 1,1 , 1; 1 , ,; , ; 1 , 1; 1 1., 1; 1 , - -- - ∆= <  =   +  = ++ - + + =  -×  × + - ×  × - - - ×  >+ × + +   bb bb bb vb b b b vb b b b D s v b vb p D s b b v b p bD s b D s b b vbp D s b b p C D s v b b D s b b p C D s v b b vbD s b b (10) В выражениях (10), (11) , 2 1 /4 () = +bb pPP - вероятность появления блока из пакетов оди наковой длины, 1() , , , D s v b 1 , , , () -1D s v b 1 , ( ,1 )+D s v b - весовые функции пакетов весом v и длиной , -1, b 1 +b соответственно. Входящая в (6) весовая функция 1 +P блока длиной -2 n bP и весом i имеет вид ( ) ( ) ( ) ( ) ( ) ( ) 1 1 min , max 0, 1 1 , , 2 , ; , , ,; ., ; , - - - = - +- -=  - -   × =  × - ×   - ×- ∑  ÷åòíîå íå÷åòíîå i n bL k i k n bL b i k i b b n bL D s i n bP D s i n bL L CC C D s k n bL LD s i k b (11) В выражении (11) =   L n b - число пакетов длиной b в пределах кодовой комбинации длиной n, а 1 ,,() - D s v n bL - весовая функция пакета весом v и длиной -. n bL Весовые функции пакетов ошибок Весовая функция пакета весом v при наличии в канале связи замираний не зависит от длины пакета x и определяется выражением ( ) ( ) ( ) ( )1 1 0 , , , ;,∞ Ω ν = ν = = Φ ∫ j v j c j j D s x D s s c p c dc (12) где ( ) j cj pc - плотность распределения вероятности нормированной интенсивности замирающего сигнала (1) . jc В выражении (12) производящая функция моментов ( ) 1 ;ΩΦ j sc разности метрик принятого ошибочно и переданного символов кода (далее - разность метрик ошибки) ( ) ( ) 1 ˆ, 1; , 0; Ω = = - = i j i j m x c m x c ii YY (13) связана с плотностью распределения ( ) 1 1;Ω w j pc случайной величины Ω1 соотношением ( ) ( ) ( ) 11 111; exp ; .∞ ΩΩ -∞ Φ = w w w ∫ jj s c s p c d (14) Линейное сложение мягких решений некогерентных квадратичных детекторов Используя результаты [9], выражение (12) в случае использования для декодирования метрики с линейным взвешиванием мягких решений некогерентных квадратичных детекторов может быть конкретизировано следующим образом: ( ) ( ) 2 2 12 ,, 2 1 214 j sn cv sDs M v ss  γ ν= - + - (15) где ( ) ( ) ( ) 0 exp ∞ =∫j c j j j M p pc p c dc (16) представляет собой производящую функцию моментов интенсивности j c сигналов (зависящую от среднего значения СШ (2) 2 2γ= sn 20 log ).=Ω b E r M N Для замираний, подчиненных закону Райса с параметром 0, ≤ <∞ K производящая функция моментов интенсивности сигналов (16) примет следующий вид [11]: ( ) ( ) ( ) ( ) 1 exp . 11 + =  + - + - jc K KpMp K p K p (17) Для замираний, подчиненных закону Накагами-m c параметром 0,5 , ≤ <∞ m выражение (16) сводится к следующему [11]: ( ) 1. -  = -  j m c pMp m (18) Для замираний, подчиненных закону Накагами-q (или Хойта) с параметром 0 1, ≤≤ q производящей функцией (16) является [11] ( ) ( ) ( ) 1/2 2 2 22 2 12 . 1 -   = - +  +  jc pq M p p q (19) Линейное сложение мягких решений некогерентных детекторов огибающей Рассмотрим теперь случай линейного взвешивания мягких решений некогерентных детекторов огибающей. Получить выражение, аналогичное (15), для указанного способа сложения затруднительно. Применим для оценки весовой функции ( ) , Dsv следующий метод. В интеграле (14) производящая функция моментов разности метрик ошибки ( ) 1 ;ΩΦ j sc является, по сути, верхней границей Чернова вероятности ошибочного приема кодового символа, которая монотонно уменьшается при увеличении мгновенного отношения СШ (3) 2 2, ξ = γ j j sn c и может быть представлена в следующем виде: ( ) ( ) ( ) ( ) 1 2 2 ;2 exp 2 , Ω Φ = ϕ γ × × - γ j j sn j sn s c C s c c (20) где () Cs - нормирующий коэффициент, () ϕ t - мультипликативная функция. Если разложить мультипликативную функцию () ϕ t в ряд ( ) 1 0 1 1 11 , ! ! ! ∞ ∞ ∞ - = = = ϕ = = = ∑ ∑ ∑ i i i i i i i i i t t t t a ia q i t i t i (21) где ( ) ( )( ) 1 2 1 01 , ,..., ,... ,2 ,..., 1 ,... , += = = + i i q q q a a i a q (22) то весовая функция (12) примет вид взвешенной суммы [12] ( ) ( ) ( ) 2 ,! ,, , , , !2 ∞ = ν = ×  γ ×-   ∑ j v sn cj iv D s v C s Aiv D c v i v i q (23) где 22 2 ,, 22 exp . 2 - -     γγ - = ×           γ ×-       j iv sn sn c j sniv j D v i v c E c v (24) Однородные полиномы Белла ( ) , ; ,Aiv q с учетом того, что ( ) 1 1,1; , = Aq q могут быть определены из следующего рекуррентного соотношения [13]: ( ) ( )1 1 1 1 ,; , 1, , 1; , 1. + + - +- = += = =  + - - >   ∑ vi i v v i i j A v i v qv C g A v i j v v q q (25) Для метрики с линейным сложением мягких решений детектора огибающей ( ) , ,; =j k i m k c Y i Y (26) производящая функция моментов разницы метрик ошибки описывается выражением [14] ( ) ( ) ( ) 1 2 2 ; exp 2 _ , , Ω  γ Φ = - ×    ×γ j sn j j sn c sc FS s FS E c s (27) где ( ) 2 1 exp 1 , 22 2 p = ++     ssFS s s erf (28) ( ) ( ) ( ) 2 00 _, exp 2 .∞ = = - - ∫ FS E p s u su u I up du (29) Используя в (29) разложение модифицированной функции Бесселя нулевого порядка ( ) 0 I up в ряд, запишем выражение (27) в виде ( ) ( ) ( ) ( ) ( ) 1 2 21 2 2 0 0 ; exp 2 exp 2 . 2! Ω ∞ ∞ + = Φ = - γ × × -- ∑ ∫ j j sn i i i i s c c FS s t u su u du i (30) «Infokommunikacionnye tehnologii» 2019, Vol. 17, No. 4, pp. 372-382 Shevchenko V.A., Pashintsev V.P. 377 Интеграл в правой части выражения (30) используется в интегральном представлении функции параболического цилиндра ,() Uas [15]: ( ) ( ) ( ) ( ) 12 2 2 0 exp 2 1 2 exp 4 , , -∞ - - = = Γ + ∫ ax x sx dx a s U a s (31) где 2 3 2. = + ai С учетом (31) выражение (30) примeт следующий вид: ( ) ( ) ( ) ( ) 1 22 ; exp 2 2 , Ω Φ= = - γ ϕ γ j j sn j sn sc C s c c (32) где ( ) ( ), = C s FS s (33) ( ) ( ) 0 2 3 2, . ! ∞ = ϕ = + ∑ i i i t t BU i s i (34) В выражении (34) ( ) ( ) ( ) ( ) ( ) 1 10 2 1 ! 2 1 ! 21 2 ! 2 1 ! 2 1 , 1. - - +- = = + = - = + = i ii i ii Bi ii iBB (35) Функция параболического цилиндра ( ) , Uas может быть найдена с использованием рекуррентного соотношения [15] ( ) ( ) ( ) 12, 3 2, 1 2, , += - - - = U n s U n s sU n s n (36) с учетом того, что ( ) ( ) ( ) ( ) ( ) 2 2 1 2, exp 4 , 1 2, 2exp 4 2 . - = - = p U s s U s s erfc s (37) Из сравнения выражений (34) и (21) можно получить формулы для определения переменных полиномов Белла (22) для метрики с линейным сложением мягких решений детекторов огибающей: ( ) ( ) 111 10 , 2 1 2, , 3 2, . --- = = - = = i i i i q ia a B U i s q a U s (38) Величина ( ) 2 2, , γ- j c sn D v i v в выражении (24) с точностью до знака аргумента является обобщенной производящей функцией моментов мгновенного отношения сигнал/шум 2 2 ξ = γ : j j sn c ( )( ) ( ) ( ) ( ) ( ) 2 2 2 0 2, , 1 exp 2 . 2 ξ ∞ φ = γ - = = ξ ξ =ξ γ ξ γ ∫ jj j i c sn i j j c j j sn j sn s D s i s p c d (39) Обобщенная производящая функция моментов мгновенного отношения сигнал/шум ( ) 2 2, , γ- j c sn D s i для распределений Райса и Накагами известна [12]. Это позволяет определить выражения для весовых функций ( ) ,ν Ds согласно (23). Для замираний по закону Райса весовая функция примет вид ( ) ( ) ( ) ( ) ( ) 2 1 , exp 2 , , , , ∞ - = +  ν= - × γ   ×-   ∑ v sn i i v iv K aDs K C s pp a B p v L A i v p q (40) где ( ) 2 1 , 2 + = γsn KK a 2 1 , 2 + = + γsn Kpv (41) ( ) ( ) ( ) 1 0 1 , 1 , , , 1. -  = -   = ii v B p v B p v ip B p v (42) Для полиномов Лаггера справедлива рекуррентная формула [16] ( ) ( ) ( ) ( ) ( ) ( ) 11 2 1 1. +- = + - - +i ii L x i x L x iL x i (43) Когда замирания подчинены распределению Накагами-m, весовая функция определится выражением ( ) ( ) ( ) ( ) 2, 2 , , , , ∞ =  ν= ×  γ  ×∑ m v sn i iv mDs C s p B p v A i vq (44) где ( ) 2 2, = + γ snpvm (45) ( ) ( ) ( ) 1 0 ,1,1 , , 1. -+-  = -   = i i B p vvmB p v ip B p v (46) При наличии в канале связи замираний Накагами-q будем иметь ( ) ( ) ( ) ( ) 2 222 22 11, 22 ,, , , , ∞ - = + ν= × γ-  ×   -  ∑ v sn i iv iv qDs C s qpa pB p a v P Aiv pa q (47) где 4 22 1 , 42 - = γsn qa q ( )2 2 22 1 , 42 + = + γsn q pv q (48) ( ) ( ) ( ) 122 0 1, , 1 ,, , , 1. -  = -  - = ii vB p a v B p a v i pa B p a v (49) Полиномы Лежандра в выражении (47) могут быть вычислены по рекуррентной формуле [15] ( ) ( ) ( ) ( ) ( ) 11 01 21 , 11 1, . +- + = - ++ = = n n n nn P x xP x P x nn P x P x x (50) Анализ полученных результатов Таким образом, разработан метод определения верхней границы вероятности ошибки на бит в «некогерентетных» каналах связи при использовании метрики с линейным сложением мягких решений демодулятора в условиях воздействия блочных замираний с распределениями Райса и Накагами. Данный метод предусматривает выполнение следующих операций. 1. Расчет вероятности ошибки на бит b P в виде совокупности выражений (5)-(12), (14) как комбинации весовых функций ( ),, νDs представляющих собой производящие функции моментов разности метрик для принятых ошибочно и переданных символов кода в пакетах ошибок. 2. В случае, когда для декодирования используется метрика с линейным сложением мягких решений квадратичного детектора, весовая функция ( ) ,ν Ds определяется выражением (15) и зависит от среднего отношения СШ (2) 2 2γ= sn 20( )log=Ω b E k n M N и производящей функции моментов интенсивности сигнала ( ). jc Mp Последняя конкретизирована для замираний по закону Райса, Накагами-m и Накагами-q в виде выражений (17)-(19), зависящих от глубины замираний (K, m, q). 3. В случае, когда для декодирования используется метрика с линейным сложением мягких решений детектора огибающей, весовая функция ( ) ,ν Ds представлена в виде взвешенной суммы (23) обобщенных производящих функций моментов мгновенного СШ возрастающего порядка с использованием полиномов Белла в качестве взвешивающих коэффициентов. Выражения для ( ) ,ν Ds конкретизированы для обобщенной производящей функции моментов мгновенного СШ при замираниях по закону Райса (40), Накагами-m (44) и Накагами-q (47). Все выражения получены в виде рекуррентных формул, не требующих значительных вычислительных затрат. Отметим, что параметры распределений K, m, q связаны однозначно с количеством замираний AF [11] ( )2 12 , 1 + = + KAF K 1 , =AF m (51) ( ) ( ) 4 22 21 , 1 + = + q AF q которое представляет собой нормированную дисперсию амплитуды сигнала ( )( ) 22 2 . -Ω = Ω j Ea AF (52) Количество замираний AF может быть использовано в качестве универсального показателя глубины замираний. К примеру, для замираний, обусловленных сцинтилляцией в ионосфере, количество замираний является квадратом индекса сцинтилляции: 2 4= AF S [20]. Разработанный метод позволяет получить зависимости вероятности ошибки на бит = b P ( ) 0, , , ,( , )b E N n k b AF LW=ψ от отношения сигнал/шум на входе приемника (Eb / N0), параметров помехоустойчивых кодов (n, k), глубины перемежения (m), определяющей длину пакетов ошибок /,= b B m глубины замираний (с учетом взаимосвязи параметров Райса K и Накагами m, q с количеством замираний AF) и закона, которому они подчинены LW, для каналов с время-селективными (блочными) замираниями. Воспользуемся полученными зависимостями для подтверждения реализуемости метода повышения помехоустойчивости каналов связи с замираниями при наличии ограничения на время задержки передачи информации tD за счет использования блочных (n, k) кодов большой длины и с высокой корректирующей способностью в сочетании с перемежением, допускающим группирование ошибок в пакеты. Отношение допустимого времени задержки передачи информации tD к времени когерентности канала связи τF определяет соотношение длины кода n и средней длины пакетов ошибок /= : b B m . = = τ D F t n n B m b (53) Из выражения (53) следует, что увеличение глубины перемежения m и уменьшение длины пакетов ошибок / = b B m требует укорочения длины кода n и наоборот. На объем алфавита кодовых символов степень группирования ошибок влияния не оказывает. Выбор 2 = mB позволяет исключить группирование ошибок. В качестве примера оценим с использованием полученных выражений помехоустойчивость каналов связи с применением двоичных () 2=M расширенных кодов Боуза - Чоудхури - Хоквингхема (БЧХ) длиной от 8 =n до 128 =n в условиях наличия релеевских замираний (когда 0, =K что соответствует 1 =m и . )1=q Распределение весов кодовых комбинаций известно для всех расширенных кодов БЧХ длиной 128≤n [17; 18]. Выберем расширенный двоичный код БЧХ с 8 =n как код минимальной длины. В канале с эталонным кодом () 8=n исключить группирования ошибок () 2= mB возможно при / 16. = nb Выражение / 16 = nb определяет также длину пакетов ошибок при использовании кодов с большей длиной. Выберем коды, которые для каждой из рассматриваемых длин n кодов (8, 32, 64, 128) обеспечивают для допустимого значения 510 -=bP наименьшее значение отношения Еb/N. Для указанных кодов на рисунке приведены зависимости вероятности ошибки на бит Pb от отношения Еb/N0, полученные согласно выражениям (5)-(12), (23), (40) применительно к каналам связи, в которых для декодирования используется метрика с линейным сложением мягких решений детекторов, огибающей при наличии блоковых ( /6 )1=nb релеевских ( 0, =K 1, =m ) 1=q замираний. Выводы Анализ графиков, представленных на рисунке, позволяет сделать следующие выводы: 1. Применение результатов разработанного метода позволяет получить искомые зависимости ( ) ( ) 0, , , , ,bb P E N n k b AF LW =ψ вероятности ошибки на бит от отношения СШ на входе приемника (Eb / N0), параметров помехоустойчивых кодов (n, k), длины пакетов ошибок / = b B m и глубины замираний (количества замираний, определяющего значения параметров Райса K и Накагами m, q) для каналов с время-селективными (блочными) замираниями с учетом ограничения на задержку передачи информации / / . τ= DF t n b 2. При заданном ограничении на задержку передачи информации tD и измеренном времени когерентности , τF глубине замираний AF и статистическому закону, которому они подчинены LW, в канале связи возможно выбрать параметры кода (n, k) и глубину перемежения (m, /, )= b B m обеспечивающие достижение допустимой вероятности ошибки 5 10-=bP при СШ на входе приемника 0/ 16 18 = …b EN дБ. Это указывает на возможность получения энергетических выигрышей 32…34 дБ по сравнению со случаем применения одиночного (неразнесенного) приема замирающих сигналов. 3. Использование «оптимальных» кодов длиной 64; =n 128 позволяет обеспечить отношение Еb / N0, требуемое для достижения вероятности ошибки на бит 5, 10-=bP на 3…4 дБ меньше по сравнению со случаем, когда используется код с восьмикратным повторением ( 8, =n . )1=k При этом учтено, что длина пакетов ошибок при увеличении длины кода растет в соответствии с соотношением / 16. = nb Для рассмотренных расширенных кодов БЧХ зависимости вероятности ошибки Pb от отношения СШ Еb / N0, определяемые согласно выражениям (5)-(12), (15) и (17) применительно к Рисунок. Зависимость вероятности ошибки Pb от отношения сигнал/шум Eb/N0 для расширенных кодов БЧХ при наличии в канале связи блоковых релеевских замираний каналам связи, в которых для декодирования используется метрика с линейным сложением мягких решений квадратичных детекторов и имеются блоковые релеевские замирания, практически не будут отличаться от приведенных на рисунке. Заключение Разработан метод расчета верхней границы вероятности ошибки на бит в «некогерентетных» каналах связи при использовании метрики с линейным сложением мягких решений демодулятора в условиях воздействия блочных замираний с распределениями Райса и Накагами. Полученные с его использованием зависимости вероятности ошибки на бит от СШ позволяют подтвердить реализуемость метода повышения помехоустойчивости каналов связи с замираниями при наличии ограничения на время задержки передачи информации за счет использования двоичных блочных (n, k) кодов большой длины и с высокой корректирующей способностью в сочетании с перемежением, допускающим группирование ошибок в пакеты.

About the authors

V. A Shevchenko

Department of Information Systems of the Ministry of Defense of the Russian Federation

Email: shevv67@mail.ru
Moscow, Russian Federation

V. P Pashintsev

North Caucasus Federal University

Email: shevv67@mail.ru
Stavropol, Russian Federation

References

Supplementary files