KKRRiT2004 A Piatyaszek [Wybrane aspekty badań transmisji wielotonowej w paśmie radiowym o szerokości 25kHz]

background image

Adam Piątyszek
Piotr Tyczka
Instytut Elektroniki i Telekomunikacji
Politechnika Poznańska, ul. Piotrowo 3a, 60-965 Poznań
[adam.piatyszek, tyczka]@et.put.poznan.pl

Warszawa, 16-18 czerwca 2004

WYBRANE ASPEKTY BADAŃ TRANSMISJI WIELOTONOWEJ

W PASMIE RADIOWYM O SZEROKOŚCI 25 kHz

*

*

Praca wykonana i sfinansowana w ramach grantu KBN nr 0 T00A 011 24

Streszczenie: W artykule przedstawiono częściowe wyniki
pracy nad projektem modemu radiowego dla transmisji w
pasmie o szerokości 25 kHz z wykorzystaniem modulacji
OFDM. Omówiono założenia projektowe dotyczące
modemu i wynikający z nich dobór parametrów oraz
scharakteryzowano wybrane elementy opracowanego
systemu symulacyjnego modemu. Przedstawiono również
uzyskane wyniki badań symulacyjnych nad jakością
transmisji modemu w wybranych kanałach z zanikami.

1. WSTĘP

Projektowanie współczesnych systemów radioko-

munikacyjnych odbywa się zwykle w ramach kilku
następujących po sobie etapów. Pierwszym z nich jest
analiza ogólnych założeń tworzonego systemu i dobór
odpowiednich parametrów. Drugim etapem jest zwykle
opracowanie modelu symulacyjnego projektowanego
systemu w celu weryfikacji przyjętych założeń realiza-
cyjnych. Kolejny etap projektowania, w przypadku re-
alizacji za pomocą procesora DSP lub układów progra-
mowalnych, to opracowanie oprogramowania poszcze-
gólnych bloków systemu i ich weryfikacja w odpowied-
nio przygotowanym środowisku testowym z wykorzy-
staniem symulatora sprzętowego i prób terenowych.
Przetestowane oprogramowanie może być podstawą
implementacji modemu.

Realizowany projekt wielotonowego modemu ra-

diowego jest przykładem tego typu zagadnienia. Głów-
nym celem projektu jest wyznaczenie parametrów
transmisyjnych projektowanego modemu oraz opraco-
wanie oprogramowania dla procesora sygnałowego,
który będzie realizował nadajnik i odbiornik tego mo-
demu.

W artykule przedstawiono dotychczasowe rezultaty

dotyczące pierwszych dwóch etapów projektowania
modemu radiowego. W punkcie drugim przeanalizowa-
no w skrócie wpływ założeń projektowych na dobór
odpowiednich technik i parametrów użytych w realizacji
systemu symulacyjnego. Punkt trzeci stanowi opis wy-
branych elementów opracowanego modelu symulacyj-
nego. Wyniki symulacyjne zostały przedstawione i omó-
wione w punkcie czwartym natomiast punkt piąty to
krótkie podsumowanie niniejszego artykułu.

2. DOBÓR PARAMETRÓW MODEMU

Głównym założeniem realizacyjnym projektowa-

nego modemu radiowego jest zastosowanie modulacji

OFDM (ang. Orthogonal Frequency Division Multi-
plexing
). Transmisja ma się odbywać w zakresie często-
tliwości radiowych od 25 MHz do 512 MHz, przy zało-
żonej szerokości pasma równej 25 kHz. Środowisko
pracy systemu to teren otwarty a maksymalny zasięg
transmisji to ok. 5-10 km. W trakcie ustalania parame-
trów należało również uwzględnić fakt, że maksymalna
prędkość poruszania się nadajnika względem odbiornika
może wynosić 100 km/h. Ponadto, w sprzyjających
warunkach system powinien pozwalać na transmisję z
przepływnością równą 64 kbit/s.

Studia literaturowe i rozważania teoretyczne wska-

zują, że jedno- lub dwudrogowy model kanału z zani-
kami Rayleigha może wystarczająco dobrze charaktery-
zować rzeczywisty kanał występujący w wyżej określo-
nych warunkach. W badaniach symulacyjnych współ-
czynniki wagowe linii opóźniającej kanału charakteryzu-
jące każdą ze ścieżek są generowane niezależnie według
metody opisanej w [1]. Przyjęto, że opóźnienie drugiej
ścieżki względem pierwszej wynosi 13,3 µs (dwa odstę-
py próbkowania sygnału w kanale), natomiast względne
tłumienie drugiej ścieżki jest równe 8,6 dB.

