25.07.2013 Views

Diskreetti matematiikka ja ekvalisointi

Diskreetti matematiikka ja ekvalisointi

Diskreetti matematiikka ja ekvalisointi

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 1<br />

DISKREETIN MATEMATIIKAN SOVELLUKSIA:<br />

KANAVA-EKVALISOINTI TIEDONSIIRROSSA<br />

Jukka Talvitie, Toni Levanen & Mikko Valkama<br />

TTY / Tietoliikennetekniikka<br />

jukka.talvitie@tut.fi, toni.levanen@tut.fi, mikko.e.valkama@tut.fi<br />

Tässä oleva esitys poh<strong>ja</strong>utuu mm. ao. kurssien sisältöön:<br />

TLT-5400 Digitaalinen siirtotekniikka<br />

TLT-5906 Digitaalisen siirtotekniikan <strong>ja</strong>tkokurssi<br />

Tarkoituksena on antaa esimerkkejä diskreetin matematiikan<br />

merkityksestä modernissa tietoliikennetekniikassa keskittyen lähinnä<br />

kanavakor<strong>ja</strong>ukseen eli <strong>ekvalisointi</strong>in.<br />

Taustaa<br />

Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 2<br />

Digitaalisen tiedonsiirron perusidea on yksinkertaisesti siirtää joukko<br />

bittejä (binäärinen lukujono) paikasta toiseen sähköisten tai<br />

sähkömagneettisten signaalien avulla. Huolimatta siirrettävän<br />

informaation diskreetistä olomuodosta, itse tiedonsiirtokanavaan<br />

lähetettävä aaltomuoto on fysikaalisista rajoitteista johtuen <strong>ja</strong>tkuvaaikainen.<br />

Oikealla tavalla toteutetun näytteenoton avulla <strong>ja</strong>tkuvaaikainen<br />

aaltomuoto saadaan kuitenkin palautettua diskreetiksi<br />

lukujonoksi vastaanottimessa ilman merkittäviä häviöitä.<br />

Ennen digitaalisten tiedonsiirtomenetelmien käyttöönottoa havaittua<br />

<strong>ja</strong>tkuva-aikaista signaalia käsiteltiin vastaanottimessa yksinomaan<br />

analogisin komponentein. Tällä tavoin signaalin hallinta vaikeutuu<br />

olennaisesti, sillä jokainen signaalille tehty prosessointitoimenpide vaatii<br />

periaatteessa oman erillisen sähköpiirin. Lisäksi analogiset<br />

sähkökomponentit ovat kalliita <strong>ja</strong> isokokoisia, minkä vuoksi niiden<br />

kaupallinen houkuttelevuus varsinkin nykymaailmassa on hyvin heikko.<br />

Tässä piileekin digitaalisen siirtotekniikan voimavara, sillä<br />

vastaanotetun <strong>ja</strong>tkuva-aikaisen signaalin näytteistetty versio on vain<br />

joukko luku<strong>ja</strong>, joiden hallinta mikropiireillä on hyvin tehokasta. Luku<strong>ja</strong><br />

voidaan helposti esimerkiksi tallentaa muistiin <strong>ja</strong> palauttaa ne sieltä<br />

myöhemmin häviöttömästi <strong>ja</strong>tkoprosessointia varten. Juuri tämän vuoksi<br />

diskreetin matematiikan tuottamat sovellukset ovat erittäin tärkeitä<br />

tietoliikennetekniikassa.<br />

Digitaalitekniikka on nykyään käytössä lähes kaikissa<br />

kulutta<strong>ja</strong>sovelluksissa kuten matkapuhelinverkoissa, yleisradio <strong>ja</strong> -TV<br />

lähetyksissä, kotien langattomissa Internet-yhteyksissä, navigaattori- <strong>ja</strong><br />

paikannuspalveluissa sekä yleisesti kodin sisäisessä tiedonsiirrossa.<br />

Jokaisessa näissä tiedonsiirtokanava aiheuttaa lähetettyyn signaalin<br />

vääristymää. Kanavaekvalisoinnilla tarkoitetaan tämän vääristymän<br />

kor<strong>ja</strong>amista <strong>ja</strong> sitä käytetään <strong>ja</strong>tkossa tämän esityksen esimerkkinä<br />

diskreetin matematiikan sovelluksesta digitaalisessa tiedonsiirrossa.


Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 3<br />

