Podręcznik

Strona: SEZAM - System Edukacyjnych Zasobów Akademickich i Multimedialnych
Kurs: Transmisja analogowa i cyfrowa
Książka: Podręcznik
Wydrukowane przez użytkownika: Gość
Data: czwartek, 3 kwietnia 2025, 06:36

1. Analogowa modulacja amplitudy

Analogowa modulacja amplitudy

1.1. Transmisja analogowa – kryteria jakości

Aby przesłać sygnał analogowy (mowa, muzyka) można wykorzystać analogowe lub cyfrowe łącze transmisyjne. To ostatnie wymaga przetworzenia sygnału analogowego na strumień binarny drogą próbkowania i kwantyzacji. Schemat analogowego i cyfrowego łącza transmisyjnego podano w Module 1, jednak zostanie on tu przypomniany dla wygody czytelnika. Modulator w systemie analogowym (Rys.1) służy do przeniesienia sygnału modulującego m(t) w odpowiednie pasmo częstotliwości, udostępniane przez kanał transmisyjny. W przypadku połączenia kablowego można zastosować transmisję w paśmie podstawowym bez modulatora, jednak i wówczas modulacja staje się niezbędna, jeśli udostępniamy kabel kilku użytkownikom. Sygnał zmodulowany s(t) zajmuje najczęściej szersze pasmo częstotliwości niż sygnał modulujący, wyjątkiem jest jednowstęgowa modulacja amplitudy. Kanał transmisyjny będziemy modelować jak filtr o paśmie B, przenoszący wszystkie składowe częstotliwościowe zawarte w sygnale s(t). Ponadto uwzględnimy szum n(t) pochodzenia termicznego. Na wyjściu kanału otrzymamy sygnał v(t) o mocy S i szum o mocy N. Stosunek mocy sygnału do szumu: SNR=\frac{S}{N} charakteryzuje jakość odebranego sygnału, a co za tym idzie, jakość kanału transmisyjnego. Demodulator ma za zadanie wydzielenie z sygnału v(t) kopii sygnału modulującego s0(t). Ze względu na szum, udaje się to tylko częściowo: na wyjściu demodulatora otrzymuje się sygnał z szumem:

m^\ast\left(t\right)=s_0\left(t\right)+n_0(t) (1)

Stosunek mocy sygnału użytecznego na wyjściu demodulatora (S0) do mocy szumu (N0) jest pewnym wskaźnikiem jakości otrzymanego sygnału: {SNR}_0=\frac{S_0}{N_0}
Aby ocenić odporność systemu transmisyjnego na szum w kanale, należy porównać SNR0 i SNR.  Duży przyrost wartości SNR w demodulatorze świadczy o dużej odporności systemu na zakłócenia. 
Innym parametrem charakteryzującym jakość systemu transmisyjnego jest wykorzystane pasmo B. Oczywiście pasmo zajęte przez sygnał zmodulowany (B) jest proporcjonalne do pasma sygnału modulującego (f_M). Oszczędny system charakteryzuje się niewielkim współczynnikiem poszerzenia pasma \frac{B}{f_M}. Niestety oba kryteria jakości systemu (odporność na zakłócenia i wąskopasmowość) są sprzeczne. Dużą odpornością na zakłócenia charakteryzują się systemy szerokopasmowe. 



Rysunek 1 Analogowe łącze transmisyjne

 

1.2. Modulacje dwuwstęgowe i jednowstęgowe AM

Modulacje dwuwstęgowe i jednowstęgowe  AM

1.3. Modulacja amplitudy – dwuwstęgowa z wytłumioną falą nośną (AM – DSB-SC double sideband suppressed carrier)

Schemat modulatora AM DSB-SC zawiera układ mnożenia sygnału modulującego m(t) (może to być sygnał mowy lub muzyki) przez sygnał fali nośnej (wzór 2 i rys.2). 

s(t)=m(t)cos(2πf0t) (2)

Rysunek 2 Schemat modulatora DSB-SC i przebieg sygnału zmodulowanego

Na tym samym rysunku pokazano sygnał DSB-SC w dziedzinie czasu. Zmniejszono częstotliwość sygnału fali nośnej żeby pokazać zmianę fazy w momencie przejścia sygnału modulującego przez zero. 
Załóżmy,że znamy widmo (transformatę Fouriera) sygnału modulującego: M(f)=F[m(t)]. Wówczas widmo sygnału zmodulowanego wyraża się wzorem:

S(f)=F[m(t)cos{(}2\pi\ f_0t)]=F[m(t)]*F[cos(2πf0t)]=\\ \quad=M(f)\ast\frac{1}{2}[δ(f-f0)+δ(f+f0)]=12M(f-f0)+12M(f+f0) (3)

Wyprowadzając wzór, wzięto pod uwagę, że transformata Fouriera iloczynu jest splotem widm, a splot z impulsem Diraca jest przesunięciem funkcji w miejsce położenia impulsu. Widma sygnałów na wejściu i wyjściu modulatora DSB-SC pokazano na Rys.3. W widmie sygnału zmodulowanego nie widać prążków widma fali nośnej, dlatego mówimy o wytłumionej fali nośnej.

Rysunek 3 Widmo sygnału modulującego i zmodulowanego w systemie DSB-SC

Pasmo częstotliwości zajęte przez sygnał zmodulowany jest dwukrotnie szersze niż pasmo sygnału modulującego. Współczynnik poszerzenia pasma \frac{B}{f_M}  ma wartość 2. 
Energie sygnału modulującego i zmodulowanego wyrażają się wzorami (4) i (5):

E_M=\int ^{f_{M}}_{-f_{M}}|M(f)|^2df (4)
E_S=\int{|S(f)|^2}df=\frac14E_M+\frac14E_M=\frac12E_M (5)

Wzór (5) jest słuszny gdy amplituda sygnału fali nośnej równa się 1. Ogólnie, energia podąża za kwadratem amplitudy i prawą stronę równania (5) należy pomnożyć przez kwadrat amplitudy sygnału fali nośnej. Pamiętajmy też o ogólnym założeniu, że obliczamy energie i moce, jakie wydzielają się w oporniku o oporności 1W, gdy napięcie i prąd mają identyczną wartość. Projektując konkretne układy nadajników, musimy przeskalować wartości energii i mocy. 
Wzór (5) można również zastosować do mocy sygnałów (uwzględniając powyższe założenia). Po modulacji moc sygnału maleje zatem do połowy: S=\frac{1}{2}P, gdzie P jest mocą sygnału modulującego. 
Na Rys.2 można zaobserwować, że tzw. obwiednia sygnału DSB-SC (można ją otrzymać za pomocą detektora szczytowego, rys.9) nie jest kopią sygnału modulującego m(t), tylko jego wartości bezwzględnej. Taki detektor nie nadaje się do odbioru sygnału DSB-SC. Należy użyć innego układu, zwanego detektorem synchronicznym (Rys.4).  

Rysunek 4 Odbiornik z detektorem synchronicznym

Sygnał użyteczny na wejściu odbiornika jest sygnałem DSB-SC: s(t)=m(t)cos(2πf0t). Ten sam sygnał pojawi się na wyjściu filtru pasmowego. Uwzględniając zależność {cos}^2{\alpha}=\frac{1}{2}+\frac{1}{2}cos{(}2\alpha), otrzymamy sygnał na wyjściu detektora synchronicznego: m(t)cos2(2πf0t)=m(t)[12+12cos(4πf0t)]. Sygnał o wysokiej częstotliwości zostanie zatrzymany w filtrze dolnopasmowym, w związku z czym na wyjściu odbiornika otrzymamy sygnał użyteczny 

s_0(t)=\frac12m(t) (6)

o mocy S_0=\frac{1}{4}P. Dwukrotne zmniejszenie amplitudy oznacza czterokrotne zmniejszenie mocy. Na rys. 5 pokazano wyniki symulacji: przebiegi czasowe i widma sygnałów w odbiorniku z detektorem synchronicznym. 

Rysunek 5 Sygnały w odbiorniku synchronicznym DSB-SC

Z modulacją DSB-SC stykamy się niemal codziennie, gdy słuchamy analogowego radia FM. Sygnał stereofoniczny dociera do nas dzięki modulacji częstotliwości, ale modulację amplitudy DSB-SC  wykorzystuje się do połączenia sygnałów kanału lewego L(t) i prawego R(t) w jeden sygnał stereofoniczny modulujący nadajnik FM. Analogowy sygnał stereofoniczny otrzymuje się w następujący sposób: 

m(t)=[L(t)+R(t)]+[L(t)-R(t)]cos(2πf_0t)+a\ cos(2π\frac{f_0}{2}t) (7)

Sygnał L(t)+R(t) wstawia się w paśmie podstawowym, natomiast L(t)-R(t)  moduluje w systemie DSB-SC się z wykorzystaniem fali nośnej o niewielkiej częstotliwości f0=38 kHz. Ponadto transmituje się sygnał harmoniczny o częstotliwości 38:2=19 kHz, zwany sygnałem pilota.  Sygnał stereofoniczny (bez pilota) pokazano na Rys. 6.

Rysunek 6 Analogowy sygnał stereofoniczny

Ze wzoru (7) wynika,  że sygnały L(t) i R(t) można rozdzielić, dokonując próbkowania sygnału stereofonicznego w odpowiednich momentach. Potrzebne są dwa układy próbkujące, jeden pobiera próbkę gdy fala nośna osiąga szczytową wartość równą 1, a drugi, gdy osiąga ona wartość równa -1. Momenty próbkowania określa się na podstawie analizy sygnału pilota.  
Widmo (kompleksowego) sygnału stereofonicznego pokazano na rys.7. Pasmo częstotliwości zajęte przez sygnały L(t) i R(t) wynosi 15 kHz.Pasmo kompleksowego sygnału stereofonicznego wynosi 38+15=53 kHz. 


Rysunek 7 Widmo analogowego sygnału stereofonicznego

1.4. Modulacja amplitudy – dwuwstęgowa z falą nośną (AM – DSB double sideband)

Wzór na modulację dwuwstęgową AM można zapisać, dodając falę nośną do sygnału DSB-SC, choć oczywiście nie jest to sposób wykorzystywany w technice. Można też interpretować ten wzór jak wynik mnożenia sygnału  m(t) z dodaną składową stałą (równą 1) przez falę nośną:

s\left(t\right)=m\left(t\right)\cos{2}\pi\ f_0t+\cos{2}\pi\ f_0t=[1+m\left(t\right)]cos2πf_0t (8)

Na Rys.8 pokazano przebieg czasowy sygnału AM-DSB. Jeśli w każdej chwili t spełniony jest warunek \left[1+m\left(t\right)\right]>0, wówczas mówimy że brak jest przemodulowań i obwiednia sygnału AM-DSB jest kopią sygnału modulującego.

Rysunek 8 Modulator AM-DSB i przebieg czasowy sygnału zmodulowanego