Rozważmy przypadek transmisji w kanale dwudro-

gowym. Zgodnie z [2], 50% pasmo koherencji można
wyznaczyć z następującego wzoru:

RMS

c

B

τ

5

1

(1)

W naszym przypadku otrzymujemy wartość ok. 46 kHz
(

τ

RMS

= 4,35

µs). Ponieważ pasmo sygnału wynosi

25 kHz można przyjąć, że użyty w symulacjach dwu-
drogowy model kanału nie modeluje kanału selektyw-
nego, nie jest też charakterystyczny dla kanału o zani-
kach typowo płaskich.

Rozważmy najbardziej niekorzystną sytuację, w

której prędkość poruszania się nadajnika względem
odbiornika jest maksymalna a częstotliwość nośna rów-
na jest 500 MHz. Otrzymujemy wtedy maksymalną
wartość częstotliwości Dopplera f

d

= 46 Hz. Czas kohe-

rencji kanału [2] można obliczyć ze wzoru:

d

c

f

T

π

16

9

(2)

i wynosi on wtedy ok. 9 ms. Minimalny odstęp próbko-
wania wynikający z szerokości pasma wynoszącej
25 kHz wynosi 40

µs. W porównaniu z oszacowanym

czasem koherencji należałoby więc zastosować możliwie
krótki odstęp modulacji symboli OFDM, aby kanał moż-

background image

można było uważać za względnie wolnozmienny.

Ze względu na wymagania dotyczące przepływno-

ści transmisji jest konieczne użycie 16 podnośnych.
Przyjęto ponadto, że przedrostek cykliczny będzie od-
powiadał ¼ czasu trwania użytecznej części symbolu
OFDM czyli 4 minimalne odstępy próbkowania. Ponie-
waż odstęp ten wynosi 40

µs, to całkowity czas trwania

symbolu OFDM jest równy 800

µs. Porównując tę war-

tość z obliczonym wcześniej czasem koherencji kanału
widzimy, że kanał jest niestety szybkozmienny.

Biorąc pod uwagę charakterystykę przyjętego mo-

delu kanału zadecydowano o konieczności użycia pew-
nych podnośnych jako pilotów w celu estymacji kanału
w odbiorniku. Ostatecznie założono, że 13 nośnych
będzie używanych do transmisji danych z wykorzysta-
niem modulacji 16-QAM, dwie podnośne będą transmi-
tować sygnały pilotów, natomiast podnośna zerowa
będzie wyłączona ze względu na działanie przetworni-
ków C/A i A/C.

3. OPIS SYSTEMU SYMULACYJNEGO

Schemat blokowy zaprojektowanego systemu

transmisyjnego z modulacją OFDM pokazano na rys. 1.
W celu zapewnienia wysokiej jakości transmisji, w roz-
patrywanym systemie zastosowano kodowanie kanało-
we. Jest to uzasadnione faktem, że kanał telekomunika-
cyjny, w którym odbywa się transmisja jest kanałem
radiowym z zanikami. Zakłócenia występujące w kanale
radiowym wprowadzają często błędy o charakterze
paczkowym, co oznacza, że błędy błędy te są skorelo-
wane. Ponieważ zdolności korekcyjne kodów kanało-
wych są ograniczone do pewnej maksymalnej liczby
błędów w każdym bloku, pojawienie się zbyt długiej
paczki błędów może spowodować niemożliwość ich
skorygowania. Aby tego uniknąć, w systemach transmi-
sji cyfrowej stosuje się przeplot, który przez rozprosze-
nie paczek błędów w czasie powoduje nadanie im bar-
dziej przypadkowego charakteru. Z tego też względu za
układem kodera kanałowego zastosowano układ
przeplotu bitowego.

Kodowanie kanałowe realizowane jest za pomocą

kodu splotowego. Wybór parametrów kodera splotowe-
go uzależniony jest przede wszystkim od wymaganej
szybkości transmisji oraz dopuszczalnej złożoności
dekodera. Podstawowymi parametrami podlegającymi
ustaleniu jest sprawność kodu R oraz długość pamięci ν.
W badaniach symulacyjnych zastosowano koder o sto-
sunkowo niskiej sprawności R = 1/2, długości pamięci
ν = 6 i wielomianach generujących g

