VYSOKÉ UČENÍ TECHNICKÉ V BRNĚ BRNO UNIVERSITY OF TECHNOLOGY
FAKULTA ELEKTROTECHNIKY A KOMUNIKAČNÍCH TECHNOLOGIÍ ÚSTAV RADIOELEKTRONIKY FACULTY OF ELECTRICAL ENGINEERING AND COMMUNICATION DEPARTMENT OF RADIO ELECTRONICS
NOVÉ METODY KANÁLOVÉHO KÓDOVÁNÍ PRO DRUŽICOVOU KOMUNIKACI NEW CHANNEL CODING METHODS FOR SATELLITE COMMUNICATION
DISERTAČNÍ PRÁCE DISSERTATION THESIS
AUTOR PRÁCE
Ing. Jaroslav Rumánek
AUTHOR
VEDOUCÍ PRÁCE SUPERVISOR
BRNO, 2009
Ing. Jiří Šebesta, Ph.D.
Abstrakt Disertační práce se zabývá problematikou nových progresivních metod kanálového kódování při přenosu datových zpráv prostřednictvím družicového transpondéru. Součástí práce je návrh systému pro přenos krátkých textových zpráv s aplikací nových turbo kódů s požadavkem na co nejnižší vysílací výkon a co nejmenší aperturu na straně uživatele. Cílem disertační práce je navrhnout systém, který bude schopen při velmi nízkých poměrech signál/šum umožnit bezchybný přenos SMS zprávy. Práce se zaměřuje především na problematiku energetické bilance, modifikaci a implementaci nových typů turbo kódů a využití jejich unikátních vlastností, na vývoj metod pro odhad chybovosti a metod pro stanovení výsledné SMS zprávy při nenulové chybovosti. Hlavním přínosem je vývoj nového typu turbo kódu, který zajišťuje optimální vlastnosti pro dané použití, dále pak vývoj nové metody pro odhad chybovosti a metody, která na základě navržené struktury rámce a teorie o turbo dekódování umožní stanovení výsledné bezchybné podoby SMS zprávy i v případě nenulové chybovosti.
Klíčová slova SMS, turbo kód, kódování, prokládání, MAP, dekódování, bitová chybovost, družice, energetická bilance, odhad bitové chybovosti, PHASE 3E, družicová komunikace
Abstract This dissertation thesis deals with new progressive channel coding methods for data transmission using satellite transponder. The design of the system for SMS transmission, in which novel turbo coding methods are applied, is discussed too. An achievement of the lowest output power and the smallest user aperture is the principal aim of the new method applications. Design of system that would be able to the error free SMS transmission by very low signal to noise ratio is analyzed in this dissertation thesis. The work is focused on energy budget, modification and implementation of new turbo code types and using unique properties, development of new bit error rate estimation methods and methods for determination of final SMS form. The main contribution is the new type of turbo code development that have optimal properties for this usage, development of new bit error rate estimation method and development of method that is able to determine final form of SMS on the basis SMS frame structure and turbo decoding theory when the bit error rate is not zero.
Keywords SMS, turbo code, coding, interleaving, MAP, decoding, bit error rate, satellite, energy budget, bit error rate estimation, PHASE 3E, satellite communication
Prohlášení Prohlašuji, že svou disertační práci na téma Nové metody kanálového kódování pro družicovou komunikaci jsem vypracoval samostatně pod vedením školitele a s použitím odborné literatury a dalších informačních zdrojů, které jsou všechny citovány v práci a uvedeny v seznamu literatury na konci práce. Jako autor uvedené disertační práce dále prohlašuji, že v souvislosti s vytvořením této disertační práce jsem neporušil autorská práva třetích osob, zejména jsem nezasáhl nedovoleným způsobem do cizích autorských práv osobnostních a jsem si plně vědom následků porušení ustanovení § 11 a následujících autorského zákona č. 121/2000 Sb., včetně možných trestněprávních důsledků vyplývajících z ustanovení § 152 trestního zákona č. 140/1961 Sb. V Brně dne 31. srpna 2009
__________________________________ autor práce
Bibliografická citace RUMÁNEK, J. Nové metody kanálového kódování pro družicovou komunikaci. Disertační práce. Brno: FEKT VUT v Brně, 2009, 88 stran.
Poděkování Děkuji vedoucímu disertační práce Ing. Jiřímu Šebestovi, Ph.D. za účinnou metodickou, pedagogickou a odbornou pomoc a další cenné rady při zpracování mé disertační práce.
V Brně dne 31. srpna 2009
__________________________________ autor práce
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Obsah 1
ÚVOD
3
1.1
Energetická bilance
4
1.2
Kapacita kanálu
4
1.3
Dopplerův posuv
5
2
PŘENOSOVÉ RÁDIOVÉ KANÁLY A DETEKCE SYMBOLŮ
7
2.1
Model kanálu a detekce symbolů
7
2.2
Rayleighův a Riceův útlumový kanál
9
3
KÓDOVÁNÍ
10
3.1 Blokové kódy 3.1.1 Konvoluční kódy 3.1.2 Zřetězené kódy
11 11 13
3.2 Turbo kódy 3.2.1 Paralelně zřetězené konvoluční (turbo) kódy 3.2.2 Maximum a posteriori a maximální pravděpodobnost 3.2.3 Dekódování Turbo kódů 3.2.4 MAP dekódovací algoritmus pro konvoluční turbo kódy
14 14 16 17 20
4
CÍLE DISERTACE
26
4.1 Analýza progresivních kanálových kódů a jejich optimalizace pro návrh systému pro přenos SMS zpráv 26 4.1.1 Výpočet energetické bilance 26 4.1.2 Analýza a modifikace Turbo kódů 27 4.1.3 Metodika pro odhad chybovosti BER 27 4.1.4 Stanovení struktury rámce SMS zprávy 27 4.1.5 Metodika pro stanovení symbolů výsledné SMS zprávy 27 4.1.6 Synchronizace 27
5
ENERGETICKÁ BILANCE
29
5.1 Výpočet energetické bilance uživatel – družice – řídící segment 5.1.1 UPLINK
29 30
5.2
34
6
Výpočet energetické bilance řídící segment – družice – uživatel
ANALÝZA TURBO KÓDŮ
36 1
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí 6.1 Výběr vhodného typu turbo kódu 6.1.1 Turbo kód pro telemetrii a družicovou komunikaci
36 36
6.2 Prokladač 6.2.1 Telemetrický prokladač 6.2.2 Pseudonáhodný prokladač 6.2.3 Prokladač s eliminací tzv. „edge efektu”
37 38 39 39
6.3
40
Implementace turbo kódu v prostředí MATLAB
6.4 Analýza kódového poměru R a délky bloku dat k 6.4.1 Kódový poměr R = 1/2 6.4.2 Kódový poměr R = 1/3 6.4.3 Kódový poměr R = 1/4 6.4.4 Kódový poměr R = 1/6 6.4.5 Zhodnocení
7
ODHAD CHYBOVOSTI ZPRÁVY
42 42 44 46 48 50
55
7.1
Odhad chybovosti BER na základě APP
55
7.2
Modifikovaná metoda pro odhad chybovosti na základě LLR
59
7.3
Srovnání metod pro odhad chybovosti BER
64
8
NÁVRH RÁMCE SMS ZPRÁVY
65
8.1 Metodika stanovení hodnot jednotlivých symbolů přenášené datové posloupnosti 8.1.1 Pravděpodobnostní metoda pro výpočet hodnot symbolů datové posloupnosti 8.1.2 Zjednodušená metoda pro výpočet hodnoty bitů 8.1.3 Výsledky metod pro stanovení hodnot symbolů 8.1.4 Zhodnocení metod pro stanovení hodnot symbolů
9
SYNCHRONIZACE
66 66 68 70 73
74
9.1
Goldovy kódy
74
9.2
Aplikace Goldových kódů
75
10
ZÁVĚRY
82
11
LITERATURA
84
2
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
1 Úvod Ve 20. století došlo k velkému průlomu bezdrátové komunikace prostřednictvím elektromagnetických vln, která je využívána pro přenos zvuku, videa a nejrůznějších datových informací. Největšího rozmachu došlo v druhé polovině 20. století. Základy moderní digitální komunikace stanovil Shannon v roce 1948, ve své fundamentální práci „Matematické základy komunikace“ [ 1 ]. Přestože od té doby uplynula dlouhá doba, jsou tyto základní informace stále aktuální. Schéma rádiového komunikačního kanálu, které navrhnul Shannon, je zobrazeno na Obr. 1.1. Toto schéma poskytuje zjednodušený pohled na způsob přenosu signálu. Reálný digitální rádiový kanál navíc obsahuje několik dalších bloků, které umožňují zlepšení přenosových vlastností [ 2 ]. rádiový kanál
šum AWGN, interference, ...
modulátor
nosná vlna
demodulátor
kodér kanálu
ochrana přenosu (FEC)
dekodér kanálu
datová kodér zdroje
zdroj signálu
komprese A/D
dekomprese D/A
přeměna typu signálu
dekodér zdroje
koncový stupeň
Obr. 1.1: Obecné schéma Shannonova rádiového komunikačního systému [ 1 ]
Disertační práce se zabývá problematikou družicové komunikace, která je velmi perspektivním odvětvím pro rozvoj mobilních komunikací. Hlavním cílem je analýza efektivních metod kanálového kódování, které umožňují spolehlivé dekódování i u signálů s velmi nízkými poměry Eb/N0. Jedná se především o turbo kódy, které mají díky navržené struktuře kodéru i dekódéru tyto unikátní vlastnosti.
3
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Při návrhu družicových spojů je potřeba klást především velký důraz na energetickou bilanci spoje. Na rozdíl od pozemních systémů je energetická bilance spoje Země – družice a zpět velmi složitou úlohou. Družicový spoj nelze jednoduše předimenzovat, protože by se to výrazně projevilo na nákladech spojených s konstrukcí družice a jejím vynesením na oběžnou dráhu, ale také změnou tepelných podmínek v těle družice [ 4 ].
1.1 Energetická bilance Energetickou bilanci spoje lze vyjádřit pomocí komunikační rovnice [ 4 ], která charakterizuje výkon nosné vlny vůči spektrální hustotě výkonu šumu v rovině konektoru přijímače. V případě, že nedochází k interferenci ani jinému rušení spoje, lze šumové pozadí signálu považovat za aditivní bílý šum s normálním rozložením (AWGN). Spektrální hustotu šumu N0 lze vyjádřit součinem k.TSRX, kde k je Boltzmanova konstanta a TSRX je ekvivalentní šumová teplota systému. Komunikační rovnice bude mít tvar (1.1)
1.2 Kapacita kanálu Kapacita kanálu C vyjadřuje množství informace, které je možné přenést kanálem za určitý čas [ 1 ]. Při přenosu působí na signál různé rušivé jevy, které zapříčiní znehodnocení signálu a s tím spojené navýšení chybovosti.
Obr. 1.2: Závislost normované přenosové kapacity (C0/B0) a poměru signál/šum (S/N) idealizovaného rádiového komunikačního systému na normované šířce pásma (B/B0) [ 1 ]
4
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Kanálovou kapacitu lze vypočítat na základě Shannon-Hartleyova vztahu (1.2) [ 1 ]. Vztah vyjadřuje maximální množství informace vyjádřené v bitech, které lze přenést systémem za jednu sekundu při nulové chybovosi BER v případě, že na systém působí gaussovský bílý šum. (1.2) Kapacitou kanálu tedy rozumíme teoreticky nejvyšší dosažitelnou bitovou rychlost idealizovanéo kanálu. V reálném systému se však k této hodnotě můžeme jen přiblížit. Požadovanou kapacitu C lze získat vhodnou kombinací parametrů B, S a N. Vyjádříme-li výkon šumu N jako součin spektrální výkonové šumové hustoty N0 a šířky pásma B tedy N = N0B a zavede-li se poměrná šířka pásma B0 = S/N0, lze vtah přepsat do normovaného tvaru (1.3) Z rovnice (1.3) je vidět, že normovaná přenosová kapacita C/B0 je funkcí jedné proměnné, kterou je poměrná šířka pásma B/B0. V grafu je zobrazena závislost poměrné šířky pásma a poměru signál/šum na poměrné šířce pásma B/B0. Z rovnosti vztahů (1.2) a (1.3) pak vyplývá (1.4) Z grafu (Obr. 1.2) je vidět, že normovaná šířka pásma je důležitým mezním bodem, kde se výkon signálu rovná výkonu šumu. Tímto bodem prochází přímka, která dělí graf na dvě části. Levá část grafu vystihuje klasické systémy, kde je úroveň signálu větší než úroveň šumu tzn. poměr signál/šum je větší než 1. V právé části jsou naopak signály, kde je signál pod úrovní šumu a hodnota poměru signál/šum je menší než 1. Přenosová kapacita se u těchto systémů blíží maximální hodnotě [ 1 ], [ 2 ].
1.3 Dopplerův posuv Při přenosu signálu Země-družice a zpět nesmíme opomenout Dopplerův posuv frekvence (DP). Vzniká vlivem rychlého pohybu družice vůči pozemskému pozorovateli a uplatňuje se v různé míře [ 4 ]. (1.5) kde
je úhel mezi vektorem družice f
je frekvence,
c
je rychlost šíření světla, je složka vektoru k pozorovateli. 5
a směrem k pozorovateli,
rychlosti
družice
do
směru
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Dopplerův posuv je odlišný pro družice na různých oběžných drahách (LEO, ICO, GEO). Na dráze GEO (GEostacionary Orbit – 35786 km) je Dopplerův posuv téměř nulový a umožňuje tak výborné využití spektra. Na nižších dráhách ICO (Intermediate Circular Orbit – 5000 15000 km) je již Dopplerův posuv nezanedbatelný. Dosahuje hodnoty cca 6 ppm ve vztahu k frekvenci zdroje signálu, který je dán větším relativním pohybem satelitu. Největší Dopplerův posuv však vzniká na drahách LEO (Low Earth Orbit < 2000 km). Na těchto drahách jsou nejmenší ztráty a lze použít menší EIRP na straně satelitu i pozemní stanice, ale kvůli velké relativní rychlosti mezi satelitem a pozorovatelem vzniká největší Dopplerův posuv frekvence (až 25 ppm – vztaženo k frekvenci zdroje signálu). Jak je vidět z rovnice (1.5) Dopplerův posuv roste přímo úměrně se vzrůstající frekvencí. To je také důvod, proč družice na drahách LEO pracují na co nejnižších kmitočtech. Důsledkem velkého Dopplerova posuvu dochází k neefektivnímu využití spektra a obtížné demodulaci signálu. V okamžiku, kdy se družice přibližuje, prudce vzrůstá směrnice tečny, kterou lze přibližně spočítat z grafu na Obr. 1.1 podle následujícího vztahu (1.6) Aby bylo možné korigovat Dopplerův posuv a zamezit tak vzniku velkých ztrát, je nutné v reálném čase predikovat aktuální polohu družice a relativní rychlost družice vůči pozemní stanici postupně přelaďovat pracovní kmitočet přijímače. K výpočtu aktuální polohy družice se používají Keplerovy elementy obsahující výchozí polohu družice a popisující její dráhu. Nejrozšířeněj ší je formát NASA tzv. dvouřádkový 2-LINE (viz. Tab. 7.3). Součástí bývá i dekódovací klíč. Tato data jsou publikována korporací NORAD [ 4 ]. Tab. 1.1: Ukázka Keplerových elementů DECODE 2-LINE ELSETS WITH THE FOLLOWING KEY: 1 AAAAAU 00 0 0 BBBBB.BBBBBBBB .CCCCCCCC 00000-0 00000-0 0 DDDZ 2 AAAAA EEE.EEEE FFF.FFFF GGGGGGG HHH.HHHH III.IIII JJ.JJJJJJJJKKKKKZ KEY: A-CATALOGNUM B-EPOCHTIME C-DECAY D-ELSETNUM E-INCLINATION FRAAN G-ECCENTRICITY H-ARGPERIGEE I-MNANOM J-MNMOTION K-ORBITNUM ZCHECKSUM TO ALL RADIO AMATEURS BT AO-07 1 07530U 74089B 07271.17170329 -.00000028 00000-0 10000-3 0 9872 2 07530 101.5057 306.1161 0012172 055.5188 304.7029 12.53572559504029
6
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
2 Přenosové rádiové kanály a detekce symbolů 2.1 Model kanálu a detekce symbolů Kanál může být charakterizován vstupní abecedou, výstupní abecedou a souborem pravděpodobností, které popisují přechody mezi vstupem a výstupem. Komunikační systémy obvykle poskytují diskrétní abecedu kanálu. Výstupy kanálu definují typ kanálu a jsou generovány demodulátorem, který rozhoduje o tom, který symbol byl přenesen. Diskrétní bezpaměťový kanál (DMC – discrete memoryless channel) je charakterizován diskrétní vstupní abecedou, diskrétní výstupní abecedou a sadou podmíněných pravděpodobností P(j|i), která rozhoduje o hodnotě přijatého j za podmínky, že bylo vysláno i. Vstupní abeceda obsahuje M symbolů (1 ≤ j ≤ M) a výstupní abeceda obsahuje Q symbolů (1 ≤ j ≤ Q). Binární symetrický kanál je speciální případ DMC, ve kterém je vstupní a výstupní znaková abeceda složena z binárních prvků (buď 1, nebo 0, a proto M = Q = 2) a podmíněné pravděpodobnosti jsou symetrické. Toto chování je určeno dvojicí pravděpodobností: (2.1) (2.2) Rovnice (2.1) popisuje pravděpodobnost, kdy bylo rozhodnuto, že byla vyslána „1“ a přijata „0“ a naopak. Proto p představuje pravděpodobnost chyby a druhá rovnice popisuje pravděpodobnost, kdy bylo rozhodnuto, že byl přijat správný symbol. Tato pravděpodobnost je dána výrazem 1 - p. Výstupem demodulátoru binárního kanálu je tedy buď logická „0“ nebo „1“. Tento způsob rozhodování je nazýván tvrdé rozhodování (hard-decision) [ 13 ]. Nedůležitější je model kanálu, kde je vstupní abeceda diskrétní, ale výstupní abeceda spojitá v rozsahu [-∞, +∞]. Příkladem je Gaussovský kanál, ve kterém je rušení způsobeno aditivním bílým šumem (AWGN). Model AWGN je používán ve velkém množství komunikačních systémů a je popsán uvedeným matematickým vztahem (rozložení hustoty pravděpodobnosti), který vyjadřuje vztah mezi hodnotou přijaté proměnné y a vyslaného symbolu uk. .
7
(2.3)
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Symbol 2 představuje disperzi AWGN. Přestože AWGN obsahuje všechny spektrální složky, je toto rozložení konečné pro filtrovaný kanál ( 2 = N0/2). Pokud budeme mít např. AWGN kanál s BPSK (Binary Phase-Shift Keying) modulací, jsou bity zobrazeny jako jednorozměrná křivka (viz. Obr. 2.1). Na základě vztahu (2.3) lze napsat dvě funkce vyjadřující hustotu rozložení podmíněné pravděpodobnosti. První z nich vyjadřuje hodnotu y za podmínky, že byla vyslána „0“ a druhá pak vyjadřuje hodnotu y za podmínky, že byla vyslána „1“. Jsou nazývány pravděpodobnostní funkce přenášených symbolů. (2.4)
(2.5) Na Obr. 2.1 jsou zobrazeny dvě pravděpodobnostní funkce. Šedá oblast vyznačuje pravděpodobnost, kdy může dojít k chybě při vyhodnocování symbolu. Levá část šedé oblasti vyjadřuje pravděpodobnost chybného příjmu v případě, že byl vyslán symbol „1“ a pravá část vyjadřuje pravděpodobnost chybného příjmu v případě, že byl vyslán symbol „0“. Z uvedeného obrázku je vidět, že body s nejvyšší hodnotou hustoty pravděpodobnosti jsou -1 pro p(y|0) a +1 pro p(y|1). Tyto hodnoty odpovídají proměnným a1 a a2 v rovnicích (2.4) a (2.5) a jsou výstupy BPSK modulátoru na straně vysílače. Jestliže je takový signál zpracováván dekodérem na straně přijímače, nazýváme tento proces měkké dekódování (soft-decision) [ 13 ]. rozhodovací úroveň p(y|0)
-4
-3
-2
p(y|1)
-1
0
1
2
3
4
Obr. 2.1: Graf funkcí hustoty rozložení podmíněné pravděpodobnosti pro AWGN kanál, BPSK modulace, s = 1
8
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
2.2 Rayleighův a Riceův útlumový kanál Nejjednodušší pravděpodobnostní model pro koeficienty kanálového FIR filtru je založen na předpokladu, že existuje velký počet statistický nezávislých odražených a rozptýlených cest s náhodnými amplitudami ve zpožďovacím okně odpovídající jednomu koeficientu. Fáze i-té cesty je 2πfc i modulo 2π. Proto fc/ i = di/λ, kde di je délka i-té cesty a λ je vlnová délka nosné. Jestliže jsou body odrazu a rozptylu relativně daleko tj. di >> λ, lze předpokládat, že je fáze pro každou cestu rovnoměrně rozdělena od 0 do 2π, a že jsou fáze různých cest nezávislé. Příspěvek každé cesty pro stanovení násobících koeficientů kanálového FIR filtru jsou (2.6) To může být modelováno jako cirkulární symetrická proměnná [ 21 ]. Každý koeficient je součtem velkých čísel těchto cirkulárních symetrických náhodných proměnných. Z toho vyplývá, že je suma mnoha malých reálných nezávislých náhodných proměnných, a proto lze na základě centrálního limitního teorému modelovat tuto proměnnou jako gaussovskou náhodnou proměnou [ 21 ]. Podobně je také fáze uniformní, je gaussián se stejným rozptylem pro fixní hodnotu . To zaručuje, že je skutečně cirkulárně symetrická . Předpokládáme tedy, že rozptyl je funkcí koeficientu a je časově nezávislý. Jestliže tedy uvažujeme Gaussovskou hustotu pravděpodobnosti, je velikost koeficientu Rayleighyova útlumu náhodná proměnná s hustotou (2.7) a hodnota
je exponenciální s hustotou rozložení (2.8)
Tento model je nazýván Rayleighův kanál, který lze použít v systémech, kde se vyskytuje mnoho cest signálu s odrazy bez přímého paprsku [ 21 ]. Existuje také velmi často používaný alternativní model, ve kterém je přímý paprsek se známou amplitudou, a proto existuje velký počet nezávislých cest. V tomto případě lze nejméně pro jednu hodnotu modelovat jako (2.9) První část výrazu odpovídá přímě cestě s uniformní fází θ a druhá část výrazu odpovídá seskupení velkého počtu odražených a rozptýlených cest nezávislých na θ. Parametr (tzv. -faktor) je poměr energie v přímě cestě k energii v odražených cestách. Čím větší je tím více je kanál deterministický. Hodnota této náhodné proměnné je nazývána jako Ricianovo rozdělení a hovoří se o Riceově kanálu [ 21 ].
9
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
3 Kódování Kódování je proces, který je využíván ve všech komunikačních systém. Lze na něj pohlížet z různých hledisek podle toho, k čemu je kódování v daném systému předurčeno [ 7 ]: kódování pro přizpůsobení informačního zdroje na kanál, kódy pro snížení nadbytečnosti, kódy pro zabezpečení proti chybovosti, kódy pro zabezpečení proti zcizení – šifrování, kódování pro zrovnoměrnění spektra přenášeného signálu – skrambling, aj. Disertační práce se zaměřuje především na problematiku kódování z hlediska zabezpečení proti chybovosti vznikající v důsledku průchodu signálu komunikačním kanálem. Obecně jsou tyto kódy nazývány kanálové kódy. Na Obr. 3.1 je zobrazeno dělení kanálových kódů. Kanálové kódy
detekční
korekční
konvoluční kódy
binární kódy
blokové kódy
binární kódy
lineární kódy
cyklické kódy
necyklické kódy
řetězcové kódy paralelní kódy
turbo kódy iterativní kódy
Obr. 3.1: Rozdělení kanálových kódů
10
nelineární kódy
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
3.1 Blokové kódy Blokové kódy jsou charakterizovány fixní délkou kódového slova n. Tento počet n prvků lze rozdělit na k informačních bitů a r = n - k zabezpečovacích (paritních) bitů. Jedná se tedy o (n, k) kódy s kódovým poměrem n/k, který vyjadřuje množství redundantních bitů přidaných za informační bity. Základní myšlenkou blokových kódů je vložení takového počtu paritních bitů, které zajistí dostatečné množství informace k pozdější detekci a opravě chyb na přijímací straně [ 7 ]. Blokové kódy lze rozdělit do dvou skupin: Lineární kódy – u těchto kódů lze odvodit libovolnou kódovou kombinaci z ostatních kódových kombinací. Je to díky tomu, že matematický základ využívá lineární algebru. Tyto kódy jsou definovány vytvářecí maticí [G] nebo v případě cyklických kódů vytvářecím mnohočlenem G(x). Dále se tyto kódy dělí na lineární systematické kódy, které mají nejprve řazeny informační bity a poté teprve zabezpečovací bity, a na ostatní lineární kódy, které mají zabezpečovací bity umístěny mezi informační bity podle určitých pravidel. Nelineární kódy – do této skupiny patří všechny kódy, které nemají vlastnosti lineárních kódů [ 9 ]. 3.1.1 Konvoluční kódy Byly poprvé uvedeny Eliasem v roce 1955 [ 16 ]. Později byly vyvinuty další účinné dekódovací postupy. Největší přínos však přinesl v roce 1967 Viterbi, jenž navrhnul a zrealizoval Viterbiho dekódovací algoritmus [ 13 ].
u
Z-1
Z-1
v
1
v
2
Z-1
Obr. 3.2: Příklad zapojení konvolučního kodéru (n, k, m) = (2, 1, 3)
Konvoluční kódování je založeno na průchodu informačních bitů posuvnými registry, kde dochází k samotnému kódování informace, jak je zobrazeno na Obr. 3.2. Z uvedeného příkladu je zřejmé, že se nejedná o kódy, které by pracovaly s kódovými slovy fixní délky. Kódování je kontinuální proces, při kterém jsou po průchodu posuvnými registry generovány výstupní bity zabezpečené posloupnosti. Důležité je si uvědomit, že v určitém čase závisí výstupní bity nejen na momentálním stavu vstupu, ale kvůli posuvným 11
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
registrům závisí i na předchozích stavech. Výstupní zabezpečená posloupnost pak obsahuje dostatečné množství informace, aby bylo možné tuto posloupnost dekódovat na přijímací straně. Každý konvoluční kodér je popsán generujícími polynomy. Polynomy z uvedeného příkladu na Obr. 3.2 mají tvar (3.1) Lze je také přepsat ve formě (3.2) Výstup je tedy založen na kovoluci vstupu u a g1, g2 (3.3) V uvedeném příkladu se jedná o nesystematický konvoluční kód, protože na výstup není dodávána originální vstupní posloupnost, ale jsou zde pouze bity zabezpečeného posloupnosti. Také zde není přítomna zpětná vazba z výstupu na vstup. Když v roce 1993 uvedl Berrou turbo kódy [ 11 ], využil při jejich realizaci rekurzivní systematické konvoluční kódy. Důvodem k použití byl především fakt, že tyto kódy dosahují mnohem lepších výsledků při nízkých poměrech signál/šum. Konvoluční kodér lze charakterizovat jako konečný stavový automat, kde obsah stavových registrů indikuje stav automatu a hodnoty výstupů závisí nejen na okamžitém stavu vstupů, ale také na momentálním stavu posuvných registrů.
u
Z-1
v
1
v
2
Z-1
Obr. 3.3: Příklad zapojení rekurzivního systematického konvolučního kodéru (n, k, m) = (2, 1, 2)
Chování kodéru lze demonstrovat pomocí stavového diagramu, jak je ukázáno na Obr. 3.4. Tento stavový diagram je nazýván trellis diagram. Konkrétně pro tento případ se může na vstupu vyskytovat 8 různých stavů, které odpovídají třístavovému kodéru. Jelikož každý interní registr může nabývat dvou možných hodnot („1“ nebo „0“), existuje 2n možných stavů. Čárkovaná čára pak odpovídá hodnotě vstupního bitu „0“ a plná čára odpovídá hodnotě vstupního bitu „1“.
12
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
výstup m = 0 (000b)
00 10 11
m = 1 (001b)
01 01
m = 2 (010b)
10
m = 3 (011b)
10 01
m = 4 (100b)
10 00
m = 5 (101b)
10 11
m = 6 (110b)
01 00
m = 7 (111b)
stavy
10 11 t=0
t=1
t=2 Vstupní bit = ‘0’ Vstupní bit = ‘1’
Obr. 3.4: Trellis diagram kódu (2, 1, 3)
3.1.2 Zřetězené kódy Tyto kódy uvedl Forney v roce 1966 [ 13 ]. Jedná se o kódy, u kterých je kaskádně řazen vnitřní a vnější kód. Představují účinnou cestu ke kompromisu mezi mírou kódového zisku a složitostí kódu. Konstrukce těchto kódů umožňuje dosažení nízké chybovosti BER. V porovnání s klasickými kódy je pomocí těchto kódů dosaženo stejného kódového zisku při menší složitosti zapojení, protože je prováděno oddělené kódování a dekódování vnějšího a vnitřního kódu. Nejpoužívanější je kombinace konvolučního kódu (vnitřní kód), který slouží pro opravu rozprostřených chyb a Reed-Solomonova kódu (vnější kód), který slouží pro opravu shluků chyb [ 13 ].
13
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
3.2 Turbo kódy Tyto kódy byly poprvé popsány v roce 1993. V současné době jsou využívány moderními komunikačními systémy především pro jejich efektivitu. Turbo kódy vykazují velký kódový zisk [ 35 ], a proto mohou zajistit nízkou bitovou chybovost BER (10 –4 až 10–7) i při malém poměru signál/šum. Nespornou výhodou je také fakt, že u turbo kódů neexistuje žádný výrazněj ší limit pro snižování chybovosti při nízkých hodnotách BER. V dekodérech je aplikován princip iterací, pomocí kterého můžeme zvýšením počtu iteračních kroků snížit chybovost. I přes řadu výhodných vlastností, mají turbo kódy své nevýhody znesnadňující jejich použití. Tou hlavní nevýhodou je vysoká latence, kterou zapříčiňuje iterační proces v dekodéru a poměrně velmi složitý způsob dekódování, jehož důvodem je složitý dekódovací algoritmus SISO (Soft Input Soft Output). Navíc nejsou přesně stanoveny optimální parametry pro kodér i dekodér a úloha pro výběr vhodných parametrů je taktéž poměrně složitá. To je zapříčiněno enormním množstvím možných kombinací a velkým počtem parametrů zahrnutých do návrhu turbo kódů. Vše závisí na konkrétní implementaci turbo kódu a kladených požadavcích. Neexistuje žádný matematický popis, kterým bychom se mohli řídit při výběru a nasledném použití turbo kódu. Jsou stanoveny určité typy algoritmů a architektury pro jednotlivá specifická použití [ 12 ], [ 14 ], [ 15 ], avšak ne vždy se dají vhodně využít. Vždy se snažíme navrhnout turbo kód, který má výhodné vlastnosti pro danou aplikaci a volíme určitý kompromis mezi složitostí turbo kódu a jeho účinností. V případě, že pro danou aplikaci nehraje roli doba, po kterou probíhá dekódovací proces, nebo máme dostatečně velký výpočetní výkon, který zajistí rychlé zpracování, je možné dosáhnout bezchybného přenosu i přes zarušený kanál. Tato myšlenka je důležitá pro realizaci systému pro přenos krátkých textových zpráv, který je uvažován. Vzhledem k tomu, že se nejedná o systém pracující v reálném čase, je možné využít turbo kódy, jako hlavní prvek pro zabezpečení proti chybovosti i při nízkých hodnotách Eb/N0. 3.2.1 Paralelně zřetězené konvoluční (turbo) kódy Jak již název napovídá, jedná se o paralelně zřetězené rekurzivní systematické konvoluční kódy (RSC). Informace je v daném okamžiku kódována prvním kodérem a zároveň po průchodu prokladačem druhým konvolučním kodérem [ 26 ]. Vstupní informace u je tedy kódována dvěma stejnými kodéry avšak v různém pořadí. Na Obr. 3.5 je zobrazen příklad zapojení turbo kodéru. Prokladač modifikuje klasickou funkci konvolučních kódů. Obvykle je prováděno kódování kontinuálně, avšak při použití prokladače probíhá kódování po blocích. Velikost prokladače definuje délku kódového slova n. V případě, že by se délka n blížila k nekonečnu, bude se chybovost při dekódování blížit nule [ 14 ]. To platí pro blokové a konvoluční kódy. Nicméně když provedeme zřetězení, jako je provedeno u turbo kódů, je možné se přiblížit kapacitě kanálu s relativně malou délkou kódového slova n (od 32 do 65536) [ 13 ]. 14
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Na Obr. 3.5 představuje symbol u vstupní nezakódovanou posloupnost, π je prokladač, u’ je nezakódovaná vstupní posloupnost po průchodu prokladačem, c1 a c2 jsou výstupní zakódované posloupnosti a s je kopie vstupní nezakódované informace, která je taktéž přenášena. Protože se přenáší i nezakódovaná vstupní posloupnost, jsou tyto kódy nazývány systematické. Výstupní posloupnost se skládá z posloupnosti s, c1 a c2. Pro příklad kodéru uvedeného na Obr. 3.5 je kódový poměr R = k/n = 1/3. Výstupní posloupnost turbo dekodéru bude mít pak tvar [s0, c10, c20, s1, c11, c21, …, sk, c1k, c2k].
s
u
Z-1
Z-1
c
1
c
2
Z-1
u' Z-1
Z-1
Z-1
Obr. 3.5: Příklad zapojení turbo kodéru
Pro RSC kodéry existují dva generující polynomy g0 a g1. Pro dopředný směr je to polynom g1 a pro zpětnou vazbu g0. Pro uvedený příklad mají polynomy následující tvar: (3.4) (3.5) Generující matice může být zapsána také v systematickém tvaru (3.6) Ve výše uvedené rovnici je I jednotková matice a P(x) je matice polynomů k x (n-k). Pro RSC kodéry je P(x) racionální lomená funkce, kde čitatel představuje dopředný směr a jmenovatel zpětnovazební směr (3.7) Výstupní posloupnost lze získat násobením vstupní posloupnosti maticí G(x). První sloupec matice vyjadřuje, že se na výstupu objevuje kopie vstupní posloupnosti. Druhý člen matice, představuje přenosovou funkci IIR filtru (kde x reprezentuje zpoždění z-1). Aby bylo možné provést násobení, lze zapsat 15
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
vstupní posloupnost v polynomickém tvaru, jak je uvedeno dále. Všechny operace sčítaní a odčítaní jsou prováděny v aritmetice modulo-2. Pro příklad předpokládejme, že chceme zakódovat vstupní posloupnost u = (11001001): u = (11001001) = 1 + x1 + x4 + x7 v=u.G v1 = (1+ x1 + x4 + x7) . 1 = 1+ x1 + x4 + x7 v2 = (1+ x1 + x4 + x7) .
systematická část
= 1 + x5 + x7 + …
v2 = (10000101…). 3.2.2 Maximum a posteriori a maximální pravděpodobnost Dekódování lze dosáhnout použitím rozhodovacího pravidla, které maximalizuje pravděpodobnost správného rozhodnutí. Existují dva druhy pravidel. Jedním z nich je maximum a posteriori (MAP) a druhým je maximální pravděpodobnost (ML – maximum likelihood) [ 15 ]. První z nich je založeno na výpočtu tzv. a posteriori pravděpodobnosti (3.8) Druhý způsob je založen na výpočtu maximální pravděpodobnosti (3.9) Dekodérem by měl pak být vybrán signál odpovídající maximu z množiny pravděpodobností, které jsou definovány uvedenými rovnicemi. Lze dokázat, že rozhodování prováděné na základě MAP kritérií je stejné jako dekódování prováděné na základě ML pravidel, jestliže jsou obě a priori pravděpodobnosti stejné resp. jestliže jsou stavy signálu {um} stejně pravděpodobné. Rovnici (3.8) lze rozšířit pomocí Bayesova pravidla (3.10) První člen ve jmenovateli odpovídá pravděpodobnosti definované rovnicí (3.9). Vezmeme-li v úvahu, že P(um) je rovna 1/M a hodnota P(y) nezávisí na přenášeném signálu, pak maximální hodnota pravděpodobnosti P(um|y) odpovídá maximální hodnotě pravděpodobnosti P(y|um). MAP kritérium je běžně používáno v dekódovacích algoritmech pro minimalizaci symbolových chyb, zatímco ML kritérium je používáno pro minimalizaci chyb kódových slov. První z nich je tedy použito pro dekódování jednotlivých symbolů, které jsou vzájemně nezávislé a druhé kritérium je použito jako sekvenční detektor, ve kterém jsou zjišťovány maximální pravděpodobnosti přijaté sekvence (3.9). ML sekvenční detektor počítá a porovnává pravděpodobnostní funkce P(y|Um), kde y je přijatý vektor a Um je jedna z možných přenesených sekvencí, která je dána podmínkou 16
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
(3.11) V AWGN kanálu, kde šum ovlivňuje nezávisle jednotlivé symboly, rozděluje ML dekodér výpočet pravděpodobnostních funkcí na jednotlivé části, které odpovídají přechodům v trellis diagramu v jednotlivých časových okamžicích (Obr. 3.4) (3.12) V rovnici (3.12) yi je i-tá část přijaté sekvence y a uim je i-tá část kódové sekvence Um. Každá i-tá pravděpodobnost je vypočítána jako rozdíl mezi i-tým přijatým symbolem kódového slova a ideálním symbolem uim. Pro dekodér s měkkým rozhodováním je běžně používána Euklidovská vzdálenost a pro dekodér s tvrdým rozhodováním je použita Hammingova vzdálenost. ML dekodér je velmi dobrý, ale jeho složitost není mnohdy akceptovatelná, protože počítá pravděpodobnosti pro všechny možnosti kódového slova (2n, kde n je délka kódového slova) a dále určuje nejlepší cestu v treliss diagramu. Celkem tedy počítá 2n x n větví. Pokud bude délka kódového slova delší než 100, není možné tuto úlohu provádět v reálném čase. Tuto úlohu velmi zjednodušuje Viterbiho algoritmus [ 23 ], který eliminuje počet možných kódových slov. V každém stavu vyhodnocuje délku jednolitých větví a na základě těchto hodnot vybere nejpravděpodobnější z nich. Tím se velmi sníží počet možných cest v trellis diagramu. 3.2.3 Dekódování Turbo kódů Viterbiho algoritmus je vhodný pro dekódování konvolučních kódů s poměrně nízkou složitostí a velkou pravidelností. Z toho důvodu jsou hojně využívány pro bezdrátovou komunikaci především v systémech vyvinutých v posledních dvaceti letech. Zisk konvolučních kódů je však vzdálen průměrně asi 4 dB od hranice kapacity kanálu [ 16 ]. Zlepšení lze dosáhnout použitím kódů s větší pamětí, což ale zapříčiní vznik složitějšího kódu. Jeho složitost se zvyšuje exponenciálně s počtem stavů kodéru, což je většinou neakceptovatelné. Zvyšují se nároky na výpočetní výkon, prodlužuje se čas dekódování a výsledkem je tedy systém s velkou latencí. Turbo kódy, uvedené v roce 1993 Berrouem [ 11 ], poskytují kódový zisk, který je vzdálen jen několik desetin od kapacity kanálu, avšak za cenu větší složitosti kodéru. Pro dekódování turbo kódů nebylo možné použít klasické algoritmy, které byly použity u konvolučních kódů, a proto bylo vyvinuto nové schéma obsahující dva SISO dekodéry propojené iterační smyčkou. Důvodem použití iterativního konvolučního dekodéru je odhad a posteriori pravděpodobnosti P(uk|y), kde uk je k-tý datový bit a y je přijatá sekvence (kódové slovo). Základní rovnice pro iterativní dekódování BPSK, kde „0“ je zobrazena jako -1 a „1“ je zobrazena jako +1. Dekodér odhaduje a posteriori poměr pro každý k-tý bit, který je dán následujícím vztahem
17
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
(3.13) Jestliže je poměr větší než 1, dekodér rozhodne, že k-tý bit je roven ‘1‘. Naopak jestliže je poměr menší než 1, dekodér vyhodnotí k-tý bit jako „0“. Úloha dekódování se velmi zjednoduší, jestliže se použije logaritmický a posteriori poměr L(uk): (3.14) V případě, že budou hodnoty datových bitů stejně pravděpodobné, lze konstatovat, že rozhodování na základě hodnoty a posteriori poměru, je stejné jako rozhodování na základě poměru pravděpodobností. V literatuře je označován vtah (3.14) jako logaritmický poměr pravděpodobností (LLR – LogLikelihood Ratio). Dekodér tedy provádí rozhodování na základě LLR (3.15) Tudíž znaménko uk je určováno pomocí tvrdého rozhodování a hodnota |L(uk)| určuje spolehlivost tohoto rozhodnutí. Rovnici LLR lze rozšířit pomocí Bayesových pravidel:
(3.16)
Dosazením rovnic (2.4) a (2.5) lze spočítat pravděpodobnost vyjádřenou výrazem (3.17) pro AWGN kanály: (3.17) kde Es je energie vysílaného signálu a N0 je rozptyl šumu. Hodnota poměru Es/N0 je známa jako poměr signál/šum. První část výrazu (3.17) reprezentuje informaci pro měkké rozhodování, která je výstupem demodulátoru a vstupuje do dekodéru. Nazýváme ji hodnota kanálu. Druhá část výrazu reprezentuje a priori informaci. (3.18) Vzhledem k tomu, že jsou obě pravděpodobnosti datových bitů stejně pravděpodobné P(uk = +1) = P(uk = -1), je hodnota a priori informace na začátku iterativního procesu nulová.
18
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
extrinsická informace
a priori informace
Le(u’)
L(u)
SISO dekodér vstupní informace
a posteriori informace
Lcy
L(u’)
Obr. 3.6: Modul SISO dekodéru
Po průchodu prvním SISO dekodérem je však hodnota a priori informace zpřesněna a již není nulová. Tato nenulová a priori informace (nyní nazývána jako extrinsická) je přiváděna na vstup druhého SISO dekodéru.
100
1 iterace
BER 10
-1
10
-2
5 iterací
10-3
10 iterací 10-4
20 iterací 10
-5
10-6 -.3
-.2
-.1
0
.1
.2
.3
.4
.5
.6
.7
.8
.9
1
Eb/N0 [dB] Obr. 3.7: Chybovost turbo kódu (R = 1/3, k = 1784b) v závislosti na počtu iterací [ 12 ]
Během průchodu druhým SISO dekodérem je její hodnota opět zpřesněna. Druhý SISO dekodér tedy poskytuje extrinsickou informaci prvnímu dekodéru. Tímto je dokončen první průchod iterační smyčkou. Na Obr. 3.6 je zobrazen modul SISO dekodéru. Na začátku iteračního procesu jsou informace vyměňující se mezi oběma SISO dekodéry nezávislé, avšak po několika iteracích jsou stále více korelované a hodnota L(uk) se zvětšuje. Nemá však význam zvyšovat počet iterací do nekonečna, protože se to pak příliš neprojeví na zvýšení účinnosti turbo 19
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
dekodéru. Situace je názorně zobrazena na Obr. 3.7, kde jsou vidět výsledky simulace turbo kódu (R = 1/3, k = 1784b) [ 12 ]. 3.2.4 MAP dekódovací algoritmus pro konvoluční turbo kódy Pro dekódování konvolučních kódů se běžně využíval a využívá Viterbiho dekódovací algoritmus [ 23 ]. Jedná se o algoritmus pro opravu chyb kódových slov. V roce 1974 byl však navržen nový typ trellis dekodéru [ 13 ]. Jedná se o dekodér využívající SISO dekódovací algoritmus, který je zaměřen na eliminaci bitových chyb. Pro dekódování každého bitu je využita tzv. a posteriori pravděpodobnost. Tento způsob dekódování byl testován na Markovově zdroji procházejícím DMC (Discrete Memoryless Channel) kanálem. Sekvence Markovova modelu lze popsat pomocí trellis diagramu (Obr. 3.8). Diagram znázorňuje kód s generujícími polynomy g0(x) = 1 + x2 + x3 a g1(x) = 1 + x1 + x3, 8 stavy a R = 1/2. výstup m = 0 (000b)
00 10 11
m = 1 (001b)
01 01
m = 2 (010b)
10
m = 3 (011b)
10 01
m = 4 (100b)
10 00
m = 5 (101b)
10 11
m = 6 (110b)
01 00
m = 7 (111b)
stavy
10 11 t=0
t=1
t=2 Vstupní bit = ‘0’ Vstupní bit = ‘1’
Obr. 3.8: Trellis diagram (8 stavů, R = 1/2) [ 13 ]
Trellis diagram reprezentuje Markovovův zdroj propojením jednotlivých stavů přechodovými pravděpodobnostmi. Jednotlivé stavy konvolučního kodéru v čase k označíme Sk = 0, 1, 2, m, …, M. Délka kodéru je dána K = log2 M + 1. Cílem kodéru je porovnat přijatou sekvenci R1N = R1, R2, …, RN a odhadnout APP stavů a přechodů Markovova procesu
20
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
(3.19) (3.20) Rovnice (3.19) reprezentuje APP stavu Sk = m. Rovnice (3.20) reprezentuje APP přechodu (Sk-1 = m’, Sk = m). Obě rovnice lze normalizovat, jestliže vezmeme v úvahu, že P{R1N} je konstanta pro přijatou sekvenci R1N (3.21) (3.22) Pravděpodobnost události, která je funkcí jednotlivých stavů lze vyjádřit jako součet příslušných λk(m). Obdobně lze také vyjádřit událost jednotlivých přechodů jako součet σk(m’, m). Aby bylo možné spočítat (3.21) a (3.22) je definována pomocná metrika (3.23) (3.24) (3.25) Na Obr. 3.9 jsou názorně zobrazeny jednotlivé symboly, které jsou využity v rovnicích (3.23), (3.24) a (3.25): αk(m) definované v rovnici (3.23) lze vidět jako pravděpodobnost bytí ve stavu m v časovém okamžiku k za předpokladu, že známe přijatou posloupnost od začátku až do okamžiku k. αk(m) tedy závisí na R1k a Sk. βk(m) definované v rovnici (3.24) lze vidět jako pravděpodobnost přijaté sekvence od okamžiku k+1 až do konce jestliže stav kodéru v čase k je m. γk(m’, m) je pravděpodobnost přijetí posloupnosti Rk a bytí ve stavu Sk daného předchozím stavem Sk-1. Reprezentuje tedy pravděpodobnost přechodu (větve) mezi předchozím a současným stavem k.
21
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
m=0
m=1
m=2
m=3
m=4
m=5
m=6
m=7
R1
Rk-1
Rk+1
Rk
Sk-1 R1 k
RN
Sk Rk+1N
k(m)
k(m)
k(m)
Obr. 3.9: Metrika závislostí jednotlivých stavů a větví v MAP algoritmu [ 13 ]
Obr. 3.9 ukazuje, že APP dekódovaného bitu dk je založena na důkladné analýze celého trellis diagramu. Kombinace výpočtů na základě pomocné metriky je optimální pro stanovení pravděpodobnosti pro jednotlivé stavy a větve v časovém okamžiku k. Aby bylo možné vypočítat hodnotu λk(m) pomocí αk(m) a βk(m), je potřeba nejdříve rozdělit druhou část výrazu uvedeného v rovnici (3.21) na dvě (pravou a levou) části vůči k a použít definici podmíněné pravděpodobnosti (3.26) První část rovnice (3.26) může být vyjádřena nezávisle na R1k na základě využití Markovovi vlastnosti, která říká, že jestliže známe Sk, pak jednotlivé události následující po čase k nezávisí na R1k (3.27) (3.28) Hodnotu σk(m’, m) lze pak vypočítat podle následujícího vztahu s využitím nadefinovaných pomocných metrik (3.29) Následující dvě rovnice (3.30) a (3.31) ilustrují výslednou podobu metrik pro výpočet pravděpodobností jednotlivých stavů 22
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
(3.30)
(3.31) Tyto rovnice lze vypočítat rekurzivně, jestliže stanovíme příslušné okrajové podmínky. αk je získána z αk-1, což naznačuje rekurzivní pohyb v dopředném směru. Naproti tomu, βk je získána z βk+1, což naznačuje rekurzivní pohyb ve zpětném směru. Okrajové podmínky pro dopřednou metriku jsou získány definicí α0 (3.32) Tyto podmínky platí, jestliže je kodér na začátku naplněn samými nulami. Tento stav odpovídá podmínce Sk = 0 s maximální pravděpodobností (tudíž α0(0) = 1). Okrajové podmínky pro zpětný směr (9.3) jsou poněkud více obtížné. Jestliže není dekodér na konci dekódování bloku vynulován, není zaručeno, že SN = 0 bude mít maximální pravděpodobnost (tudíž βN(0) = 1). (3.33) Ukončení trellis dekodéru je velmi důležité, avšak zvyšuje náročnost při dekódování turbo kódů pomocí MAP algoritmu. Ukončení je nutné provést z důvodu existence zpětnovazební smyčky. Pokud bychom tento krok neprovedli, není zaručeno, že budou na konci dekódování všechny stavy dekodéru v nule. V případě, že ponecháme dekodér bez ukončení, zapříčiníme tím snížení kódového zisku. Rovnice (9.3) popisuje metriku přechodů. (3.34)
Na základě popisu, který uvedl Berrou lze rovnici přepsat do tvaru (3.35)
kde q(.|.) může nabývat buď hodnoty 0, nebo 1. π(Sk = m, Sk-1 = m’) je pravděpodobnost přechodů mezi stavy a je porovnávána se statistikou vstupu kodéru. Tento výraz odpovídá a priori informaci. Při první iteraci je její hodnota
23
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
nulová a a priori P{dk = 0} = 1/2).
hodnoty
jsou
stejně
pravděpodobné
(P{dk = 1} =
Z toho důvodu bude q(.|.) vždy rovno 0, a proto lze rovnici (3.35) rozdělit na přechod do 0 a do 1. (3.36) (3.37) V případě RSC kodéru jsou hodnoty Rk rozděleny na dvě nezávislé části ,
(3.38)
kde je systematická část přijaté posloupnosti v čase k a je kódovaná část v čase t. Pravděpodobnost přechodu je nezávislá na stavech Sk a Sk-1, protože systematická část závisí pouze na nezakódovaném vstupu. Proto lze rovnice (3.36) a (3.37) přepsat do tvaru (3.39)
(3.40) APP pravděpodobnost lze pak spočítat jako (3.41) Podmíněnou pravděpodobnost lze upravit tak, aby byla do výpočtu zahrnuta metrika uvedená v předchozím textu. Dále je také poměr APP přepsán do logaritmického tvaru, aby bylo možné získat požadovaný měkký (soft) výstup.
(3.42)
Na základě rovnic (3.39) a (3.40) lze považovat systematickou část za nezávislou na stavech Sk a Sk-1 lze výraz rozložit na jednotlivé části (3.43) Pravá část výrazu v rovnici (3.43) representuje extrinsickou informaci, která je příspěvkem MAP dekodéru ke zlepšení LLR pro bit dk v další iteraci. Extrinsická informace efektivně procházející dalším dekodérem je získána odečtením a priori informace od druhého členu rovnice (9.3), jelikož a priori 24
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
informace byla přijata od druhého dekodéru v předchozí polovině iterační smyčky. Extrinsická informace generovaná dekodérem by měla být nezávislá na a priori informaci přijaté tímto dekodérem. (3.44)
25
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
4 Cíle disertace 4.1 Analýza progresivních kanálových kódů a jejich optimalizace pro návrh systému pro přenos SMS zpráv Cílem disertace je analýza efektivních progresivních kanálových kódů pro datovou družicovou komunikaci s aplikačním využitím pro systém přenosu krátkých textových zpráv s využitím družice PHASE 3E na vysoce eliptické dráze HEO (Highly Elliptical Orbit). Systém by měl být navržen tak, aby poskytovaná služba byla dostupná s co nejmenšími nároky kladenými na uživatele tzn. co nejnižší vysílací výkon, co nejmenší apertura. Systém pro přenos krátkých textových zpráv bude složen z následujících částí (Obr. 4.1): počítač s obslužným softwarem na straně uživatele, vysílací anténa, družice, přijímací anténa, server s obslužným softwarem. Družice PHASE 3E
Uživatel
PC
Server
kódovací dekódovací software
software pro zpracování SMS zpráv
Obr. 4.1: Koncepce systému pro přenos krátkých textových zpráv
4.1.1 Výpočet energetické bilance Aby bylo možné systém provozovat je potřeba v první části provést výpočet energetické bilance, na jejímž základě budou stanoveny parametry jednotlivých částí komunikačního řetězce. S ohledem na využití družicové komunikace, jsou hlavním problémem velké ztráty a rušení na přenosové cestě. Systém musí být navržen tak, aby byly jednotlivé části optimálně využity tzn. nesmí být předimenzovány ani poddimenzovány. Při návrhu je potřeba počítat
26
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
s jistými omezeními. Těmi hlavními jsou parametry družice, které jsou pevně dány [ 55 ]. 4.1.2 Analýza a modifikace Turbo kódů Hlavní myšlenkou pro úspěšnou realizaci systému je implementace vhodného druhu kódování, které umožní pracovat se signály s velmi nízkými poměry Eb/N0. S příchodem turbo kódů se objevily nové možnosti a díky jejich unikátním vlastnostem se začaly využívat i pro tento účel. Dokáží pracovat s velmi slabými signály a dosahují poměrně velkého kódového zisku narozdíl od ostatních typů běžně používaných kódů. Cílem disertační práce je výběr vhodného turbo kódu a jeho následná modifikace s cílem optimalizovat jeho vlastnosti pro dané použití. Turbo kód je potřeba analyzovat z hlediska volby vhodného kódového poměru, délky bloku dat a typu prokladače. Tyto tři parametry mají stěžejní vliv na efektivitu turbo kódu. 4.1.3 Metodika pro odhad chybovosti BER Aby bylo možné v reálném systému provádět výsledné zpracování zprávy, je zapotřebí stanovit metodiku pro odhad chybovosti BER. Finální zpracování zprávy se odvíjí od hodnoty dosažené chybovosti, a proto je tato úloha důležitou součástí systému. Bez znalosti hodnoty chybovosti by nebylo možné stanovit zda je výsledná podoba SMS zprávy relevantní či nikoliv. Současně je nutná znalost dosažené chybovosti BER i pro vlastní proces dekódování turbokódem kódovaných dat. Důvodem je iterační technika dekódování, kdy algoritmus dekodéru musí na základě stanoveného odhadu BER ukončit dekódování po příslušné iteraci (při překročení vhodně stanovené prahové hodnoty BER). 4.1.4 Stanovení struktury rámce SMS zprávy Struktura rámce SMS zprávy je taktéž nedílnou součástí systému a vychází z předpokladu, že je stanovena metoda pro odhad chybovosti. Velikost rámce musí vycházet z vlastností vybraného turbo kodéru a jeho struktura bude navržena s ohledem na možnost získání bezchybné SMS zprávy i při nenulové chybovosti, které bude pravděpodobně dosaženo. 4.1.5 Metodika pro stanovení symbolů výsledné SMS zprávy Vzhledem k předpokládanému dosažení nenulové chybovosti je nutná aplikace metodiky pro stanovení symbolů výsledné SMS zprávy. Tato metodika bude vycházet ze stanovené struktury SMS zprávy a vlastností turbo kódů. Především bude vycházet z teoretických poznatků popisující dekódování turbo kódů. 4.1.6 Synchronizace Předpokladem navrženého systému je jeho použití při velmi nízkých poměrech Eb/N0. Aby bylo možné aplikovat navržený turbo kód, který pracuje se striktně definovaným blokem dat, metody pro odhad chybovosti a stanovení výsledného podoby SMS zprávy, je zapotřebí provádět rámcovou synchronizaci. 27
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Pro tuto úlohu musí být zvolena vhodná metoda, která dokáže pracovat při již zmiňovaných nízkých hodnotách poměru Eb/N0 a je pomocí ní možné stanovit začátek rámce SMS zprávy. Rovněž je nutné zajistit synchronizaci symbolovou a u koherentní demodulace i synchronizaci fáze nosné. Tato část není v práci dále podrobně rozebírána, neboť je synchronizace symbolová a fáze nosné úzce svázána s použitou modulační technikou.
28
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
5 Energetická bilance Stěžejní úlohou při návrhu systému pro přenos krátkých textových zpráv je výpočet energetické bilance. Jedná se o analýzu jednotlivých částí družicového spoje a stanovení patřičných parametrů, které mají nezanedbatelný vliv na přenos signálu. Cílem disertační práce je vytvoření takového systému, který bude klást minimální požadavky na uživatelskou část. Řídící segment pak bude navržen s dostatečným výkonem, aby mohl systém spolehlivě fungovat.
5.1 Výpočet energetické bilance uživatel – družice – řídící segment Prvotní výpočet energetické bilance ve směru (uživatel – družice – řídící segment) je prováděn pro parametry uvedené v tabulce. Na základě výpočtu budou stanoveny kritické části systému a volba vhodných parametrů, pro dosažení optimálního návrhu družicového spoje. Tab. 5.1: Parametry systému UPLINK
DOWNLINK
Fv = 435 MHz
Fd = 2,3 GHz
D = 0,35 m
D=2m
Gtmax = 3 dB R = 50 000 km Lftx = 0,5 dB Parametry družice Grmax = 6 dB
Gtmax = 6 dB
Konverzní zisk = 150 dB
Lftx = 0,5 dB
F = 1 dB Lfrx = 0,5 dB Ostatní parametry B = 3 kHz 3dB
= 2°
29
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
5.1.1 UPLINK Výpočet EIRP pozemní stanice (5.1) Přídavné ztráty (5.2) Ztráty šířením pro vzestupnou trasu (5.3) (5.4) Šumová teplota antény (5.5) (5.6) Jakostní číslo přijímacího systému G/T (5.7) Výpočet poměru C/N0 a C/N (5.8) Úroveň signálu přijatého na družici (5.9) Úroveň šumu přijatého na družici (5.10) Konverzní zisk družice 150dB Výkon užitečného signálu a šumu na straně vysílací antény družice (5.11) Obdobně jako pro vzestupnou část byl proveden výpočet i pro sestupnou část. Výsledky jsou graficky vyneseny na obrázcích Obr. 5.1 a Obr. 5.2. Z uvedených výsledků je patrné, že pro zvolené parametry nevycházejí hodnoty příliš dobře. Cílem je, aby bylo dosaženo vhodného poměru signál/šum pro vysílací výkon cca 10 W. Z grafů je zřejmé, že teprve při hodnotě vysílací výkonu 27 W se dostáváme do stavu, kdy poměr signál/šum je roven 1 (tzn. 0 dB). 30
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Je tedy nutné změnit ty parametry, u kterých je to možné, a zlepšit tak úrovně signálu na přijímací straně. Za zcela fixní je nutné brát parametry družice. Parametry, které můžeme měnit, jsou: výkon vysílací stanice, anténa uživatele a na straně serveru. Vzhledem k tomu, že se snažím o dosažení co nejmenšího vysílacího výkonu, nebudeme příliš měnit hodnotu vysílacího výkonu ani velikost antény. Ta má být kompaktní a dobře přenositelná. Výkon přijatého signálu a šumu -168 -170 -172
P, N [dB(W)]
-174 -176 -178 -180 -182 -184 -186
0
5
10
15
20 25 30 Vysílaný výkon Ptx [W]
35
40
45
50
Obr. 5.1: Výkon přijatého signálu a šumu pro průměr paraboly D = 2 m
Poměr S/N na přijímací straně 5
Poměr S/N [dB]
0
-5
-10
-15
0
5
10
15
20 25 30 35 Vysílaný výkon Ptx [W]
40
45
Obr. 5.2: Poměr S/N na straně serveru pro průměr paraboly D = 2 m
31
50
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Jediný parametrem, který lze měnit, je zisk antény (paraboly) na straně serveru. Zde je možnost zvětšením průměru paraboly dosáhnout většího zisku a získat tak na straně serveru signál s lepším poměrem signál/šum. I zde však jsou jistá omezení. Je to především prostor pro instalaci velké paraboly a také cena. Tab. 5.2: Úroveň signálu přijatého na straně serveru v závislosti na průměru paraboly (D [m]) PRX [dB(W)] PTX [W]
D=1m
D=2m
D=3m
D=4 m
2
-188,97
-182,95
-179,44
-176,95
4
-185,96
-179,94
-176,43
-173,94
6
-184,20
-178,18
-174,67
-172,18
8
-182,95
-176,93
-173,42
-170,93
10
-181,98
-175,96
-172,45
-169,96
12
-181,19
-175,17
-171,66
-169,17
14
-180,52
-174,50
-170,99
-168,50
16
-179,94
-173,92
-170,41
-167,92
18
-179,43
-173,41
-169,90
-167,41
20
-178,97
-172,95
-169,44
-166,95
22
-178,56
-172,54
-169,02
-166,54
24
-178,18
-172,16
-168,65
-166,16
26
-177,83
-171,81
-168,30
-165,81
28
-177,51
-171,49
-167,98
-165,49
30
-177,21
-171,19
-167,68
-165,19
32
-176,93
-170,91
-167,40
-164,91
34
-176,67
-170,65
-167,13
-164,65
36
-176,42
-170,40
-166,89
-164,40
38
-176,18
-170,17
-166,65
-164,16
40
-175,96
-169,94
-166,43
-163,94
42
-175,75
-169,73
-166,22
-163,73
44
-175,55
-169,53
-166,01
-163,53
46
-175,35
-169,34
-165,82
-163,33
48
-175,17
-169,15
-165,64
-163,15
50
-174,99
-168,97
-165,46
-162,97
Z uvedené tabluky a následného grafického zobrazení, je vidět vliv průměru (efektivní plochy) paraboly na hodnotu signálu přijatého na straně serveru. Předchozí výpočty byly provedeny pro parabolu o průměru 2 m. Z uvedených grafů vyplývá, že velikost paraboly má velký vliv na zlepšení parametrů signálu. Zvětšením průměru na D = 3 m dojde ke zvýšení zisku paraboly a poměr signál/šum je roven 1 již při hodnotě vysílacího výkonu 12 W. Čímž dosáhneme snížení hodnoty vysílacího výkonu o více než polovinu. Jednotlivé hodnoty, při kterých úroveň signálu vzroste nad úroveň šumu, jsou 32
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
vyznačeny v Tab. 5.2. Tato hodnota by měla být postačující vhledem k tomu, že budou implementovány turbo kódy. Pomocí nich by měl být systém schopen se signály, které nemají příliš dobrý poměr signál/šum. Výkon přijatého signálu a šumu -160
D=4m D=3m D=2m
-165
P, N [dB(W)]
-170
D=1m
-175 -180 -185 -190 -195
0
5
10
15
20 25 30 35 Vysílaný výkon Ptx [W]
40
45
50
Obr. 5.3: Velikost přijatého signálu a šumu v závislosti na průměru paraboly (D)
Poměr S/N
Poměr S/N [dB] na přijímací straně
10
D=4m D=3m D=2m
5 0
D=1m -5 -10 -15 -20 -25
0
5
10
15
20 25 30 35 Vysílaný výkon Ptx [W]
40
45
Obr. 5.4: Hodnota poměru S/N v závislosti na průměru paraboly (D)
33
50
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
5.2 Výpočet energetické bilance řídící segment – družice – uživatel Vhledem k tomu, že systém musí umožnit komunikaci v obou směrech, byl proveden výpočet energetické bilance pro zadané parametry i v opačném směru. Jedná se o obdobný výpočet jako v předchozím případě s tím rozdílem, že se uživatel dostal do pozice příjmu a omezujícím faktorem není vysílací část, ale naopak přijímací část. Je nutné zajistit dostatečný vysílací výkon, aby bylo možné přijmout signál pomocí malé patchové antény a dekódovat ho pomocí obslužného softwaru na straně uživatele. Byly provedeny výpočty pro různé průměry paraboly. Výsledky jsou uvedeny v Tab. 5.3, ze které vyplývá, že vhodným kompromisem je volba paraboly o průměru D = 3 m. Jak je vidět při přenosu signálu ze serveru k uživateli jsou výkonové úrovně signálu poněkud horší, než tomu bylo v opačném směru. S vhodnou velikostí paraboly a využitím turbo kódů by však měl systém pracovat bez větších komplikací. Tab. 5.3: Úroveň signálu přijatého na straně uživatele v závislosti na průměru paraboly (D [m]) PRX [dB(W)] PTX [W]
D=1m
D=2m
D=3m
D=4m
2
-190,32
-184,43
-181,04
-178,67
4
-187,31
-181,42
-178,03
-175,66
6
-185,55
-179,66
-176,27
-173,90
8
-184,30
-178,41
-175,02
-172,65
10
-183,33
-177,44
-174,05
-171,68
12
-182,54
-176,65
-173,26
-170,89
14
-181,87
-175,98
-172,59
-170,22
16
-181,29
-175,40
-172,01
-169,64
18
-180,78
-174,89
-171,50
-169,13
20
-180,32
-174,43
-171,04
-168,67
22
-179,90
-174,01
-170,62
-168,26
24
-179,53
-173,64
-170,25
-167,88
28
-178,86
-172,97
-169,58
-167,21
30
-178,56
-172,67
-169,28
-166,91
32
-178,28
-172,39
-169,00
-166,63
34
-178,01
-172,12
-168,73
-166,37
36
-177,77
-171,88
-168,49
-166,12
38
-177,53
-171,64
-168,25
-165,88
40
-177,31
-171,42
-168,03
-165,66
42
-177,10
-171,21
-167,82
-165,45
44
-176,89
-171,00
-167,61
-165,25
46
-176,70
-170,81
-167,42
-165,05
48
-176,52
-170,63
-167,24
-164,87
50
-176,34
-170,45
-167,06
-164,69
34
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Výkon přijatého signálu a šumu -160
D=4m D=3m D=2m
-165
P, N [dB(W)]
-170 -175
D=1m
-180 -185 -190 -195
0
5
10
15
20 25 30 35 Vysílaný výkon Ptx [W]
40
Obr. 5.5: Velikost přijatého signálu a šumu na straně uživatele v závislosti na průměru paraboly (D)
35
45
50
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
6 Analýza turbo kódů 6.1 Výběr vhodného typu turbo kódu Pro každou aplikaci je zapotřebí vybrat vhodný typ turbo kódu. Jednotlivé turbo kódy jsou navrženy tak, aby co nejvíce splňovaly požadavky systému, ve kterém mají pracovat. Příkladem mohou být turbo kódy pro magnetický či optický záznam dat, vysokorychlostní komunikační systémy, satelitní komunikaci, telemetrii aj. Každý systém klade jiné požadavky na vlastnosti turbo kódu. Pro přenos krátkých textových zpráv prostřednictvím družicové komunikace byl zvolen jako základ turbo kód určený pro telemetrii a družicovou komunikaci [ 12 ], [ 32 ]. Tento typ turbo kódu je určen pro přenos krátkých datových bloků a je schopen pracovat s nízkými poměry Eb/N0. Lze tedy předpokládat, že modifikací tohoto turbo kódu, lze získat kód, který bude vhodný pro přenos krátkých textových zpráv (SMS) s využitím družicové komunikace. 6.1.1 Turbo kód pro telemetrii a družicovou komunikaci Turbo kodér se skládá ze dvou konvolučních kodérů, kdy na jeden je signál přiváděn přímo a na druhý přes prokládací člen. Výstupy kodérů jsou multiplexovány se vstupním neupraveným signálem. Výsledkem je zabezpečená posloupnost bitů, která je přenášena komunikačním systémem. Pro experimenty je využit turbo kód, jehož základem je kodér na Obr. 6.1 [ 12 ], [ 32 ]. Vstupní rámec toho kódu se skládá z následujících částí: rámec primární hlavičky: 48 bitů, které se liší podle typu použití datový rámec: (k – 64) datových bitů, rámec pro opravu chyb: 16 paritních kontrolních bitů pro opravu chyb generovaných podle polynomu 1 + D5 + D12 + D16. Vstupní rámec kódu by měl mít délku k = 1784, 3568, 7136 nebo 8920 bitů. Struktura turbo kodéru zobrazená na Obr. 6.1 se skládá z paralelního zapojení dvou shodných lineárních, systematických, rekurzivních konvolučních kodérů E1 a E2. Tvary jejich polynomů v dopředném směru jsou (1 + D + D3 + D4), (1 + D2 + D4) a (1 + D + D2 + D3 + D4). Zpětnovazební polynom má pak tvar (1 + D3 + D4).
36
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
E1 informaèní bity
u1,i
vstupní rámec (k bitù)
D
D
D
D p11,i
prokladaè
p12,i
p13,i u1,i
permutovaný informaèní rámec
D
D
D
D p21,i
v dy pøipojen
p22,i
nikdy nepøipojen
p23,i 1/2
støídavì pøipojen
1/3
1/4
1/6
E2
Obr. 6.1: Blokové schéma turbo kodéru pro telemetrii a družicovou komunikaci [ 12 ]
6.2 Prokladač Prokládání (interleaving) dat je technika hojně se vyskytující v komunikačních systém, jako prostředek ke zvýšení odolnosti proti chybovosti. V komunikačním kanálu se často vyskytují tzv. shluky chyb. Aby bylo možné tyto chyby účinně eliminovat, je potřeba tyto chyby rozprostřít. Tuto úlohu zastává právě prokladač. S příchodem turbo kódů vzrostla také důležitost prokladačů, jelikož jejich výkon závisí na kódových sekvencích, které se objevují na výstupu turbo kodéru. Úkolem prokladače je generování takových posloupností, aby jejich korelace byla co nejmenší. Účinnost prokladače závisí na zvolených parametrech. Kvalita prokladače tedy úzce souvisí s celkovou efektivitou turbo kódů a je jedním z hlavních faktorů, díky nimž se přibližujeme Shannonovu limitu [ 64 ]. Prokladač je tedy tím prvkem, který má hlavní vliv na rozdělení Hammingovy váhy kódové sekvence. Jak je již obecně známo, Hammingova váha je dána počtem „1“ v kódové sekvenci. S tím úzce souvisí pojem minimální Hammingova vzdálenost, což je počet „1“, ve kterých se liší dvě kódová slova. Čím menší je Hammingova vzdálenost, tím větší je pravděpodobnost chybně dekódovaného slova [ 65 ]. Minimální Hammingova vzdálenost je jedním z nejdůležitějších parametrů, který určuje kódový zisk při vysokých poměrech signál/šum (> 2 dB). V případě malých poměrů signál/šum dochází při dekódování ke zvýšení chybovosti zapříčiněné vysokou úrovní šumu. Jestliže snížíme úroveň šumu (zvýšíme 37
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
poměr Eb/N0) dojde ke zmenšení chybovosti díky povaze vstupních datových sekvencí. Obecně řečeno, malá Hammingova váha vstupní posloupnosti zapříčiňuje malou Hammingovu váhu posloupnosti výstupní. Abychom zamezili vzniku takových posloupností, obsahují turbo kódy prokladač, který záměnou pozice bitů umožní generování kódových sekvencí s vyšší Hammingovou vahou. Výstupní posloupnost je složena z bitů u, y1 a y2 a tudíž také Hammingova váha je dána třemi dílčími hodnotami w(u), w(y1) a w(y2). Použitím dobrého prokladače pak dosáhneme výslednou posloupnost s relativně velkou Hammingovou vahou a dosáhneme tak větší odolnost proti chybovosti, kterou s sebou přináší turbo kódy. Pro experimenty jsem aplikoval tři typy prokladačů, které jsem zvolil na základě jejich parametrů. u
u
y1
RSC 1
Prokladač
y2
RSC 2
Obr. 6.2: Obecné schéma zapojení turbo kodéru
6.2.1 Telemetrický prokladač Jedná se o prokladač standardně používaný pro telemetrii a družicovou komunikaci [ 12 ]. Blok prokladače má délku k. Je charakterizován zákonem permutace π: Zk Zk. Vstupní bity konvolučního kodéru E1 prokládány nejsou, zatímco bity vstupující do E2 jsou prokládány podle . Permutační zákon prokladače π: Zk podle následujícího algoritmu:
Zk byl navrhnut Berrouem a pracuje
Je dáno k, kde k = k1k2 (k1 = 8, k2 = k/8). Pak pro každé s od 0 do (k – 1) se vypočítá matematických operací:
38
(s) podle následujících
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
kde pq nabývá jedné z osmi následujících hodnot: p1 = 31; p 2 = 37; p 3 = 43; p4 = 47; p5 = 53; p 6 = 59; p7 = 61; p 8 = 67 [ 12 ]. 6.2.2 Pseudonáhodný prokladač Princip pseudonáhodného prokladače spočívá v generování pseudonáhodné posloupnosti, na jejímž základě dochází prokládání dat. Výhodou prokladače tohoto typu je jednoduchost, rychlost a účinnost. 6.2.3 Prokladač s eliminací tzv. „edge efektu” Tento prokladač byl původně aplikován v turbo kódech s krátkými bloky dat. Na základě poznatků z [ 64 ] je tento prokladač použit pro experimenty s delšími bloky dat (k = 1024, 4096, 8464 a 9216 bitů). Princip prokladače s periodou p a bloky dat m x m je popsán v následujících krocích: 1) Data jsou zapisována po sloupcích a následně vyčítána po řádcích z matice matice velikosti m x m. 2) Aplikace operace modulo p, na jejímž základě jsou data rozdělena do p skupin: Skupina 0, Skupina 1, Skupina2, … , Skupina p-1 kde jsou data ve Skupiněi složena z těch, které jsou na pozici k = i mod p, i = 0, 1, 2, … , p – 1. 3) Následně jsou data z jednotlivých skupin umísťována po blocích délky p (viz. Obr. 6.3c): Skupina 0 Skupina 1 . . . Skupina p-1
Blok0,1 , Blok0,2 , ... Blok1,1 , Blok1,2 , … . . . . . . . Blokp-1,1 , Blokp-1,2 , …
4) Proložení dat z bodu 3 na základě následující sekvence Blok0,1 , Blok 2,1 , ... , Blokp-1,1 Blok0,2 , Blok 2,2 , ... , Blokp-1,2 5) Je provedena inverze dat z bodu 4, z důvodu zamezení tzv. „edge“ efektu.
39
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
0
8
16 24 32 40 48 56
0
8
16 24 32 40 48 56
1
9
17 25 33 41 49 57
1
9
17 25 33 41 49 57
2
10 18 26 34 42 50 58
2
10 18 26 34 42 50 58
3
11 19 27 35 43 51 59
3
11 19 27 35 43 51 59
4
12 20 28 36 44 52 60
4
12 20 28 36 44 52 60
5
13 21 29 37 45 53 61
5
13 21 29 37 45 53 61
6
14 22 30 38 46 54 62
6
14 22 30 38 46 54 62
7
15 23 31 39 47 55 63
7
15 23 31 39 47 55 63
a) Data zapsaná po sloupcích 0
24 48
1
9
8
32 56 16 40
33 57 17 41
49 10 18 42
b) Data rozdělená do p skupin
3
2
63 55 31
25
7
47 23 62 39
15 54 46 22 61 38 14 53
26 50 11
30
6
45 37 13 52 29
34 58 19 27 51 12 35 59
44 21 60 36 28
20 43
59 35 12 51 27 19 58 34
5
4
28 36 60 21 44
29 52 13 37 45
6
30
11 50 26
53 14 38 61 22 46 54 15
25
39 62 23 47
40 16 56 32
7
31 65 63
c) Data po prokladu na základě bodu 4
2
3
4
5
43 20
42 18 10 49
41 17 57 33 8
9
1
48 24
0
d) Data po inversi na základě bodu 5
Obr. 6.3: Prokladač s p = 3, délka bloku 8 x 8 bitů
6.3 Implementace turbo kódu v prostředí MATLAB Pro zvolený turbo kód byly vytvořeny čtyři základní modely v prostředí MATLAB – Simulink pro poměry R = 1/2 , R = 1/3, R =1/4 a R = 1/6. Na Obr. 6.4 je vidět celkové zapojení kodéru i dekodéru. Pro experimenty je zvolen AWGN kanál, který se svými vlastnostmi blíží prostředí, kterým se šíří signál při družicové komunikaci. Pro analýzy byly implementovány jednotlivé typy prokladačů a jejich výsledky jsou prezentovány v následující části. Pro každý graf je zobrazena černá křivka, vyjadřující chybovost pro nezakódovaný signál (NC – no coded).
40
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Obr. 6.4: Celkové schéma modelu komunikačního řetězce s turbo kodérem a dekodérem (R = 1/6)
41
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
6.4 Analýza kódového poměru R a délky bloku dat k 6.4.1 Kódový poměr R = 1/2 pseudonáhodný k = 1784
telemetrický k = 1784
speciální k = 1024
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 0,0
0,5
1,0
1,5
2,0
2,5
Eb/N0
Obr. 6.5: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/2
pseudonáhodný k = 3568
telemetrický k = 3568
speciální k = 4096
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 0,0
0,5
1,0
1,5
2,0 Eb/N0 [dB] 2,5
Obr. 6.6: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/2
42
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Z uvedených grafů pro kódový poměr R = 1/2 je vidět, že minimálních chybovostí je dosaženo při hodnotě Eb/N0 = 1,6 dB. Minimální chybovostí, jak bude uváděno i dále v textu, je myšlen bod, kdy se začíná projevovat pouze zbytková chybovost a její hodnota již klesá jen velmi pomalu. Z průběhu grafů je dále patrné, že celkově se pro kódový poměr R = 1/2 nejlépe chová turbo kód s pseudonáhodným prokladačem. Ostatní dva typy prokladačů dosahují v tomto případě poněkud horších výsledků. pseudonáhodný k = 7136
telemetrický k = 7136
speciální k = 8464
NC
BER[-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 0,0
0,5
1,0
1,5
2,0 Eb/N0 [dB] 2,5
Obr. 6.7: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/2
pseudonáhodný k = 8920
telemetrický k = 8920
speciální k = 9216
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 0,0
0,5
1,0
1,5
2,0 Eb/N0 [dB] 2,5
Obr. 6.8: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/2
43
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
6.4.2 Kódový poměr R = 1/3 pseudonáhodný k = 1784
telemetrický k = 1784
speciální k = 1024
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04 -0,5
0,0
0,5
1,0
1,5 Eb/N0 [dB] 2,0
Obr. 6.9: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/3
pseudonáhodný k = 3568
telemetrický k = 3568
speciální k = 4096
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,5
0,0
0,5
1,0
1,5 Eb/N0 [dB] 2,0
Obr. 6.10: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/3
Pro kódový poměr R = 1/3 jsou výsledky pro všechny tři typy implementovaných prokladačů obdobné. Nepatrně horších výsledků dosahuje turbo kód se speciálním prokladačem. Tento turbo kód sice dosahuje menší chybovosti (pro k = 4096 a 8464 bitů) avšak při hodnotách Eb/N0 cca o 0,2 dB vyšších než ostatní dva typy. 44
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Na rozdíl od předchozích výsledků pro kódový poměr R = 1/2 se zde zmenšila hodnota Eb/N0, při které klesne chybovost na minimum. pseudonáhodný k = 7136
telemetrický k = 7136
speciální k = 8464
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,5
0,0
0,5
1,0
1,5 Eb/N0 [dB] 2,0
Obr. 6.11: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/3
pseudonáhodný k = 8920
telemetrický k = 8920
speciální k = 9216
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,5
0,0
0,5
1,0
1,5 Eb/N0 [dB] 2,0
Obr. 6.12: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/3
45
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
6.4.3 Kódový poměr R = 1/4 pseudonáhodný k = 1784
telemetrický k = 1784
speciální k = 1024
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04 -0,5
0,0
0,5
1,0
1,5 Eb/N0 [dB] 2,0
Obr. 6.13: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/4
pseudonáhodný k = 3568
telemetrický k = 3568
speciální k = 4096
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,5
0,0
0,5
1,0
1,5 Eb/N0 [dB] 2,0
Obr. 6.14: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/4
Pro kódový poměr R = 1/4 jsou výsledky velmi různorodé. Pro k = 1784 a 8920 bitů dosahuje nejlepších výsledků turbo kód s telemetrickým 46
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
prokladačem. Pro k = 7136 bitů (turbo kód s pseudonáhodným a telemetrickým prokladačem) a k = 8464 bitů (turbo kód se speciálním prokladačem) jsou výsledky téměř totožné. Velmi dobrých hodnot však dosahuje turbo kód se speciálním prokladačem pro k = 4096 bitů, který v celém průběhu dosahuje nejmenší chybovosti. pseudonáhodný k = 7136
telemetrický k = 7136
speciální k = 8464
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,5
0,0
0,5
1,0
Eb/N0 [dB]
1,5
Obr. 6.15: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/4
pseudonáhodný k = 8920
telemetrický k = 8920
speciální k = 9216
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,5
0,0
0,5
1,0
Eb/N0 [dB]
Obr. 6.16: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/4
47
1,5
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
6.4.4 Kódový poměr R = 1/6 pseudonáhodný k = 1784
telemetrický k = 1784
speciální k = 1024
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04 -0,5
0,0
0,5
1,0
Eb/N0 [dB]
1,5
Obr. 6.17: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/6
pseudonáhodný k = 3568
telemetrický k = 3568
speciální k = 4096
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,5
0,0
0,5
Eb/N0 [dB]
1,0
Obr. 6.18: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/6
Pro kódový poměr R = 1/6 dosahuje nejlepších výsledků turbo kód s telemetrickým prokladačem. S kódovým poměrem R = 1/6 je dosaženo minimálních chybovostí pro velmi nízké hodnoty Eb/N0 a jednotlivé
48
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
charakteristiky poměrně prudce klesají. Je však zřejmé, že těchto hodnot je dosaženo za cenu velkého navýšení počtu paritních bitů. pseudonáhodný k = 7136
telemetrický k = 7136
speciální k = 8464
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,5
0,0
0,5
Eb/N0 [dB]
1,0
Obr. 6.19: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/6
pseudonáhodný k = 8920
telemetrický k = 8920
speciální k = 9216
NC
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,5
0,0
0,5
Eb/N0 [dB]
Obr. 6.20: Graf závislosti bitové chybovosti BER na poměru Eb/N0, kódový poměr R = 1/6
49
1,0
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
6.4.5 Zhodnocení Jednotlivé simulace byly provedeny pro již zmíněné kódové poměry R = 1/2, 1/3, 1/4 a 1/6. Tyto kódové poměry byly zvoleny, protože jsou standardně používány s vybraným typem turbo kódu. V závislosti na výběru kódového poměru se odvíjí mimo jiné také náročnost na výpočetní výkon a tudíž i dobu zpracování.
Eb/N0 [dB]
Eb/N0 pro BER 10E-2
Eb/N0 pro BER 10E-3
1,80 1,60 1,40 1,20 1,00 0,80 0,60 0,40 0,20 0,00
BER [-]
Obr. 6.21: Porovnání hodnot chybovosti pro Eb/N0 = 10-2 a Eb/N0 = 10-3 (R = 1/2)
BER [ - ]
2,5E-04 2,0E-04 1,5E-04 1,0E-04 5,0E-05 0,0E+00
Eb/N0 [ dB ]
Obr. 6.22: Porovnání hodnot chybovosti (R = 1/2)
50
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Z uvedených výsledků je patrné, že se vzrůstajícím počtem paritních bitů klesá hodnota poměru Eb/N0, při které chybovost klesne na přijatelnou úroveň. Dalším faktorem, který má neméně důležitý vliv na chybovost je délka bloku dat k. Na Obr. 6.21 až Obr. 6.28 je provedeno srovnání turbo kodérů s jednotlivými typy prokladačů. Obr. 6.21, Obr. 6.23, Obr. 6.25 a Obr. 6.27 znázorňují hodnoty Eb/N0, při kterých je dosaženo chybovosti BER = , resp. BER = . Obr. 6.22, Obr. 6.24, Obr. 6.26 a Obr. 6.28 znázorňují hodnoty Eb/N0, při kterých je dosažena minimální chybovost. Hodnocení výsledků je provedeno s ohledem na účel použití. Dále je brán v úvahu fakt, že při přenosu SMS zprávy postačuje chybovost BER = . Tato úvaha se vztahuje k délce SMS zprávy. Přeneseme-li např. zprávu o délce 100 znaků (tj. 800b) je pravděpodobné, že bude zpráva obsahovat pouze jeden chybný znak, což plně postačuje. Na druhou stranu je potřeba zajistit, aby byly bezchybně přeneseny identifikační části zprávy jako adresa odesilatele a adresáta. Tato část je řešena v kapitole 8, zabývající se návrhem rámce SMS zprávy.
Eb/N0 [dB]
Eb/N0 pro BER 10E-2
Eb/N0 pro BER 10E-3
1,00 0,90 0,80 0,70 0,60 0,50 0,40 0,30 0,20 0,10 0,00
BER [-]
Obr. 6.23: Porovnání hodnot chybovosti pro Eb/N0 = 10-2 a Eb/N0 = 10-3 (R = 1/3)
Z Obr. 6.22 je vidět, že pro kódový poměr R = 1/2 je nejmenší chybovosti dosaženo při implementaci speciálního prokladače s délkou bloku k = 4096 bitů pro eliminaci tzv. „edge efektu“. Chybovost klesá na hodnotu BER = při poměru Eb/N0 = 1,6 dB. Provedeme-li však srovnání na základě grafu (viz. Obr. 6.21) lze říci, že ani jeden z typů implementovaných prokladačů nepřináší lepší výsledky než ostatní. Na Obr. 6.24 jsou zobrazeny hodnoty chybovostí pro kódový poměr R = 1/3. Zde vycházejí nejmenší hodnoty chybovosti pro varianty s nejdelšími bloky dat k = 7136 a 9216 bitů. Velmi dobré úrovně chybovosti dosáhne také turbo kodér se speciálním prokladačem. Nevýhodou je však, že chybovost klesne na přijatelnou úroveň až při hodnotě Eb/N0 = 1,0 dB, což je o 0,2 dB více 51
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
BER [ - ]
než u ostatních variant. Z grafu (viz. Obr. 6.23) je patrné, že nejmenší chybovosti je dosaženy s nejdelšími bloky dat. Vezmeme-li však v úvahu, že při délce bloku k = 3568 bitů, resp. k = 4096 bitů, přenášíme přibližně poloviční množství bitů, jsou výsledky pro tuto variantu poměrně dobré. 4,5E-04 4,0E-04 3,5E-04 3,0E-04 2,5E-04 2,0E-04 1,5E-04 1,0E-04 5,0E-05 0,0E+00
Eb/N0 [ dB ]
Obr. 6.24: Porovnání hodnot chybovosti (R = 1/3)
Varianta s kódovým poměrem R = 1/4 (Obr. 6.26) přináší opět nejmenší hodnoty chybovosti při použití prokladačů s nejdelšími bloky dat. V případě, že srovnáváme výsledky z Obr. 6.25 je vidět, že implementací telemetrického a speciálního prokladače s délkou bloku dat k = 3568 bitů, resp. k = 4096 bitů dosáhneme velmi nízké chybovosti. Eb/N0 pro BER 10E-2
Eb/N0 pro BER 10E-3
1,20 Eb/N0 [dB]
1,00 0,80 0,60 0,40 0,20 0,00
BER [-]
Obr. 6.25: Porovnání hodnot chybovosti pro Eb/N0 = 10-2 a Eb/N0 = 10-3 (R = 1/4) 52
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
BER [ - ]
Tato hodnota je srovnatelná i s turbo kódy s delšími bloky dat. V případě uvedených dvou variant přináší nejlepší výsledky turbo varianta, se speciálním prokladačem (k = 4096 b). Při poměru Eb/N0 = 0,8 dB klesne chybovost na hodnotu BER = . 9,0E-04 8,0E-04 7,0E-04 6,0E-04 5,0E-04 4,0E-04 3,0E-04 2,0E-04 1,0E-04 0,0E+00
Eb/N0 [ dB ]
Obr. 6.26: Porovnání hodnot chybovosti (R = 1/4)
Pro kódový poměr R = 1/6 (Obr. 6.28) dosáhneme lepších výsledků než v předchozím případě. Těchto hodnot je však dosaženo za cenu velkého navýšení počtu paritních bitů, a proto se jeví jako efektivnější varianta s kódovým poměrem R = 1/4. Eb/N0 pro BER 10E-2
Eb/N0 pro BER 10E-3
0,60
Eb/N0 [dB]
0,50 0,40 0,30 0,20 0,10 0,00
BER [-]
Obr. 6.27: Porovnání hodnot chybovosti pro Eb/N0 = 10-2 a Eb/N0 = 10-3 (R = 1/6)
53
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
BER [ - ]
6,0E-04 5,0E-04 4,0E-04 3,0E-04 2,0E-04 1,0E-04 0,0E+00
Eb/N0 [ dB ]
Obr. 6.28: Porovnání hodnot chybovosti (R = 1/6)
Celkově tedy vychází nejlépe varianta s kódovým poměrem R = 1/4, kdy je implementován speciální prokladač s eliminací tzv. „edge efektu“ s délkou bloku dat k = 4096 bitů (viz. Obr. 6.14). Tato varianta je vybrána s ohledem na dosaženou hodnotu chybovosti (BER = při Eb/N0 = 0,8 dB), složitost zapojení turbo kodéru, délku bloku dat a účel použití.
54
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
7 Odhad chybovosti zprávy V kapitole 6.4 byla provedena analýza turbo kódu a byla vybrána varianta, která přináší optimální výsledky pro přenos SMS zprávy. Jedná se o turbo kód s implementovaným prokladačem pro eliminaci tzv. „edge efektu“ s poměrem R = 1/4. Programové prostředí MATLAB Simulink umožňuje přesně stanovit chybovost BER a chování implementovaného turbo kódu. V reálném systém to však není možné, protože na přijímací straně není známa vyslaná posloupnost. Proto je potřeba se zaměřit na problematiku odhadu chybovosti BER. Jednotlivé metody pro odhad chybovosti jsou popsány v [ 66 ], [ 68 ] a [ 72 ]. Je zapotřebí vybrat a modifikovat vhodnou metodu, na jejímž základě bude možné stanovit chybovost přijaté zprávy s ohledem na akceptovatelnou výpočetní náročnost.
7.1 Odhad chybovosti BER na základě APP V systému, kde je použit dekodér se „soft“ výstupem, je výstupní hodnotou logaritmický poměr pravděpodobností (LLR). Uvažujme systém na Obr. 7.1, kde je i-tý přenesený bit, ni je gaussovský šum s rozptylem a yi je přijatý signál. Dále předpokládejme dokonalý odhad amplitudy signálu na straně přijímače a uvažujme relativní chybu odhadu 2, kterou představuje faktor . Rozložení pravděpodobnosti LLRi lze pak vyjádřit vztahem (7.1) [ 66 ]. (7.1) kde
označíme normální rozdělení průměru m s rozptylem
. Na základě uvedených nezakódovaného přenosu rovna
hypotéz
je
chybovost
a
BERNC
(7.2) kde
.
(7.3)
Chybovost BER je dána poměrem . Parametry α a však nelze odhadnout z pLLR(LLRi), jelikož se překrývají části křivek Gaussovského rozložení. Navíc statistické momenty rozdělení LLR závisí na a . Chybovost BER tedy nemusí splňovat vztah (7.4). (7.4)
55
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
2
ni
x
A
0,
+
x
2
yi
Odhad LLR
LLRi
xˆ
A
Obr. 7.1: Model přenosové cesty
Pomocí metody uvedené v [ 66 ] lze chybovost odhadnout použitím a posteriori pravděpodobností (APP), které se objevují na výstupu turbo dekodéru. V případě, že je uvažována přenosová cesta zobrazená na Obr. 7.1 a délka rámce k, dostaneme na výstupu a posteriori pravděpodobnosti APPi(+A) a APPi(-A) pro jednotlivé přenesené bity xi. Odhad pravděpodobnosti lze pak matematicky popsat pomocí vztahu (7.5) (7.5) Ekvivalentní systém pro odhad chybovosti byl navržen také v [ 69 ] a [ 71 ] a je popsán rovnicí (7.6). (7.6) Rovnice (7.5) a (7.6) vyjadřují to samé, jelikož platí (7.7) a (7.8). (7.7) (7.8) Na Obr. 7.2, Obr. 7.3 a Obr. 7.4 jsou uvedeny grafy závislosti chybovosti BER na poměru Eb/N0. Modrá křivka představuje průběh originálního průběhu chybovosti a červená křivka průběh odhadu chybovosti na základě popsané metody.
56
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
BER_original
BER_odhad
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04 -0,3
-0,1
0,1
0,3
0,5
0,7
0,9
1,1
1,3
1,5
Eb/N0 [dB] Obr. 7.2: Odhad chybovosti BER na základě APP (Turbo kód se speciálním prokladačem, R = 1/4, k = 1024b)
BER_oroginal
BER_odhad
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,3
-0,1
0,1
0,3
0,5
0,7
0,9
Obr. 7.3: Odhad chybovosti BER na základě APP (Turbo kód se speciálním prokladačem, R = 1/4, k = 4096b)
57
1,1 1,3 Eb/N0 [dB]
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
BER_original
BER_odhad
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,3
-0,1
0,1
0,3
0,5
0,7
0,9
1,1 1,3 Eb/N0 [dB]
Obr. 7.4: Odhad chybovosti BER na základě APP (Turbo kód se speciálním prokladačem, R = 1/4, k = 8464b)
Z uvedených grafů (viz. Obr. 7.2, Obr. 7.3, Obr. 7.4) je patrné, že metoda odhadu chybovosti BER na základě a posteriori pravděpodobnosti (APP) podává poměrné dobré výsledky. Na Obr. 7.5 je provedené procentuální srovnání chybovosti pro délky bloku dat k = 1024, 4096 a 8464 bitů. I když chyba není příliš velká, největší odchylku vykazuje odhad chybovosti pro k = 4096 bitů dat turbo kódu, který je vybrán pro přenos SMS zpráv, což není příliš vhodné, protože vybraný turbo kód pracuje právě s touto délkou bloku.
Odchylka [%]
100,0
Odchylka (k = 1024b) [%]
Odchylka (k = 4096b) [%]
Odchylka (k = 8464b) [%]
90,0 80,0 70,0 60,0 50,0 40,0 30,0 20,0 10,0 0,0 -0,2
0,0
0,2
0,4
0,6
0,8
Obr. 7.5: Zhodnocení odhadu chybovosti na základě APP
58
1,0
1,2 Eb/N0 [dB]
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
7.2 Modifikovaná metoda pro odhad chybovosti na základě LLR Tato metoda je navržena na základě poznatků z [ 75 ]. Pro jednotlivé bity přenášené posloupnosti lze vypočítat pravděpodobnost na základě vztahu (9.3), kde (n) je logaritmický poměr pravděpodobností. (7.9)
Vyhledávací tabulka 1 1 N
n
Vyhledávací tabulka 2
Zpětný prokladač
Dekodér
Obr. 7.6: Blokové schéma systému pro odhad chybovosti na základě LLR [ 75 ]
Princip metody popsané v [ 75 ] zachycuje blokové schéma zobrazené na Obr. 7.6. Zpětný prokladač rozdělí data na systematickou část a paritní posloupnosti. Na základě hodnot systematické posloupnosti je provedeno tzv. „tvrdé“ rozhodování. Tyto získané hodnoty jsou použity pro stanovení polohy přepínače. Na základě polohy přepínače je následně vybrána vyhledávací tabulka, podle níž je stanovena hodnota příspěvku pravděpodobnosti P(n) (viz. Tab. 7.3). Tab. 7.1: Vyhledávací tabulka pro LOG-MAP dekodér (přepínače v poloze log. „1“ – tzn. HARD hodnota > 0) SOFT hodnoty xSi Od
Do
Příspěvek pravděpodobnosti P(n)
-∞
-2,708005
0
-2,708005
-1,46634
0,125
-1,46634
-0,78846
0,250
-0,78846
-0,25131
0,375
-0,25131
0,251314
0,500
0,251314
0,788457
0,625
0,788457
1,466337
0,750
1,466337
2,70805
0,875
2,70805
∞
1 59
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Na základě experimentů bylo zjištěno, že popsaná metoda nepřináší příliš dobré výsledky pro odhad chybovosti zvoleného turbo kódu, který pracuje při velmi nízkých poměrech Eb/N0, a proto byla navržena modifikace uvedené metody pro odhad chybovosti BER. Modifikovaná metoda zobrazená na Obr. 7.7, je založena na principu popsaném v předchozím textu. Není zde však počítána celková pravděpodobnost na základě součtu příspěvků pravděpodobností jednotlivých bitů, ale pro každý bit je na základě sady rozhodovacích pravidel (viz. Tab. 7.3) stanoveno, zda byl daný bit dekódován správně či špatně. HARD hodnoty (xHi)
SOFT hodnoty (xSi)
Zpětný prokladač
Sada rozhodovacích pravidel
Dekodér
Obr. 7.7: Blokové schéma modifikované metody pro odhad chybovosti BER
Tab. 7.2: Sada rozhodovacích pravidel HARD hodnoty (xHi)
SOFT hodnoty (xSi)
Výstupní bity (yi)
xHi > 0
xSi ≥ -t1
1
xHi > 0
xSi < -t1
0
xHi ≤ 0
xSi ≤ t1
0
xHi ≤ 0
xSi > t1
1
xHi > -t2 & xHi < t2
xSi > -t2 & xSi < t2
chyba
t1 – rozhodovací úroveň 1, t2 – rozhodovací úroveň 2
Jak je uvedeno v Tab. 7.3 princip určení chybovosti je poměrně jednoduchý, avšak je zapotřebí správně stanovit hodnoty rozhodovacích úrovní (t1 a t2). Hodnota rozhodovací úrovně (t1 = 2,70805) je stanovena na základě Tab. 7.1. Je to mezní hodnota, která s velkou pravděpodobností určí, zda byl bit dekódován správně či špatně. Poněkud složitější situace nastává při určení hodnoty rozhodovací úrovně (t2). V případě, že hodnota xHi a xSi spadá do intervalu vymezeným rozhodovací úrovní t2, je velmi pravděpodobné, že se jedná o chybně dekódovaný bit. Tato skutečnost vyplývá z Tab. 7.1. Hodnota byla stanovena na základě experimentálních měření chybovosti BER pro různé turbo kódy. Stanovení hodnoty t2 je důležitou úlohou, avšak přesnost této hodnoty není kritická, jelikož se jedná o odhad chybovosti. 60
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
BER_original
BER_odhad (t2 = 0,5)
BER_odhad (t2 = 1,5)
BER_odhad (t2 = 2)
BER_odhad (t2 = 1)
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04 -0,3
-0,1
0,1
0,3
0,5
0,7
0,9
1,1
1,3 1,5 Eb/N0 [dB]
Obr. 7.8: Odhad chybovosti BER na základě LLR (Turbo kód se speciálním prokladačem, R = 1/4, k = 1024b)
BER_original
BER_odhad (t2 = 0,5)
BER_odhad (t2 = 1,5)
BER_odhad (t2 = 2)
BER_odhad (t2 = 1)
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,3
-0,1
0,1
0,3
0,5
0,7
0,9
1,1
Obr. 7.9: Odhad chybovosti BER na základě LLR (Turbo kód se speciálním prokladačem, R = 1/4, k = 4096b)
61
1,3 1,5 Eb/N0 [dB]
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
BER_original
BER_odhad (t2 = 0,5)
BER_odhad (t2 = 1,5)
BER_odhad (t2 = 2)
BER_odhad (t2 = 1)
BER [-]
1,0E+00
1,0E-01
1,0E-02
1,0E-03
1,0E-04
1,0E-05 -0,3
-0,1
0,1
0,3
0,5
0,7
0,9
1,1
1,3 1,5 Eb/N0 [dB]
Obr. 7.10: Odhad chybovosti BER na základě LLR (Turbo kód se speciálním prokladačem, R = 1/4, k = 8464b)
Na Obr. 7.8, Obr. 7.9 a Obr. 7.10 jsou uvedeny výsledky simulací odhadu chybovosti s využitím modifikované metody pro tři různé varianty turbo kódu. V každém grafu je zobrazena křivka znázorňující originální hodnoty (černá) a čtyři křivky odhadu chybovosti pro různé hodnoty rozhodovací úrovně (t2 = 0,5; 1; 1,5; 2). Pro jednotlivé délky bloku dat k = 1024, 4096 a 8464 bitů je na Obr. 7.11, Obr. 7.12 a Obr. 7.13 provedené procentuální zhodnocení výsledků odhadu chybovosti. Z uvedených grafů vyplývá, že pro hodnoty t2 = 0,5 a t2 = 2 je odhad chybovosti nepřesný a dosahuje chyby až 80 %. Velmi dobrých výsledků je však dosaženo pro hodnoty t2 = 1 a t2 = 1,5. Pro tyto dvě hodnoty t2 se odhad chybovosti velmi přibližuje originálnímu průběhu (viz. Obr. 7.8, Obr. 7.9 a Obr. 7.10). Pro zvolený turbo kód se speciálním prokladačem (R = 1/4, k = 4096 b), který je použit pro přenos SMS zpráv, dosahuje nejlepších výsledků metoda pro t2 = 1. Chyba nepřesahuje 35 % a především v oblasti hodnot chybovosti BER = 10-2 a 10-3, která je stěžejní pro správnou funkčnost systému, je velmi malá odchylka od originálního průběhu (viz. Obr. 7.12).
62
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Odchylka (t2 = 0,5) [%]
Odchylka (t2 = 1,0) [%]
Odchylka (t2 = 1,5) [%]
Odchylka (t2 = 2,0) [%]
Odchylka [%]
100,0 80,0 60,0 40,0 20,0 0,0
BER [-] Obr. 7.11: Zhodnocení odhadu chybovosti na základě LLR (k = 1024b)
Odchylka (t2 = 0,5) [%]
Odchylka (t2 = 1,0) [%]
Odchylka (t2 = 1,5) [%]
Odchylka (t2 = 2,0) [%]
Odchylka [%]
100,0 80,0 60,0 40,0 20,0 0,0
BER [-] Obr. 7.12: Zhodnocení odhadu chybovosti na základě LLR (k = 4096b)
Odchylka (t2 = 0,5) [%]
Odchylka (t2 = 1,0) [%]
Odchylka (t2 = 1,5) [%]
Odchylka (t2 = 2,0) [%]
Odchylka [%]
100,0 80,0 60,0 40,0 20,0 0,0
BER [-] Obr. 7.13: Zhodnocení odhadu chybovosti na základě LLR (k = 8464b)
63
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
7.3 Srovnání metod pro odhad chybovosti BER Pro odhad chybovosti byly aplikovány dvě metody. První metoda popsaná v [ 66 ] využívá pro odhad chybovosti a posteriori pravděpodobnost (APP) a druhá modifikovaná metoda vycházející z [ 75 ], pak využívá pro výpočet odhadu LLR hodnoty. Podrobné výsledky jsou uvedeny v kapitole 7.1 a 7.2. Obě metody vykazují velmi dobré výsledky a jsou použitelné pro odhad chybovosti zvoleného turbo kódu. Při bližším pohledu na výsledné hodnoty obou metod však vyplývá, že lepších výsledků dosahuje modifikovaná metoda. V Tab. 7.3 jsou uvedeny výsledné hodnoty obou metod. Modifikovaná metoda dosahuje větší přesnosti při hodnotě rozhodovací úrovně t2 = 1, která má stěžejní vliv na přesnost odhadu chybovosti. Maximální odchylka je 34,6 % a to pouze pro krajní hodnotu BER = 8,3.10-5, kde se již vyskytuje zbytková chybovost. V oblasti, kde chybovost nabývá hodnot BER = 10-2 ÷ 10-3, je odchylka maximálně 17,5 %. Jak již bylo zmíněno v předchozím textu, výsledný odhad chybovosti závisí na hodnotě rozhodovací úrovně t2. Přesnost parametru t2 však není kritická a pro dosažení kvalitního odhadu chybovosti se může pohybovat v rozmezí 1 ≤ t2 ≤ 1,5. Tab. 7.3: Srovnání výsledků metod pro odhad chybovosti
BER [-]
Odchylka [%] Metoda s využitím APP
Odchylka [%] Metoda s využitím LLR (t2 = 1)
2,37E-01
16,4
11,3
1,99E-01
24,8
0,4
9,52E-02
31,5
18,5
1,92E-02
39,4
14,0
1,90E-04
11,2
8,3
1,21E-04
31,2
17,5
1,07E-04
33,2
15,4
8,30E-05
42,2
34,6
Odchylka *%+ Metoda s využitím APP
Odchylka *%+ Metoda s využitím LLR (t2 = 1)
Odchylka [%]
100,0 80,0 60,0 40,0 20,0 0,0
BER [-] Obr. 7.14: Grafické srovnání výsledků metod pro odhad chybovosti 64
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
8 Návrh rámce SMS zprávy Návrh struktury rámce SMS zprávy je důležitou součástí systému pro přenos SMS zpráv. Délka rámce se odvíjí od délky bloku dat k zvoleného turbo kódu. Jedná se o turbo kód s implementovaným speciálním prokladačem, délkou bloku k = 4096 bitů a kódovým poměrem R = 1/4. Rámec zprávy je navržen tak, aby i v případě výskytu chyb došlo k jednoznačné identifikaci zprávy, tj. aby bylo možné zajistit bezchybný příjem adresy příjemce i odesilatele. Dále je nutné stanovit metodu pro odhad chybovosti BER přenesené zprávy. Odhadem chybovosti se zabývá kapitola 7, kde jsou navrženy, a porovnány jednotlivé metody. adresa příjemce
adresa odesilatele
délka SMS zprávy
Informační bity
12B
12B
8B
480B
512B
Obr. 8.1: Rámec SMS zprávy
Na Obr. 8.1 je zobrazena navržená struktura rámce o velikosti 512B, tj. 4096b. Rámec je doplněn o počáteční a koncovou část, která jsou zapotřebí pro synchronizaci. Tato část se však nevztahuje k procesu samotného kódování. Informační část rámce se skládá z následujících částí: Adresa příjemce (12B) – slouží pro identifikaci příjemce. Adresa má délku 2B, avšak pro eliminaci chyb a úspěšnou identifikaci příjemce je zopakována šestkrát. Adresa odesilatele (12B) – slouží pro identifikace odesilatele. Struktura zprávy je stejná jako v případě adresy příjemce. Délka SMS zprávy (8B) – slouží pro stanovení délky zprávy. Každá SMS zpráva může mít různý počet znaků (max. 160 znaků). Na základě této informace je SMS zpráva rozdělena na n bloků, což je využito pro eliminaci chyb. Informační bity (480B) – obsahuje text SMS zprávy o délce max. 160 znaků.
65
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
8.1 Metodika stanovení hodnot jednotlivých symbolů přenášené datové posloupnosti Rámec SMS zprávy je složen ze čtyř výše uvedených částí. Kromě informačních bitů textové části zprávy jsou důležitou součástí rámce adresa příjemce, adresa odesilatele a délka SMS zprávy. Každá část je nezbytná pro bezchybný přenos a správné nalezení adresáta, proto je navržen mechanismus pro výpočet hodnot jednotlivých symbolů. 8.1.1 Pravděpodobnostní metoda pro výpočet hodnot symbolů datové posloupnosti Metoda využívá k výpočtu hodnoty xHi, které se vyskytují na vstupu dekodéru a „soft“ hodnoty xSi vyskytující se na výstupu dekodéru. Bitové posloupnosti jsou rozděleny na čtyři části (adresa odesilatele, adresa příjemce, délka SMS zprávy, SMS zpráva). Výsledná hodnota bitů přenášené zprávy je stanovena na základě následujících kroků: transformace hodnot (xTi) – transformace je prováděna na základě hodnoty bitů přijaté posloupnosti na vstupu dekodéru. Pravidla pro transformaci popisuje (8.1). (8.1) výpočet pravděpodobnosti (pi) – výpočet pravděpodobnosti je založen na rovnicích (8.2) a (8.3) [ 75 ]. S výpočtem úzce souvisí hodnoty bitů na vstupu dekodéru. Jestliže je hodnota bude pravděpodobnost vypočítána pomocí (8.2). Pokud je však hodnota je proveden výpočet pravděpodobnosti na základě rovnice (8.3). (8.2) (8.3) transformace hodnot pravděpodobnosti (pTi) – pro vypočítané hodnoty pravděpodobnosti z předchozího kroku je provedena transformace hodnot na základě (8.4). (8.4) výpočet výsledné hodnoty bitu (yj) – výsledná hodnota bitů je vypočítána pomocí (9.1). Hodnota k vyjadřuje počet opakování části zpracovávaného úseku rámce SMS zprávy. (8.5) 66
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
HARD hodnoty (xHi) xi
Zpětný prokladač
SOFT hodnoty (xSi) Dekodér
Výpočet pravděpodobnosti
Transformace: xHi > 0 => xTi = 1 xHi < 0 => xTi = -1
yj
Aritmetické operace
Transformace: pi: <0;1> => pTi: <-1;1>
Obr. 8.2: Princip metody pro určení hodnot bitů přenášené datové posloupnosti
8.1.1.1
Příklad
Pro uvedený příklad je uvažována posloupnost x = 11011B, která se čtyřikrát opakuje. posloupnost bitů na vstupu dekodéru (xHi) XH1
X H2
X H3
X H4
X H5
X H6
X H7
X H8
X H9
XH10
2,5535
2,4246
1,7874
-0,8916
4,6683
0,6377
-0,4168
-3,1659
2,9367
0,0390
XH11
XH12
XH13
XH14
XH55
XH16
XH17
XH18
XH19
XH20
1,3889
0,4740
0,2559
0,3112
4,1475
-1,1247
0,6230
-1,5507
0,0069
2,6374
posloupnost bitů na výstupu dekodéru (xSi) XS1
X S2
X S3
X S4
X S5
X S6
X S7
X S8
X S9
XS10
3,500
1,000
-5,250
-4,000
1,750
2,250
-0,750
-5,625
6,125
7,500
XS11
XS12
XS13
XS14
XS55
XS16
XS17
XS18
XS19
XS20
-1,250
2,375
-3,125
-1,500
5,125
1,750
1,875
-4,250
2,500
2,250
transformace hodnot (xTi) XT1
X T2
X T3
X T4
X T5
X T6
X T7
X T8
X T9
XT10
1
1
1
-1
1
1
-1
-1
1
1
XT11
XT12
XT13
XT14
XT55
XT16
XT17
XT18
XT19
XT20
1
1
1
1
1
-1
1
-1
1
1
67
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
výpočet pravděpodobnosti (pi) p1
p2
p3
p4
p5
p6
p7
p8
p9
p10
0,97
0,73
0,01
0,98
0,85
0,90
0,68
1,00
1,00
1,00
p11
p12
p13
p14
p55
p16
p17
p18
p19
p20
0,22
0,91
0,04
0,18
0,99
0,15
0,87
0,99
0,92
0,90
transformace hodnot pravděpodobnosti (pTi) pT1
pT2
pT3
pT4
pT5
pT6
pT7
pT8
pT9
pT10
0,94
0,46
-0,99
0,96
0,70
0,81
0,36
0,99
1,00
1,00
pT11
pT12
pT13
pT14
pT55
pT16
pT17
pT18
pT19
pT20
-0,55
0,83
-0,92
-0,64
0,99
-0,70
0,73
0,97
0,85
0,81
výpočet výsledné hodnoty bitu (yj) y1
y2
y3
y4
y5
y6
y7
y8
y9
y10
0,94
0,46
-0,99
-0,96
0,70
0,81
-0,36
-0,99
1,00
1,00
y11
y12
y13
y14
y55
y16
y17
y18
y19
y20
-0,55
0,83
-0,92
-0,64
0,99
0,70
0,73
-0,97
0,85
0,81
určení hodnoty bitů na základě tzv. tvrdého rozhodování y1,6,11,16 = 1,9
y2,7,12,17 = 1,67
y3,8,13,18 = -3,87
y4,9,14,19 = 0,24
y5,10,15,20 = 3,5
1
1
0
1
1
8.1.2 Zjednodušená metoda pro výpočet hodnoty bitů Zjednodušená metoda pro stanovení symbolů využívá ke své funkci „soft“ hodnot, které jsou výstupem turbo dekodéru. Na Obr. 8.3 je zobrazen vývojový diagram popisující princip celého systému. Po přijetí rámce SMS zprávy je provedeno rozdělení na čtyři části (adresa příjemce, adresa odesilatele, délka SMS zprávy a informační část SMS zprávy). Adresa příjemce i odesilatele jsou rozděleny na 6 x 2B. Část obsahující informaci o délce SMS zprávy je rozdělena na 8 x 1B. Pro každý bit jednotlivých částí je vypočten aritmetický průměr, na jehož základě je následně stanovena hodnota bitu tzv. tvrdým rozhodováním. Zpracování adresy příjemce a odesilatele tímto krokem končí. Část zprávy obsahující informaci o délce informační části SMS zprávy je následně použita pro rozdělení této části na n stejně dlouhých bloků. Každá zpráva může mít délku max. 160 znaků, tj. 160B, z čehož vyplývá, že se informační část zprávy může opakovat maximálně 4krát. Pro každý bit je opět spočítána hodnota aritmetického průměru a pomocí tvrdého rozhodování je stanovena hodnota jednotlivých bitů. Výhodou toho způsobu zpracování zprávy je, že lze využít situace, kdy je SMS zpráva kratší než 160 znaků. V takovém případě se může opakovat SMS zpráva více než 4krát, a tudíž lze pro stanovení hodnoty jednotlivých bitů využít více nezávislých hodnot přenášené datové posloupnosti. 68
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Start
Rozdělení rámce na 4 části
Rozdělení
Rozdělení
Rozdělení
Adresa příjemce 6x 2B
Adresa odesilatele 6x 2B
Délka zprávy 10x 1B
SMS zprávy
Výpočet aritmetického průměru pro jednotlivé n-tice bitů
Výpočet aritmetického průměru pro jednotlivé n-tice bitů
Výpočet aritmetického průměru pro jednotlivé n-tice bitů
Výpočet aritmetického průměru pro jednotlivé n-tice bitů
Tvrdé rozhodování (HARD decision)
Tvrdé rozhodování (HARD decision)
Tvrdé rozhodování (HARD decision)
Adresa příjemce
Adresa odesilatele
Délka zprávy
Rozdělení
Tvrdé rozhodování (HARD decision) Informační část SMS zprávy
Výsledná SMS zpráva
Konec
Obr. 8.3: Blokové schéma zjednodušené metody pro zpracování přijaté SMS zprávy a určení hodnot jednotlivých bitů datové posloupnosti
8.1.2.1
Příklad
Pro uvedený příklad je uvažována posloupnost x = 11011B, která se čtyřikrát opakuje. posloupnost „soft“ hodnot na výstupu dekodéru XS1
X S2
X S3
X S4
X S5
X S6
X S7
X S8
X S9
XS10
3,500
1,000
-5,250
-4,000
1,750
2,250
-0,750
-5,625
6,125
7,500
XS11
XS12
XS13
XS14
XS55
XS16
XS17
XS18
XS19
XS20
-1,250
2,375
-3,125
-1,500
5,125
1,750
1,875
-4,250
2,500
2,250
69
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
rozdělení posloupnosti na čtyři čtveřice bitů XS1
X S2
X S3
X S4
X S5
3,500
1,000
-5,250
-4,000
1,750
X S6
X S7
X S8
X S9
XS10
2,250
-0,750
-5,625
6,125
7,500
XS11
XS12
XS13
XS14
XS15
-1,250
2,375
-3,125
-1,500
5,125
XS16
XS17
XS18
XS19
XS20
1,750
1,875
-4,250
2,500
2,250
výpočet aritmetického průměru hodnot pro jednotlivé čtveřice bitů (8.6) určení hodnoty bitů na základě tzv. tvrdého rozhodování sa1
= 1,56
sa2
1
= 1,13
sa3
1
= -4,56
sa4
= 0,78
0
sa5
1
= 4,16 1
8.1.3 Výsledky metod pro stanovení hodnot symbolů Na Obr. 8.4 až Obr. 8.10 jsou zobrazeny chybovosti BERsms v závislosti na počtu iterací. Jednotlivé grafy zobrazují hodnoty chybovosti pro Eb/N0 = –0,2 ÷ +1 dB. Obecně platí na základě uvedených výsledků, že pomocí metod pro stanovení hodnot bitů výsledné SMS zprávy je dosaženo menší chybovosti pro jednotlivé poměry Eb/N0 v porovnání s chybovostí dosažené bez využití uvedených metod. Je však zapotřebí brát v úvahu, že je těchto hodnot dosaženo za cenu navýšení redundance. Modrou barvou jsou znázorněny hodnoty chybovosti bez využití metod pro stanovení hodnot bitů, zelenou barvou jsou znázorněny hodnoty chybovosti získané pomocí pravděpodobnostní metody a červenou barvou jsou zobrazeny hodnoty chybovosti pro zjednodušenou metodu.
BER [-]
Chybovost bez využití navržených metod
Zjednodušená metoda
Pravděpodobnostní metoda
3,0E-02 2,5E-02 2,0E-02 1,5E-02 1,0E-02 5,0E-03 0,0E+00 1
2
3
4
5
6
7
8
9
10
Počet iterací [-]
Obr. 8.4: Chybovost výsledné SMS zprávy (Eb/N0 = –0,2dB)
70
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
BER [-]
Chybovost bez využití navržených metod
Zjednodušená metoda
Pravděpodobnostní metoda
2,5E-02 2,0E-02 1,5E-02 1,0E-02 5,0E-03 0,0E+00 1
2
3
4
5
6
7
8
9
10
Počet iterací [-]
Obr. 8.5: Chybovost výsledné SMS zprávy (Eb/N0 = 0dB)
BER [-]
Chybovost bez využití navržených metod
Zjednodušená metoda
2,5E-02 2,0E-02 1,5E-02 1,0E-02 5,0E-03 0,0E+00 1
2
3
4
5
6
7
8
9
10
Počet iterací [-]
Obr. 8.6: Chybovost výsledné SMS zprávy (Eb/N0 = 0,2dB)
BER [-]
Chybovost bez využití navržených metod
Zjednodušená metoda
Pravděpodobnostní metoda
2,0E-02 1,8E-02 1,6E-02 1,4E-02 1,2E-02 1,0E-02 8,0E-03 6,0E-03 4,0E-03 2,0E-03 0,0E+00 1
2
3
4
5
6
7
8
9
10
Počet iterací [-]
Obr. 8.7: Chybovost výsledné SMS zprávy (Eb/N0 = 0,4dB)
71
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
BER [-]
Chybovost bez využití navržených metod
Zjednodušená metoda
Pravděpodobnostní metoda
1,8E-02 1,6E-02 1,4E-02 1,2E-02 1,0E-02 8,0E-03 6,0E-03 4,0E-03 2,0E-03 0,0E+00 1
2
3
4
5
6
7
8
9
10
Počet iterací [-]
Obr. 8.8: Chybovost výsledné SMS zprávy (Eb/N0 = 0,6dB)
BER [-]
Chybovost bez využití navržených metod
Zjednodušená metoda
Pravděpodobnostní metoda
1,6E-02 1,4E-02 1,2E-02 1,0E-02 8,0E-03 6,0E-03 4,0E-03 2,0E-03 0,0E+00 1
2
3
4
5
6
7
8
9
10
Počet iterací [-]
Obr. 8.9: Chybovost výsledné SMS zprávy (Eb/N0 = 0,8dB)
BER [-]
Chybovost bez využití navržených metod
Zjednodušená metoda
Pravděpodobnostní metoda
1,2E-02 1,0E-02 8,0E-03 6,0E-03 4,0E-03 2,0E-03 0,0E+00 1
2
3
4
5
6
7
8
9
10
Počet iterací [-]
Obr. 8.10: Chybovost výsledné SMS zprávy (Eb/N0 = 1dB)
72
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Z uvedených průběhů je na první pohled zřetelné snížení chybovosti na více než poloviční hodnoty. Cílem navržených metod je však získání bezchybné SMS zprávy, a proto jsou nejdůležitějšími výsledky hodnoty zobrazené na Obr. 8.8 a Obr. 8.9 a Obr. 8.10. Pro poměr Eb/N0 = 0,6 dB je dosaženo nulové chybovosti při počtu 10 iterací s využitím pravděpodobnostní metody. Pomocí zjednodušené metody nelze pro tento poměr Eb/N0 dosáhnout nulové chybovosti. Při poměru Eb/N0 = 0,8 dB je dosaženo nulové chybovosti již pro obě metody, avšak i v tomto případě podává pravděpodobnostní metoda lepší výsledky, protože chybovost klesne na nulovou hodnotu již po 8 iteracích. Pro poměr Eb/N0 = 1 dB je opět lepší pravděpodobnostní metoda, protože je výsledná chybovost nulová již po 4 iteracích. 8.1.4 Zhodnocení metod pro stanovení hodnot symbolů Obě metody jsou založeny na navržené struktuře rámce SMS zprávy (Obr. 8.1). Při určování jednotlivých symbolů využívají obě metody opakující se části zprávy. První metoda využívá výpočet pravděpodobností a následnou transformaci hodnot. Druhá zjednodušená metoda je založena na výpočtu aritmetického průměru „soft“ hodnot na výstupu dekodéru a následném tvrdém rozhodování. Obě metody přináší velmi dobré výsledky a umožňují snížení chybovosti výsledné zprávy při stejném nebo dokonce nižším poměru Eb/N0 v porovnání s chybovostí, kterou lze docílit bez použití navržené struktury rámce. Ve výsledném hodnocení vykazuje pravděpodobnostní metoda lepších výsledků než zjednodušená metoda (viz. Obr. 8.4 až Obr. 8.10).
73
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
9 Synchronizace Proces synchronizace fáze nosné, symbolu a rámce významně ovlivňuje celý detekční proces a nedostatečná synchronizace je příčinou zvýšení chybovosti [ 80 ]. Také při použití turbo kódů je velmi důležitou součástí systému, z čehož vyplývá, že lze jejich unikátních vlastností využít pouze v případě, že je tato úloha provedena velmi precizně. Vlastní problematika synchronizace symbolové a nosné je natolik rozsáhlá a je úzce svázána s aplikovanou modulační technikou, že ji nelze v této práci pojmout a tato stať bude zaměřena pouze na synchronizaci rámcovou. Stěžejní požadavek na přijímací systém je funkčnost i při velmi nízkých poměrech Eb/N0, pro které je nasazení turbo kódu ještě úspěšné. To však naráží na problém spolehlivé rámcové synchronizace. Cílem je tedy nalézt takový přístup, kdy i při tak nízkých poměrech Eb/N0 lze rámcovou synchronizaci spolehlivě zajistit. Typickými systémy, které úspěšně pracují i se signály pod úrovní šumu, jsou systémy CDMA, které se staly inspiraci pro řešení tohoto problému. Pro zajištění rámcové synchronizace je využita vhodná rozprostírací posloupnost aplikovaná v CDMA systémech jako unikátní synchronizační slovo. V uvažovaném systému jsou navrženy posloupnosti generované pomocí Goldových kódů, které se vyznačují unikátními synchronizačními vlastnostmi [ 78 ].
9.1 Goldovy kódy Tyto kódy byly navrženy v roce 1967 [ 78 ] a jsou pojmenovány po jejich tvůrci Dr. Robertu Goldovi. Jedná se o binární sekvence s velmi dobrými korelačními vlastnostmi [ 79 ] (viz. Obr. 9.1). Sadu n Goldových sekvencí lze získat z tzv. preferovaného páru o délce L = 2n – 1 operací modulo-2. Korelace Goldových kódů nabývá tří hodnot {–1, –t(m) , t(m) – 2 }, kde (9.1) Goldovy kódy jsou pro své velmi dobré korelační vlastnosti používány jako rozprostírací posloupnosti v systémech CDMA a v satelitní navigaci GPS.
74
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
300 250 200
R [-]
150 100 50 0 -50 0
200
400
600 n [-]
800
1000
1200
Obr. 9.1: Autokorelační průběh funkce Goldova kódu o délce L = 1023b
9.2 Aplikace Goldových kódů Použití Goldových kódů pro časovou synchronizaci vychází z korelačních vlastností těchto kódů uvedených v předchozí kapitole. Blokové schéma systému pro synchronizaci je zobrazeno na Obr. 9.3. Princip systému spočívá v umístění Goldovy posloupnosti před rámec zakódované SMS zprávy (viz. Obr. 9.2). Synchronizační Goldova posloupnost
Zakódovaná SMS zpráva
511b
4096b
Obr. 9.2: Struktura rámce se synchronizační Goldovou posloupností L = 511b
Přijatá posloupnost je vzorkována takovou frekvencí, aby pro každý bit bylo získáno osm vzorků. Pomocí korelace navzorkované bitové posloupnosti a kopie Goldovy posloupnosti se systém snaží nalézt začátek rámce SMS zprávy. Předpokladem pro správnou funkci systému je možnost přesného synchronizace i pro nízké poměry Eb/N0. Na Obr. 9.4, Obr. 9.5 a Obr. 9.6 jsou uvedeny výsledky korelace při aplikaci Goldovy posloupnosti o délce 127b, 511b a 1023b při poměru Eb/N0 = 0dB. Z grafických průběhů je patrné, že použití Goldových kódů přináší velmi dobré výsledky i při velmi nízkém poměru Eb/N0 a lze je použít pro synchronizaci.
75
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Při detailnějším pohledu na grafy uvedené na Obr. 9.4, Obr. 9.5 a Obr. 9.6 však zjistíme, že při pouhém odečtení maximální hodnoty korelační funkce dochází k jisté chybě při určení hodnoty vzorku. Mnohem přesnější určení hodnoty vzorku dosáhneme na základě výpočtu z dvojic vzorků, které se vyskytují symetricky vlevo a vpravo od vrcholu. Na uvedených grafech není chyba příliš patrná, avšak pro nižší poměry Eb/N0 je chyba relativně velká (viz. Obr. 9.7). Vzorkovač Korelace
Hladání začátku rámce SMS zprávy
Zpracování dekódovaných dat
Výsledná SMS zpráva
Generátor Goldovy posloupnosti
Turbo dekodér
Obr. 9.3: Blokové schéma časové synchronizace
500
450
400
400 R [-]
350
300
300
250
R [-]
200
200
2030
2035
2040 n [-]
2045
2050
100
0
-100 0
500
1000
1500
2000 2500 n [-]
3000
3500
4000
4500
Obr. 9.4: Výsledky korelace přenášené bitové posloupnosti a Goldova kódu (127b), Eb/N0 = 0 dB 76
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
2000
2500 1800
1600 R [-]
2000
1400
1500
1200
R [-]
1000
1000
8175
8180
8185 n [-]
8190
8195
500
0
-500 0
2000
4000
6000
8000 10000 12000 14000 16000 18000 n [-]
Obr. 9.5: Výsledky korelace přenášené bitové posloupnosti a Goldova kódu (511b), Eb/N0 = 0 dB
4000
4500
3800 3600
4000
3400 3200
R [-]
3500
3000 2800
3000
2600 2400
2500
2200
R [-]
2000
2000
1.6365
1.637
1.6375
1.638 n [-]
1.6385 4
x 10
1500 1000 500 0 -500 0
0.5
1
1.5
2 n [-]
2.5
3
3.5 4
x 10
Obr. 9.6: Výsledky korelace přenášené bitové posloupnosti a Goldova kódu (1023b), Eb/N0 = 0 dB
77
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
2200
2000
2000 1800
R [-]
1600
1500
1400
1200
1000
1000 8180
R [-]
8175
8185 n [-]
8190
8195
500
0
-500 0
2000
4000
6000
8000 10000 12000 14000 16000 18000 n [-]
R [-]
Obr. 9.7: Výsledky korelace přenášené bitové posloupnosti a Goldova kódu (511b), Eb/N0 = -10 dB
2200
(xv, yv)
2000 1800 1600 1400
(x0, y0) (x1, y1)
1200 1000 8175
8180
8185
8190
n [-]
8195
Obr. 9.8: Aproximace vrcholu a korelační funkce a určení chyby na základě podobnosti trojúhelníků
Chybu lze eliminovat výpočtem, na jehož základě lze přesně stanovit hodnotu vzorku, ve kterém nabývá korelační funkce maxima. Pro stanovení této 78
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
hodnoty lze vycházet z podobnosti trojúhelníků, jestliže provedeme aproximaci vrcholu korelační funkce (viz. Obr. 9.8). Rovnice (9.2) vychází z uvedené podobnosti trojúhelníků a na jejím základě lze provést odvození pro výpočet hodnoty vzorku, ve kterém nabývá korelační funkce maxima (9.3). (9.2) (9.3) Pro uvedenou metodu byly provedeny simulace a výsledky jsou zobrazeny v grafech na Obr. 9.10 a Obr. 9.11. Grafy znázorňují chybu časové synchronizace (Δδ/Tb) v závislosti na hodnotě poměru Eb/N0. Simulace byly provedeny pro hodnoty Eb/N0 = -10 ÷ 1 dB. Chyba časové synchronizace se pohybuje pro Goldův kód L = 127b v rozmezí 0,0269 až 16,27 bitů, L = 511b v rozmezí 0,0136 až 0,0506 bitu a L = 1023b v rozmezí 0,01 až 0,0346 bitu. Nejdůležitější částí křivky je však oblast pro Eb/N0 = 0 ÷ 1 dB, ve které pracují navržené metody systému. Chyba v této části nabývá hodnot pro L = 127b v rozmezí 0,0269 až 0,0298 bitu, L = 511b v rozmezí 0,0136 až 0,0157 bitu a L = 1023b v rozmezí 0,01 až 0,0111 bitu. Z toho plyne, že pro všechny uvedené typy Goldova kódu je chyba natolik malá, že nebude mít zásadní vliv na stanovení hodnoty bitu, a tudíž plně postačuje i kratší verze Goldova kódu o délce L = 127b. Uvedené výsledky se vztahuji k pulzům typu REC, u kterých je předpokládán ideální obdélníkový průběh. Velikost odchylky bude samozřejmě záviset na tvaru impulzů symbolů vyskytujících se v reálném systému, i tak by však měl být synchronizace pomocí Goldových kódů spolehlivá, protože je chyba časové synchronizace velmi malá. Uvedené posloupnosti je možné využít také pro synchronizaci symbolovou včetně možného vkládání kratších Goldových posloupností do vlastního informačního rámce SMS zprávy pro udržení časové synchronizace.
79
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Δδ/Tb [-]
0,06
0,05
0,04
0,03
0,02
0,01
0,00 -5
-4
-3
-2
-1
0
1
Eb/N0 [dB]
Δδ/Tb [-]
18 16 14 12 10 8 6 4 2 0 -11
-10
-9
-8
-7
-6
-5
-4
-3
-2
-1
0
1
Eb/N0 [dB] Obr. 9.9: Chyba časové synchronizace v závislosti na poměru Eb/N0 Goldova posloupnost (127b)
80
2
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Δδ/Tb [-]
0,06
0,05
0,04
0,03
0,02
0,01
0,00 -11
-10
-9
-8
-7
-6
-5
-4
-3
-2
-1
0
1
2
Eb/N0 [dB] Obr. 9.10: Chyba časové synchronizace v závislosti na poměru Eb/N0 Goldova posloupnost (511b)
Δδ/Tb [-]
0,040 0,035 0,030 0,025 0,020 0,015 0,010 0,005 0,000 -11
-10
-9
-8
-7
-6
-5
-4
-3
-2
-1
0
1
Eb/N0 [dB] Obr. 9.11: Chyba časové synchronizace v závislosti na poměru Eb/N0 Goldova posloupnost (1023b)
81
2
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
10 Závěry Disertační práce se zabývá problematikou nových progresivních metod kanálového kódování při přenosu datových zpráv prostřednictvím družicového transpondéru. Součástí práce bylo navrhnout systém, který bude schopen při velmi nízkých poměrech signál/šum umožnit bezchybný přenos SMS zpráv. Práce se zaměřuje především na problematiku energetické bilance, modifikaci a implementaci nových typů turbo kódů a využití jejich unikátních vlastností, na vývoj metod pro odhad chybovosti a metod pro stanovení výsledné SMS zprávy při nenulové chybovosti. Pro uvedený systém byl proveden výpočet energetické bilance systému ve směru uživatel – družice – řídící segment a zpět s ohledem na fakt, že uživatel bude moci přistupovat k systému pomocí malé „patchové“ antény. Jak je uvedeno v [ 4 ], družicový systém nelze jednoduše předimenzovat, a proto byl brán zřetel na optimální využití jednotlivých částí sytému. Výsledkem energetické bilance je stanovení průměru paraboly D = 3 m na straně řídícího segmentu pro uvažovaný vysílací výkon PTX = 12 W na straně uživatele. Hlavním cílem disertace byl výběr a modifikace vhodného typu turbo kódu, který je schopen pro relativné malé bloky dat zajistit co nejnižší chybovost BER při velmi nízkých poměrech Eb/N0. Pro systém byl zvolen turbo kód často využívaný pro telemetrii a družicovou komunikaci [ 12 ]. Pro zvolený typ turbo kódu pak byly prováděny modifikace v podobě změny implementovaného prokladače – telemetrický, pseudonáhodný a prokladač s eliminací tzv. „edge efektu“. Byly navrženy varianty pro kódový poměr R = 1/2, R = 1/3, R = 1/4 a R = 1/6, které jsou standardně využívány i v jiných komunikačních systémech. Pro uvedené varianty byl turbo kodér testován s délkou bloku dat k = 1024, k = 1784, k = 3568, k = 4096, k = 7136, k = 8464 a k = 9216 bitů v závislosti na typu použitého prokladače. Nejlepších výsledků dosáhl turbo kód s implementovaným prokladačem pro eliminaci „edge efektu“, s kódovým poměrem R = 1/4 a délkou bloku dat k = 4096 bitů. Tento kód dosáhl, v porovnání s ostatními typy, nejnižší chybovosti BER = při Eb/N0 = 0,8 dB. Kódy s kódovým poměrem R = 1/6 sice dosahují nižší chybovosti avšak za cenu poměrně velkého navýšení redundance, a proto byla vybrána uvedená variantu turbo kódu. Výsledky jsou stanoveny na základě vlastností, které mají uvedené typy kódů v oblasti, kde chybovost klesá na hodnoty BER = 10-2 ÷10-3. Tato hodnota chybovosti je zapotřebí pro následné zpracování zprávy. Další část práce je zaměřena na návrh metod pro odhad chybovosti BER. Jedná se o nedílnou součást systému, na jejímž základě lze stanovit, kdy klesne chybovost na dostatečně nízkou úroveň, aby bylo možné dosáhnout bezchybného přenosu SMS zprávy po výsledném zpracování. Pro odhad chybovosti byla implementována metoda popsaná v [ 66 ], která je založena na 82
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
APP pravděpodobnostech a modifikovaná metoda, jejímž základem je popis v [ 75 ]. Obě metody dosahují velmi dobrých výsledků, avšak s využitím modifikované metody je dosaženo při správné volbě parametru t 1 lepších výsledků. Pro modifikovanou metodu je maximální odchylka 17,5 % pro oblast, kde chybovost nabývá hodnot BER = 10-2 ÷ 10-3, která je, jak již bylo zmiňováno, pro zpracování nejdůležitější. Dalším tématem, které je zpracováno v disertační práci je návrh rámce SMS zprávy a metodika stanovení hodnoty bitů datové posloupnosti. Rámec zprávy se skládá z adresy příjemce (6 x 2B), adresy odesilatele (6 x 2B), délky zprávy (8 x 1B) a samotného textu zprávy (480B). Rámec zprávy je navržen s ohledem na konečné zpracování, kdy je zapotřebí při nenulové chybovosti získat výslednou SMS zprávu. Pro tento účel byly navrženy dvě metody. První z nich je pravděpodobnostní metoda a druhou je zjednodušená metoda. Obě metody využívají navrženou strukturu rámce a opakující se bitové sekvence. Pravděpodobnostní metoda je založena na výpočtu pravděpodobností a následné transformaci hodnot, kdežto zjednodušená metoda využívá ke stanovení hodnot bitů pouze aritmetický průměr a tvrdé rozhodování. Obě metody přináší velmi dobré výsledky, avšak pravděpodobnostní metoda umožní získat bezchybnou SMS zprávu již při poměru Eb/N0 = 0,6dB a při poměrech Eb/N0 = 0,8dB a Eb/N0 = 0,8dB je dosaženo nulové chybovosti při menším počtu iterací. Posledním tématem disertační práce je časová synchronizace, která je velmi důležitá pro správnou funkci navržených metod a především pro efektivní využití turbo kódů. Jak je uvedeno v [ 79 ] a [ 80 ] nepřesná synchronizace vede ke zvýšení chybovosti. Proto je důležité tuto úloha zvládnout velmi precizně. Pro synchronizaci jsou navrženy Goldovy kódy, které mají velmi dobré korelační vlastnosti (viz. Obr. 9.1). Princip navržené metody spočívá v umístění Goldovy posloupnosti před rámec SMS zprávy, generování kopie posloupnosti na straně příjmu a výpočtu jejich korelace. Byly simulovány tři typy Goldových kódů s délkou L = 127b, L = 511b a L = 1023b. Pro poměr Eb/N0 = 0 ÷ 1 dB nabývá chyba časové synchronizace (Δδ/Tb) hodnot pro L = 127b v rozmezí 0,0269 až 0,0298 bitu, L = 511b v rozmezí 0,0136 až 0,0157 bitu a L = 1023b v rozmezí 0,01 až 0,0111 bitu. Tyto výsledky jsou velmi dobré a plně postačuje použití Goldova kódu s délkou L = 127b. Uvedené posloupnosti je možné využít také pro synchronizaci symbolovou, včetně možného vkládání kratších Goldových posloupností do vlastního informačního rámce SMS zprávy pro udržení synchronizace. Dalším téma, které se nabízí pro řešení, je problematika již zmiňované symbolové synchronizace a také synchronizace nosné. Tato problematika je úzce svázána s použitou modulační technikou, která však již není předmětem této disertační práce.
83
11 Literatura [1] [2] [3] [4] [5] [6] [7] [8] [9] [ 10 ] [ 11 ]
[ 12 ] [ 13 ] [ 14 ] [ 15 ]
[ 16 ] [ 17 ] [ 18 ] [ 19 ] [ 20 ] [ 21 ] [ 22 ] [ 23 ]
DOBEŠ, J. ŽALUD, V. Moderní radiotechnika. Praha: BEN, 2000, 768 s. ISBN 80-7300-1322. ŽALUD, V. Moderní radioelektronika. 1. vyd. Praha : BEN, 2004, 656 s. ISBN 80-86056-473. ŽALUD, V., Radioelektronika. Praha: ČVUT, 1993. ISBN 80-01-01042-2 KASAL, M. Směrové a družicové spoje. Vysoké učení technické v Brně. 2. vyd. Brno, 2005. 107 s. ISBN 80-214-2496-6. FREEMAN, R. L. Fundamentals of Telecommunications. Wiley-IEEE Press, 2005, ISBN 0471710458. ADÁMEK, J., Kódování. Praha, SNTL, 1989. NĚMEC, K., Datová komunikace. Skripta [online]. 2003 [cit. 2008-05-08] HANUS, S., Rádiové a mobilní komunikace. Skripta [online]. 2003 [cit. 2008-04-23]. HANUS, S. Bezdrátové a mobilní komunikace. 1. dotisk vyd. Brno : VUT, 2005. 134 s. ISBN 80-214-1833-8. JIROUŠEK, R., IVÁNEK, J., MÁŠA, P., TOUŠEK, J., VANĚK, N. Principy digitální komunikace. LEDA. Voznice, 2006. 309 s. ISBN 80-7335-084-1. BERROU, C., GLAVIEX, A., THITIJSHIMA P. Near Shannon Limit Error Correcting Coding And Decoding: Turbo Codes. IEEE International Conference on Communications, Geneva, Switzerland, May 1993, Vol. 2/3, p. 1064-1071. SRIPIMANWAT, K., Turbo Code Applications. Netherlands, Springer, 2005. ISBN 10 14020-3686-8. GIULIETTI, A., BOUGARD, B., PERRE, L., Turbo Codes: Desirable and Designable. United States of America, Kluwer Academic Publisher Group, 2004. ISBN 1-4020-7660-6. VUCETIC, B., Turbo Codes. Principles and applications. Netherlands, Springer, 2000, 386 s. ISBN 0-7923-7868-7. SOLEYMANI, M., GAO, Y., VILAIPORNSAWAI, U. Turbo Coding for Satelite and Wireless Communications. Netherlands, Kluwer Academic Publisher, 2002, 214 s. ISBN 1-40207197-3 BURR, A.., Modulation and Coding for Wireless Communications. England, Pearson Education, 2001, 360 s. ISBN 0-201-39857-5. BIGLIERI, E.., Coding for Wireless Channels. United States of America, Springers Science and Business Mediam Inc., 2005, 427 s. ISBN 1-4020-8083-2. GOLDSMITH, A.., Wireless Communications. United States of America, Cambridge University Press, 2005, 644 s. ISBN 0-521-83716-2. VITERBI, A. J., OMURA, J. K., Principles of Digital Communiation and Coding. United States of America, McGraw-Hill, Inc., 1979, 560s., ISBN 0-07-067516-3. GUTHERY, B. S., CRONIN, M. J., Coding for Wireless Channels. United States of America, McGraw-Hill, Inc., 2002, 306 s. ISBN 0-07-137540-6. VISWANATH, T., TSE, D., Fundamentals of Wireless Communication. United States of America, Cambridge University Press, 2005, 564 s. ISBN 978-0-521-84527-4. ZIEMER, R., PETERSON, R., Digital communication. New Jersey, Prentice-Hall, 1992. ISBN 0-13-896481-5. HAGENAUER, J.; HOEHER, P., A Viterbi algorithm with soft-decision outputs and its applications, Global Telecommunications Conference, 1989, and Exhibition.
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
[ 24 ] [ 25 ]
[ 26 ]
[ 27 ] [ 28 ] [ 29 ] [ 30 ] [ 31 ] [ 32 ]
[ 33 ] [ 34 ]
[ 35 ]
[ 36 ] [ 37 ] [ 38 ] [ 39 ] [ 40 ] [ 41 ] [ 42 ] [ 43 ] [ 44 ] [ 45 ] [ 46 ]
Communications Technology for the 1990s and Beyond. GLOBECOM '89., IEEE , vol., no., pp. 1680-1686 vol.3, 27-30 Nov 1989. VUCETIC, B., YONGHUI Li, PEREZ, L.C., FAN JIANG, Recent Advances in Turbo Code Design and Theory, Proceedings of the IEEE , vol.95, no.6, pp.1323-1344, June 2007. HAGENAUER, J, ROBERTSON, P, PAPKE, L. Iterative (”turbo”) decoding of systematic convolutional codes with the map and the sova algorithms. In ITG Tagung Codierung fur Quelle, Kanal und Ubertragung. October 1994. s. 21-29. LEE, L.-N., HAMMONS, A.R., Jr., FENG-WEN Sun, EROZ, M., Application and standardization of turbo codes in third-generation high-speed wireless data services, Vehicular Technology, IEEE Transactions on , vol.49, no.6, pp.2198-2207, Nov 2000. VLČEK, K., Turbo-kódy a rádiový přenos dat. Sdělovací technika, 2000, č.8, s. 24-26. HRDINA, Z., VEJRAŽKA, F. Digitální radiová komunikace. Skripta ČVUT v Praze. Praha, vydavatelství ČVUT, 1995. ISBN 80-01-01059-7. SÝKORA, J., Digitální rádiová komunikace II. Skripta ČVUT v Praze. Praha,vydavatelství ČVUT, 1995. ISBN 80-01-01374-X. BERROU, C, GLAVIEUX, A. Near optimum error correcting coding and decoding: Turbocodes. IEEE Trans. Communication, vol 44, no. 10, 1996, p. 1261 -1270. MASWY, J.L., Foundation and methods of channel encoding. Proc. Int. Conf Information Theory and Systems, vol. 65, NTG-Fachberichte, Berlin, 1978, p. 148- 157 HAGENAUER, J.; OFFER, E.; PAPKE, L., Iterative decoding of binary block and convolutional codes, Information Theory, IEEE Transactions on , vol.42, no.2, pp.429-445, Mar 1996. MUDER, D.J., Minimal trellises for block codes, Information Theory, IEEE Transactions on , vol.34, no.5, pp.1049-1053, Sep 1988. Calzolari, G.-P., Vassallo, E., Habinc, S., CCSDS telemetry channel coding: the turbo coding option, New Technologies, New Standards (Ref. No. 1998/519), IEE 5th CCSDS Workshop , vol., no., pp.5/1-5/6, 9 Nov 1998. BENEDETTO, S., MONTORSI, G., Unveiling turbo codes: some results on parallel concatenated coding schemes, Information Theory, IEEE Transactions on , vol.42, no.2, pp.409-428, Mar 1996. BENEDETTO, S., MONTORSI, G., Design of parallel concatenated convolutional codes, Communications, IEEE Transactions on , vol.44, no.5, pp.591-600, May 1996. Lossles Data Compression. Blue Book [online]. 1997 [cit. 2008-09-10]. Dostupný z WWW:
. TM Space Data Link Protokol. Blue Book [online]. 2003 [cit. 2008-09-10]. Dostupný z WWW: . TM Space Packet Protocol. Blue Book [online]. 2003 [cit. 2008-09-10]. Dostupný z WWW: . TM Communications Operation Procedure-1. Blue Book [online]. 2003 [cit. 2008-09-10]. Dostupný z WWW: . Time Code Formats. Blue Book [online]. 2002 [cit. 2008-09-10]. Dostupný z WWW: . Orbit Data Message. Blue Book [online]. 2004 [cit. 2008-09-10]. Dostupný z WWW: . Space Communications Protocol Specification (SCPS)—Network Protocol (SCPS-NP). Blue Book [online]. 1999 [cit. 2008-09-10]. Dostupný z WWW: . Space Communications Protocol Specification (SCPS)—Security Protocol (SCPS-SP). Blue Book [online]. 1999 [cit. 2008-09-10]. Dostupný z WWW: . Space Communications Protocol Specification (SCPS)—Transport Protocol. Blue Book [online]. 1999 [cit. 2008-09-10]. Dostupný z WWW: . Telemetry Summary of Concept and Rationale. Green Book [online]. 1987 [cit. 2008-09-10]. Dostupný z WWW: .
85
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí [ 47 ] Overview of Space Communications Protocols. Green Book [online]. 2007 [cit. 2008-09-10]. Dostupný z WWW: . [ 48 ] TM Synchronization and Channel Coding--Summary of Concept and Rationale. Silver book. [online]. 2006 [cit. 2008-09-10]. Dostupný z WWW: . [ 49 ] Telemetry Channel Coding. [online]. 2002 [cit. 2008-09-10]. Dostupný z WWW: . [ 50 ] Packet Telemetry. Silver book. [online]. 2000 [cit. 2008-09-10]. Dostupný z WWW: . [ 51 ] Packet Telemetry Services. Silver book. [online]. 1996 [cit. 2008-09-10]. Dostupný z WWW: . [ 52 ] Packet Telemetry Service Specification. Silver book. [online]. 2001 [cit. 2008-09-10]. Dostupný z WWW: . [ 53 ] Flexible Serially Concatenated Convolutional Turbo Codes with Near-Shannon Bound Performance for Telemetry Applications Orange book.. [online]. 2001 [cit. 2008-09-10]. Dostupný z WWW: . [ 54 ] EUTELSAT [online]. 2005 [cit. 30. 5. 2007]. Dostupný z WWW: <www.eutelsat.com>. [ 55 ] AMSAT - The Radio Amateur Satellite Corporation [online]. 2004 , 2008 [cit. 2008-02-10]. Dostupný z WWW: . [ 56 ] WALKE, B. Satelite-UMTS – Specification of Protocols and Traffic Performance Analysis. Aachener Beitrager zur Mobil- und Telekommunikation, 2005, ISBN 3-86130-176-8. [ 57 ] HENG-PING, X., SHENG, P., XUN, Z., MENG-TIAN, R., Puncturing Period on the Performance of Punctured Turbo Codes. Wuhan University of Natural Sciences, vol. 7, no. 3, 2002, 1007-1202. [ 58 ] ROWITCH, D.N.; MILSTEIN, L.B., On the performance of hybrid FEC/ARQ systems using rate compatible punctured turbo (RCPT) codes, Communications, IEEE Transactions on , vol.48, no.6, pp.948-959, Jun 2000. [ 59 ] ZHOU; G., LIN, T-S.; Wang; W., LINDSEY, W.C.; LAI, D.; CHEN, E.; SANTORU, J., On the concatenation of turbo codes and Reed-Solomon codes, Communications, 2003. ICC '03. IEEE International Conference on , vol.3, no., pp. 2134-2138 vol.3, 11-15 May 2003. [ 60 ] FERRARI, M.; OSNATO, F.; SITI, M.; VALLE, S.; BELLINI, S., Performance of concatenated Reed-Solomon and turbo codes with non ideal interleaving, Global Telecommunications Conference, 2001. GLOBECOM '01. IEEE , vol.2, no., pp.911-915 vol.2, 2001. [ 61 ] LIU; C.Y. TANG; H., LIN; S., FOSSORIER, M.P.C., An interactive concatenated turbo coding system, Vehicular Technology, IEEE Transactions on , vol.51, no.5, pp. 998-1010, Sept. 2002. [ 62 ] FERRARI, M., BELLINI, S., Low BER Turbo Codes for Sattelite Communications. [online]. 2001 [cit. 2009-04-02]. Dostupný z WWW: < http://conferences.esa.int/01C14>. [ 63 ] TEPE, K.E.; ANDERSON, J.B., Turbo coding behavior in Rayleigh fading channels without perfect interleaving, Military Communications Conference, 2001. MILCOM 2001. Communications for Network-Centric Operations: Creating the Information Force. IEEE , vol.2, no., pp. 1157-1164 vol.2, 2001. [ 64 ] WANG, D.; KOBAYASHI, H., On design of interleavers with practical size for turbo code s, Communications, 2000. ICC 2000. 2000 IEEE International Conference on , vol.2, no., pp.618-622 vol.2, 2000. [ 65 ] ASOODEH, S.; REZZAZE, A., Design of Optimal Period Interleaver in Turbo Codes, Information Theory, 2007. CWIT '07. 10th Canadian Workshop on , vol., no., pp.37-40, 6-8 June 2007. [ 66 ] STRINATI, E.C.; SIMOENS, S.; BOUTROS, J., Error Rate Estimation Based on Soft Output Decoding and its Application to Turbo Coding, Wireless Communications and Networking Conference, 2007.WCNC 2007. IEEE , vol., no., pp.72-76, 11-15 March 2007. [ 67 ] STRINATI, E.C.; SIMOENS, S.; BOUTROS, J., New error prediction techniques for turbo coded OFDM systems and impact on adaptive modulation and coding, Personal, Indoor and Mobile Radio Communications, 2005. PIMRC 2005. IEEE 16th International Symposium on , vol.2, no., pp.1116-1119 Vol. 2, 11-14 Sept. 2005. 86
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí [ 68 ] LETZEPIS, N.; GRANT, A., Bit Error Rate Estimation for Turbo Decoding, Communications, IEEE Transactions on , vol.57, no.3, pp.585-590, March 2009. [ 69 ] LOELIGER, H.-A., A posteriori probabilities and performance evaluation of trellis codes, Information Theory, 1994. Proceedings., 1994 IEEE International Symposium on , vol., no., pp.335-, 27 Jun-1 Jul 1994. [ 70 ] BAHL, L., COCKE, J., JELINEK, F., RAVIV, J., Optimal decoding of linear codes for minimizing symbol error rate (Corresp.), Information Theory, IEEE Transactions on , vol.20, no.2, pp. 284-287, Mar 1974. [ 71 ] LAND, I., HOEHER, P. A., Log-Likelihood Values and Monte Carlo Simulation – some fundamentals result. International symposium on Turbo Codes & Related Topics, Brest, France, 2000. [ 72 ] ZHAI, F., FAIR, I. J., New error detection techniques and stopping criteria for turbo decoding, Electrical and Computer Engineering, 2000 Canadian Conference on , vol.1, no., pp.58-62 vol.1, 2000. [ 73 ] El GAMAL, H., HAMMONS, A.R., Jr., Analyzing the turbo decoder using the Gaussian approximation, Information Theory, IEEE Transactions on , vol.47, no.2, pp.671-686, Feb 2001. [ 74 ] FORNEY, G. D., Maximum-Likelihood Sequence Estimation of Digital Sequences in the Presence of Intersymbol Interference, Information Theory, IEEE Transaction, vol.18, no.3, pp.363-378, May 1972. [ 75 ] Ericsson Inc., Plano, United States. Method and Devices for Calculating Bit Error Rate of Received Signal. Inventor: MOELKER, D.-J., Patent No.: US 2007/0162788 A1, Public date: Jul. 12, 2007. [ 76 ] VAN STRALEN, N.A., TAYLOR, D.P., Joint synchronization and demodulation of trellis codes, Communications, IEEE Transactions on , vol.42, no.6, pp.2268-2275, Jun 1994. [ 77 ] HERZET, C., WYMEERSCH, H., MOENECLAEY, M., VENDENDORPE, L., On MaximumLikelihood Timing Synchronization, Communications, IEEE Transactions on , vol.55, no.6, pp.1116-1119, June 2007. [ 78 ] GOLD, R., Optimal binary sequences for spread spectrum multiplexing (Corresp.), Information Theory, IEEE Transactions on , vol.13, no.4, pp. 619-621, Oct 1967. [ 79 ] FAZEL, K., KAISER, S. Multi-Carrier and Spread Spectrum Systems. England, John Wiley & Sons Ltd., p.298, ISBN 0-470-84899-5. [ 80 ] Vitthaladevuni, P.K., Alouini, M.-S., Effect of imperfect phase and timing synchronization on the bit-error rate performance of PSK modulations, Communications, IEEE Transactions on , vol.53, no.7, pp. 1096-1099, July 2005. [ 81 ] HRDINA, Z., VEJRAŽKA, F., Digitální rádiová komunikace. Praha, Vydavatelství ČVUT, 1998, ISBN 80-01-01820-2. [ 82 ] RUMÁNEK, J. Družicové komunikační systémy. Elektrorevue - Internetový časopis (http://www.elektrorevue.cz), 2007, roč. 1, č. 37, s. 1-13. ISSN: 1213-1539. [ 83 ] RUMÁNEK, J. Družicové komunikační systémy. Sdělovací technika, 2008, č. 2, s. 4-7, 0036-9942. [ 84 ] RUMÁNEK, J. Systémy pro zpracování signálů prostřednictvím zvukové karty. In Sborník příspěvků konference ZVŮLE 2007. Brno: Vysoké učení technické v Brně, Fakulta elektrotechniky a kumunikačních technologií, 2007. s. 164-167. ISBN 978-80-214-34684. [ 85 ] RUMÁNEK, J. Nové modulační techniky pro datovou družicovou komunikace s nízkou přenosovou rychlostí. In Seminář o řešení projektu GA ČR 102/03/H109 . Brno: Vysoké učení technické v Brně, Fakulta elektrotechniky a komunikačních technologií, Ústav radioelektroniky, 2007. s. 38 (1 s.). ISBN 978-80-214-3444-8. [ 86 ] RUMÁNEK, J.; ŠEBESTA, J. Nové techniky pro datovou družicovou komunikaci . In Pokročilé metody, struktury a komponenty elektronické bezdrátové komunikace. 2008. s. 1-3. ISBN 978-80-214-3753-1.
87
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí [ 87 ] RUMÁNEK, J.; ŠEBESTA, J. Generátor PSK31. In Kondor 2008. Trencianske Jastrabie, Slovakia: Institute of Radio Electronics, Faculty of Electrical Engineering and Communication, 2008. s. 161-163. ISBN 978-80-214-3394-6. [ 88 ] RUMÁNEK, J. Družicové komunikační systémy. Sdělovací technika, 2008, roč. 1, č. 2, s. 37. ISSN: 0036-9942. [ 89 ] RUMÁNEK, J. Družicové komunikační systémy. Elektrorevue - Internetový časopis (http://www.elektrorevue.cz), 2007, roč. 1, č. 37, s. 1-13. ISSN: 1213-1539. [ 90 ] RUMÁNEK, J. Systémy pro zpracování signálů prostřednictvím zvukové karty . In Sborník příspěvků konference ZVŮLE 2007. Brno: Vysoké učení technické v Brně, Fakulta elektrotechniky a kumunikačních technologií, 2007. s. 164-167. ISBN 978-80-214-34684. [ 91 ] RUMÁNEK, J. Nové modulační techniky pro datovou družicovou komunikace s nízko u přenosovou rychlostí. In Seminář o řešení projektu GA ČR 102/03/H109 . Brno: Vysoké učení technické v Brně, Fakulta elektrotechniky a komunikačních technologií, Ústav radioelektroniky, 2007. s. 38 (1 s.). ISBN 978-80-214-3444-8. [ 92 ] RUMÁNEK, J.; ŠEBESTA, J. The SMS Satellite System Design. In Proceedings of the 19th International conference Radioelektronika 2009. Božetěchova 133, Brno, 612 00: MJ Servis, 2009. s. 259-261. ISBN 978-80-214-3865-1.
88
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Seznam zkratek AMSAT
Radio amatérská satelitní organizace (The Radio Amateur Satellite Corporation)
APP
A posteriori pravděpodobnost (A posteriori probability)
AWGN
Aditivní bílý šum s normálním rozložením (Aditive White Gaussian Noise)
BER
Bitová chybovost (Bit Error Rate)
BPSK
Dvoustavové fázové klíčování (Binary Phase Shift Keying)
CDMA
Kódový multiplex (Code Division Multiple Access)
DMC
Diskrétní bezpaměťový kanál (Discrete Memoryless Channel)
DP
Dopplerův posuv frekvence
EIRP
Ekvivalentní isotropický vyzářený výkon (Equivalent Isotropically Radiated Power)
FEC
Dopředná korekce chyb (Forward Error Correction)
GEO
Geostacionární eliptická dráha (Geostacionar Earth Orbit)
GPS
Globální navigační systém (Global Positioning System)
HEO
Vysoce eliptická dráha (Highly Elliptical Orbit)
ICO
Střední kruhová dráha (Intermediate Circular Orbit)
LEO
Nízká kruhová dráha (Low Earth Orbit)
LLR
Logaritmický poměr pravděpodobností (Logaritmic Likelihood Ratio)
MAP
Maximum a posteriori
ML
Maximální pravděpodobnost (Maximum Likelihood)
NASA
Národní úřad pro letectví a kosmonautiku (National Aeronautics and Space Administration)
NORAD
Severoamerické velitelství protivzdušné obrany (North American Aerospace Defense Command)
RSC
Rekurzivní systematický konvoluční kodér (Recursive Systematic Convolution Code)
SISO
Měkký vstup – měkký výstup (Soft Input Soft Output)
SMS
Krátká textová zpráva (Short Message Service)
89
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
Seznam symbolů úhel mezi vektorem družice
a směrem k pozorovateli
pole zpětné stavové metriky (dk)
log. pravděpodobnost pro bit dk, MAP dekódovací algoritmus
2
disperze kanálu
C
výkon nosné vlny
c
rychlost šíření světla
c1, c2
kódovaná posloupnost
D
průměr paraboly
dk
k-tý informační bit
e
Eulerovo číslo
Eb
energie jednoho bitu
f
frekvence
F
šumové číslo
G
vytvářecí matice
g
generující polynom
G(x)
vytvářecí mnohočlen
GR
zisk přijímací antény
GT
zisk vysílací antény
I
jednotková matice
k
délka bloku dat
k
počet informačních bitů
L(uk)
logaritmický a posteriori poměr pravděpodobností
L0
útlum kanálu
LT
přídavné ztráty
M
počet symbolů
N
výkon šumu
n
délka posloupnosti
N0
spektrální hustota šumu
P
pravděpodobnost 90
Nové techniky pro datovou družicovou komunikaci s nízkou přenosovou rychlostí
PRX
úroveň přijatého signálu
PT
výkon vysílače
Q
počet stavů
r
počet paritních bitů
R
kódový poměr
R1N
přijatá sekvence (pozice 1 až N)
S
výkon signálu
s
systematická část
Sk
stavy kodéru v k-té pozici
Ts
šumová teplota antény
TSRX
ekvivalentní šumová teplota systému vztažená ke vstupu přijímače
u
vstupní posloupnost
uk
k-tý dekódovaný bit
Um
jedna z možných kódových sekvencí
y
přijatá bitová sekvence
α
pole dopředné stavové metriky
αk(m’,m) pravděpodobnost přechodu mezi stavy m‘, m v k-té pozici γk(m)
pravděpodobnost bytí ve stavu m na pozici k
π
Ludolfovo číslo
91
Curriculum Vitae Name:
Jaroslav RUMÁNEK
Born:
August 6th 1982 in Vsetín
Contact:
[email protected]
Phone:
+420 774 545 915
Website:
www.foto-rumanek.cz
Education 2001 – 06
Technical University of Brno / Department of Radio Electronics Pre-graduate study of Radio Electronics State exam passed in June 2006 Diploma thesis Elaboration of Algorithm for Computer Aided Optimization of Band-Pass Networks in June 2006
2006 – 09
Technical University of Brno / Department of Radio Electronics Ph.D. study of Electronics
Experience 1/08 – 12/08 Development of software for Zoning measurement system (C++ Builder). 1/01 – 8/09 Development of websites (graphic design, PHP, HTML). Languages English, German Other activities Programming (C++ Builder, Matlab, Matlab Simulink, PHP, HTML, MySQL, XML, Pascal, Assembler) Photography Graphic design