Wówczas jako demodulatora można użyć bardzo prostego układu detektora szczytowego z diodą i obwodem RC (Rys.9). 
Widmo sygnału AM-DSB naszkicowano na Rys.10. Jest to widmo sygnału DSB-SC, do którego dodano prążki widma fali nośnej (porównaj z Rys.3.  Transmisja fali nośnej wymaga dodatkowej mocy nadajnika. Ta dodatkowa moc nie niesie żadnej informacji – jest ona zawarta we wstęgach bocznych. Jej wielką zaletą jest wykorzystanie prostego detektora, niewymagającego generowania fali nośnej po stronie odbiorczej. Taki detektor nazywamy niekoherentnym. 

Rysunek 9 Detektor obwiedni w układzie detektora szczytowego

Rysunek 10 Widmo sygnału AM-DSB

 

 

1.5. Modulacja amplitudy – jednowstęgowa z wytłumioną falą nośną (AM – SSB-SC single sideband suppressed carrier)

Informacja osygnale modulującym jest ukryta we wstęgach bocznych widma. Ze względu na symetrię, zawierają one tę samą informację. Sygnał zmodulowany może więc zawierać jedną wstęgę, drugą można wytłumić. Zmniejsza się w ten sposób pasmo kanału wykorzystywanego do transmisji (Rys.11). Współczynnik poszerzenia pasma \frac{B}{f_M}  ma wartość 1.
Gdy zachowana zostaje wstęga górna, mówimy o górnowstęgowej modulacji SSB (USSB – Upper Single Sideband). W przeciwnym przypadku mamy do czynienia z modulacją dolnowstęgową (LSSB – Lower Single Sideband). 
Do odbioru sygnałów AM-SSB służy detektor synchroniczny (Rys.4). Na Rys.12 pokazano wyniki symulacji: przebiegi czasowe i widma sygnałów w różnych punktach nadajnika i odbiornika. Proszę zwrócić uwagę na fakt, że prawa część widma sygnału modulującego staje się górną wstęgą sygnału USSB, po czym w odbiorniku, na wyjściu detektora synchronicznego, niejako wraca na swoje miejsce. To samo można powiedzieć o lewej części widma, która jest lustrzanym odbiciem części prawej. 
Ze względu na oszczędność pasma, modulacja SSB jest chętnie wykorzystywana przez amatorów – krótkofalowców. 

Rysunek 11 Generowanie sygnału AM-SSB metodą filtracyjną

 

Rysunek 12 Przebiegi czasowe i widma sygnałów w modulatorze i demodulatorze AM-SSB

Aby zbadać odporność modulacji AM-SSB na szum w kanale transmisyjnym, porównajmy stosunek mocy sygnału do szumu na wyjściu kanału i na wyjściu odbiornika. Załóżmy, że sygnał modulujący ma moc P, amplituda sygnału fali nośnej jest równa 1 a gęstość mocy szumu w kanale wynosi hW/Hz.  Wówczas na wyjściu modulatora dwuwstęgowego (DSB-SC) otrzymujemy sygnał o mocy S_{DSB}=\frac{1}{2}\ P  (p.1.2.1). Po stłumieniu jednej wstęgi moc zmniejsza się dwukrotnie:   S=\frac{1}{4}\ P.  Sygnał SSB o tej mocy pojawia się na wejściu odbiornika. Detektor synchroniczny to układ mnożący przez falę nośną, podobnie jak modulator dwuwstęgowy (Rys.3). Na jego wyjściu pojawiają się dwie przesunięte kopie widmowe sygnału wejściowego o amplitudzie zmniejszonej do połowy, a więc o mocy zmniejszonej 4 razy. Całkowita moc sygnału na wyjściu detektora synchronicznego wynosi \frac{1}{4}S+\frac{1}{4}S=\frac{1}{2}S. Połowa tej mocy znajduje się w zakresie wysokich częstotliwości i zostanie wytłumiona w dolnopasmowym filtrze (Rys.4). W paśmie podstawowym pozostanie połowa mocy: S_0=\frac{1}{2}\cdot\frac{1}{2}S\ =\frac{1}{4}S . Jest to moc sygnału użytecznego na wyjściu odbiornika. 
Podobnie przeprowadzimy analizę mocy szumu. Pasmo zajęte przez sygnał SSB wynosi fM, a więc szum na wyjściu kanału ma moc : N=\eta f_M. W detektorze synchronicznym pojawią się dwie niepokrywające się kopie widmowe szumu, podobnie jak to było z sygnałem SSB (Rys.12). Na wyjściu detektora synchronicznego moc szumu zmniejszy się dwukrotnie: \frac{1}{4}N+\frac{1}{4}N=\frac{1}{2}N, z czego na pasmo podstawowe (niskich częstotliwości) przypada połowa:N_0=\frac{1}{2}\cdot\frac{1}{2}N=\frac{1}{4}N. Moc szumu na wyjściu odbiornika  zmalała więc 4 razy w porównaniu z mocą szumu na jego wejściu. To samo stało się z mocą sygnału użytecznego: również zmalała czterokrotnie. Oznacza to, że stosunek mocy sygnału użytecznego do mocy szumu nie zmienił się po przejściu przez odbiornik:

(SNR_0=\frac{S_0}{N_0}=\frac{S}{N}=SNR) (9)

Nie jest to wysoka odporność na szum w kanale (są modulacje oferujące znaczne powiększenie SNR w odbiorniku, jednak wymagają one szerokiego pasma).  Przy współczynniku poszerzenia pasma równym 1 nie jest możliwe uzyskanie większej odporności na zakłócenia , tzn. {SNR}_0\le SNR.  Uzyskana równość (wzór 9) jest największą wartością SNR_0, jaką można uzyskać. Wynika to z twierdzenia Shannona o przepustowości kanału x.


 

 

1.6. Modulacja częstotliwości i modulacja fazy

Modulacja częstotliwości i modulacja fazy

1.7. Modulacja kąta i jej podstawowe parametry

Modulacja częstotliwości (FM – Frequency Modulation) i modulacja fazy (PM = Phase Modulation) należą do tzw. modulacji kąta. Kopia sygnału modulującego jest transmitowana w fazie lub w częstotliwości sygnału zmodulowanego, natomiast nie jest możliwe jej odzyskanie z amplitudy tego sygnału. Na rys.13 pokazano modulator kątowy i naszkicowano możliwy sygnał wyjściowy modulatora. Rysunek ma charakter poglądowy. 

 

Rysunek 13 Modulacja kąta

Sygnał modulujący m(t) może wpływać na fazę chwilową Q(t) bezpośrednio lub przez całkę. W tym pierwszym przypadku mówimy o modulacji fazy

\mathrm{\Theta}(t)=2\pi\ f_0t+D_pm(t) (10)

W tym drugim przypadku mamy do czynienia z modulacją częstotliwości:

\mathrm{\Theta}(t)=2\pi\ f_0t+D_f\int_{t0}^tm(λ)dλ (11)

Współczynniki Dp i Df możemy traktować jako wzmocnienie toru akustycznego w nadajniku, można też je wykorzystać jako amplitudę sygnału modulującego (wówczas przyjmujemy, że m(t) ma największą wartość równą 1). 
Częstotliwość chwilowa sygnału to szybkość zmian fazy tego sygnału

f(t)=\frac{1}{2\pi}\frac{∂}{∂t}Θ(t) (12)

We wzorze dzielimy pochodną fazy przez 2p, aby uzyskać częstotliwość w Hz, w przeciwnym wypadku mielibyśmy prędkość kątową w radianach ma sekundę.  Różniczkując wzory (10) i (11), otrzymujemy wzory na częstotliwość chwilową sygnału PM i FM. 

PM:\qquad f(t)=f_0+\frac{D_p}{2\pi}\frac{∂}{∂t}m(t) (13)
FM:\qquad f(t)=f_0+\frac{D_f}{2\pi}m(t) (14)

Sygnał modulujący wpływa więc bezpośrednio na częstotliwość chwilową w modulacji FM. 
Dla harmonicznego sygnału modulującego określa się parametry sygnału zmodulowanego, takie jak indeks modulacji i dewiacja częstotliwości  Weźmy sygnał harmoniczny o amplitudzie równej 1 i częstotliwości fm:\ m\left(t\right)=\cos{(2\pi f_mt)}. Jeśli chcemy zmienić amplitudę tego sygnału, to możemy wykorzystać współczynniki wzmocnienia Dp i Df. Po podstawieniu do wzoru (10), otrzymujemy fazę chwilową sygnału PM: 

\mathrm{\Theta}(t)=2\pi\ f_0t+D_pcos(2πf_mt) (15)

Największa wartość odchyłki fazy chwilowej sygnału zmodulowanego  od fazy fali nośnej nosi nazwę indeksu modulacji i jest oznaczana przez b. W modulacji fazy, indeks modulacji wynosi \beta=D_p, czyli zależy od amplitudy sygnału modulującego. Podstawiając m\left(t\right)=\cos{(2\pi f_mt)} do (13), otrzymujemy wzór na częstotliwość chwilową sygnału PM:

f(t)=f_0+\frac{D_p}{2\pi}2πf_m(-sin(2πf_mt))=f_0-D_pf_msin(2πf_mt) (16)

Największa wartość odchyłki częstotliwości chwilowej sygnału zmodulowanego  od częstotliwości fali nośnej nosi nazwę dewiacji częstotliwości. Dewiacja częstotliwości określa w pierwszym przybliżeniu pasmo częstotliwości zajęte przez sygnał zmodulowany.  Dla modulacji PM wynosi ona \mathrm{\Delta F}=D_pf_m i zależy zarówno od amplitudy jak i częstotliwości sygnału modulującego. Ponadto słuszna jest zależność:

\beta=\frac{\mathrm{\Delta F}}{f_m} (17)

Obliczmy teraz indeks modulacji i dewiację częstotliwości dla modulacji FM.  Faza chwilowa wyraża się wzorem (11), co po podstawieniu sygnału harmonicznego wynosi

\mathrm{\Theta}(t)=2\pi f_0t+D_f\int_0^tcos(2πf_mλ)dλ=2πf_0t+\frac{D_f}{2πf_m} sin(2πf_mt) (18)

Największa odchyłka fazy, czyli indeks modulacji wynosi więc \beta=\frac{D_f}{2\pi f_m}. Z kolei dewiację częstotliwości otrzymujemy ze wzoru (14), co po podstawieniu sygnału harmonicznego wynosi

f(t)=f_0+\frac{D_f}{2\pi}cos(2πf_mt) (19)

Dewiacja częstotliwości zależy od amplitudy sygnału modulującegoi nie zależy od częstotliwości: \mathrm{\Delta F}=\frac{D_f}{2\pi}. Nie oznacza to, że częstotliwość sygnału modulującego nie wpływa w ogóle na pasmo częstotliwości sygnału FM. Wpływ ten istnieje, ale jest znacznie mniejszy niż w przypadku modulacji PM.  M.in. z tego względu analogowa modulacja częstotliwości jest bardziej popularna niż analogowa modulacja fazy. W radiofonii komercyjnej na falach ultrakrótkich (częstotliwości rzędu 100 MHz) wykorzystuje się modulację FM. Z tego względu w dalszej części podręcznika zajmiemy się modulacją FM. 
Zauważmy, że dla modulacji FM również spełniona jest zależność (17): \beta=\frac{\mathrm{\Delta F}}{f_m}.
 

1.8. Widmo i pasmo sygnału FM

Faza sygnału FM jest opisana wzorem (11), a jego przebieg czasowy wynosi

s\left(t\right)=A\ cos{[}2\pi f_0t+D_f\int_{t_0}^{t}m\left(\lambda\right)d\lambda] (20)

Gdy odchyłka fazy chwilowej D_f\int_{t_0}^{t}{m(\lambda)d\lambda} jest duża, np. przekracza 1 radian, mówimy o szerokopasmowej modulacji częstotliwości: (wideband frequency modulation – WBFM).  Na rys.14 pokazano wyniki symulacji nadajnika FM. Sygnał modulujący jest sygnałem harmonicznym o amplitudzie m: m(t)=m\ cos(2\pi f_mt). Można zaobserwować powiększanie się odchyłki częstotliwości chwilowej wraz ze wzrostem amplitudy sygnału, zgodnie ze wzorem (14): \ f\left(t\right)=f_0+\frac{D_f}{2\pi}\ m(t)

Rysunek 14 Sygnały FM dla  m(t)=m cos(2pfmt)

W ogólnym przypadku obliczenie widma lub gęstości mocy sygnału WBFM nie jest sprawą prostą. Przybliżonego oszacowania gęstości mocy, a co za tym idzie, pasma sygnału FM można dokonać metodą aproksymacji quasi-statycznej. Przyjmuje się wówczas, że sygnał modulujący m(t) jest wolnozmiennym sygnałem przypadkowym o znanej gęstości prawdopodobieństwa wartości chwilowych (próbek) p_m(x).

Rysunek 15 Aproksymacja quasi-statyczna: gęstość prawdopodobieństwa wartości chwilowych sygnału modulującego, gęstość prawdopodobieństwa wartości częstotliwości chwilowej i oszacowana gęstość mocy sygnału FM

Zgodnie ze wzorem (14) wartości chwilowej sygnału modulującego równej x, odpowiada częstotliwość chwilowa sygnału FM, równa f=f_0+\frac{D_f}{2\pi}x.  Gęstość prawdopodobieństwa zmiennej f pokazano na rys.15.  Następnie przyjmujemy, że moc sygnału FM podąża za częstotliwością chwilową. Jest to przybliżenie, jednak tym dokładniejsze, im wolniej zmienia się częstotliwość chwilowa f (stąd nazwa aproksymacja quasi-statyczna). Przy tym założeniu gęstość mocy sygnału FM jest kopią gęstości prawdopodobieństwa częstotliwości chwilowej. Na rys.15 narysowano dwustronną gęstość mocy. 
Znając krańcowe wartości częstotliwości chwilowej, można oszacować pasmo sygnału FM. Wynosi ono 

B=\frac{D_f}{\pi}x_{max}=2\ ∆F (21)

Gdzie xmax jest amplitudą sygnału modulującego, a ∆F=\frac {D_f}{2π}x_{max} jest dewiacją częstotliwości, czyli największą jej odchyłką od częstotliwości fali nośnej.     

Dokładne obliczenie widma sygnału FM jest możliwe dla harmonicznego sygnału modulującego: m(t)=cos(2\pi f_mt).  Wówczas przebieg czasowy sygnału FM o amplitudzie fali nośnej A=1, wynosi:

s\left(t\right)=cos{[}2\pi f_0t+D_f\int_{0}^{t}{cos{(}2\pi f_m\lambda)}d\lambda]=\\=cos{[}2\pi f_0t+\frac{D_f}{2\pi f_m}sin{(}2\pi f_mt)]=cos[2πf_0t+β\ sin(2πf_mt)]=\\=cos{[}2\pi f_0t]cos[β\ sin(2πf_mt)]-sin[2πf_0t]sin[β\ sin⁡(2πf_mt)] (22)

Wzór (22) przedstawia w istocie dwie modulacje AM – dwuwstęgowe z wytłumioną nośną. Sygnałami modulującymi są cos(bsin(2pfmt)) i sin(bsin(2pfmt)). Są to funkcje okresowe, ich okres wynosi 1/fm. Na rys.16-18 pokazano przebiegi tych sygnałów w obrębie jednego okresu dla wybranych wartości indeksu modulacji b. Gdy indeks modulacji b wynosi 1 zauważa się drugą harmoniczną sygnału cos(bsin(2pfmt)) i pierwszą harmoniczną sygnału sin(bsin(2pfmt)) – rys.16. Przy b=3 widać wyraźnie wpływ czwartej harmonicznej sygnału cos(bsin(2pfmt)) i trzeciej harmonicznej sygnału sin(bsin(2pfmt)) – rys.17. Przy b=10 oba sygnały mają liczne składowe -odpowiednio parzyste i nieparzyste harmoniczne – rys.18.

   Rysunek 16 Sygnały cos[bsin(2pt)] (po lewej) i sin[bsin(2pt))] (po prawej), b=1

 