0

= 171

8

, g

1

= 133

8

.

Kod ten jest powszechnie uznawany za standard w tele-

komunikacji [3] i znalazł zastosowanie między innymi w
naziemnym i satelitarnym segmencie systemu telewizji
cyfrowej DVB oraz w standardzie bezprzewodowej sieci
lokalnej Hiperlan/2. Układ kodera tego kodu splotowego
pokazano na rys. 2. Liczba stanów kodera wynosi
2

ν

= 64.

T

T

T

T

T

T

wyj cie 0

wyjcie 1

wejcie

Rys. 2. Koder splotowy o parametrach: R = 1/2, v = 6,

g

0

= 171

8

, g

1

= 133

8

Z uwagi na prostotę realizacyjną i przesyłanie da-

nych w postaci bloków o skończonej długości zastoso-
wano układ przeplotu blokowego. Układ taki można
przedstawić w postaci macierzy elementów pamięci
zawierającej I wierszy oraz J kolumn. W rozpatrywanym
systemie dane zapisywane są do macierzy „wiersz po
wierszu”, a odczytywane „kolumna po kolumnie”.
Optymalny dobór wielkości i organizacji macierzy prze-
plotu stanowił przedmiot badań symulacyjnych i ich
wyniki zostały opisane w punkcie czwartym. W każdym
przypadku długość wiersza macierzy przeplotu stanowi-
ła wielokrotność 13 bitów odpowiadających w przypad-
ku modulacji 16-QAM ¼ symbolu OFDM. Liczba wier-
szy wynikała z kolei z długości wiersza i całkowitej
długości bloku kodowego.

Bity z wyjścia układu przeplotu podawane są na

układ modulatora 16-QAM, a następnie symbole wyj-
ściowe danych z modulatora są multipleksowane z sym-
bolami pilotów. Zestaw 13 symboli danych i 2 pilotów
jest podawany na wejście bloku szybkiej odwrotnej
transformacji Fouriera (IFFT) w celu utworzenia symbo-
lu OFDM. Ze względu na konieczność interpolacji pró-
bek symbolu OFDM zdecydowano się na wykonanie tej
operacji za pomocą IFFT o niestandardowej strukturze i
rozmiarze (N = 96). Taki zabieg eliminuje konieczność
zaprojektowania dodatkowego filtru interpolującego i
dzięki temu uniknięto dodatkowych zniekształceń pró-
bek transmitowanych. Po interpolacji i utworzeniu sym-
bolu OFDM następuje dodanie przedrostka cyklicznego
(CP) i tak przygotowane dane są transmitowane poprzez
kanał radiowy.

Na wejściu odbiornika znajduje się układ estymacji

wstawianie

16/96

IFFT

96/16

usuwanie

usuwanie

pilotów

preambuł

FFT

preambuł

korektor

estymator

kanału

Viterbiego

dekoder

bitowy

rozplot

obliczanie

metryk

koder

splotowy

przeplot

bitowy

modulator

16−QAM

wstawianie

pilotów

+CP

−CP

AWGN

offsetu

korekcja

z zanikami

kanał

Rys. 1. Schemat blokowy opracowanego systemu symulacyjnego

background image

przesunięcia częstotliwości (ang. frequency offset
estimation
) i jego korekcji. W układzie tym realizowana
jest także synchronizacja w dziedzinie czasu, jednakże
wyniki przedstawione w punkcie czwartym zostały wy-
generowane przy założeniu idealnej synchronizacji cza-
sowej. Następnie usuwany jest przedrostek cykliczny a
próbki symbolu OFDM są decymowane i demodulowa-
ne za pomocą algorytmu FFT. W kolejnym kroku, układ
estymatora kanału wyznacza aktualną wartość transmi-
tancji kanału dla każdej z używanych nośnych. Wyko-
rzystuje on w tym celu interpolację liniową estymat
transmitancji kanału obliczonych przy pomocy znanych
symboli pilotów. Informacja z estymatora kanału służy
do korekcji odebranych danych w układzie korektora,
zgodnie z algorytmem ZF (ang. Zero Forcing).

Skorygowane przez korektor odebrane symbole