Tiedonsiirtojärjestelmän rakenne <strong>ja</strong> siirtokanavan vaikutus<br />

Kantataajuisen (nollataajuuden ympäristössä si<strong>ja</strong>itsevan) digitaalisen<br />

siirtojärjestelmän periaatteellinen lohkokaavio on esitetty alla olevassa<br />

kuvassa.<br />

Ensimmäisessä vaiheessa lähetettävä bittijono muutetaan symboleiksi.<br />

Tässä joukko peräkkäisiä bittejä muutetaan tietyksi symboliaakkoston<br />

määräämäksi lukuarvoksi. Alla on esimerkin vuoksi esitetty 4-tasoisen<br />

reaalisen symboliaakkoston rakenne:<br />

Esimerkki 4-tasoisesta<br />

symboliaakkostosta<br />

Bittiyhdistelmä Symboli<br />

00 -3<br />

01 -1<br />

10 1<br />

11 3<br />

Pääasiallinen tarkoitus bittien kuvaamisessa symboleiksi on kasvattaa<br />

järjestelmän taajuuskaistan käytön tehokkuutta, sillä<br />

signaalikohinasuhteen <strong>ja</strong> kanavan vaikutusten lisäksi digitaalisen<br />

siirtojärjestelmän tiedonsiirtokapasiteettiin vaikuttaa ainoastaan käytetty<br />

taajuuskaistanleveys. Tämä puolestaan määräytyy siitä, kuinka monta<br />

symbolia (diskreettiä lukuarvoa) aikayksikköä kohden halutaan siirtää.<br />

Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 4<br />

Jos esimerkiksi käytettävää taajuuskaistaa on 100 Hz, niin teoriassa<br />

100 symbolia voidaan siirtää sekunnissa. Tällöin siirrettäessä bittejä<br />

”sellaisenaan” (binäärinen symboliaakkosto), saavutettava bittinopeus<br />

on 100 bittiä/s, kun taas edellä esitetyn 4-tasoisen symboliaakkoston<br />

avulla saavutetaan samalla kaistalla bittinopeus 200 bittiä/s. Mitä<br />

suurempi symboliaakkosto on, sitä enemmän bittejä yhtä lähetettyä<br />

symbolia kohden voidaan siirtää. Toisaalta, olettaen lähetysteho<br />

kiinnitetyksi, suuremman symboliaakkoston huonona puolena on<br />

suurempi herkkyys kohinan <strong>ja</strong> häiriöiden aiheuttamille virheille.<br />

Diskreetit symbolit muunnetaan <strong>ja</strong>tkuva-aikaiseksi signaaliksi<br />

lähetinsuodattimen avulla. Tässä valittua (<strong>ja</strong>tkuva-aikaista)<br />

pulssimuotoa g(t) painotetaan diskreeteillä symbolien arvoilla, jolloin<br />

lähetettävä signaali on muotoa<br />

•<br />

Â<br />

S() t = A g( t -kT)<br />

k =-•<br />

k<br />

missä T on kahden peräkkäisen symbolin välinen aika. Alla olevassa<br />

kuvassa on esimerkki käytettäessä pulssimuotona kanttipulssia<br />

(huomaa kuvassa myös yleiset kriteerit käytettävälle pulssimuodolle):<br />

Kanttipulssin si<strong>ja</strong>sta käytännöllisempi ratkaisu on käyttää ”pyöreämpiä”<br />

pulssimuoto<strong>ja</strong>, kuten ns. nostettu<strong>ja</strong> kosinipulsse<strong>ja</strong>. Erityisen olennaista<br />

on kuitenkin ymmärtää edellä esitettyjen pulssimuotokriteerien merkitys.<br />

Koska näytehetkellä pulssi saa aina arvon yksi, tietyllä<br />

symbolia<strong>ja</strong>nhetkellä signaalista S(t) otettu näytteen arvo on aina samalla<br />

a<strong>ja</strong>nhetkellä lähetetyn symbolin Ak arvo. Lisäksi koska pulssin arvo saa<br />

aina arvon nolla muilla symbolia<strong>ja</strong>nhetkillä, eri pulssit/symbolit eivät<br />

häiritse toisiaan.


Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 5<br />

Lähetetty <strong>ja</strong>tkuva-aikainen signaali S(t) muokkautuu kanavassa <strong>ja</strong><br />

kohina summautuu vääristyneen aaltomuodon päälle. Tämän jälkeen<br />