Rysunek 17 Sygnały cos[bsin(2pt)] (po lewej) i sin[bsin(2pt))] (po prawej), b=3

 

Rysunek 18 Sygnały cos[bsin(2pt)] (po lewej) i sin[bsin(2pt))] (po prawej), b=10

Amplitudy składowych harmonicznych omawianych sygnałów opisane są funkcjami Bessela, których argumentem jest indeks modulacji b. Położenie prążków widma obu sygnałów pokazano na rys.19, a na rys.20 – położenie prążków widma sygnału FM, dla harmonicznego sygnału modulującego o częstotliwości fm. Modulacja amplitudy opisana wzorem (22) przesuwa widmo sygnału modulującego w prawo i w lewo o częstotliwość sygnału modulującego.  Jeśli amplituda fali mośnej wynosi 1, wówczas amplituda składowych widmowych zmniejsza się dwukrotnie  (Rys.3). Na rys.20 pokazano prawostronną część widma, nie opisano również fazy składowych harmonicznych.

 

Rysunek 19 Położenie prążków widmowych sygnału cos(bsin(2pfmt)) (u góry) i sygnału sin(bsin(2pfmt))  (u dołu)

 

Rysunek 20 Położenie prążków widmowych sygnału FM dla harmonicznego sygnału modulującego

Z wykresów funkcji Bessela (rys.21) wynika, że dla małych wartości indeksu modulacji x=b, tylko nieliczne z nich osiągają wartości na tyle różne od zera, żeby można było mówić o ich wpływie na widmo i pasmo sygnału. I tak, dla wartości bliskich zeru, widmo składa się z fali nośnej (amplituda J0) i dwóch prążków bocznych (amplituda J1). Jest to tzw. Wąskopasmowa modulacja częstotliwości. 

Rysunek 21 Rodzina funkcji Bessela argumentu x

Zwiększenie indeksu modulacji o 1 powoduje “oderwanie się” kolejnej funkcji Bessela od wartości bliskich zeru i pojawienie się dwóch prążków w widmie sygnału FM, co skutkuje poszerzeniem pasma o podwojoną częstotliwość sygnału modulującego. Ta obserwacja stoi u podstaw tzw. prawa Carsona. Prawo Carsona określa szerokość pasma, w którym zawarte jest 98% mocy sygnału FM:

B=2(\beta+1)f_m=2(\frac{\mathrm{\Delta F}}{f_m}+1)f_m=2(\mathrm{\Delta F}+f_m) (23)

W komercyjnej radiofonii FM obowiązuje norma ustalająca wartość dewiacji częstotliwości równą DF=75kHz. Pasmo sygnału modulującego wynosi  fM=53 kHz  (jest nim sygnał stereofoniczny, rys.7). Obliczamy stąd pasmo B=256 kHz dla transmisji stereofonicznej w pasmie UKF
Jeżeli częstotliwość sygnału modulującego jest niewielka, wówczas pasmo jest równe podwojonej dewiacji częstotliwości, co już otrzymaliśmy metodą aproksymacji quasi-stacjonarnej (wzór 21).

 

1.9. Odbiór sygnału FM

Odbiornikiem sygnału FM może być każdy dyskryminator częstotliwości, tj. układ przetwarzający częstotliwość na sygnał napięciowy lub prądowy. Jest bardzo wiele takich układów, np. zliczanie przejść przez zero w jednostce czasu. Odbiornik z dyskryminatorem częstotliwości pokazano na rys.22. Oprócz dyskryminatora zawiera on filtr pasmowy na wejściu i dolnopasmowy na wyjściu. Ma to na celu ograniczenie mocy szumu na wyjściu odbiornika.

Rysunek 22 Odbiornik FM z dyskryminatorem częstotliwości

Częstotliwość chwilowa sygnału FM jest opisana wzorem (14): f\left(t\right)=f_0+\frac{D_f}{2\pi}m(t), jednak dyskryminator jest tak skalibrowany, że nie odtwarza stałego sygnału odpowiadającego częstotliwości fali nośnej, tylko sygnał proporcjonalny do m(t). 
Modulacja częstotliwości jest bardziej odporna na szumy w kanale, niż modulacja amplitudy czy fazy. Szum w kanale w podobny sposób „atakuje” amplitudę i fazę fali nośnej przenikając, jako szum o płaskiej gęstości mocy, na wyjście odbiornika. W modulacji częstotliwości szum też przenika na wyjście odbiornika, jednak jego gęstość mocy nie jest stała, tylko narasta w funkcji częstotliwości. Ma ona kształt paraboli przechodzącej przez punkt (0,0). Oczywiście  pasmo szumu jest ograniczone filtrem wyjściowym, przepuszczającym sygnał tylko w paśmie sygnału modulującego (od 0 do częstotliwości  fM). W zakresie niskich częstotliwości gęstość mocy szumu na wyjściu odbiornika FM jest mniejsza niż na wyjściu odbiornika PM czy AM. Mniejsza jest też całkowita moc szumu na wyjściu.
Dla wyjaśnienia tego problemu porównajmy wzory (10) i (11). Wynika z nich, że modulator fazy, jeśli podać na jego wejście scałkowany sygnał modulujący, działa jak modulator częstotliwości. Jeśli do demodulacji zastosujemy dyskryminator fazy, to na jego wyjściu otrzymamy scałkowany sygnał mowy czy muzyki, trzeba więc go będzie zróżniczkować. Wraz z nim przez układ różniczkujący przejdzie szum z dyskryminatora fazy. Układ różniczkujący jest filtrem górnopasmowym, którego charakterystyka częstotliwościowa rośnie w funkcji częstotliwości. W ten sposób na wyjściu odbiornika FM pojawia się szum gęstości mocy rosnącej w funkcji częstotliwości.
Dla sygnału modulującego m(t)=cos(2\pi f_mt) można otrzymać następującą zależność, charakteryzującą odporność na szumy w kanale:

SNR_0=\frac{3}{2}\ \frac{D_f^2}{4\pi^2f_M^2}SNR=\frac{3}{2}\beta^2SNR (24)

Można ten wzór przepisać, wyrażając wartości SNR w decybelach:

SNR_0[dB]=10\ log_{10}\frac32+20\log_{10}β+SNR[dB] (25)

Jakość sygnału na wyjściu demodulatora, wyrażona przez SNR0,  zależy od indeksu modulacji b.  Jeśli porównamy ten wzór z zależnością SNR_0=SNR , charakteryzującą modulację AM-SSB, widzimy znaczną przewagę modulacji FM. SNR0  może znacznie przewyższać stosunek mocy sygnału do szumu na wyjściu kanału (SNR). Płacimy jednak za to poszerzeniem pasma, gdyż pasmo zależy od indeksu modulacji b: wzór (23).
Wzory (24) i (25) są słuszne dla niezbyt zaszumionych kanałów. Przy bardzo dużym szumie, częstotliwość chwilowa sumy sygnału i szumu nie podąża za częstotliwością chwilową sygnału i na wyjściu odbiornika pojawia się silny szum impulsowy, uniemożliwiając odbiór. Ten tzw. efekt progowy pokazano na rys.23.

Rysunek 23 Efekt progowy w modulacji FM

 

Porównanie modulacji AM i FM można w skrócie podsumować następująco:

AM: 

  • wąskie pasmo (zwłaszcza modulacja jednowstęgowa SSB)

  • mała odporność na szumy w kanale

FM:

  • szersze pasmo sygnału zmodulowanego 

  • większa odporność na szumy w kanale

 

2. Kody transmisyjne i modulacje cyfrowe

Kody transmisyjne i modulacje cyfrowe

2.1. Jakość systemu transmisji cyfrowej

Schemat blokowy cyfrowego łącza transmisyjnego był pokazany w Module 1, tym niemniej został tu powtórzony (Rys.24). 
W procesie próbkowania, kwantyzacji i kodowania sygnał modulujący zostaje przetworzony na ciąg bitów. Do ich transmisji wykorzystuje się szereg elementarnych sygnałów, zwanych symbolami transmisyjnymi. Zbiór symboli tworzy kod transmisyjny lub modulację cyfrową. To pierwsze określenie jest stosowane najczęściej do opisu transmisji w paśmie podstawowym, w zakresie niskich częstotliwości. Do opisu transmisji w zakresie wysokich częstotliwości, gdy wykorzystywana jest fala nośna, używamy pojęcia modulacji cyfrowej. Nie zawsze jest to regułą, czasem mówimy o modulacji w paśmie podstawowym. 

Rysunek 24 Cyfrowe łącze transmisyjne

