CC BY
50
УДК 621.391
С.В. ЗАЙЦЕВ*
МЕТОД ПІДВИЩЕННЯ ДОСТОВІРНОСТІ ПЕРЕДАЧІ ІНФОРМАЦІЇ В РОЗПОДІЛЕНИХ ПРОГРАМОВАНИХ БЕЗПРОВІДНИХ СИСТЕМАХ З ТУРБОКОДАМИ
Чернігівський державний технологічний університет, Чернігів, Україна
Анотація. У статті запропоновано новий метод підвищення достовірності передачі інформації в розподілених програмованих безпровідних системах, сутність якого полягає в застосуванні для підвищення надійності прийняття рішень про декодований біт додаткової інформації, отриманої шляхом кодування жорстких рішень з виходу кожного компонентного декодера турбокоду. Ключові слова: разподілені програмовані безпровідні системи, турбокоди.
Аннотация. В статье предложен новый метод повышения достоверности передачи информации в распределенных программируемых беспроводных системах, сущность которого заключается в применении для повышения надежности принятия решений про декодированный бит дополнительной информации, полученной путем кодирования жестких решений с выхода каждого компонентного декодера турбокода.
Ключевые слова: распределенные программируемые беспроводные системы, турбокоды.
Abstract. This paper proposes a new method for increasing the reliability of data transmission for distributed programmable wireless systems, the essence of which is to apply to increase the reliability of decision-making about the decoded bits of additional information obtained by encoding the tough decisions from the output of each component of turbo code decoder.
Keywords: distributed programmable wireless systems, turbo codes.
1. Вступ
Значне розширення за останнє десятиліття частотної смуги мереж передачі даних привело до появи нового класу безпровідних систем - розподілених програмованих безпровідних систем, які базуються на принципах SDR (software defined radio). Розподілені програмовані безпровідні системи містять декілька обчислювальних вузлів, об’єднаних високошвидкіс-ним інтерфейсом передачі даних [1]. Зазначені системи використовують апаратні засоби для виконання функцій під керуванням програмного забезпечення. Розроблювачі цих засобів покладають великі надії на відкриту архітектуру SDR-мереж наступного покоління, де обладнання безпровідних мереж можуть випускати одні виробники, а функції й режими роботи буде визначати програмне забезпечення інших розроблювачів.
Розподілені програмовані безпровідні системи дозволять підвищити якість обслуговування за допомогою високої адаптивності, спрощують розробку програмної частини обладнання й створюють нові можливості при реалізації нових стандартів [2].
Одним із режимів роботи розподілених програмованих безпровідних систем є режим роботи в умовах впливу організованих завад, які за певних умов здатні істотно погіршувати характеристики достовірності передачі інформації [3]. Одним із перспективних напрямів протидії організованим завадам є ефективні сигнально-кодові конструкції на основі завадостійких турбокодів (ТК) з характеристиками завадостійкості, близькими до межі К. Шеннона [4, 5].
Виникає завдання підвищення достовірності передачі інформації в розподілених програмованих безпровідних системах з турбокодами.
Метою роботи є розробка методу підвищення достовірності передачі інформації в розподілених програмованих безпровідних системах з турбокодами.
© Зайцев С.В., 2013
ISSN 1028-9763. Математичні машини і системи, 2013, № 3
2. Виклад основного матеріалу
Розглянемо структурні схеми кодера (передавальна частина) та декодера (приймальна частина) турбокоду.
На рис. 1, 2 показані структурні схеми двокомпонентного кодера та ітеративного декодера ТК відповідно при паралельному з'єднанні двох рекурсивних систематичних зго-рточних кодів (РСЗК).
Схема кодера турбокоду використовує РСЗК зі швидкістю 1/ п виду (1, ё1/ & о, , &п-1 / & о), де & о - поліноміальний генератор зворотного зв'язку, а g1,..., gn_1 -
поліноміальні генератори прямих зв'язків. Кожен РСЗК виконує кодування інформаційної послідовності за своєю діаграмою, структура якої залежить від поліноміальних генераторів РСЗК [6, 7].
Рис. 1. Структурна схема двокомпонентного кодера ТК
Послідовність на виході кодера ТК має вигляд: X = (XС, XП), де XС = и - система- • П —П1 П2 . „ .
тичний вихід кодера, а X = (X , X ) - перевірочний вихід кодера ТК. При цьому
П1 П11 П1у . „ . „пт- ч П2 П21 П2у . „
X = (X ,..., X ) - перевірочний вихід РСЗК 1, X = (X ,..., X ) - перевірочний вихід РСЗК 2, V - загальна кількість перевірочних символів кожного РСЗК кодера ТК.
Демодульована послідовність символів подається на декодери 1 і 2 (рис. 2):
—1 —С1 —П1 —П1 /І7П11 Т7шу
У = (ЬСУ , ЬСУ ) - для декодера 1, де У = (У ,..., У ), Ьс _ параметр канальної на. . . . —2 —С2 —П2 —П2 —П21 — П2у
дійності . Відповідно У = (ЬСУ ,ЬСУ ) - для декодера 2, де У = (У ,...,У ).
— С1 —С —02 . . . ...
У = У , У - послідовності систематичних символів з урахуванням відповідної операції перемеження (П).
Розглядається такт роботи в момент часу ,.
Дійсне представлення “м'якого” рішення або логарифмічне відношення функцій правдоподібності (ЛВФП) поза декодером визначається виразом [5, 6]
ДА і Уі) =1п—+17 +1прт1—+1) = La(х() + Ь(у, І ), (1)
р( у, 1 = -1) р( х, = _1)
де Ь(у, | х,) - ЛВФП у,, яке одержується шляхом виміру у, на виході каналу при чергуванні умов, що може бути переданий х, =+1 або х, =-1, а Ьа (х,) - апріорне ЛВФП біта даних х,. Для спрощення позначень рівняння (1) може бути переписане таким чином [5, 6]:
Д(хі ) = ДС (у, ) + Да (х, ). (2)
Тут Ьс (у,) означає, що член ЛВФП виходить у результаті канальних вимірів, зроблених у приймачі. Для систематичних кодів ЛВФП на виході декодера дорівнює [5, 6]
Ь( х,) = Ь (х) + Ье (х). (3)
Ітерація 1
Рис. 2. Структурна схема двокомпонентного декодера ТК
У цьому виразі Ь'(х() - ЛВФП поза демодулятором (на вході декодера), а Ье (х,) -“зовнішнє” ЛВФП, що представляє зовнішню інформацію, яка випливає з процесу декодування. Вихідне ЛВФП декодера прийме вид
Ь( х) = Ьс (Уі ) + Ьа (хі) + Ье (хі). (4)
Знак Ь(х() є твердим рішенням про символ х,, а модуль | Ь(х,) | - ступенем надійності (правдоподібності) цього рішення.
Декодер і у відповідності зі своїм алгоритмом виробляє “м'які” рішення про декодовані символи (вихідне ЛВФП), які складаються з трьох частин [5, 6]:
Е(хС) = Ес • ус1 + Еа(Xе) + Ее(Xе), (5)
де хС - систематичний символ кодера ТК.
При цьому “зовнішня” інформація декодера і про символ ххс, що є апріорною для декодера і +1 (з урахуванням операції перемежіння), прийме вид [5, 6]
Ее (хС) = Д+ЧхС) = Е (хС) - Еа(хС) - Ьс • уС1. (6)
Декодер і +1, одержавши апріорні відомості про інформаційні символи, робить аналогічні обчислення, визначаючи свою “зовнішню” інформацію про символ ххс [6]:
С(хС) = їа(хС) = Е+1(хС) - Е^) - Ьс • уС2, (7)
яка надходить на вхід декодера і + 2 наступної ітерації декодування.
Після виконання необхідної кількості ітерацій або у випадку примусової зупинки ітеративної процедури декодування виносяться рішення про декодовані символи:
хС = Г1, якщ° П (хС) > 0 (8)
х [0, якщо Е (хС ) < 0
Як відомо, декодування символів ТК відбувається за діаграмою відповідного РСЗК. При цьому обчислюються перехідна рекурсія, пряма рекурсія, зворотна рекурсія, ЛВФП на виході декодера та параметр “зовнішньої” інформації [6].
Рішення про декодований переданий біт приймається за інформацією про стан каналу, прийняті інформаційні та перевірочні біти, а також за апріорною інформацією, яка
отримується за результатами декодування попереднім декодером. Надійність прийняття рішення про декодований біт інформації можна підвищити за рахунок використання додаткової інформації при розрахунку ЛВФП про переданий біт.
Для підвищення надійності прийняття рішення про декодований біт пропонується використовувати при розрахунку ЛВФП наступної ітерації декодування інформацію, отриману шляхом повторного кодування жорстких рішень з виходу кожного компонентного декодера ТК. Структурна схема, яка реалізує цей алгоритм, показана на рис. 3.
Отримаємо математичну модель ітеративного декодування ТК з урахуванням додаткової інформації при розрахунку ЛВФП про переданий біт.
Введемо такі позначення:
N - розмір кадру переданих символів; к - поточний індекс часу, к є{1, N};
ск = {С, сі,..., сі} - кодовий символ, вироблений РСЗК, скт є (-1,1), т є (1, і);
ґ 1,5 2, р 3, р і, р ч І 1,5 2, р 3, р і,р І , л .ч ^ .
хк = (хк ,хк ,хк ,...хк Xіхк ,хк ,хк ,•••,хк ]є (-ЛА) модульований символ. А
- канальний коефіцієнт, для каналу з адитивним білим гаусівським шумом (АБГШ) А = 1;
Ук = (у*5, УкР, УкР,... уїр) - прийнятий символ; у1 = (у1, у2,... yN ) - один кадр прийнятого символу.
Для двох подій А і В їх спільна ймовірність позначається як Р(А, В) і дорівнює:
Р(А, В)= Р(А)Р(В / А), (9)
Р( А, В / С ) = Р( А / С )Р(В / А, С), (10)
P(A) = £ P(A, B), (11)
B
P(A / B )= P(A), (12)
p(ab ) = p(a)p(b ). (13)
У 1974 році був запропонований алгоритм BCJR (Bahl-Cocke-Jelinek-Raviv) або алгоритм декодування по максимуму апостеріорної ймовірності MAP (maximum a posteriori probabilities), який здійснює розрахунок апостеріорної ймовірності кожного декодованого символу [4, 7].
Логарифмічне відношення функцій правдоподібності (ЛВФП) L(uk) для випадкової двійкової змінної uk визначається таким чином:
P(uk = 1/yN)
'P(Uk = 0/yN)
Рішення може бути прийняте по знаку L(uk), тобто
~k = sign[L(uk )]. (15)
З урахуванням (9) і (11), вираз (14) перепишеться таким чином:
Р( = 1 ^ w р( N) EP(^k-1 = ^ uk = 1,У?)/Р(У?)
W ^ = З)/Ptf)) = loS[i Р(„- = ,,, = 0, )/Р( )]. (16)
s
Ітерація i
Рис. З. Структурна схема модифікованого двокомпонентного декодера ТК
Використовуючи пару (як_І5 ик), sk , перехід по решітчастій діаграмі буде однозначно визначений. Тобто сума по У для спільної ймовірності в чисельнику рівносильна підсумовуванню по и + для переходів з попереднього стану = я' в поточний 5 = ^, викликаних інформаційним символом и{ = 1 і, подібно в знаменнику, підсумовування по и для переходів з попереднього стану в поточний, викликаних інформаційним символом
и = 0.
Прийнятий символ може бути розбитий на три частини. Перша частина містить спостереження перед моментом часу к , друга частина - поточне спостереження, третя частина - спостереження після моменту к :
/НУ-1, У> , Ум}. (17)
Підставивши (17) в (16), отримаємо
X Р^к _1 = ^ ик =1 У1к _\ Ук, уГ+1)/р( у1 )
Р(ик) = і08[^ ^ ^ к_1 N \ /р/ ], (18)
X Р(їк_1 = Я, ик = 0, У1 , Ук, Уш ) /р(У1 )
и
де X ( ) - сума по усіх можливих переходах (як_1, sk) у момент к , викликаних інформа-
+
и
ційним символом и( = 1, і X ( ) - сума по всіх можливих переходах (як_1, sk) у момент к ,
и
викликаних інформаційним символом и{ = 0 .
При використанні формули (9) спільна ймовірність в (18) перетвориться в умовну ймовірність.
X р(як _1 = , У1к-1 )Р(ик = 1 Ук, Ук+1 / sk _1 = , У1к-1) /р( у1 )
Р(ик ) 1об[ Хр(^к _1 = ^ ук-1) р(ик = 0, Ук, ук+1 / ^ _1 = ^ У1к-1)/ р(у1)]
и-
X р( Як -1 = ^> У1к -1) р(ик =1 Ук> ук+1/ sk -1 = я')/ р( у?)
= 1ост[ и +___________________________________________________]
“ X р(*‘-1 = ^ ук-1) р(“к = 0, Ук, ук+1 / Як-1 = Я')/р( уГ )‘. (19)
и-
Використовуючи формулу (10), отримаємо
р(ик = и Ук, УкГ+1/ Як-1 = ^ = р(УкГ+1/ Як-1 = ^ ик = и Ук)р(ик = и Ук/ Як-1 = ^ =
= Р(УкГ+1 /Як = Я)р(ик = u, Ук /Як-1 = Я,). (20)
Нехай ак(я) А р(як = ^У1к), Рк-1(Я,) А р(УГ /Як-1 = Ю, і
7к (я', я) А р(ик = и,ук /як-1 = я') . Підставивши (20) в (19), маємо
X а к-1(Я')р к (я)7к (Я'>Я) / р( У? )
Ц"к) = і08[ X а к-1 (Я' )Ь к Мї к (*'. Я)/р( уГ )]. (21)
и -
Обчислення а к (^, р к-Д^ ї к Я).
7к (5,5) - це умовна ймовірність того, що існує перехід із стану ,ук-1 у стан 8к , викликаний символом ик. Використавши формулу (10), 7к (5 ,5) можна переписати таким чином:
7к(^ 5) А Р(ик = u, Ук / 5к-1 = 5) = р(Ук / 5к-1 = ^ ик = и)р(ик = и / 5к-1 = ^
= Р(Ук/ ск = с)р(ик = и). (22)
ак (5) - спільна ймовірність у стані 5 для моменту часу к . Використовуючи формулу (11) і (9), вираз для ак (5) запишеться таким чином:
ак (5)АР(5к = ^ У1к) = Е Р(>*к-1 = 5', ик = и
(5)АР(Як = ^ У1к) = Е Р(8к-1 = 5', ик = u, У1к-1, Ук) =
5 '/5к =5
= Е Р(5к-1 = 5', У1к-1 )Р(ик = u, Ук/ 5к-1 =5', У1к-1) =
5'/ 5к =5
= Е Р(5к-1 = 5', У1к-1)Р(ик = и Ук / 5к-1 = 5') =
5 '/ 5к =5
= Е а-1(5')7 к5)=Е а-1(5')7 к s), (23)
к-1V / і к V ’ / к -
5'/5к =5 5'
де X ( ) - сума по усіх можливих станах sk_l, які закінчуються у стані 5к = 5. Початкова
5'/5к = 5
Г1 , якщо 5 = 1 умова: а0(5) = <{ .
І 0, якщо 5 Ф 1
Ьк-1 (•? ) - умовна ймовірність у стані У для моменту часу к -1. Вона може бути за-рекурсивно:
Рк-1(5,)АР(Ук/ 5к-1=5,) = Е Р(ик =u,Ук/ 5к-1=5,) =
5/5к-1=5'
= Е Р(УкЛ+1/ 5к-1 = 5', ик = u, Ук)Р(ик = u, Ук/ 5к-1 = 5') =
5 /5к-1=5 '
= Е рСу<"+1/ 5 = *)р(и<- = и у<-/ 5<-1 = 5') =
= X Ьк (5)Ук (5', 5), (24)
5 / 5к_1 =5'
де X ( ) - сума по усіх можливих станах 5, які починаються зі стану 5к _1 = 5і. Початкова
5 / Як _1 =5'
Г1 , якщо 5 = 1
умова: р! = ^ ^ , .
[0, якщо 5 Ф 1
З метою захисту від надмірності даних при апаратній реалізації а к (5) і Р к _1(5|) мають бути нормалізовані відносно у!1, тобто ак (5) і Рк_1(5|) мають бути розділені на
р( уN), де
р(у!) = р( ук-1) р(у! / ук-1). (25)
*/ -Ч-1 =5
З урахуванням нормалізації, вираз (21) запишеться у такому вигляді:
X а к _1(5 ,)Тк(5', 5)рк(5) / р( уЦ)
Ь(ик) = 1ое[^-----------------------------—] =
X ак_1(5 ,)Тк(5', 5)рк(5)/ р(ух )
и_
X ак_1(5 ,)тк (5', 5)рк(5) а (5,) Р (5)
и+______________________ У" ак ~1 (5 ) у („ І „) рк (5)
=1ов[ р(ук-1)р(уN / ук-1) = 1ов[^ р(у1к_1)Ук( , ) р(у! / уГ) -
X ак _1(5 ,)їк(5', 5)рк о) у ак ■1(5') у (5 - 5)_______РкС*)____
и_ к_^, к_к £Пук")Ук( ,)р(у!/уГ)
р( уГ 1) р( у! / уГ 1)
Розділивши ак_!(5,) А ру5) , ~к (5) А р(уРк /^к_1) , отримаємо
X а к _1(5,)т к5)~к(5)
1(ик) = Ь§[у ~ ^ ,• ш^]. (27)
X а к _1(5)У к(5,5)рк(5)
и _
Так як
р(у!) = X р(5к = 5, ук) = X а к (5), (28)
розділимо а к (5) на р(у1к ). В результаті маємо
а(5) а ом = р(*к=5, ук > = а(5) =
р( у,к) X р(5- = 5 у,") X ак(5) X а к _1 (5,)у к (5',5) X а к _1 (5,)у к 5 5)
5'
XX а к _ 1( 5,)у к (s', 5) XX а к _ 1( 5,)у к (s', 5)
= XX а к .2(і')ук .1(5', 5)р( уN / ук-1).
Розділивши Р к _1(5|) на (31) і використовуючи рівняння (24), маємо Р к _1(5'):
рк _1(5= Рк _1( 5
р(ук_1 /ук_2) XXак_2(5’)Ук_1(5',5)р(уN /ук-1)
(29)
р(у1) можна представити як декомпозицію
р( у! ) = р( ук_1, уN)=р( ук_1)р( уN / уГ) = р(у!_1 / ук_2) р( ук_2). (зо)
З урахуванням (23), (28) і (30), отримаємо
р( ук _2) = р( ук -) = XX а к _2(5 ')у к _(5', 5) =
р(у1 ) 5 5' р(у1 )
(31)
5 5
5 5
X ~* (л)У * (л\ л)
XX а * _
2(л')У*-1(л\ л) XX а * - 2(ОУ*-1(^', Л)'
З рівняння (22) видно, що ук (5 ,5) складається з двох частин.
Нехай р(ск = _1) = р_ и р(ск = 1) = р+ . Використовуючи р(с[ = _1) + р(с\ = 1) = 1 і р(ик ) = р(с1), визначимо
Е(и*) А 1об р(и* ^ = 1ов^_4- =І0^ Р(С* 1
Р(и* = 0)
Р(с* =-*) 1 - Р(с* = 1)
Р Р
Ьа (с1 ) = Іов ~^ = 1оёт-^,
(33)
Р(с1 =-1) = Р-
Р 7
V/ Р+
-+%
Р /
і / Р+ (
Р
Р. (
Р.
1
"У*
1 +
Р
Р+
,-1
Р
1 +
Р
Р+
"Vі
(34)
(35)
Об'єднуючи (34) і (35) разом, отримаємо
Еа ( 4 )
Р (СІ) = ■
Еа (С1 ),1
Е (С* ),1
1 + Є
-Е (С* )
Є
= Лкє 2
(36)
де А - функція від 4 (ск ) і Ак =------1- .
1 + Є_^ (Ск )
Розглядаючи РСЗК з кодовою швидкістю 1/ q (для кожного вхідного біта РСЗК виробляє q біт), ймовірність отримання символу ук за умови, що був переданий символ ск , р(У = ук / ск ), може бути апроксимована таким чином:
р(У = ук /ск) = р(У = ук /ик)»1ітр(ук <У <ук +А/и = ик) =
А®0
Ііт р(у* /и* )А = Ііш р(у* / х* )А,
А®0
А®0
(37)
де р(у* / х*) - функція щільності розподілу ймовірності символу у* за умови передачі символу х*, А - довільне мале позитивне число, яким можна знехтувати при обчисленні. Так як Р(у* /и*) пропорційна р(у* / х*), можна записати таке:
Р( у* / ик) = Р( уї, у*2*’р,/3р,-/С , С*2, ) = Р( уї, у**’р, у*’р,■•■/х^, х*2, р, х*р ,-••) =
*
.
2
Є
*
*
2
2
1=2
Так як використовується канал без пам'яті,
І ■’і’1-1')’ -X(>-’ір -хір Г (>і‘ № )2 -X (уУ)2+ х?-' Г
Р(у* / и*) ¥ Є 2—5 ^ 2— = є 2—' — "2 ‘L
'=2
2—;
Є
— ".X —
а« '=2 °и
, р . '1’У I р
у*^’. XXі у*р+х* іх*р
= В*Є
— п
— п
2—
г ]+(*г)2
- Xі—2-2— '', р і, р '', р
" у * = у* + х *
(у*,’ )2+(4,' )2 (у* 1 +х
де Вк = Є
Об'єднуючи (36) і (38), вираз для ук(5 ,5) запишеться як
* ,р
, і, р '‘>У I і ,р
у*'# .X іу* +Х* іх*
А*Є 2
Ук(5, 5) = р( ук/ ск = с) р(ик = и) = ВкЄ
Ак і Вк можуть бути опущені при обчисленні. Для каналу з АБГШ
(39)
N р-іїп
(40)
де А02 - відношення енергії біта до спектральної щільності потужності шуму, Ес - енергія
кодового символу, Еь
_ р ‘ 1
енергія некодованого біта, ЕС = ЯЕь
-2 = N0 = _і_
п 2 2-Я-Л? 2-А
2 , р А —, припустимо, що Ес = 1.
0 "0
Замінюючи вираз для аП в виразі (39), отримаємо
У* (’', Л) = В*-Є
1’ і’ 9 Іу'>р+г*'>р\- г'^-р
- $ ,уІу +х* )х* — ^ -2
—п і—2 —п
1’ 1’ 9 їй'#+г*'р- гl’P
уі -хі +^\ук-+Г±_ГЇ_
-4*-Є
2 _
4*- В*-Є
—
+£
-п
-є
Еа (сі ) С1
----------С*
2 * _
= 4*‘в*ехР
у*’-4-СІ + £(у*,р + х-*р)-4-с*
l=2
2
2
*
2
+ “ -Еа (С*)- С*
= 4*-в*-ехр
1 т ' 1\ „1 , 4- 4 - Л0 1 ,,1,’ _1
“-Еа(Сі)- С* +
2 р
у* С*
" 9 (
ехр X
_ l=’ і
4-4-А2 1
р
Т-(у*-р + р)-С*
2
2-4(с*)-С* + Ес-і-у*”-С*
ехр
-і-у- + хГ)-
X і Ес- 2- (у*,р+хТ
(41)
де
к
4 - A - h2
Lc A --------— , ge(s',s) A exp
Остаточно,
p
У (Lc- 1 - (vr + Xf) - cl j
L(uk ) = log[
У а k-i(s,)y і s)p-(s)
u+_____________________________
У а і-i(s,)y і s)p-(s)
],
(42)
У ак-1(s')ук(s',s) r
P , 4 "lv ’ ' ~ [і, якщо s = 1
а k (s) = *„ P---------, а0(s) = l
У У аі-і(s 1 її(s',s) І0, якщо s *1
p--i(s,) =
У Pi (s)yі s)
_____s___________________
У У аі-’(s,)yі-iC^s)
, PN (s) =
1, якщо s =1 I 0, якщо s Ф 1
g(s', s) ^ exp
i- La (ck) -ck + Lc-’-y.,s-ck
exp
У f Lc-(vp + X*lp)-cl
i =2
Використовуючи правило Байєса, ЛВФП для MAP декодера може бути записано як
L(ui)=iog
P(uk = 1/yN)
P(uk = 0/yN)
= log
P( yN / Щ = 1)
P(yN / ui = 0)
+ log
P(ui =k) P(ui = 0)
(43)
ик - випадкова змінна, яка з рівною ймовірністю набуває значень «1» чи «0», тобто р(ик = 1) = р(ик = 0) у випадку ФМ-2. У турбодекодері декодер 1 приймає послідовність
І, s *2,p
/15 *2Р * кР \ • • / |1 5 * ,2 р * \к Р \
виду ук = (ук , у к у к >, а декодер 2 - послідовність ук = (у к , у к, — у к ). Декодер 1 приймає інформацію з декодера 2, яка є для нього апріорною. З (27) і (41) маємо
_ Xак-1(5 'Уїк(5\ 5)Рк(5)
Е(ик) = І08
У а і-i(s ')Yі(s', s)p і(s)
= log-
Уа і-i(s,) -l~k (s)-exp u+ ’ La (cl)-ck + ’ Lc - vk‘-cl exp У [ Lc-(vr + X^)- ck j
Уа--1(s,) -f~k (s)-exp u~ 1 La (c'i> ck + ’ Lc - y'k‘-c't exp У ^ Lc-(y,p + X-Г)- ck j
(44)
Так як Ьа(с*) и Ьс - у1*’ не залежать від ик і стану решітчастої діаграми РСЗК, вони можуть бути виключені з виразу (44). З урахуванням цього виразу, (44) перепишеться таким чином:
Xа*-і(’')- Ь* (’)- У*. (’’, ’)
Ци*) = Ц(с* ) + Ьс - у'*']+ Іов £ ~ ~ , , , (45)
'Уа-i(s')- Рк (s)- у- (s', s)’
u
s s
ss
u
u
де уе(з',з) А ехр
£ [ і,-і-(у,р+Лр )-
У
Iе(сі) - апріорна інформація про біт и;,
розрахована попереднім декодером, іс - - канальна інформація про прийнятий сим-
вол, третя частина цього виразу представляє зовнішню інформацію, яка може бути передана на наступний декодер як апріорна.
Розглянемо особливості обчислення вихідного ЛВФП для декодера і, використовуючи алгоритм декодування Мар.
Позначимо попередній стан діаграми 8І-Х = з', а поточний - = з . На рис. 4 пока-
заний приклад діаграми для РСЗК виду (1, 7/5). Для кожного переходу діаграми обчислюється перехідна рекурсія уі (з', з) зі стану / в стан з за формулою
У(^ з) ~ ехР
-(і (Xе) +1, у)+1с -£ (у1Ґ + х‘Ґ) - X-;"
і=2
(46)
де хгс, хір, і є (2, д) - відповідно систематичний символ кодера ТК і перевірочні символи РСЗК до проходження каналу з шумами; ус, у‘(,р, і є (2, д) - систематичний символ кодера ТК і перевірочні символи РСЗК після проходження каналу з шумами; Іа (х^) - апріорна інформація і -го декодера; і, - параметр канальної “надійності”; д - кількість перевірочних символів РСЗК, х і , і є (2, д) -перевірочні символи РСЗК, які використовуються як додаткова інформація при розрахунку ЛВФП про переданий біт.
Для кожного стану діаграми визначається пряма рекурсія а{(з) (починаючи з початку блока - при прямому обчисленні) за формулою:
а і (*) = £ а і-і (з 0 • у і (з ^ з), де а і-і (з')
з'
- нормована пряма рекурсія і -1 такту. Підсумовування виконується по
/
всіх попередніх станах з , для яких існує перехід у поточний стан з .
5 =
5 =
з =
Рис. 4. Діаграма РСЗК виду (1, 7/5)
Одержавши прямі рекурсії для всіх станів даного такту, обчислюється параметр нормалізації, що визначається сумою даних рекурсій:
X а* (5) = XX ~ *-і (5') • у 5).
5 5 /
З урахуванням параметра нормалізації, вираз прямої рекурсії прийме вид
Xа *-і(5 ')• У*(5', 5)
~ *(5) = ^ ^ ~ , ^-----Т~' 7 .
XX ~ *-і( 5)^ у *(^ 5)
(47)
(48)
Для кожного стану діаграми визначається зворотна рекурсія Р?-1(з/), починаючи з
кінця блока при зворотному обчисленні, за формулою:
с
;
з з
Pt-iOO = Z ~t(s)^ gt(s/, s),
де Ь (5) - нормована зворотна рекурсія такту *. Підсумовування відбувається по всіх по... , точних станах 5 , в які існує перехід з попереднього стану 5 .
З урахуванням параметра нормалізації, отриманого при обчисленні прямої рекурсії,
вираз зворотної рекурсії прийме вид
Z p t(s )^ gt(s', s)
s__________________________
ZZ a t-i( s ')• gt(s', s).
(50)
ЛВФП на виході декодера i записується таким виразом:
Z a t-i(s,)- ~t(s)- gt (s', s)
(s',s )
L (xC) = log ut=+.
Z a t-i(s,)- pt(s)^ gt (s\s)
(s',s )
ut =-1
Z a t-i( s ')• ~t(s )exp
= log
.• xC ■L (xC)+Lc y)
exp
L.^ ’• Z У•p + x• t•p )• x;•p
2 i = 2
Z a t-i( s ')• ~t(s )exp
• (La (xC ) + Lc • yC )
exp
L
,1,^ (,,;p + x Г) •
t = 2
Таким чином, перевірочні символи x t , i є (2, q), t є (1, N), які вироблені кодером
ітерації декодування i -1, передаються до відповідного декодера ітерації i і використовуються для підвищення надійності прийняття рішень при розрахунку ЛВФП про переданий біт.
Для підтвердження ефективності запропонованого методу підвищення достовірності передачі інформації проведено імітаційне моделювання характеристик
достовірності передачі інформації в каналі з адитивним білим гаусівським шумом та організованими завадами.
Використовуючи програмний продукт Borland C++ Builder 6.0, моделювались такі елементи системи передачі інформації: джерело псевдовипадкової послідовності,
двокомпонентний кодер ТК, модулятор ФМ-2, канал з адитивним білим гаусівським шумом та організованими завадами, демодулятор ФМ-
2, модифікований згідно з запропонованим методом двокомпонентний декодер ТК з алгоритмом декодування Map, модуль розрахунку середньої ймовірності бітової помилки декодування.
На рис. 5 показано результати моделювання характеристик достовірності передачі інформації при впливі в каналі адитивного
h’, дБ
Рис. 5. Результати моделювання характеристик достовірності передачі інформації
ss
і
C
г, p
x
x
2
білого гаусівського шуму (відношення сигнал/шум h0’ = Еъ / N0 = 12 дБ) та організованої широкополосної шумової завади (з різними значеннями відношення сигнал-завада h2’ = Еъ / Nj) без використання та з використанням розробленого методу. При моделюванні використовувався ТК з поліноміальними генераторами g0 = 7, g. = 5 , псевдовипадковим перемежувачем, швидкістю кодування R = 1/3, кількістю біт у блоці N = 1000, 8 ітераціями декодування, алгоритмом декодування Map.
Результати моделювання свідчать, що додаткове використання перевірочних символів x*'tp, і є (2, q), t є (1, N), які виробляються кодерами на кожній ітерації декодування та передаються на декодери наступної ітерації при розрахунку ЛВФП про передані біти, сприяють підвищенню достовірності передачі інформації за рахунок підвищення надійності прийняття рішень. Так, при відношенні сигнал-завада в каналі h’ = 0,5 дБ використання запропонованого методу зменшує середню ймовірність бітової помилки декодування з PB = 6 10-3 до PB = 10-3, тобто кількість помилок зменшується в 6 разів. При зміні відношення сигнал-завада в каналі від h’ = 0 дБ до h’ = 2 дБ кількість помилок зменшується в 4...6 разів.
З. Висновки
У статті запропоновано метод підвищення достовірності передачі інформації в розподілених програмованих безпровідних системах, сутність якого полягає в застосуванні для підвищення надійності прийняття рішень про декодований біт додаткової інформації, отриманої шляхом кодування жорстких рішень з виходу кожного компонентного декодера турбокоду.
Отримані нові аналітичні співвідношення в алгоритмі декодування турбокодів по максимуму апостеріорної ймовірності, які враховують додаткову інформацію при розрахунку логарифмічних відношень функцій правдоподібності про передані біти.
Результати моделювання характеристик достовірності передачі інформації в розподілених програмованих безпровідних системах свідчать, що запропонований метод дозволяє зменшити кількість помилок при декодуванні в 4...6 разів.
СПИСОК ЛІТЕРАТУРИ
1. Сорохтин Е.М. Распределенные программно-определяемые радиосистемы / Е.М. Сорохтин, С.А. Минеев // Вестник Нижегородского университета им. Н.И. Лобачевского. - 2010. - № 5 (2). -С. 383 - 388.
2. Burns P. Software Defined Radio for 3G / Burns P. - Boston: Artech House, 2003. - 279 p.
3. Помехозащищенность систем радиосвязи с расширением спектра сигналов методом псевдослучайной перестройки рабочей частоты / [В .И. Борисов, В.М. Зинчук, А.Е. Лимарев и др.]. - М.: Радио и связь, 2000. - 384 c.
4. Berrou C. Near Shannon limit error-correcting coding and decoding: turbo-codes / C. Berrou, A. Gla-vieux, P. Thitimajshima // Proc. Int. Conf. оп Commun., ICC-93. - Geneva, Switzerland, 1993. - May. -P.1064 - 1070.
5. Скляр Б. Цифровая связь. Теоретические основы и практическое применение / Скляр Б. - [2-е изд.]. - М.: Вильямс, 2003. - 1104 с.
6. Woodard J. Comparative Study of Turbo Decoding Techniques: An Overview / J. Woodard, L. Hanzo // IEEE Transactions on Vehicular Technology. - 2000. - Vol. 49, N 6. - P. 2208 - 2232.
7. Berrou C. Near optimum error correcting coding and decoding: turbo-codes / C. Berrou, A. Glavieux // IEEE Trans. on Commun. - 1996. - Vol. 44, N 10. - P. 1261 - 1271.
Стаття надійшла до редакції 22.03.2013