vastaanottosuodatuksessa käytössä olevan taajuuskaistan ulkopuolinen<br />

kohina suodatetaan pois <strong>ja</strong> signaalista otetaan näytteitä<br />

symbolia<strong>ja</strong>n<strong>ja</strong>kson T välein.<br />

Kuten jo mainittiin, näytteistys poimii saapuvasta signaalista vain tietyt<br />

a<strong>ja</strong>nhetket: kT (k=-…). Toisin sanoen vastaanotettua signaalia<br />

tarkastellaan vain symbolia<strong>ja</strong>nhetkillä, minkä vuoksi edellä kuvatusta<br />

lohkokaaviosta voidaan kehittää efektiivinen diskreettiaikainen malli:<br />

Tässä lohkokaaviossa ei enää esiinny <strong>ja</strong>tkuva-aikaisia signaaleita.<br />

Analyyttisissa tarkasteluissa tämä malli on kuitenkin riittävä, mikäli<br />

voidaan olettaa, että pulssimuoto toteuttaa edellä mainitut kriteerit <strong>ja</strong><br />

näytteenotto onnistuvat ideaalisesti. Nyt vastaanotettua näytejonoa<br />

voidaan kuvata diskreetillä tavalla seuraavasti:<br />

Rk = Ak * pk + Zk<br />

missä pk on kanavan (näytteistetty) diskreettiaikainen impulssivaste <strong>ja</strong> Zk<br />

summautuvaa kohinaa. Merkki * kuvaa konvoluutiota, jonka avulla<br />

saadaan laskettua lineaaristen järjestelmien input-output-vasteita (kts.<br />

Sitikka-materiaali ”Integroinnin sovelluksia tiedonsiirtotekniikassa”).<br />

Jos kanavan impulssivaste ulottuu useamman symboliaikavälin T<br />

alueelle, niin kanavassa syntyy symbolien välistä keskinäisvaikutusta:<br />

ISI (Inter Symbol Interference). Esimerkiksi kanava, jonka diskreetti<br />

impulssivaste on muotoa pk=-0.4(k+1)+(k)+0.8(k-1) ((t) on<br />

yksikköimpulssifunktio) voidaan esittää graafisesti seuraavaan tapaan:<br />

Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 6<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

-0.2<br />

Kanavan impulssivaste p k<br />

-0.4<br />

-1 -0.8 -0.6 -0.4 -0.2 0<br />

k<br />

0.2 0.4 0.6 0.8 1<br />

Kyseisen kanavan tapauksessa vastaanotettu näytejono voidaan nyt<br />

kirjoittaa muodossa<br />

R = A * p + Z = - 0.4A + A + 0.8A<br />

+ Z<br />

k k k k k- 1 k k+ 1 k<br />

Näin ollen jokainen lähetetty symboli nähdään vastaanottimessa<br />

painotettuna kolmen peräkkäisen symbolin summana. Taajuustasossa<br />

tämä näkyy spektrin vääristymisenä, jossa eri taajuudet vaimenevat<br />

toisiinsa nähden eri tavalla.<br />

Ekvalisaattorin tehtävänä on pienentää tätä kanavassa syntyvää ISI:ä.<br />

Ongelmana on siis löytää sellainen diskreettiaikainen suodatin (ts.<br />

suodattimen kertoimet/tapit …,ck-1, ck, ck+1,…), jonka avulla ISI<br />

minimoituu. Aikatasossa tämä tarkoittaa järjestelmän kokonaisvasteen<br />

(kanava+ekvalisaattori) pakottamista lähelle yksikköimpulssin vastetta,<br />

kun taas taajuustasossa järjestelmän tuottama amplitudispektri pyritään<br />

saamaan mahdollisimman tasaiseksi amplitudiarvon 1 ympäristöön.<br />

Erilaisia <strong>ekvalisointi</strong>menetelmiä löytyy kir<strong>ja</strong>llisuudesta useita, joista ns.<br />

zero-forcing-menetelmä esitellään seuraavilla sivuilla.


Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 7<br />

Lineaarinen <strong>ekvalisointi</strong> (zero-forcing-menetelmä)<br />

Ehkä yksinkertaisin <strong>ja</strong> intuitiivisin menetelmä ISI:n poistamiseen on ns.<br />

zero-forcing-ekvalisaattori. Tässä lähtökohtana on etsiä sellainen<br />

diskreettiaikainen ra<strong>ja</strong>tun kestoinen (Finite Impulse Response)<br />