Jeśli wykorzystujemy dwa symbole, mamy do czynienia z kodem binarnym (modulacją binarną). Wówczas symbole transmisyjne będziemy oznaczać przez s0(t) i s1(t), przypisując odpowiednim symbolom wartości logiczne 0 i 1.. Mając M symboli, będziemy je numerować od 1 do M i będziemy mówić o modulacji M-wartościowej. W modulacji M-wartościowej jeden symbol może transmitować \log_2{M} bitów. Np. symbolom modulacji 4-wartościowej przypisujemy pary bitów: 00, 01, 10, 11.
Ciąg transmitowanych symboli opisują następujące parametry:

  • Szybkość modulacji, czyli liczba symboli na sekundę. Jeśli czas trwania symbolu wynosi T, to szybkość modulacji wynosi R_m=\frac{1}{T} bodów (jednostka bod [Bd] pochodzi id nazwiska francuskiego badacza E.Baudota). 
  • Szybkość transmisji binarnej (przepływność binarna)R_b [bit/s] zależy od szybkości modulacji i od wartościowości modulacji:
    R_b=lR_m,\qquad l={log}_2{M} (26)
    Gdzie l jest liczbą bitów przenoszona w jednym symbolu.
  • Energia symbolu E_i=\int_{0}^{T}{s_i^2\left(t\right)dt} i średnia energia E_s=\sum_{i=1}^{M}{P_iE_i} (Pi jest prawdopodobieństwem nadania i-go symbolu).  
  • Moc średnia sygnału (mierzona zwykle na wyjściu kanału) S=\frac{E_s}{T}=E_sR_m
  • Energia na bit E_b=\frac{S}{R_b} , czyli energia wydatkowana na przeniesienie jednego bitu.
  • Pasmo częstotliwości B jest proporcjonalne do szybkości modulacji: W typowych warunkach B\approx R_m=\frac{1}{T} (patrz rys.8 Moduł 1). Najmniejszą wartością, wynikająca z twierdzenia Nyquista, jest B=R_m/2
  • Efektywność widmowa jest miarą wykorzystania pasma częstotliwości. Mówi ona, ile bitów można przesłać na sekundę w paśmie 1 Hz:
    E_f=\frac{R_b}{B}  (27)
  • Stopa błędów binarnych BER (bit error rate), czyli prawdopodobieństwo przekłamania jednego bitu Pe. Zależy ona nie tylko od jakości odebranego sygnału (SNR), ale od szybkości transmisji Rb, oraz od zastosowanej modulacji.  
  • Stosunek mocy sygnału do szumu na wyjściu kanału SNR=S/N , wyrażany najczęściej w decybelach [dB]. Często stosowana jest inna miara jakości sygnału na wyjściu kanału: stosunek energii na bit (Eb) do gęstości mocy szumu na wyjściu kanału (h). W transmisji binarnej, gdy każdy symbol przenosi 1 bit, R_b=R_m\approx B  i obie miary są zbliżone:
    SNR=\frac{S}{N}=\frac{E_bR_b}{\eta B}\approx\frac{E_bR_m}{\eta B}=\frac{E_b}{\eta} (28)
  • Jeśli używamy transmisji cyfrowej do przenoszenia dźwięku lub obrazu, to odporność systemu transmisyjnego na szumy w kanale możemy przedstawić jako funkcję {SNR}_0=f\left(SNR\right) , podobnie jak w transmisji analogowej (p.1.1)

Wiele z wymienionych wyżej parametrów stanowi kryterium jakości systemu transmisyjnego. Pożądana byłaby duża szybkość transmisji bez błędów binarnych, wąskie pasmo (ze względu na oszczędność), wysoka efektywność widmowa i odporność na szumy. Niestety postulaty te są sprzeczne. Graniczne wartości niektórych z wymienionych parametrów będą podane na końcu tego modułu.

2.2. Kody transmisyjne

Na rys.8 w Module 1 naszkicowano przebiegi czasowe popularnych kodów. Strumienie binarne w układach cyfrowych przesyła się na bliskie odległości, wówczas stosuje się binarny kod unipolarny:  

s_1\left(t\right)\neq0,\qquad s_0\left(t\right)=0\qquad 0\le t (29)

Symbol s1(t) może mieć kształt prostokątny, może stanowić półokres funkcji sinus, itp. Wadą tego kodu jest niezerowa wartość średnia, co oznacza przepuszczanie prądu stałego przez łącze kablowe.  Nie ma tej wady binarny kod bipolarny:

s_1\left(t\right)\neq0,\qquad s_0\left(t\right)=-s_1(t)\qquad 0\le t (30)

Do prawidłowego odbioru symboli niezbędna jest synchronizacja elementowa – zlokalizowanie położenia symboli na osi czasu.  Oba wymienione kody nie zapewniają synchronizacji gdy nadawany jest np. długi ciąg jednakowych symboli o kształcie prostokątnym.  Z tego względu stosuje się kody z powrotem do zera (return to zero, RZ), do transmisji wykorzystując połowę szczeliny czasowej. Prowadzi to jednak do podwojenia szerokości pasma. Pewnym kompromisem jest kod AMI (alternate mark inversion) – rys.25.

Rysunek 25 Kod transmisyjny AMI i zmodyfikowany kod AMI

W kodzie AMI wartości logicznej „0” odpowiada brak sygnału, a wartości „1” impuls, na przemian dodatni i ujemny. Zagrożeniem dla utrzymania synchronizacji elementowej byłaby długa sekwencja zer, z tego względu zmodyfikowano kod, przypisując niektórym zerom logicznym impulsy wstawione niezgodnie z zasadą zmiany znaku. W ten sposób, bez skrócenia impulsów (a więc bez poszerzenia pasma) udało się utrzymać synchronizację elementową. 

2.3. Modulacje cyfrowe

Aby wykorzystać wyższe częstotliwości, np. do transmisji bezprzewodowej, potrzebna jest fala nośna o odpowiedniej częstotliwości f0. Najprostszą modulację otrzymamy, mnożąc sygnał kodu unipolarnego przez nośną. Jest to cyfrowa modulacja amplitudy zwana OOK (on-off-keying). Na rys.26 naszkicowano przebieg czasowy i gęstość mocy sygnału OOK. 

Rysunek 26 Modulacja OOK w dziedzinie czasu i częstotliwości, wg [skrypt]

Widmo sygnału OOK jest przesuniętym do częstotliwości f0 widmem sygnału kodu unipolarnego. Patrz rys.3 przedstawiający efekt mnożenia sygnału m(t) przez falę nośną. To samo odnosi się do gęstości mocy. Sygnał kodu unipolarnego posiada składową stałą, więc w jego widmie istnieje prążek na częstotliwości równej zero. Po pomnożeniu przez falę nośną prążek przemieszcza się w prawo i w lewo o częstotliwość f0 (na rys. 26 pokazano tylko prawą część gęstości mocy). Koncentracja mocy na częstotliwości f0 nie przenosi informacji, lecz wymaga zwiększenia mocy nadajnika. Brak tej koncentracji mocy poprawiłby parametry modulacji , w szczególności obniżyłaby się wartość SNR na wyjściu kanału przy tej samej jakości sygnału na wyjściu odbiornika (czyli przy tej samej wartości SNR0). Większy stosunek SNR0/SNR (lub większa różnica, jeśli wielkości te wyrażono w dB) świadczyłaby o lepszej odporności na szumy w kanale. 
Prążek fali nośnej zniknie, jeśli w modulatorze (rys.26) użyjemy kodu bipolarnego zamiast unipolarnego. Sygnał kodu bipolarnego ma zerową składową stałą. Otrzymamy w ten sposób binarną modulację fazy (BPSK – binary phase shift keying), zwaną też niekiedy dwuwartościową modulacją amplitudy. Symbolami transmisyjnymi są tu sygnały harmoniczne o przeciwnej fazie (i co za tym idzie, przeciwnej polaryzacji), np.

s_1\left(t\right)=\cos{2\pi f_0t},\qquad s_0\left(t\right)=-\cos{2\pi f_0t}\qquad 0\le t (31)

Popularną modulacją cyfrową jest różnicowa modulacja fazy (DPSK -differential phase shift keying). Tutaj wykorzystuje się tę samą parę sygnałów, co w PSK. Informacja jest jednak kodowana w różnicy faz pomiędzy kolejnymi sygnałami. Np. zmiana fazy o 180 stopni oznacza logiczne „1”, a utrzymanie tej samej fazy  logiczne „0” (konwencja może też być odwrotna). Przykład podano na rys.27.
Symbole opisane równaniem

s_i\left(t\right)=A_i\cos{(2\pi f_0t+\varphi_i)} (32)

czyli sygnały o tej samej częstotliwości, różniące się amplitudą o fazą początkową, są przedstawiane w postaci punktu na płaszczyźnie zespolonej – rys.28.  Zbiór M takich punktów tworzy konstelację opisującą daną modulację. 

 

Rysunek 27 Dwuwartościowa modulacja PSK i DPSK wg [skrypt]

 

Rysunek 28 Punkt konstelacji przedstawiający sygnał (32)

Na rys.29 pokazano konstelacje 4 – 8-wartościowych modulacji PSK. Umożliwiają one transmisję odpowiednio 2 i 3 bitów w jednym symbolu. 4-wartościowa modulacja PSK jest jednocześnie 4- wartościową QAM (quadrature amplitudę modulation), gdyż można ją otrzymać z dwóch 2-wartościowych modulacji amplitudy: jedna wykorzystuje nośną kosinusoidalną a druga sinusoidalną. Na rys.30 pokazano konstelacje 16- i 32- wartościowej modulacji QAM, oraz jedną z licznych kombinacji AM-PM (modulacja amplitudy i fazy).  

Rysunek 29 Modulacja PSK 4-  i 8-wartościowa

 

Rysunek 30 Modulacje 16-wartościowe QAM i AM-PM oraz 32-wartościowa QAM

Modulacje wielowartościowe umożliwiają lepsze wykorzystanie pasma kanału. Pasmo zależy od czasu trwania symbolu (T). Nieskończenie długi symbol opisany wzorem (32) miałby nieskończenie wąskie pasmo. Im częściej zmieniają się parametry sygnału (amplituda, faza), tym szersze jest pasmo. Pasmo jest więc proporcjonalne do szybkości modulacji R_m=\frac{1}{T} . W modulacji 2-wartościowej każdy symbol przenosi jeden bit, a w M-wartościowej  l=\log_2{M}  bitów. Rośnie zatem efektywność widmowa:

\frac{R_b}{B}=\frac{R_ml}{B}=\frac{2Bl}{B}=2l (33)

We wzorze (33) podstawiono R_m=2B, jest to wartość graniczna, w praktyce szybkość modulacji jest nieco mniejsza.
Innym rodzajem modulacji jest kluczowanie częstotliwości (FSK-frequency shift keying). W binarnej FSK występują dwie częstotliwości fali nośnej – jednej z nich jest przypisana wartość logiczna 1 , a drugiej 0. 
Nośnych może być więcej, najczęściej tworzą one równoodległy na szereg prążków widmowych. Tak powstaje OFDM (Orthogonal Frequency Division Multiplexing). [1]
Strumień binarny jest rozdzielany między N modulatorówamplitudy pracujących z falami nośnymi o różnych częstotliwościach (rys.31). Widma sygnałów zmodulowanych pokrywają się częściowo, tym niemniej sygnały można rozdzielić po stronie odbiorczej, dzięki ich wzajemnej ortogonalności. 


Rysunek 31 Modulator OFDM (wg [1])

Sygnał OFDM można wygenerować, stosując odwrotną dyskretną transformatę Fouriera (IDFT) – rys.32. Wybierając N symboli do jednoczesnej transmisji z wykorzystaniem N modulatorów, tworzymy pożądane widmo DFT. Należy pamiętać o „lustrzanym” odbiciu i o tym, że próbka o indeksie zerowym i próbka w pozycji środkowej są liczbami rzeczywistymi, a więc nie można in nadać pożądanej fazy. Potrzeba zatem 2N+2 próbek widma. Po zastosowaniu IDFT  otrzymuje się sygnał OFDM w dziedzinie czasu. Modulator amplitudy pozwoli na przeniesienie go w dowolny zakres częstotliwości.

Rysunek 32 Widmo DFT sygnału OFDM

3. Odbiór optymalny sygnałów cyfrowych

Odbiór optymalny sygnałów cyfrowych

3.1. Optymalna reguła decyzyjna w transmisji binarnej

W kanale transmisyjnym występuje wiele czynników utrudniających odbiór: szum termiczny, ograniczone pasmo, tłumienie, zaniki, zakłócenia impulsowe, zniekształcenia nieliniowe. Załóżmy, że występuje jedynie szum pochodzenia termicznego o gęstości mocy h W/Hz (dwustronnie h/2 W/Hz). Kod transmisyjny jest binarny, a pasmo kanału równe pasmu sygnału zmodulowanego (rys.33). Wówczas można przyjąć, że symbole transmisyjne pojawiają się na wyjściu kanału w swym pierwotnym kształcie: s_0^\prime=s_0,\quad s_1^\prime=s_1. Moc szumu na wyjściu kanału jest równa N=hB. Zakładamy również, że zapewniona jest synchronizacja elementowa, tzn. znane są momenty rozpoczęcia i zakończenia kolejnych symboli. 

Rysunek 33 Uproszczony model kanału transmisyjnego

Zaszumiony sygnał v(t), pojawiający się na wyjściu kanału, powinien być zmierzony w taki sposób, aby wynik pomiaru zawierał całą informację o odebranym symbolu i szumie Ze względu na ograniczone pasmo (B), właściwą metodą pomiarową jest próbkowanie w częstotliwością 2B. W czasie trwania jednego symbolu zgromadzimy w ten sposób N=2BT próbek – rys.34.

 