16-QAM podlegają następnie przetwarzaniu w celu
uzyskania na ich podstawie miękkodecyzyjnych metryk
wejściowych dla dekodera Viterbiego kodu splotowego.
W celu uzyskania miękkich metryk dla poszczególnych
bitów zawartych w transmitowanych symbolach
16-QAM, stosowana jest suboptymalna metoda przed-
stawiona w pracy [4]. Rozpatrując i-ty symbol odebrany
r

i

oraz j-ty bit w tym symbolu, przyrosty metryk

λ

k

i,j

(k = 0, 1) dla każdego bitu b

i,j

= k obliczane są

następująco:

2

,

2

,

,

|

|

)

|

(

ln

k

j

i

i

i

j

i

i

k

j

i

x

y

h

k

b

r

p

=

=

λ

(3)

gdzie y

i

jest skorygowanym przez korektor sygnałem

odebranym a x

0

i,j

i x

1

i,j

są sygnałami z konstelacji

modulacji 16-QAM, które posiadają na j-tej pozycji
bloku odwzorowanego w sygnał z tej konstelacji,
odpowiednio, bit “0” i bit “1” i są jednocześnie
najbliższe sygnałowi y

i

w sensie kwadratu odległości

Euklidesa. Wielkość h

i

jest próbką odpowiedzi

impulsowej kanału, dostarczaną przez estymator kanału.

Wyznaczone na podstawie równania (3) przyrosty

metryk

λ

k

i,j

podlegają następnie operacji rozplotu w celu

przywrócenia uporządkowania zgodnego z wyjściem
kodera splotowego w nadajniku. Miękkodecyzyjny
dekoder Viterbiego w oparciu o przyrosty metryk

λ

k

i,j

dokonuje dekodowania sekwencji kodowej i dostarcza
na swoim wyjściu estymaty nadanych bitów bloku
danych.

4. REZULTATY SYMULACYJNE

W celu weryfikacji przyjętych założeń realizacyj-

nych modemu, wykonanych zostało szereg badań symu-
lacyjnych. W symulacjach uwzględniono zarówno przy-
padek kanału z jedną jak i z dwoma ścieżkami. W każ-
dym przypadku w symulacjach przesyłano 10 mln. bitów
co zapewnia wystarczającą wiarygodność rezultatów na
poziomie bitowej stopy błędów (BER) wynoszącej 10

-5

.

Na rys. 3 i 4 zostały przedstawione rezultaty badań nad
relacją pomiędzy liczbą kolumn i wierszy w układzie
przeplotu bitowego. W badaniach tych przyjęto, że dłu-
gość bloku danych wynosi 52 000 bitów co odpowiada
transmisji 1000 symboli OFDM.

Na rys. 3 pokazano rezultaty dla jednodrogowego

modelu kanału z zanikami Rayleigha i maksymalnej
częstotliwości Dopplera f

d

= 10 Hz. Jak widać, najlepsze

rezultaty przynosi zastosowanie układu przeplotu posia-
dającego macierz o 26 kolumnach i 2000 wierszy
(26 x 2000 bitów). Nieco gorszą jakość transmisji (od-
powiednio ok. 1,5 dB i 2,5 dB straty przy BER = 10

-5

)

uzyskuje się dla przeplotów o macierzy 52 x 1000 bitów
oraz 13 x 4000 bitów. Zwiększanie długości wiersza
macierzy przeplotu powyżej 52 bitów przynosi
natomiast już bardzo wyraźną, kilkudecybelową stratę.

Eb/N0

BER

10 dB

15 dB

20 dB

25 dB

2.E-05

3.E-05

5.E-05

7.E-05

1.E-04

2.E-04

3.E-04

5.E-04

7.E-04

1.E-03

2.E-03

3.E-03

5.E-03

7.E-03

2.E-02

3.E-02

5.E-02

7.E-02

1.E-01

5.E-01

26 x 2000

52 x 1000

104 x 500

260 x 200

520 x 100

13 x 4000

Rys. 3. Wpływ wymiarów macierzy przeplotu na bitową

stopę błędów w kanale jednodrogowym (fd = 10 Hz)

Eb/N0

BER

10 dB

15 dB

20 dB

2.E-05

3.E-05

5.E-05

7.E-05

1.E-04

2.E-04

3.E-04

5.E-04

7.E-04

1.E-03

2.E-03

3.E-03

5.E-03

7.E-03

2.E-02

3.E-02

5.E-02

7.E-02

1.E-01

5.E-01

52 x 1000; fd=2,5Hz

52 x 1000; fd=40Hz

104 x 500; fd=2,5Hz

104 x 500; fd=40Hz

26 x 2000; fd=2,5Hz

26 x 2000; fd=40Hz

13 x 4000; fd=2,5Hz

13 x 4000; fd=40Hz

Rys. 4. Bitowa stopa błędów w zależności od wymiarów

macierzy przeplotu i maksymalnej częstotliwości

Dopplera w kanale jednodrogowym

Z kolei na rys. 4 przedstawiono wyniki symulacji w

tym samym kanale dla dwóch innych maksymalnych
częstotliwości Dopplera: f

d

= 2,5 Hz i f

d

= 40 Hz. W

przypadku wyższej częstotliwości Dopplera obserwuje-
my lepszą jakość transmisji a wpływ relacji pomiędzy
liczbą kolumn i wierszy macierzy przeplotu jest znacznie
mniejszy – krzywe stopy błędów dla systemów o prze-
plotach 26 x 2000 bitów, 52 x 1000 bitów i 104 x 500
bitów dla BER = 10

-5

praktycznie się pokrywają. Dla

sytuacji f

d

= 2,5 Hz najlepsze wyniki uzyskuje się dla

przeplotu o macierzy 13 x 4000 bitów i pogorszenie
jakości transmisji pogłębia się wyraźnie wraz ze wzro-
stem liczby kolumn w macierzy przeplotu.

Na podstawie przytoczonych wyników można

background image

stwierdzić, że przeplot o organizacji 52 x 1000 lub
26 x 2000 bitów wydaje się być najlepszym rozwiąza-
niem dla rozważanego systemu.

Porównanie jakości transmisji systemu z kodowa-

niem i bez kodowania w kanale jednodrogowym dla
różnych wartości maksymalnej częstotliwości Dopplera
f

d

ilustruje rys. 5. W przypadku systemów z kodowaniem

długość bloku danych wynosi 5200 bitów. Na rysunku
tym umieszczono także krzywą dla systemu niekodowa-
nego działającego w kanale AWGN. Zysk wynikający z
kodowania uzależniony jest od częstotliwości Dopplera i
dla BER

= 10

-

3 wynosi od kilku decybeli (dla

f

d

= 2,5 Hz) do ponad 15 dB (dla f

d

= 40 Hz). Dla

niższych stóp błędów zysk ten jest nawet większy. Z
rys. 5 widać także, że krzywe stopy błędów dla systemu
niekodowanego leżą blisko siebie natomiast jakość
transmisji systemu z kodowaniem, z uwagi na
zastosowany przeplot, silnie zależy od częstotliwości
Dopplera.

Rys. 6. Wpływ długości bloku kodowego na jakość

transmisji w kanale dwudrogowym w zależności od

maksymalnej wartości częstotliwości Dopplera

Eb/N0

BER

10 dB

15 dB

20 dB

25 dB

2.E-05

3.E-05

5.E-05

7.E-05

1.E-04

2.E-04

3.E-04

5.E-04

7.E-04

1.E-03

2.E-03

3.E-03

5.E-03

7.E-03

2.E-02

3.E-02

5.E-02

7.E-02

1.E-01

5.E-01

fd=2,5Hz; 5,2kb

fd=10Hz; 5,2kb

fd=40Hz; 5,2kb

fd=2,5Hz; 26kb

fd=10Hz; 26kb

fd=40Hz; 26kb

fd=2,5Hz; 52kb

fd=10Hz; 52kb

fd=40Hz; 52kb

Eb/N0

BER

5 dB

10 dB

15 dB

20 dB

25 dB

2.E-05

3.E-05

5.E-05

7.E-05

1.E-04

2.E-04

3.E-04

5.E-04

7.E-04

1.E-03

2.E-03

3.E-03

5.E-03

7.E-03

2.E-02

3.E-02

5.E-02

7.E-02

1.E-01

5.E-01

niekod.; AWGN

niekod.; fd=2,5Hz

niekod.; fd=10Hz

niekod.; fd=40Hz

kod.; fd=2,5Hz

kod.; fd=40Hz

kod.; fd=10Hz

Rys. 5. Porównanie jakości transmisji dla systemów z

kodowaniem i bez kodowania w kanale jednodrogowym

Na rys. 6 przedstawiono porównanie jakości

transmisji systemów z kodowaniem w kanale
dwudrogowym dla różnych wartości częstotliwości
Dopplera oraz długości bloku danych (odpowiadających
wielkości przeplotu bitowego). Analiza wyników
pozwala stwierdzić, że im niższa jest częstotliwość
Dopplera f

d

tym silniej jakość transmisji zależy od

długości bloku danych. W szczególności, w najlepszym
przypadku (dla f

d

= 40 Hz) zmniejszenie długości bloku

danych o połowę (z 52 000 bitów na 26 000 bitów) nie
przyniosło praktycznie żadnej straty w jakości
transmisji.

5. PODSUMOWANIE

W artykule opisano projekt wielotonowego

modemu radiowego dla transmisji w pasmie o szerokości
25 kHz. Przedstawiono założenia dotyczące modemu i
dobór jego parametrów. Następnie zaprezentowano
opracowany system symulacyjny oraz wyniki badań
przeprowadzonych z jego wykorzystaniem. Na
podstawie otrzymanych wyników symulacji można
stwierdzić, że zaproponowane kodowanie protekcyjne

wraz z przeplotem bitowym przynoszą istotną poprawę
jakości transmisji w porównaniu do systemu
niekodowanego. Zysk kodowania silnie zależy od
wartości maksymalnej częstotliwości Dopplera i może
przekraczać nawet 15 dB przy stopie błędów na
poziomie 10

-3

. Wyniki symulacji pozwoliły tez ustalić

pragmatyczne

wymiary

układu

przeplotu

przy

określonym założeniu dotyczącym długości bloku
danych. Następnym etapem prac nad modem będzie
opracowanie oprogramowania poszczególnych bloków
systemu, przeznaczonych do implementacji za pomocą
procesora sygnałowego.

Autorzy dziękują osobom współuczestniczącym w

projekcie badawczym KBN, w szczególności prof. K.
Wesołowskiemu, dr. H. Boguckiej, dr. Z. Długaszew-
skiemu oraz mgr. A. Langowskiemu za udział w ustale-
niu założeń transmisyjnych, wyznaczeniu wielu parame-
trów modemu i zaproponowanie konkretnych rozwiązań.

LITERATURA

[1] Y. R. Zheng and C. Xiao, “Simulation Models

With Correct Statistical Properties for Rayleigh
Fading Channels,” IEEE Trans. Commun., Vol. 51,
No. 6, 2003, 920-928

[2] T. S. Rappaport, Wireless Communications

Principles & Practice, Prentice Hall, Upper Saddle
River, New Jersey, 1996

[3] S. G. Wilson, Digital Modulation and Coding,

Prentice Hall, Upper Saddle River, New Jersey,
1996

[4] T. May, H. Rohling and V. Engels, “Performance

analysis of Viterbi decoding for 64-DAPSK and
64-QAM modulated OFDM signals,” IEEE Trans.
Commun.
, Vol. 46, No. 2, 1998, 182-190


Wyszukiwarka

Podobne podstrony:
Andrzej STRUŻYNA Sztuka gier wideo Wybrane aspekty badań nad estetycznością i artyzmem gier kompute
Kowal J Wybrane aspekty badań jakościowych w systemach informacyjnych
EKONOMICZNE I SPOŁECZNE ASPEKTY BADAŃ PRZESIEWOWYCH WYBRANYCH NOWOTWORÓW W POLSCE”
KKRRiT2002 K Wesołowki [Transmisja wielotonowa i jej zastosowania w systemach radiowyc]
Cw 3 patologie wybrane aspekty
Wybrane aspekty gospodarki finansowej NBP 17072009 2
bóle krzyża (FLAGGERMUS), Wybrane aspekty diagnostyki lędźwiowo
rozrod wyk 2011 01 19, Wybrane aspekty rozrodu małych przeżuwaczy
Gazprom jako narzędzie polityki energetycznej Rosji w Europie wybrane aspekty
17 Wybrane aspekty metabolizmu ostatnie osiągnięcia naukiid 17370 ppt
WYBRANE ASPEKTY BEZPIECZEŃSTWA DANYCH BANKOWYCH
Charakterystyka wybranych aspektów percepcji słuchowej u dzieci dyslektycznych
Filiciak EURO w internecie wybrane aspekty
Wybrane aspekty starzenia sie
Aktywa przedsiębiorstwa Wybrane aspekty

więcej podobnych podstron