suodatin, joka pyrkii poistamaan ISI:n kokonaan. Tämän ongelman<br />

analysointi on huomattavan paljon helpompaa suorittaa taajuustasossa,<br />

sillä siellä konvoluutio voidaan esittää kertolaskuna. Muuntamalla<br />

aikatason vasteet z-muunnoksella taajuustason siirtofunktioksi,<br />

vastaanotetut näytteet voidaan esittää seuraavasti:<br />

Rz () = PzAz ()() + Nz ()<br />

missä P(z), A(z) <strong>ja</strong> N(z) ovat kanavan impulssivasteen, symbolijonon <strong>ja</strong><br />

kohinan z-muunnokset. Ekvalisoidut näytteet saadaan tällöin z-tasossa<br />

ilmaistua<br />

Qz () = Cz ()( PzAz ()() + Nz ())<br />

= CzPzAz () ()() + CzNz () ()<br />

Tästä nähdään suoraan, että mikäli kanavan vaikutus halutaan<br />

vastaanotetusta signaalista kokonaan poistaa, tulee ekvalisaattorin<br />

siirtofunktio valita siten, että C(z)P(z)=1. Toisin sanoen ekvalisaattorin<br />

siirtofunktio on muotoa<br />

Cz () = 1/ Pz ()<br />

Termi zero-forcing juontaakin juurensa nimenomaan tästä<br />

lähestymistavasta, jossa ISI ikään kuin ”pakotetaan nollaan”.<br />

Ekvalisoidut symbolit voidaan nyt esittää muodossa<br />

Qz () = CzPzAz () ()() + CzNz () ()<br />

1 1<br />

= PzAz ()() + Nz ()<br />

Pz () Pz ()<br />

Nz ()<br />

= Az () +<br />

Pz ()<br />

Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 8<br />

Tarkastellaan seuraavaksi esimerkkitilannetta, jossa edellä esitettyyn 3tappiseen<br />

kanavaan lähetetään symbolit 1, -1, 1 <strong>ja</strong> 1. Alla olevassa<br />

kuvassa on esitetty lähetetyn sekvenssin aikatason esitykset<br />

siirtojärjestelmän eri vaiheissa.<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

-0.2<br />

-0.4<br />

-0.6<br />

-0.8<br />

-1<br />

Lähetetty sekvenssi<br />

-1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 4<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

-0.2<br />

-0.4<br />

-0.6<br />

-0.8<br />

-1<br />

Vastaanotettu sekvenssi<br />

-1 -0.5 0 0.5 1 1.5<br />

k<br />

2 2.5 3 3.5 4<br />

1.2<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

-0.2<br />

Kanavan impulssivaste p(k)<br />

-0.4<br />

-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1<br />

Systeemin kokonaisvaste<br />

-0.2<br />

-5 -4 -3 -2 -1 0<br />

k<br />

1 2 3 4 5<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

-0.2<br />

-0.4<br />

2<br />

1.5<br />

1<br />

0.5<br />

0<br />

-0.5<br />

Kanavan ulostulo<br />

-1<br />

-1 -0.5 0 0.5 1 1.5<br />

k<br />

2 2.5 3 3.5 4<br />

Ekvalisaattorin vaste<br />

-0.6<br />

-4 -3 -2 -1 0 1 2 3 4<br />

Kanavan ulostulossa nähdään selvästi, että alkuperäinen sekvenssi on<br />

voimakkaasti vääristynyt ISI:n vaikutuksesta. Itse asiassa, jos<br />

päätökset/arvaukset lähetetystä sekvenssistä tehtäisiin tähän signaaliin<br />

perustuen, symbolien ilmaisussa tapahtuisi luultavimmin virhe, sillä 3.<br />

lähetetty symboli on lähempänä -1:ä kuin 1:ä. Käytetty ekvalisaattori on<br />

tässä 9-tappinen eli sen impulssivaste on muotoa c-4, c-3,…c0,…, c3, c4.<br />

Ekvalisaattorin kertoimet on laskettu zero-forcing-menetelmään<br />

perustuen siten, että C(z)=1/P(z). Systeemin kokonaisvaste ( pk ck<br />

* tai<br />

z-tasossa P(z)C(z)) on kuvan perusteella hyvin lähellä<br />

yksikköimpulssifunktiota, joten ISI on selvästi pienentynyt. ISI saadaan<br />

poistettua sitä tarkemmin, mitä enemmän ekvalisaattoriin sisällytetään<br />