Rysunek 34 Pomiar odebranego symbolu metodą próbkowania

Mając wynik pomiaru, czyli wektor  \mathbf{v}=[v_1,\ v_2,\ldots,v_N]t,  (t – transpozycja), należy teraz opracować algorytm decyzyjny, który każdemu możliwemu wektorowi v przypisywałby decyzję: 0 lub 1. Oznacza to, że N-wymiarową przestrzeń należy podzielić na dwa zbiory decyzyjne: V0 i V1.  Na rys.35 pokazano podział przestrzeni 2-wymiarowej.  Punkty s_0^\prime,\ s_1^\prime to wyniki pomiaru symboli bez szumu, a p_0(v),\ p_1(v) to gęstości prawdopodobieństwa otrzymania wyniku pomiaru równego v przy nadaniu s0 i s1 Oznaczmy jeszcze przez P0 i P1 prawdopodobieństwa nadania obu symboli (oczywiście P0+P1=1).

Rysunek 35 Algorytm decyzyjny w przestrzeni 2-wymiarowej

Błąd wystąpi w dwóch przypadkach: podjęcia decyzji 1 gdy nadano 0 i podjęcia decyzji 0 gdy nadano 1. Prawdopodobieństwo błędnego odbioru wynosi:

P_e=P_0P\left(1\middle|0\right)+P_1P(0|1) (34)

Prawdopodobieństwa warunkowe są równe P\left(1\middle|0\right)=\int_{V_1}^{\ }{p_0\left(v\right)dv},   P\left(0\middle|1\right)=\int_{V_0}^{\ }{p_1\left(v\right)dv}

Ze względu na fakt, że \int_{V_0}p_1(v)dv+\int_{V_1}p_1(v)dv=1 , to ostatnie można przepisać w postaci \int_{V_0}p_1(v)dv=1-\int_{V_1}p_1(v)dv. Ostatecznie (34) można zapisać następująco:

P_e=P_0\int_{V_1}p_0(v)dv+P_1 [1-\int_{V_1}p_1(v)dv]=P_1+\int_{V_1}[P_0p_0(v)-P_1p_1(v)]dv (35)

Minimalizując Pe, mamy wpływ jedynie na podział przestrzeni obserwacji sygnału na dwa zbiory decyzyjne, możemy jedynie zdefiniować V1. Optymalna reguła decyzyjna polega na zakwalifikowaniu do V1 wszystkich obserwacji v, dla których funkcja podcałkowa jest ujemna. Całka we wzorze (35) przyjmie wówczas wartość ujemną o jak największej wartości bezwzględnej.  Reasumując, optymalna reguła decyzyjna jest następująca:

P_0p_0\left(v\right)

​​​​P_0p_0(v)\geq P_1p_1(v)⇒v∈V_0

(36)

Najczęściej oba symbole pojawiają się z tym samym prawdopodobieństwem: P_0=P_1=\frac{1}{2}. Wówczas na decyzję ma wpływ jedynie wartość gęstości prawdopodobieństwa: gdy p_1\left(v\right)>p_0\left(v\right) podejmujemy decyzję „1” i odwrotnie.
Należy zauważyć, że kształt funkcji gęstości prawdopodobieństwa p_1\left(v\right) i p_0\left(v\right) jest identyczny, gdyż zależy jedynie od szumu (jest to dwuwymiarowa zmienna o rozkładzie gaussowskim).  Warunek p_1\left(v\right)>p_0\left(v\right) będzie spełniony w punkcie v, który leży bliżej centralnego punktu s_1^\prime niż punktu s_0^\prime. Zamiast obliczać wartości p_1\left(v\right) i p_0\left(v\right), można porównywać odległości ||v-s_0^\prime||,\quad||v-s_1^\prime||. Chodzi o odległości euklidesowe ||x||=\sqrt{\sum_{i=1}^{N}x_i^2}, gdyż gaussowska gęstość prawdopodobieństwa jest stała w punktach jednakowo odległych od punktu centralnego (rys.35). Ostatecznie optymalna reguła decyzyjna w przypadku gdy  P_0=P_1=\frac{1}{2} jest następująca:

\left|\left|v-s_1^\prime\right|\right|

||v-s_1^\prime||\geq||v-s_0^\prime||\Rightarrow v\in V_0

(37)

Oczywiście można też porównywać kwadraty odległości. Zasada minimum odległości opisana wzorem (37) jest łatwa w stosowaniu i intuicyjnie prosta, trzeba jednak pamiętać o warunku P_0=P_1=\frac{1}{2}

3.2. Realizacje odbiornika optymalnego

Zasadę optymalnego odbioru w transmisji binarnej (wzór 37) można przepisać w taki sposób, aby odnosiła się do sygnałów ciągłych:

\int_0^T[v(t)-s_1(t)]^2dt

\int_0^T[v(t)-s_1(t)]^2dt\geq\int_0^T[v(t)-s_0(t)]^2dt⇒"0"

(38)

Pierwszy wiersz (38) można przepisać w następującej postaci:

\int_{0}^{T}{v^2(t)dt-2}\int_{0}^{T}{v(t)s_1(t)dt+\int_{0}^{T}{s_1^2(t)dt\

a następnie uprościć i podzielić stronami przez (-2):

\int_{0}^{T}{v(t)s_1(t)dt-\frac{1}{2}\int_{0}^{T}{s_1^2(t)dt}} >\int_{0}^{T}{v(t)s_0(t)dt-\frac{1}{2}\int_{0}^{T}{s_0^2(t)dt}}\Rightarrow"1" (39)

Wzór (39) przedstawia zasadę działania odbiornika optymalnego w układzie korelatora – rys.36. 

Rysunek 36 Odbiornik optymalny kodu binarnego w układzie korelatora

Istotnie, całka \int_{0}^{T}{v(t)s_1(t)dt} przedstawia korelację sygnału odebranego z jednym z symboli.  Z kolei E_1=\int_{0}^{T}{s_1^2(t)dt\ } jest energią symbolu.
Wzór (39) można przepisać w nieco inny sposób:

\int_0^Tv(t)[s_1(t)-s_0(t)]dt>\frac12E_1-\frac12E_0\ ⇒\ "1" (40)

Jeśli korelacja sygnału odebranego v(t) z sygnałem różnicowym p(t)=s_1(t)-s_0(t) przekracza próg decyzyjny \frac{1}{2}E_1-\frac{1}{2}E_0 , wówczas podejmujemy decyzję, że odebrany sygnał jest zaszumionym sygnałem s_1(t). W przeciwnym wypadku uznajemy ten sygnał za s_0(t). Schemat zmodyfikowanego korelatora pokazano na rys.37. Korelację oznaczono przez y. 

 

Rysunek 37 Korelator z progiem decyzyjnym

Przy okazji warto zapoznać się z popularnym w teorii sygnałów pojęciem filtru dopasowanego. Układ obliczania korelacji (mnożenie i całkowanie na rys.37) można zastąpić filtrem dopasowanym do sygnału p(t). 
Z wyjścia filtru pobieramy próbkę pod koniec szczeliny czasowej przeznaczonej na odbiór symbolu – rys.38. Jest ona równa korelacji y z rys.37, co wynika z następujących przekształceń.

Rysunek 38 Odbiornik optymalny z filtrem dopasowanym

Odpowiedź impulsowa filtru wynosi h(t) = p(T-t). Sygnał wyjściowy filtru jest splotem sygnału wejściowego z odpowiedzią impulsową:

y(t)=v(t)\ast h(t)=\int_{-\infty}^{\infty}{v(\tau)h(t-\tau)d\tau}=\int_{-\infty}^{\infty}{v(\tau)p(T-t+\tau)d\tau}

Próbka sygnału y(t) pobrana w momencie t=T wynosi

y=y(t)\left|\begin{matrix}\ \\t=T\\\end{matrix}\right.=\int_{-\infty}^{\infty}{v(\tau)p(\tau)d\tau}=\int_{0}^{T}{v(\tau)p(\tau)d\tau} (41)

Kończy to dowód, że obliczanie korelacji można zastąpić filtracją dopasowaną.
Odbiornik przetwarzający próbki sygnałów (wzór 37), jak i odbiorniki przetwarzające sygnały ciągłe: wzór 38, rys. 26, rys.37, rys.38, realizują tę samą zasadę minimum odległości i są sobie równoważne. W szczególności oferują tę samą stopę błędów. Wynosi ona 

P_e=\frac{1}{2}erfc{(}\sqrt{\frac{E_{10}}{4\eta}\ \ \ }) (42)

gdzie E_{10}=\int_{0}^{T}{p^2(t)dt} jest energią sygnału różnicy p(t)=s_1(t)-s_0(t), h jest gęstością mocy szumu na wyjściu kanału, a erfc{(}x)=\frac{2}{\sqrt\pi}\int_{x}^{\infty}{exp{(}-u^2)du} jest tzw. komplementarną funkcją błędu. Jest ona silnie malejąca (rys.39), dlatego dążymy do jak największej wartości E10. Jest to oczywiste, gdyż energia E10 jest miarą różnicy między symbolami, która powinna być jak największa, jeśli mamy te sygnały rozpoznać w obecności szumu.  Pamiętajmy poza tym, że wzór (42) obowiązuje dla transmisji binarnej przy założeniu, że oba symbole pojawiają się z jednakowym prawdopodobieństwem.

Rysunek 39 Komplementarna funkcja błędu

3.3. Kod unipolarny i kod bipolarny

Dla każdej pary symboli  s_1\left(t\right),\ \ s_0(t)  można zastosować odbiornik optymalny. Prawdopodobieństwo błędnego odbioru symbolu Pe (42) będzie tym mniejsze, im większa jest różnica sygnałów, której miarą jest energia E10. Energię tę można zwiększyć przez zwiększenie amplitud obu symboli, wymagałoby to jednak zwiększenia mocy nadajnika. Moc nadajnika możemy obliczyć, dzieląc średnią energię symbolu transmisyjnego przez czas jego trwania:

S=\frac{E_1P_1+E_0P_0}{T} (43)

Tę samą moc zmierzymy na wyjściu kanału, gdyż zakładamy że kanał nie tłumi sygnału (rys.33). Jeżeli P_0=P_1=\frac{1}{2} , wówczas S=\frac{1}{2T}(E_1+E_0).  Problemem jest taki dobór symboli  s_1\left(t\right),\ \ s_0(t)  aby osiągnąć maksimum E10 przy danej mocy sygnału S.  Nietrudno wykazać, że taka para symboli powinna różnić się znakiem, sam kształt symbolu nie ma znaczenia:

s_1\left(t\right)=-s_0(t) (44)

Jest to bipolarny kod transmisyjny. Oba sygnały tego kodu mają jednakową energię E_1=E_0=ST. Sygnał różnicowy p(t) wynosi p(t)=s_1(t)-s_0(t)=2s_1(t), a jego energia jest 4 razy większa niż E1 (pamiętamy że energia jest całką z sygnału podniesionego do kwadratu):E_{10}=4E_1=4ST. Podstawiając E10 do wzoru (42), otrzymuje się wzór na stopę błędów kodu bipolarnego: 

P_e=\frac{1}{2}erfc{(}\sqrt{\frac{E_{10}}{4\eta}\ \ })=\frac{1}{2}erfc{(}\sqrt{\frac{ST}{\eta}}\ ) (45)

Stopa błędów maleje, gdy zwiększa się moc sygnału S i czas jego trwania T. W kodzie binarnym transmituje się 1/T symboli na sekundę i tyle samo bitów na sekundę. Oznacza to, że powolna transmisja jest bardziej odporna na szumy w kanale. Wzór (45) jest słuszny dla każdej pary sygnałów spełniających warunek (44): może to być para impulsów prostokątnych, ale także dwuwartościowa modulacja PSK (31).
Zajmijmy się teraz kodem unipolarnym, w którym jeden z symboli jest sygnałem równym zero:  s_0\left(t\right)=0.
Średnia moc sygnału wynosi  S=\frac{E_1P_1+E_0P_0}{T}=\frac{E_1P_1}{T}=\frac{E_1}{2T}  Sygnał różnicowy p\left(t\right)=s_1\left(t\right)-s_0\left(t\right)=s_1\left(t\right), w związku z tym ) E_{10}=E_1=2ST. Podstawiając E10 do wzoru (42), otrzymuje się wzór na stopę błędów kodu unipolarnego:

P_e=\frac{1}{2}erfc{(}\sqrt{\frac{E_{10}}{4\eta}}\ \ )=\frac{1}{2}erfc{(}\sqrt{\frac{ST}{2\eta}\ \ }) (46)

Porównując ze wzorem (45), zauważamy, że przy tej samej mocy nadajnika, szybkości transmisji i szumie w kanale, we wzorze (46) mamy mniejszy argument funkcji erfc, a więc większą wartość tej funkcji i większe prawdopodobieństwo błędnego odbioru symbolu. Reasumując, kod bipolarny jest bardziej odporny na szumy w kanale– rys,40.  Wzór (46) określa stopę błędów dla każdej modulacji binarnej, w której jeden z dwóch sygnałów jest równy zeru, np. modulacji OOK (rys.26).  

Rysunek 40 Stopa błędów kodu unipolarnego i bipolarnego

Pamiętamy, że porównywaliśmy kod unipolarny i bipolarny w warunkach, gdy oba sygnały były nadawane z jednakowym prawdopodobieństwem. Jedynie w tych warunkach obowiązuje zasada minimum odległości (37) i wzór (42) jest słuszny.  
Naruszenie tego warunku wydawałoby się korzystne w przypadku kodu unipolarnego. Jeżeli sygnał s_0\left(t\right)=0 pojawia się częściej niż s_1\left(t\right)\neq0, czyli P_1, wówczas średnia moc sygnału unipolarnego wynosi  S=\frac{E_1P_1+E_0P_0}{T}=\frac{E_1P_1}{T}.  Zmniejszając P1, można zwiększać energię E1, zachowując stałą moc nadajnika S. Większa energia to większa amplituda symbolu s_1\left(t\right), co sprawia, że łatwiej ten symbol odróżnić od symbolu zerowego.  Można się zatem spodziewać zmniejszenia stopy błędów Pe. Czytelnik będzie się mógł o tym przekonać, wykonując ćwiczenie symulacyjne (Ćw.4 – transmisja cyfrowa – kody transmisyjne).
Dlaczego nie stosujemy w praktyce tego sposobu zwiększania odporności na szumy w kanale? Odpowiedź podsuwa teoria informacji – naruszając równowagę prawdopodobieństw P1 i P0, zmniejszamy ilość transmitowanej informacji. 
Ilość informacji przenoszona przez jeden sygnał jest odwrotnie proporcjonalna do prawdopodobieństwa nadania tego sygnału. Jest to intuicyjnie jasne: rzadko występujące zdarzenie, gdy już wystąpi, dostarcza nam dużej porcji informacji. Definiując ilość informacji używa się też logarytmu, gdyż przy zaistnieniu dwóch niezależnych zdarzeń, ich prawdopodobieństwa się mnoży, a ilości informacji dodaje. Przy podstawie logarytmu równej 2, wyrażamy ilość informacji w bitach.  I tak, sygnał  s_0(t) przenosi {log}_2{\frac{1}{P_0}} bitów, a sygnał s_1(t) przenosi {log}_2{\frac{1}{P_1}} bitów.  Średnia ilość informacji przenoszona przez jeden symbol jest nazywana entropią i wyraża się wzorem:

H=P_0{log}_2{\frac{1}{P_0}}+P_1{log}_2{\frac{1}{P_1}}\quad \frac{bit}{symbol}] (47)