tappe<strong>ja</strong>. Tappien lukumäärä kasvattaa kuitenkin järjestelmän<br />

laskennallista taakkaa, minkä vuoksi niiden määrää joudutaan aina<br />

tapauskohtaisesti rajoittamaan.<br />

Aikatason tarkastelun lisäksi tilannetta voidaan havainnollistaa myös<br />

taajuustasossa. Taajuustason esitys saadaan ottamalla Fourier-


Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 9<br />

muunnokset (kts. Sitikka-materiaali ”Integroinnin sovelluksia<br />

tiedonsiirtotekniikassa”) edellisessä kuvassa esitetyistä aikatason<br />

vasteista. Alla oleviin kuviin on piirretty kanavan, ekvalisaattorin <strong>ja</strong><br />

järjestelmän kokonaisvasteen amplitudispektrit.<br />

Amplitudivaste<br />

Taajuusvasteet kanavalle <strong>ja</strong> ekvalisaattorille<br />

1.6<br />

1.4<br />

1.2<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

Kanavan vaste<br />

Ekvalisaattorin vaste<br />

0<br />

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1<br />

normalsoitu taajuus<br />

Amplitudivaste<br />

1.6<br />

1.4<br />

1.2<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

Systeemin kokonaisvaste<br />

0<br />

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1<br />

normalsoitu taajuus<br />

Kuvista havaitaan, että ekvalisaattorin vaste on käänteinen verrattuna<br />

kanavan vasteeseen. Kanavassa vaimentuneita taajuuksia voimistetaan<br />

<strong>ja</strong> päinvastoin, minkä vaikutuksesta kokonaisvaste ”tasoittuu”<br />

amplitudiarvon 1 ympäristöön.<br />

Vaikka zero-forcing-menetelmä intuitiivisesti vaikuttaakin erittäin<br />

järkevältä, siihen liittyy kuitenkin ikävä kohinan voimistumisilmiö.<br />

Ekvalisaattori voimistaa signaalia niillä taajuuksilla, joissa kanavan<br />

vaimennus on suuri. Tällöin, itse hyötysignaalin lisäksi, ekvalisaattori<br />

voimistaa tarpeettomasti myös kohinaa. Tämä nähdään selkeästi jo<br />

edellä esitetystä kaavasta (ekvalisoidun sekvenssin z-muunnos):<br />

Nz ()<br />

Qz () = Az () +<br />

Pz ()<br />

Kyseisen ilmiön vuoksi ISI:n poistoon käytetäänkin usein muunlaisia<br />

menetelmiä, kuten esimerkiksi MSE-ekvalisaattoria (Minimum Square<br />

Error), jossa ISI:ä ei aluperin pyritäkään kohinan tehosta riippuen<br />

poistamaan aivan kokonaan. Eräs toinen mahdollisuus on käyttää ns.<br />

sekvenssi-ilmaisua, jossa kanavan aiheuttama ISI otetaan huomioon<br />

itse symbolien ilmaisuprosessissa (esim. Viterbi-algoritmi), jolloin<br />

varsinainen <strong>ekvalisointi</strong> ennen ilmaisua jää tarpeettomaksi.<br />

Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 10<br />

Tuntemattoman <strong>ja</strong>/tai muuttuvan kanavan hallinta<br />

Zero-forcing-menetelmässä ekvalisaattorin kertoimet määritetään<br />

suoraan kanavan vasteen perusteella. Yleensä kanavan vaste ei ole<br />

kuitenkaan tunnettu, minkä vuoksi zero-forcing-menetelmän käyttö<br />

sellaisenaan ei ole mahdollista. Lisäksi esimerkiksi langattomassa<br />

tiedonsiirrossa kanavan vaste muuttuu <strong>ja</strong>tkuvasti, minkä vuoksi<br />

ekvalisaattorin tappikertoimia joudutaan <strong>ja</strong>tkuvasti päivittämään. Tähän<br />

tarkoitukseen adaptiivinen <strong>ekvalisointi</strong> tarjoaa laskennallisesti tehokkaat<br />

työkalut. Alla olevassa kuvassa on esitetty lineaarisen adaptiivisen<br />

ekvalisoinnin periaatteellinen lohkodiagrammi.<br />

Kohinan summauksen jälkeen vastaanotin havaitsee kanavassa<br />

vääristyneen symbolisekvenssin Rk. Perusideana adaptiivisessa<br />

ekvalisoinnissa on minimoida virhe ekvalisoidun sekvenssin Qk <strong>ja</strong><br />

oikean sekvenssin Ak välillä. Käytännössä tämä onnistuu esimerkiksi<br />

ennalta määrätyn, vastaanottimessa tunnetun, pilottisignaalin avulla.<br />

Näin ollen tiedetään minkälainen ekvalisoidun sekvenssin Qk tulisi olla<br />

<strong>ja</strong> siten ekvalisoinnista aiheutuva virhe voidaan laskea suoraan<br />

pilottisignaalin <strong>ja</strong> Qk:n erotuksella. Pilottisignaalin si<strong>ja</strong>sta vastaavana<br />

referenssisignaalina voidaan käyttää myös ekvalisoidusta signaalista<br />

tehtyjä symbolipäätöksiä, mikäli voidaan olettaa, että päätökset ovat<br />

enimmäkseen oikeita (esim. yli 90% todennäköisyydellä).<br />

Minimoitaessa keskimääräistä (neliö)virhettä päädytään usein melko<br />

suuriin lineaarisiin yhtälöryhmiin, joiden ratkaisussa tarvittava matriisi-


Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 11<br />

inverssi saattaa olla laskennallisesti ottaen epäkäytännöllinen. Jos<br />

edelleen pidetään mielessä, että kanava muuttuu koko a<strong>ja</strong>n, tappe<strong>ja</strong><br />

joudutaan <strong>ja</strong>tkuvasti päivittämään <strong>ja</strong> laskettavien matriisi-inverssien<br />

määrä aikayksikköä kohden vain kasvaa. Tähän ongelmaan tehokkaan<br />

ratkaisun antavat erilaiset iteratiiviset laskentamenetelmät, joista LMSalgoritmi<br />

(Least Mean Squares) on eräs tunnetuimmista. Tässä<br />

ekvalisaattorin tappe<strong>ja</strong> päivitetään jokaisen ekvalisaattoriin saapuvan<br />

näytteen perusteella, minkä vuoksi <strong>ja</strong>tkuvat matriisi-inverssit jäävät<br />

tarpeettomiksi <strong>ja</strong> laskennallinen taakka kevenee.<br />

LMS-algoritmi perustuu ns. gradienttialgoritmiin, jossa neliövirhe<br />

minimoidaan pyrkimällä kohti gradientin nollakohtaa (minimikohta)<br />

iteratiivisesti. Tämä tapahtuu siirtymällä aina pieni askel kerrallaan kohti<br />

negatiivisen gradientin suuntaan kunnes nollakohta saavutetaan.<br />

Alkuperäinen gradientti-algoritmi tarvitsee toimiakseen tiedot<br />

vastaanotetun signaalin tilastollisista ominaisuuksista: näytteiden<br />

autokorrelaatiomatriisi <strong>ja</strong> ristikorrelaatio referenssisymbolin välillä.<br />

Nämä ovat kuitenkin yleensä tuntemattomia. LMS-algoritmi eroaa<br />

gradienttialgoritmista juuri tässä mielessä, sillä se arvioi nämä<br />

tilastolliset ominaisuudet ”hetkellisesti” perustuen suoraan<br />

vastaanotettuihin näytteisiin. LMS-algoritmin voidaan toiminta voidaan<br />

lyhyesti esittää seuraavalla iteratiivisella prosessilla:<br />

*<br />

k+ 1 = k+ bEkk<br />

c c r , missä<br />

T<br />

k = k -cr<br />

k k<br />

E A<br />

c<br />

r<br />

k -L<br />

0 L<br />

(laskettu virhe)<br />

= [ c ,..., c ,..., c ] (ekvalisaattorin tapit a<strong>ja</strong>nhetkellä k)<br />

= [ R ,..., R ,..., R ] (vastaanotettu sekvenssi)<br />

k k+ L k k-L Tässä on ns. askelparametri, joka määrittää kuinka suuri askel kohti<br />

arvioitua negatiivisen gradientin suuntaan otetaan. Jos on liian suuri,<br />

algoritmi muuttuu epästabiiliksi (”hajoaa käsiin”), <strong>ja</strong> toisaalta, jos on<br />

liian pieni, algoritmi ei ehdi konvergoitua pilottisekvenssin aikana. Alla<br />

olevassa kuvassa on esitetty LMS algoritmin toiminnallinen rakenne<br />