Biorąc pod uwagę, że P_0+P_1=1, otrzymujemy entropię jako funkcję P1. Wykres tej funkcji (rys.40) został otrzymany przez C. Shannona w połowie XX w. Przesyłamy 1 bit informacji na symbol transmisyjny, gdy oba sygnały pojawiają się z jednakowym prawdopodobieństwem. Gdy naruszymy tę równowagę, zmniejszamy ilość transmitowanej informacji. 

Rysunek 41 Entropia w transmisji binarnej

 

3.4. Odbiorniki modulacji wielowartościowej

Modulacje wielowartościowe zapewniają lepsze wykorzystanie pasma kanału. Pasmo ogranicza liczbę symboli przesyłanych w ciągu sekundy (baud rate), natomiast z wartościowością związana jest liczba bitów zakodowanych w jednym symbolu (\log_2{M} bitów na symbol w modulacji M-wartościowej). 
Jeżeli symbole s_1\left(t\right),s_2\left(t\right),\ldots\ s_M\left(t\right) pojawiają się z jednakowym prawdopodobieństwem (P_i=\frac{1}{M}), wówczas obowiązuje zasada minimum odległości: poszukuje się symbolu najbliższego (w sensie odległości euklidesowej) odebranemu sygnałowi v(t). Dla M=4 – wartościowej modulacji i sygnałów reprezentowanych przez N=2 próbki pokazano to na rys. 42. Granice obszarów decyzyjnych przebiegają w równej odległości od najbliższych punktów reprezentujących niezaszumione symbole.

Rysunek 42 Zasada minimum odległości – modulacja 4-wartościowa

Zasada minimum odległości prowadzi również do innej równoważnej struktury odbiornika optymalnego – do korelatora. Przetwarzając sygnały ciągłe, wygodniej jest sformułować tę zasadę w następujący sposób:

min{|}|v-s_i||^2\ \rightarrow\ \int_0^T[v(t)-s_i(t)]^2dt (48)

Minimalizację  wyrażenia \int_{0}^{T}{v^2\left(t\right)dt}-2\int_{0}^{T}{v\left(t\right)s_i\left(t\right)dt}+\int_{0}^{T}[s_i(t)]^2dt można zastąpić maksymalizacją wyrażenia o przeciwnym znaku. Po usunięciu pierwszej całki (nie zależy ona od symbolu) i podzieleniu przez 2 otrzymujemy regułę decyzyjną:

max{\{{}\int_{0}^{T}{v\left(t\right)s_i\left(t\right)dt}-\frac{1}{2}E_i}|\} (49)

Jej implementacja jest identyczna jak dla modulacji binarnej – porównaj rys. 36 i rys.43.

Rysunek 43 Odbiornik optymalny w układzie korelatora

Powiększenie wartościowości prowadzi do szybszej transmisji, jednak czyni modulację bardziej wrażliwą na szumy w kanale. Można się o tym przekonać, wykonując ćwiczenie 5: transmisja cyfrowa – modulacje cyfrowe. 

Wartościowość modulacji może sięgać kilkudziesięciu (rys.30) a nawet kilkuset, szczególnie w przypadku modulacji QAM i AM-PM. Odbiornik z rys.43 musiałby wówczas zawierać wielką liczbę elementarnych układów korelujących. Dla sygnałów typu PSK, QAM, AM-PM, opisanych wzorem (32): s_i\left(t\right)=A_i\cos{(2\pi f_0t+\varphi_i)} liczbę korelatorów można zredukować do dwóch. Taki odbiornik nazywa się odbiornikiem kwadraturowym (rys.44). Jest on równoważny układowi korelatora w tym sensie, że dla tych samych sygnałów wejściowych zwraca te same decyzje.  Odbiornik kwadraturowy oblicza korelacje sygnału wejściowego z ortogonalną parą sygnałów cos i (-sin).  Gdyby nie było szumu, to otrzymana para korelacji (x,y) stanowiłaby współrzędne punktu konstelacji, opisującego odebrany sygnał: \left(x,y\right)=(A_i\cos{\varphi_i},A_i\sin{\varphi_i}) – rys.28. Pokażmy, że x=A_i\cos{\varphi_i} .


Rysunek 44 Odbiornik kwadraturowy

Na wejście podajemy i-ty symbol bez szumu:  v(i)=A_icos{(}2\pi f_0t+φi). Po pomnożeniu przez  \frac{2}{T}cos{(}2\pi f_0t) otrzymamy, stosując wzór  \cos{\alpha\ \cos{\beta}}=\frac{1}{2}[\cos{\left(\alpha-\beta\right)+\cos{\left(\alpha+\beta\right)}}] :

v(t)\frac{2}{T}cos{(}2\pi f_0t)=\frac{2A_i}T{cos(2πf_0t+φi)cos(2πf_0t)}=\\=\frac{2A_i}{T}\frac{1}{2}{{cos\varphi}_i{+}cos{(}4\pi f_0t+φi)}=\frac {A_i}{T}cosφi+\frac{A_i}{T}cos(4πf_0t+φi)  

Całkując sygnał o wartości stałej :  A_i/Tcos{(}\varphi_i), otrzymamy A_icos{(}\varphi_i). Jeśli w przedziale czasu 0-T znajduje się całkowita liczba okresów sygnału cos{(}4\pi f_0t) , wówczas całka z tego sygnału jest równa zeru. Ostatecznie otrzymujemy:

x=\int_0^Tv(t)\frac2Tcos(2πf_0t)dt=\int_0^T\frac{A_i}{T}cos(φi)dt+0=A_icosφi  

W podobny sposób można pokazać, że y=A_i{sin\varphi}_i{\ }. Gdy na wejście podamy symbol z szumem, otrzymamy parę korelacji x, y , czyli zmierzone z błędem współrzędne punktu konstelacji odebranego symbolu. Należy teraz obliczyć odległość euklidesową od punktu (x,y) do wszystkich punktów konstelacji stosowanej modulacji, wybrać najmniejszą odległość i podjąć decyzję. 
Na rys.45 pokazano wyniki pomiaru (x,y) dla 25000 symboli modulacji 16QAM. Przy silnym szumie(SNR=6 dB)  odbiornik optymalny nie działa prawidłowo – stopa błędów przekracza 0.1. 

Rysunek 45 Odbiór 25000 symboli 16QAM z wykorzystaniem odbiornika kwadraturowego

Stopa błędów odbiornika optymalnego modulacji wielowartościowej zależy od wartościowości modulacji (in większa wartościowość, tym trudniej wybrać właściwy symbol), konstelacji punktów i szumu. Konstelacja pokazana na rys 45 dość równomiernie wypełnia płaszczyznę, można jednak znaleźć lepszą (patrz ćwiczenie symulacyjne nr 5). Porównując modulacje, należy tez wziąć pod uwagę pasmo i efektywność widmową. Na rys.46 porównano szereg modulacji, notując efektywność widmową i SNR na wyjściu kanału przy ustalonej stopie błędów 0.0001.Obniżenie dopuszczalnej stopy błędów w niewielkim stopniu wpływa na wyniki. SNR wyrażono jako stosunek energii przypadającej na 1 bit do gęstości mocy szumu na wyjściu kanału (28).  Zbadano następujące modulacje:

  • Modulacje amplitudy AM (zwane też ASK – Amplitude Shift Keying). Wszystkie punkty konstelacji leżą na prostej przechodzącej przez początek układu współrzędnych. Dopuszcza się rozmieszczanie punktów po obu stronach początku układu. 
  • Modulacje fazy (kluczowanie fazy - PSK – Phase Shift Keying). Wszystkie symbole mają tę samą amplitudę, punkty konstelacji leżą więc na okręgu. 
  • QAM (Quadrature Amplitude Modulation)–rys.30.

Modulacje wielowartościowe oferują lepszą efektywność widmową, ale wymagają lepszego kanału (wyższe wartości SNR). QAM jest bardziej odporna na szumy w kanale niż ASK i PSK, dzięki wykorzystaniu różnych amplitud i faz swoich konstelacjach. Nie jest możliwe jednoczesne osiągnięcie wysokiej efektywności widmowej i odporności na szum w kanale. Na lewo od krzywej oznaczonej jako granica Shannona nie istnieje żadna modulacja. Jeśli stosunek energii na bit do gęstości mocy szumu jest niższy niż \frac{1}{\sqrt2} (czyli -1.505 dB), wówczas transmisja z niską stopą błędów w ogóle nie jest możliwa. O granicznych parametrach transmisji będzie mowa w końcowej części tego modułu.

Rysunek 46 Porównanie modulacji cyfrowych

 

3.5. Modulacja DPSK

Optymalne odbiorniki modulacji PSK są odbiornikami koherentnymi – wymagają dokładnego odtworzenia sygnałów referencyjnych (w odbiorniku kwadraturowym są to sygnały \cos{(2\pi f_0t)} i -\sin{(2\pi f_0t)} zsynchronizowane fazowo z sygnałem v(t) docierającym do odbiornika). Synchronizację najczęściej osiąga się przed rozpoczęciem transmisji danych. W przypadku jej utraty można ją odzyskać śledząc ciąg odebranych symboli. Załóżmy, że w dwuwartościowej (binarnej) modulacji PSK odbierane są symbole \pm\sin{(2\pi f_0t)}, rys.47.  Synchronizator może zacząć śledzić jeden z tych dwóch sygnałów. W przypadku gdy jest to sygnał niewłaściwy, wszystkie decyzje odbiornika będą błędne. Aby temu zapobiec, koduje się kolejne bity nie bezpośrednio w fazie, ale w różnicy faz dwóch sąsiednich sygnałów. Np. zachowanie tej samej fazy początkowej oznacza logiczne 1, a przesunięcie fazy o 180 stopni - logiczne 0. W ten sposób zdefiniowaliśmy dwuwartościową różnicową modulację fazy (DPSK – Differential Phase Shift Keying).  W modulacji DPSK, nawet przy błędnej fazie sygnału referencyjnego, jesteśmy w stanie wykryć, czy w odbieranym sygnale wystąpiło, czy nie wystąpiło przesunięcie fazy. 

Rysunek 47 Sygnał binarnej modulacji PSK, właściwy sygnał referencyjny i sygnał referencyjny przesunięty o 180 stopni

Odbiornik DPSK może być zrealizowany jako odbiornik koherentny lub niekoherentny. Odbiornik koherentny odtwarza sygnał referencyjny (właściwy lub przesunięty), rozpoznaje kolejne sygnały \pm\sin{(2\pi f_0t)} i porównuje dwie kolejne decyzje: jeśli są identyczne to wyprowadza wartość logiczna 1 a jeśli przeciwne – wyprowadza 0. W pierwszym etapie działa tak, jak odbiornik PSK. Gdyby nie było szumu w kanale, to rozpoznawałby wszystkie sygnały właściwie lub wszystkie niewłaściwie. Porównywanie dwóch kolejnych decyzji pozwoliłoby na bezbłędny odbiór w systemie DPSK. Gdy w kanale jest szum, wówczas rozpoznawanie kolejnych sygnałów odbywa się z błędem, którego prawdopodobieństwo podano we wzorze (45): P_e^{PSK}=\frac{1}{2}erfc{(}\sqrt{\frac{ST}{\eta}} ).  Pojedynczy, izolowany błąd doprowadzi do powstania dwóch błędów w systemie DPSK, ponieważ błędna decyzja będzie porównana z poprzednią i z następną.  Dokładny wzór na stopę błędów koherentnego odbiornika DPSK otrzymamy, biorąc pod uwagę dwa zdarzenia, które do tego błędu prowadzą:

  • Pierwszy z porównywanych sygnałów jest rozpoznany dobrze, a drugi źle. Prawdopodobieństwo takiego zdarzenia wynosi  (1-P_e^{PSK})P_e^{PSK}
  • Pierwszy z porównywanych sygnałów jest rozpoznany źle, a drugi dobrze. Prawdopodobieństwo takiego zdarzenia wynosi  P_e^{PSK}(1-P_e^{PSK})

Pozostałe przypadki, tzn. dwa sygnały rozpoznane źle i dwa sygnały rozpoznane dobrze – nie prowadzą do błędu w systemie DPSK. 
Ostatecznie otrzymujemy prawdopodobieństwo błędu koherentnego odbiornika DPSK:

P_e^{DPSK}=2P_e^{PSK}\left(1-P_e^{PSK}\right)\approx2P_e^{PSK}=erfc\sqrt{\frac{ST}{\eta}} (50)

Odbiornik niekoherentny mierzy bezpośrednio przesunięcie fazowe dwóch sąsiednich sygnałów. Może być zrealizowany jako odbiornik kwadraturowy. Błąd fazowy b w odtworzeniu sygnałów referencyjnych powoduje rotację odebranych punktów (x,y), co nie przeszkadza w prawidłowym pomiarze kąta Df między kolejnymi sygnałami – rys.48.Bezpośredni pomiar kąta jest jednak obarczony większym błędem niż błąd rozpoznawania jednego z dwóch sygnałów. Skutkiem tego jest większa stopa błędów odbiornika niekoherentnego niż koherentnego:

P_e=\frac{1}{2}e^{-\frac{E_b}{\eta}}=\frac{1}{2}e^{-\frac{ST}{\eta}} (51)

 

Rysunek 48 Bezpośredni pomiar kąta w odbiorniku kwadraturowym DPSK

Ćwiczenie symulacyjne nr 5 pozwoli Państwu na porównanie modulacji PSK, oraz DPSK z odbiorem koherentnym i niekoherentnym. 

4. Graniczne parametry transmisyjne

W p.2.1 zostały sformułowane kryteria jakości transmisji cyfrowej. Pożądane byłoby wąskie pasmo zmodulowanego sygnału, a jednocześnie duża szybkość modulacji (liczba symboli transmitowanych w ciągu sekundy). Niestety są to wymagania sprzeczne: symbole o krótkim czasie trwania mają szersze pasmo (rys.8, Moduł 1). Pożądana byłaby duża przepływność binarna (liczba bitów na sekundę) i niska stopa błędów. Tu też jest sprzeczność - ze wzoru 45 dla kodu bipolarnego wynika, że krótki czas przeznaczony na transmisję bitu (T)  prowadzi do wzrostu prawdopodobieństwa błędnego odbioru. Cenna byłaby duża efektywność widmowa, ale jest ona osiągalna tylko przy wysokim stosunku mocy sygnału do szumu w kanale (rys.46). Poza tym, dążymy do jak największej odporności na szumy w kanale, czyli jak najwyższej jakości odtworzonej po stronie odbiorczej mowy czy muzyki (maksimum SNR0) przy danym stosunku mocy sygnału do szumu na wyjściu kanału (SNR). 
Rozpatrzmy dokładniej wymienione problemy, zaczynając od szybkości modulacji (liczba symboli na sekundę). 
Niech nadawane symbole mają postać krótkich impulsów (rys.49).  Przesyłana informacja jest zakodowana w ich amplitudach: x_0,\ x_1,x_2\ldots  Jest to tzw. PAM (pulse amplitude modulation). Impulsy są nadawane co T sekund (w praktyce będą to mikro- lub nanosekundy). Stosując odpowiedni filtr, można tym impulsom nadać dowolny kształt, pod warunkiem nieprzekraczania czasu trwania równego T. 

Rysunek 49 Sygnał nadawany (PAM)

Kanał transmisyjny jest dolnopasmowy o paśmie B [Hz]. Modelem kanału jest filtr idealny o charakterystyce częstotliwościowej podanej na rys.50. 

 

Rysunek 50 Charakterystyka częstotliwościowa kanału

Odpowiedź impulsowa kanału jest odwrotną transformatą Fouriera charakterystyki częstotliwościowej (patrz zad.4, p.8.1, Moduł1). 

 

Rysunek 51 Odpowiedź impulsowa kanału

Na rys.52 pokazano odpowiedź kanału na 3 impulsy z rys.49 (zachowano ten sam kolor impulsu i odpowiedzi kanału).  Oczywiście na wyjściu kanału obserwujemy sumę tych sygnałów i w ogólnym przypadku nie można odczytać amplitud x_0,\ x_1,x_2\ldots  Jest to tzw. interferencja międzysymbolowa (ISI – Intersymbol Interference). 

 

Rysunek 52 Odpowiedź kanału na sygnał PAM

Interferencję międzysymbolową można zneutralizować, jeśli wykorzysta się regularne przejścia przez zero odpowiedzi impulsowej kanału (rys.51). Nadając kolejne impulsy w momentach  0,\frac{1}{2B},\frac{2}{2B},\frac{3}{2B}\ldots , umożliwiamy idealną  demultipleksację, gdyż w każdym z tych momentów nie zeruje się reakcja kanału tylko na jeden impuls (rys.53). 

 

Rysunek 53 Odpowiedź kanału na sygnał PAM gdy T=1/(2B)

Nadając impulsy co 1/(2B) sekund, transmitujemy 2B symboli na sekundę. W ten sposób otrzymujemy Twierdzenie Nyquista o szybkości modulacji:  W kanale o paśmie B [Hz] można przesyłać do  2B impulsów (symboli) na sekundę bez interferencji międzysymbolowej. 
W twierdzeniu Nyquista nie mówi się nic o szumie w kanale i błędach w transmisji. Gdy nie ma szumu w kanale, wówczas można transmitować (oczywiście w teorii, nie w praktyce) impulsy o nieskończonej liczbie różnych amplitud. Wówczas nawet jeden impuls może przenieść nieskończoną ilość informacji. 
Szum w kanale istnieje zawsze, w związku z tym rozpatrzmy problem jego wpływu na stopę błędów w transmisji.  Weźmy dla przykładu wzór na stopę błędów w kodzie bipolarnym (45): P_e=\frac{1}{2}erfc{(}\sqrt{\frac{ST}{\eta}}\ ). Czy można otrzymać stopę błędów dowolnie bliską zeru? Teoretycznie jest to możliwe, gdy \frac{ST}{\eta}\rightarrow\infty, a więc nieograniczonej mocy nadajnika lub nieskończenie długim symbolu, a więc szybkości transmisji dążącej do zera. W swoich pracach opublikowanych w połowie XX w. Claude Shannon wykazał, że można jednak uczynić stopę błędów dowolnie bliską zeru, jeśli rozpoznawać będziemy całą przesłaną wiadomość a nie pojedyncze symbole. Ideę Shannona realizuje się w praktyce, tworząc kody zabezpieczające przed błędami, wymagające opóźnionej decyzji [2]. 
Jakie warunki muszą być spełnione, aby udało się zredukowanie stopy błędów do wartości bliskiej zeru? 
Załóżmy że nadajemy pewną wiadomość składającą się z dużej liczby (k) impulsów: s_1,s_2,\ldots,s_k.  Potrzebny nam będzie nadajnik o mocy S=E(s_i^2). Kanał ma pasmo B i charakterystykę częstotliwościową jak na rys.50. 
Aby zneutralizować interferencje międzysymbolowe, impulsy nadajemy co 1/(2B), czyli przesyłamy 2B impulsów na sekundę. Na wyjściu kanału obserwujemy sygnał użyteczny s(t) jak na rys.53 i szum n(t) o mocy N=\eta B. Ich suma to sygnał v(t)=s(t)+n(t), który podajemy na wejście odbiornika. Pasmo tego sygnału wynosi B, więc można przeprowadzić jego pomiar, pobierając 2B próbek na sekundę (twierdzenie o próbkowaniu). Jednocześnie mamy 2B impulsów na sekundę, więc wystarczy jedna próbka na nadany impuls. Cała wiadomość zawiera zatem k próbek, a jej transmisja trwa k/(2B) sekund. 
Tak więc wiadomość możemy przedstawić w postaci k-wymiarowego wektora \overset{-}s=s_1,s_2,\ldots,s_k, a szum w postaci wektora  \overset{-}n=n_1,n_2,\ldots,n_k. Sygnał odebrany z szumem to wektor \overset{-}v=\overset{-}s+\overset{-}n
Typowa wiadomość wykorzystuje całą moc nadajnika, czyli  S\approx \frac{1}{k}\sum\nolimits ^{k}_{i=1} (s_{i} )^{2} =\frac{1}{k}|| \overset{-}{s}|| ^{2}. Stąd kwadrat normy wektora \overset{-}{s} można oszacować jako || \overset{-}{s}|| ^{2}\approx kS. Tak więc wektory wiadomości leżą w odległości zbliżonej do \sqrt{kS} od początku układu współrzędnych.
Moc szumu można oszacować jako  N\approx\frac{1}{k}\sum_{i=1}^{k}{(n_i)^2}=\frac{1}{k}||\overset{-}{n}||^{2}. Oszacowanie jest tym dokładniejsze, im większe k. Stąd  || \overset{-}{n}||^{2}=\sum_{i=1}^{k}(n_i)^2≈kN.  Wektory szumu leżą przy powierzchni kuli o promieniu \sqrt{kN}. Rysunek 54 to rysunek poglądowy przedstawiający wektory wiadomości, szumu i sygnału odebranego. 

Rysunek 54 Wektory wiadomości, szumu i sygnału odebranego oraz ich normy