tarkasteltaessa ekvalisaattorin l:nen tapin [ck]l päivitysprosessia:<br />

Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 12<br />

Tutkitaan seuraavaksi tilannetta, jossa käytetään LMS-algoritmia tutun<br />

3-tappisen kanavan (kts. sivu 5) <strong>ekvalisointi</strong>in. Määritetään<br />

vastaanotetun sekvenssin signaali-kohinatehosuhteeksi 20dB <strong>ja</strong><br />

käytetään algoritmissa 9-tappista suodatinta sekä askelpituutta =0.002.<br />

Alla olevaan kuvaan vasemmalla on havainnollistettu absoluuttisen<br />

virheen |Ek| käyttäytyminen iteraatioiden edetessä. Kohinasta johtuen<br />

virheen käyrä ei ole tasainen mutta konvergoituu selkeästi nollan<br />

ympäristöön. Oikean puoleisessa kuvassa ekvalisaattorin eri tappien<br />

(siis 9 kpl yhteensä) lukuarvot on esitetty iteraatioiden funktiona.<br />

3.5<br />

3<br />

2.5<br />

2<br />

1.5<br />

1<br />

0.5<br />

Ekvalisaattorin virhe<br />

0<br />

0 100 200 300 400 500 600 700 800 900 1000<br />

iteraationumero<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

-0.2<br />

-0.4<br />

Ekvalisaattorin tappikertoimet<br />

-0.6<br />

0 100 200 300 400 500 600 700 800 900 1000<br />

iteraationumero<br />

Jatketaan vielä esimerkkiä <strong>ja</strong> tarkastellaan ekvalisoitu<strong>ja</strong> symboleita Qk iteraatioiden edetessä. Käytetään lähetykseen satunnaisia symboleita,<br />

jotka valitaan joukosta Ak Œ{ -3, - 1,1,3} . Alla oleviin kuviin on<br />

havainnollistettu tilannetta ilman <strong>ekvalisointi</strong>a (vasemmalla) <strong>ja</strong><br />

ekvalisoinnin kanssa (oikealla).


Näytearvo<br />

6<br />

4<br />

2<br />

0<br />

-2<br />

-4<br />

-6<br />

Ekvalisoimattomat näytteet a<strong>ja</strong>n funktiona<br />

1000 2000 3000 4000 5000 6000 7000 8000 9000<br />

k<br />

Diskreetin matematiikan excursio: kanava-<strong>ekvalisointi</strong> tiedonsiirrossa / 13<br />

Näytearvo<br />

6<br />

4<br />

2<br />

0<br />

-2<br />

-4<br />

-6<br />

Ekvalisoidut näytteet a<strong>ja</strong>n funktiona<br />

500 1000 1500 2000 2500 3000 3500 4000<br />

Vasemmanpuoleisesta kuvasta on selvää, ettei symboli-ilmaisin pysty<br />

ilman <strong>ekvalisointi</strong>a tuottamaan luotettavia symbolipäätöksiä.<br />

Oikeanpuoleisessa kuvassa LMS-algoritmin avulla ilmaisimelle tulevat<br />

näytearvot vaikuttavat jo huomattavasti luotettavimmilta. Alussa on<br />

havaittavissa algoritmin ns. oppimis<strong>ja</strong>kso, jonka aikana ekvalisaattorin<br />

tapit konvergoituvat ”oikeisiin” arvoihinsa. Kuten ehkä havaita<br />

saattaa(?), alkuarvoina ekvalisaattorin kaikille tapeille on tässä annettu<br />

nolla-arvo.<br />

Lopuksi on vielä syytä huomauttaa, että vaikka kaikki edellä annetut<br />

esimerkit poh<strong>ja</strong>utuvatkin reaalilukuihin, samat metodit ovat yhtä lailla<br />

käytettävissä myös kompleksilukujen tapauksessa. Itse asiassa<br />

useimmat modernit tiedonsiirtojärjestelmät käyttävät nimenomaan<br />

kompleksiarvoista symboliaakkostoa. Tämän vuoksi ekvalisaattorin<br />

toteutuksen täytyy myös olla kompleksinen (kts. sitikka-materiaali<br />

”Kompleksiluvut <strong>ja</strong> radiosignaalit”). Käytännössä kompleksiset<br />

suodattimet voidaan kuitenkin aina toteuttaa käyttämällä kahta<br />

rinnakkaista reaaliarvoista suodatinta.<br />

k

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!