Istotne znaczenie ma norma wektora szumu w przestrzeni wysokowymiarowej. Zmienna losowa  ||\overset{-}n||^2=\sum^k _{i=1}(n_i)^2 ma rozkład chi kwadrat o k stopniach swobody, gdyż próbki szumu to nieskorelowane zmienne o rozkładzie gaussowskim [2]. Na rys.55 pokazano rozkłady znormalizowanej zmiennej  \frac{1}{k}||\overset {-}n||^2 dla różnych wartości k. Wariancja szumu wynisi 1. Norma wielowymiarowych wektorów szumu dąży do wartości stałej. 

Rysunek 55 Gęstość prawdopodobieństwa zmiennej  \frac{1}{k}\sum_{i=1}^{k}n_i^2

Ta obserwacja ma znaczenie kluczowe. Wskazuje ona, że nadając pewną wiadomość  \overset{-} s, odbierzemy wektor leżący przy powierzchni kuli o promieniu \sqrt{kN} i środku w punkcie  \overset{-} s – rys.56.

Rysunek 56 Wynik pomiaru wiadomości leży na powierzchni otaczającej ją kuli

W odbiorze stosujemy zasadę minimum odległości. Jeśli zatem transmitowane wiadomości będą leżały dostatecznie daleko od siebie, to otaczające je kule szumu nie zetkną się. Każdy wynik pomiaru leżący na powierzchni kuli wskaże nam środek kuli jako najbliższą wiadomość i będzie to decyzja nieobciążona błędem. Tak więc możliwa jest bezbłędna transmisja wiadomości pod warunkiem że będą długie (duża wartość k).
Pozostaje problem, ile takich niezachodzących na siebie kul szumu można rozmieścić. Jesteśmy ograniczeni mocą nadajnika.  Sygnał odebrany v(t)=s(t)+n(t) ma moc średnią S+N, a więc wektory sygnałów odebranych \overset {-} v=\overset {-} s+\overset {-} n leżą w odległości około \sqrt{kS+kN} od początku układu współrzędnych. Wewnątrz kuli o promieniu  R=\sqrt{kS+kN} musimy rozmieścić kule szumu o promieniu r=\sqrt{kN}. Ile takich kul można rozmieścić? Tyle, ile wynosi ich stosunek objętości, a więc  [\frac{R}{r}]^k=[\frac{S+N}{N}]^{k/2}=[1+\frac SN]^{k/2}. Tyle różnych wiadomości można odebrać bezbłędnie. Oznacza to bezbłędne przesłanie {log}_2{\ }[1+\frac{S}{N}]^{k/2}=\frac k2log2 [1+\frac SN] bitów.
Transmisja całej wiadomości trwa \frac{k}{2B} sekund, a więc w ciągu sekundy wysyłamy (bez interferencji międzysymbolowych) C=\frac{2B}{k}\ \ \frac{k}{2}{log}_2{\left[1+\frac{S}{N}\right]=B\ }{log}_2{\left[1+\frac{S}{N}\right]} bitów. W sposób uproszczony pokazaliśmy uzasadnienie bardzo ważnego twierdzenia w telekomunikacji:
Tw. Shannona o przepustowości kanału
W kanale o szerokości pasma B, przy stosunku mocy sygnału do szumu na wyjściu kanału równym \frac{S}{N}=\frac{S}{\eta B} można przesyłać z dowolnie niską stopą błędów 

C=B\ log_2[1+\frac SN]=B\ log_2[1+\frac {S}{ηB}]\ \ bit/s (52)

Tę największą szybkość transmisji nazywa się przepustowością kanału (channel capacity). Im większy stosunek mocy sygnału do mocy szumu, tym większa przepustowość. Kanał bez szumu ma przepustowość nieskończoną, gdyż jesteśmy w stanie bezbłędnie rozpoznać nieskończoną liczbę różnych sygnałów.

Gdy  \frac{S}{N}\rightarrow\infty  to  C\rightarrow\infty (53)

Ciekawszym problemem jest wpływ pasma na przepustowość kanału. Pasmo B występuje we wzorze (52) dwukrotnie. Zmniejszając pasmo do zera, zmniejszamy szum na wyjściu kanału (N=\eta B), ale maleje również liczba transmitowanych symboli (twierdzenie Nyquista). Stosując regułę de l’Hospitala można obliczyć granicę (pozostawiam to czytelnikowi jako ćwiczenie): 

Gdy  B\rightarrow0  to  C\rightarrow0 (54)

Tak więc zmniejszanie pasma nie jest korzystne z punktu widzenia szybkości transmisji.  Z drugiej strony, poszerzenie pasma zapewnia większą przepustowość kanału mimo większej mocy szumu:

Gdy  B\rightarrow\infty\   to  C=B\ log_2[1+\frac S{ηB}]=\frac Sηlog_2(1+\frac S{ηB})^\frac{ηB}{S}=\frac Sηlog_2\ e

gdyż \lim \limits _{x\rightarrow\infty} (1+\frac{1}{x})^x=e

(55)

Pomimo że noc szumu dąży do nieskończoności (N=\eta B\rightarrow\infty) , kanał jest zdolny do bezbłędnej transmisji z maksymalną przepływnością określoną wzorem (55). 
Porównajmy wąskopasmową (pasmo B1 Hz) i szerokopasmową (pasmo B2 Hz) transmisję w kanale, charakteryzującym się gęstością mocy szumu W/Hz  - rys.57. Mamy do dyspozycji nadajnik o mocy S watów i dwie możliwości:  

  1. Transmisja wąskopasmowa: wykorzystujemy wąskie pasmo B1 Hz aby otrzymać wysoki stosunek mocy sygnału do szumu  SNR=\frac S{ηB_1}
  2. Transmisja szerokopasmowa: wykorzystujemy szersze pasmo B2 >> B1 Hz

 

Rysunek 57 Wykorzystanie kanału w transmisji wąskopasmowej (po lewej) i szerokopasmowej (po prawej)

Rysunek 57 sugeruje, że bardziej odporna na szum jest transmisja wąskopasmowa, gdyż gęstość mocy sygnału jest większa niż gęstość mocy szumu. Z kolei w transmisji szerokopasmowej gęstość mocy sygnału znajduje się poniżej gęstości mocy szumu. Wynika to z faktu, że oba prostokąty na rys.57 mają to samo pole równe mocy sygnału S. 
Wrażenie to jest jednak fałszywe – większą szybkość transmisji otrzymujemy w przypadku szerokopasmowym. Proszę się o tym przekonać przeglądając rozwiązane zadanie 11,.p.5.1

Ważnym parametrem transmisyjnym jest efektywność widmowa. Mówi ona o tym, ile bitów na sekundę  można przesłać w kanale o szerokości pasma 1 Hz. Przypomnijmy sobie oznaczenia:

  • Rb – szybkość transmisji [bit/s],   
  • Rb/B – efektywność widmowa [bit/s/Hz]
  • Eb– energia zużyta na transmisję 1 bitu [J=Ws]

Załóżmy, że osiągnęliśmy graniczną szybkość transmisji równą przepustowości kanału:
R_b=C=B\ {log}_2{\left[1+\frac{S}{N}\right]}=B\ {log}_2{\left[1+\frac{E_bR_b}{\eta B}\right]}. Przekształcając wzór Shannona (52) wykorzystaliśmy fakt, .ze moc sygnału S jest to energia na sekundę, czyli energia na bit (Eb) razy liczba bitów na sekundę (Rb). Przekształcając dalej, otrzymuje się: \frac{R_b}{B}={log}_2{\left[1+\frac{E_bR_b}{\eta B}\right]}  a w końcu 

 \frac{E_b}{\eta}=\frac{B}{R_b}\left[2^\frac{R_b}{B}-1\right] (56)

Wzór (56) określa najmniejszą wartość SNR (wyrażoną jako stosunek energii na bit i gęstości mocy szumu) jaka jest niezbędna do zapewnienia bezbłędnej transmisji z efektywnością widmową Rb/B. Jest to granica Shannona przedstawiona na Rys. 46. Gdy pasmo kanału rośnie do nieskończoności, wówczas efektywność widmowa maleje do zera.  Graniczna wartość SNR wynosi w tej sytuacji ln(2), czyli około -1.6 dB: 

Gdy   B→∞    to   \frac {E_b}{η} → ln(2)    (57)

Przy niższych wartościach SNR nie może być mowy o bezbłędnej transmisji. 

Z twierdzeniem Shannona o przepustowości kanału wiąże się graniczna odporność na zakłócenia. Przypominamy że chodzi o zależność  wartości SNR na wyjściu odbiornika od wartości SNR na wyjściu kanału: SNR0=f(SNR). Pożądana byłaby jak najwyższa jakość sygnału na wyjściu odbiornika (np. mowy czy muzyki) przy niskiej jakości sygnału na jego wejściu: {SNR}_0\gg SNR. Niestety, istnieje graniczna wartość odporności na zakłócenia, która wyraża się wzorem:

 SNR_0\le\left[1+\frac{f_M}{B}SNR\right]^\frac{B}{f_M}-1    (58)

Parametrem jest tu współczynnik poszerzenia pasma B/fM czyli pasmo sygnału zmodulowanego przez pasmo sygnału modulującego (np. sygnału mowy). We wzorze (58) stosunek mocy sygnału do szumu w kanale jest zdefiniowany następująco: \ SNR=\frac{S}{N}=\frac{S}{\eta\ f_M}. Jako moc szumu w kanale przyjęto moc szumu przypadającą na pasmo sygnału przed modulacją. 
Rodzinę funkcji opisaną wzorem (58) pokazano na rys.58. Wąskopasmowe modulacje nieposzerzające pasma (B=fM) mają niską odporność na szum w kanale (SNR0=SNR). Im szersze pasmo, tym większa odporność. Nie wzrasta ona jednak do nieskończoności, nawet przy nieograniczonym paśmie istnieje krzywa graniczna – rys.58.  Ograniczenie SNR sygnału wyjściowego (np. muzyki transmitowanej bezprzewodowo) odnosi się do każdej transmisji – nawet analogowej. Potwierdza to analiza porównawcza modulacji AM i FM. Szerokopasmowa modulacja FM jest bardziej odporna na szumy w kanale. 

Rysunek 58 Graniczna odporność na zakłócenia

Wzór (58) można wyprowadzić z twierdzenia Shannona o przepustowości kanału. 
Graniczna szybkość transmisji w kanale:

C=B\ {log}_2{(}1+\frac{S}{\eta B})=B\ {log}_2{(}1+\frac{f_M}{B}\frac{S}{\eta f_M})=B\ {log}_2{(}1+\frac{f_M}{B}SNR)  

Szybkość transmisji na wyjściu odbiornika: C_0=f_M{log}_2{(}1+SNR_0).
Liczba bitów na wyjściu odbiornika nie może być większa niż na jego wejściu, gdyż odbiornik nie jest źródłem informacji, on tylko przetwarza otrzymaną informację: C_0\le C. Przepisując tę nierówność otrzymujemy kolejno:

{log}_2{(}1+SNR_0)\le\frac{B}{f_M}{log}_2{(}1+\frac{f_M}{B}SNR) 

{log}_2{(}1+SNR_0)\le{log}_2{(}1+\frac{f_M}{B}SNR)^\frac{B}{f_M} 

1+SNR_0\le(1+\frac{f_M}{B}SNR)^\frac{B}{f_M} 

SNR_0\le(1+\frac{f_M}{B}SNR)^\frac{B}{f_M}-1 

 

Ostatnia postać jest identyczna ze wzorem (58).

 

5. Biblioteka

Z Modułem 2 są związane dwa ćwiczenia symulacyjne: Transmisja cyfrowa – kody transmisyjne i Transmisja cyfrowa – modulacje cyfrowe. Pierwsze z nich zapoznaje studenta z kodem transmisyjnym unipolarnym i bipolarnym. Można porównać odporność obu kodów na szum w kanale w sytuacji, gdy oba symbole nadawane są z tym samym prawdopodobieństwem i różnym prawdopodobieństwem. Drugie ćwiczenie wprowadza w binarne i wielowartościowe modulacje oparte na fali nośnej: modulacje PSK, QAM, DPSK. Zapoznaje studenta z działaniem odbiornika kwadraturowego, pozwala na konstruowanie własnych modulacji. 
Laboratorium i materiały pomocnicze opisano w Module 4.