889


E1. Narysuj schemat blokowy ilustrujący współczesne implementacje (realizacje fizyczne) systemów cyfrowego przetwarzania sygnałów (CPS)

Konwersja analogowo-cyfrowa:

0x01 graphic
0x01 graphic
0x01 graphic

(zegar) (okres próbkowania)

„Sample and hold” 0x01 graphic

gdzie: 0x01 graphic

Koncepcyjna reprezentacja systemu konwersji analogowo cyfrowej:

0x01 graphic

konwerter właściwy (digitizer)

E2. Omów poznane klasyfikacje sygnałów. Podaj przykłady.

Sygnały:

- deterministyczne(czas ściśle określony w dowolnej chwili)

- stochastyczne (nieokreślony i można określić jedynie za pomocą parametrów statystycznych)

Podstawowe kryteria klasyfikacji sygnałów wiążą się z ich cechami czasowymi. Dotyczą zakresu, ograniczoności i ciągłości zbiorów wartości argumentu i zbiorów wartości funkcji.

Sygnały:(rysunki po kolei jak nazwy)

- ciągłe z czasem ciągłym (analogowe),

- ciągłe z czasem dyskretnym (dyskretne),

- dyskretne z czasem ciągłym (skwantowane),

- dyskretne z czasem dyskretnym (cyfrowe).

0x01 graphic

E3.Jakie są zalety i wady CPS w porównaniu z techniką analogową?

ZALETY:

-gwarantowana dokładność; zależy ona tylko od liczby bitów (liczba przyjmująca dwie wartości zero i jeden)

- perfekcyjna powtarzalność; nie ma zmian wywołanych tolerancja elementów (np. nagranie cyfrowe można kopiować wielokrotnie bez jakiegokolwiek pogorszenia jakości sygnału)

- nie występuje dryft temperatury, ani skutki starzenia

- postęp w technologii półprzewodników pozwala na uzyskanie zwiększonej niezawodności, mniejszych wymiarów, niższych kosztów, niższego poboru mocy i większej szybkości

- większa elastyczność, systemy cps można oprogramować i przeprogramowywać, by uzyskać rózne funkcje bez modyfikacji sprzętu

- lepsze charakterystyki, można realizować funkcje nieosiągalne w technice analogowej(np. liniowe charakterystyki fazowe, zespolone algorytmy filtracji)

- w pewnych przypadkach informacja może być dostępna tylko w postaci cyfrowej i jedyną możliwością jest jej przetwarzanie cyfrowe

Wady cyfrowego przetwarzania sygnałów (CPS)

układy cyfrowe nie są jeszcze znane),

częstotliwości (jednak zakres częstotliwości ciągle rośnie),

przetworników AC i CA, które mogą być skomplikowane i drogie,

silnych (np. sterujących głośnikiem) jest bardzo trudne;

wymagane jest stosowanie odpowiednio przedwzmacniaczy i

układów mocy w takich przypadkach,

wymaga układów CPS o znacznie szerszym paśmie w porównaniu

do układów analogowych.

E4.Objaśnij różnice w przetwarzaniu sygnału analogowego filtrem analogowym, a systemem z filtrem cyfrowym.

(coś tego mało?? )

cyfrowy: istnieje koherentna granica pulsacji widma

analogowy: granica pulsacji jest rozmyta (dąży do zera)

E5. Objaśnij zjawisko aliasingu na przykładzie próbkowania sinusoidy.

Dla częstotliwości próbkowania mniejszej od dwukrotnej największej częstotliwości (pulsacji) sygnału (0x01 graphic
) nie jesteśmy w stanie na wyjściu uzyskać takiego sygnału podanego na wejście. Dzieje się tak wskutek nakładania się na siebie widma sygnału próbkowanego (aliasing). Na przykładzie sinusoidy mamy:

0x01 graphic

Transformata Fouriera 0x01 graphic

0x01 graphic

jego widmo:

0x01 graphic
później znajduje się wykres widma sygnału spróbkowanego 0x01 graphic
z częstotliwością próbkowania (pulsacja próbkowania) 0x01 graphic

0x01 graphic

a więc sygnału uzyskamy na wyjściu (zrekonstruowany):

0x01 graphic

bez aliasingu:

0x01 graphic

sygnał na wyjściu (w pełni odtworzony):

0x01 graphic

Obserwujemy więc zmianę w pulsacji (częstotliwości) sygnału zrekonstruowanego w przypadku wystąpienia zjawiska aliasingu.

E6. Napisz twierdzenie Nyquista o próbkowaniu równomiernym. Co nazywamy częstotliwością Nyquista, szybkością Nyquista i przedziałem Nyquista?

Tw. Niech 0x01 graphic
będzie sygnałem o ograniczonym widmie, tzn.:

0x01 graphic
tzn.dla0x01 graphic
0x01 graphic

Wówczas 0x01 graphic
jest jednocześnie określone przez ciąg swoich próbek:

0x01 graphic

jeśli:

0x01 graphic
0x01 graphic
częstotliwość próbkowania

lub inaczej:

0x01 graphic

czyli jeśli częstotliwość próbkowania jest co najmniej dwukrotnie większe od maksymalnej częstotliwości w widmie sygnału 0x01 graphic
.

Minimalne częstotliwości próbkowania, przy której można jeszcze odtworzyć sygnał 0x01 graphic
bez zniekształceń, tzn częstotliwość 0x01 graphic
nazywa się szybkością Nyquista. Odpowiadający jej maksymalny okres próbkowania 0x01 graphic
nazywa się przedziałem Nyquista. Częstotliwość Nyquista to 0x01 graphic
.

E7. Narysuj schemat blokowy przetwornika C/D traktując przetwarzanie dwustopniowo> Napisz wzory na przebiegi czasowy x(t) i jego widmo X(jΩ) na wyjsciu pierwszego bloku.

C/D przetwornik czasu ciągłego na dyskretny

0x01 graphic

0x08 graphic
-czas ciagly

0x08 graphic

0x08 graphic

Przedstawiony zostal schemat blokowy reprezentujący idealny konwerter czasu ciągłego na czas dyskretny, czyli idealny układ próbkujący. W ogólności operacja probkowania nie jest odwracalna. Aby zapewnic odtwarzalność sygnal Xc(t) musi spełniać zalozenie Fp>2FM

Reprezentacja operacji probkowania, jako procesu dwustopniowego:

- mnożenie przez ciag impulsow Diraca

- zamiana sekwencji impulsow na ciag dyskretny.

0x08 graphic
Dziedzina czasu 0x01 graphic
-dystrybucja grzebieniowa

Z właściwości próbkującej 0x01 graphic

0x01 graphic
, bo => 0x01 graphic

Dziedzina częstotliwości - znajdujemy widmo xp(t)

Oznaczamy xc(jΩ)=0x01 graphic
, Ω[rad/s] - pulsacja analogowa

0x01 graphic

0x01 graphic

Wniosek: Widmo sygnalu xp(t) jest superpozycja okresowa powtarzających się widm sygnalu ciągłego xC(t). Kolejne składniki Xp(jΩ) sa poprzesuwane o całkowita wielokrotność pulsacji próbkowania Ωp.

E8. Zilistruj przypadki konwersji C/D: bez aliasingu i z aliasingiem, poslugujac sie widmem.
Przyjmujemy, ze 0x01 graphic
dla 0x01 graphic
tzn. Najwieksza niezerowa skladowa w widmie0x01 graphic
jest skladowa pulsacji 0x01 graphic

0x08 graphic

0x08 graphic
W wyniku konwersji C/D widmo 0x01 graphic
zostaje pomnozone przez dystrybucje grzebieniowa δ(t).

0x08 graphic
0x08 graphic
0x01 graphic

0x08 graphic

0x08 graphic
0x08 graphic

Otrzymujemy 0x01 graphic

  1. bez aliasingu

0x08 graphic

  1. z aliasingiem (widma składowe nie są rozłączne (nakładają się) -> zniekształcenie aliasowe)

0x08 graphic

T

E9. Objaśnij w jaki sposób można odtworzyc sygnal ciągły xC(t) na podstawie sygnalu xp(t) o reprezentacji w postaci delty Diraca?

Można odtworzyc xC(t) na podstawie xp(t) za pomocą idealnego filtru dolnoprzepustowego, takiego, ze:

0x08 graphic
Charakterystyka czestotliwosciowa filtru dolonoprzepustowego

0x08 graphic
<- filtr rekonstrukcyjny

0x08 graphic
Filtr rekonstrukcyjny

Wówczas w układzie:

0x08 graphic

Widmo sygnalu wyjściowego ma ta sama postac, jak widmo sygnalu xC(t)

0x01 graphic

Musi być spełnione twierdzenie o probkowaniu.

E10. Do czego sluzy formula interpolacyjna Shannona i jak wyglada? Podaj przykłady interpolacji graficznej.

Formula interpolacji Shannona stosowana jest do ostworzenia xC(t) na podstawie x [n].

Przyjmujemy, ze Ωgp/2=Π/T

0x08 graphic
Wtedy: Hd(jΩ)=T dla |Ω|<Π/T; 0 dla|Ω|>Π/T

Odp. Impulsowa z def. F-1

0x01 graphic

Wówczas:

0x01 graphic

FORMULA INTERPOLACYJNA SHANNONA

Ponieważ hd(0)=1 oraz hd(nT)=0 dla n=+-1,+-2,.. jeżeli x[n]= xc(nT) wówczas: X0(nT)= XC(nT) dla wszystkich całkowitych n, niezależnie o T (probki x0 i xC sa takie same) Caly przebieg xC jest skonstruowany na bazie funkcji sinc()

0x01 graphic
gdzie 0x01 graphic
; sin(0)=1

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic

0x01 graphic
0x01 graphic

Idealny system (X0(t)=XC(t)) odtwarzający sygn. Ciągły.

0x01 graphic

E11. Na czym polega idea przetworniak D/C (schemat)? Co oznacza skrót?

0x08 graphic

0x08 graphic

D/C - przetwornik czasu dyskretnego na ciagly

E12. Narysuj schematy blokowe:

a) idealnego dyskretnego-czasowego systemu przetwarzania sygnałów analogowych

b) systemu cyfrowego przetwarzania sygnałów analogowych.

0x08 graphic

0x08 graphic
Ad a)

0x08 graphic

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
Ad b)

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic

AAF -anti aliasin filter, SiH - sample and hold, A/D - analog to digital converter

D/A - digital to analog converter RF - reconstruction filter

E13. Zdefiniuj i narysuj charakterystykę częstotliwościową analogowego filtru antyaliasingowego (AAF).

AAF - stosuje się aby wymusić warunek, żeby częstotliwość próbkowania była co najmniej 2-krotnie większa od maksymalnej częstotliwości w widmie sygnału analogowego.

Charakterystyka :
0x01 graphic

Filtr praktyczny ma pasmo przejściowe od Ωpas i pasmo zaporowe od Ωstop, przy czym :

Ωpas <Ωstop <= Ωo/2

Ωp - pulsacja próbkowania

Ωpas - pulsacja graniczna pasma przepustowego AAF tj. najwieksza pulsacja jaka ma być wiarygodnie zachowana w sygnale sprobkowanym.

Ωstop - pulsacja graniczna pasma zaporowego

E14. Objaśnij podstawowe cechy analogowego filtru Butterwortha jako filtru AAF. Na jakiej podstawie określa się rząd tego filtru (liczbę biegunów transmitancji) ?

* Antyaliasingowy filtr dolnoprzepustowy Butterwortha jest filtrem maksymalnie płaskim, o charakterystyce amplitudowej malejącej monotonicznie w paśmie zaporowym.

* Rząd filtru można określić na podstawie charakterystyki amplitudowej podniesionej do kwadratu :

0x01 graphic
0x01 graphic
, gdzie N - rząd filtru

* Maksymalnie plaska charakterystyka amplitudowa dla Ω=0 wynika stad iż (2N-1) pochodnych z funkcji 0x01 graphic
w tymże punkcie (Ω=0) jest równych zeru.

* Wzmocnienie opisywanego filtru wyraża się wzorem :

G(Ω)= 10 log 0x01 graphic
= 20 log 0x01 graphic
, z którego łatwo odczytać ze dla :

- Ω=0 => G(0) = 0

- Ω= Ωc => G(Ωc) = 10 log(1/2) 0x01 graphic
-3dB

(dlatego właśnie pulsacje odcięcia nazywamy 3-decybelową)

* Róznica poziomów tłumienia (w [dB]) dla Ωpas i Ωo okreslona jest przez :

10 log 0x01 graphic
0x01 graphic
0x01 graphic
, gdzie :

Ωo = Ωp - Ωpas

Ωp - minimalna częstotliwość próbkowania

Ωpas - pulsacja graniczna pasma przepustowego

Ωc- pulsacja odciecia

E15. W jakim celu stosuje się nadprobkowanie sygnału analogowego? Na czym ono polega ? Objaśnij posługując się rysunkiem widma sygnalu.

Nadprobkowanie sygnalu analogowego stosuje się w celu zmniejszenia rzedu filtru (filtry analogowe wysokiego rzedu sa bardzo drogie). Idea polega na tym, ze o rzedzie filtru decyduje szerokość pasma przejściowego. Rzad N jest tym mniejszy im szersze pasmo przejściowe, a wiec im wieksza częstotliwość próbkowania Ωp.

0x01 graphic

Przy dwukrotnym nadprobkowaniu sygnalu, pasmo przejściowe filtru jest ponad 3-krotnie szersze. Te wymagania spelnia filtr o mniejszym rzedzie.

E16. Narysuj konfiguracje fizyczna konwersji analogowo-cyfrowej z układem „samlple & hold” i konwerterem „walsciwym” A/O. Objaśnij dzialanie układu ZOH i zilustruj odpowiednim rysunkiem.

0x01 graphic

Na wstepie układu S & H sygnal zostaje przemnożony z sygnalem „grzebieniowym” - ciagiem impulsow Dirac'a. Układ ZOH pobiera w takt wystepowania impulsow probki i podtrzymuje je na swym wyjsciu do nastepnego impulsu. 0x01 graphic

E17. Koncepcyjna reprezentacja systemu konwencji analogowo-cyfrowej z próbkowaniem i

kwantowaniem i kodowaniem binarnym.

0x01 graphic

Kod offsetowy binarny - symbole sa przypisywane w kolejności numerycznej rozpoczynając od najmniej znaczącego poziomu kwantowania.

Kod U2 - najbardziej na lewo położony bit (MSB) jest bitem znaku a pozostale bity reprezentuja liczby całkowite lub ulamki.

Liczba poziomow kwantowania : 0x01 graphic

Krok kwantowania : 0x01 graphic

B - liczba bitow w „B+1 bitowym” slowie, rozdzielczość konwertera.

E18. Kwantyzator

system nieliniowy transformujący próbkę wejsciowa x[n] w próbkę x^[n] nalezaca do skończonego zbioru scisle określonych wartości. Operacje otrzymywania skwantyzowanej probki zapisujemy : x^[n]=Q(x[n]). Operacja kwantyzacji nie jest odwracalna.

E19. Założenia statystycznego modelu błędu kwantowania.

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

x[n] x[n]=x[n]+e[n]

e[n]

Założenia statystycznego modelu błędu kwantowania:

  1. Ciąg błędów e[n] jest dyskretnym, stacjonarnym procesem stochastycznym.

  2. Ciąg błędów jest nieskorelowany z ciągiem x[n].

  3. Zmienne losowe procesu błędu są nieskorelowane tzn. ciąg błędu jest białym szumem.

  4. Rozkład prawdopodobieństwa procesu błędu jest równomierny w zakresie błędu kwantowania.

0x08 graphic

P[e]

0x08 graphic
0x08 graphic
0x08 graphic
1/Δ

0x08 graphic
-Δ/2 Δ/2

Czyli: e[n] to biały szum o równomiernym rozkładzie prawdopodobieństwa, wartość średnia to 0, wariancja (moc szumu):

0x01 graphic

E20. Zdefiniuj współczynnik SNR.

SNR - Signal to Noise Ratio - stosunek mocy sygnału do mocy szumu - miara degradacji sygnału przez szum addytywny.

Dla (B+1) - bitowego kwantyzatora bipolarnego:

0x01 graphic
[dB]

σe2 - moc szumu

σx2 - moc sygnału do kwantyzacji

Xm - zakres kwantowania

σx - wartość średniokwadratowa amplitudy sygnału (RMS)

Model szumowy praktycznie działa dla 8 lub więcej bitów przy pełnym wysterowaniu. Obliczenia SNR są prawdziwe jedynie wtedy, gdy starannie dopasujemy zakres sygnału wejściowego do pełnego zakresu konwertera A/D

E21. W pewnym systemie CPS a) sinosoidalny sygnał o amplitudzie Xm; b) sygnał o gaussowskim rozkładzie amplitudy próbkuje się z częstotliwością Fp=20kHz i kwantuje za pomocą kwantyzatora bipolarnego w pełnym zakresie 2Xm=2V. Określ liczbę bitów kwantyzatora taką, aby średniokwadratowy błąd kwantyzacji σe<1μV. Oblicz przepływność bitową i zakres dynamiczny kwantyzatora.

Wariancja błędu (szumu) kwantyzacji: 0x01 graphic
, gdzie krok kwantyzacji 0x01 graphic
. Stąd średniokwadratowy błąd kwantyzacji: 0x01 graphic
. Z tego wzoru B wynosi:

0x01 graphic

Opowiada to liczbie 2B+1=220=1048576 - poziomów kwantyzacji.

Dla takiego B:

0x01 graphic

Przepływność bitowa wynosi zatem:

0x01 graphic

Zakres dynamiczny (SNR) kwantyzacji sterowanego sygnałem gaussowskim wynosi około:

0x01 graphic
a wysterowanego w pełni sinusoidą:

0x01 graphic
.

E22. Jaką minimalną częstotliwość próbkowania Fp należy zastosować do sygnału x(t)=cos(2*pi*Fo*t) by można go było odtworzyć jednoznacznie na podstawie ciągu próbek x[n] pobranych równomiernie co T sekund. Jaki zachodzi związek między Fp a T.

Fp > 2*Fo <= wynika z tw.Nyquista bo sygnal jest sinusoidalny czyli w widmie posiada deltę Diraca, więc wzór jest bez znaku równości

T = 1/Fp <= wiadomo

E23. Jak definiujemy sygnał?

Def. encyklopedyczna:

Przebieg fizyczny (elektryczny, radioelektryczny lub optyczny) słuzący do przenoszenia informacji między dwoma odbiornikami.

Klasyfikacja sygnałów:

1.

2.

0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

E24. Ogólny podział systemów przetwarzania sygnału.

Klasyfikacja systemów dyskretnych:

jeśli x1[n]=x2[n] dla 0x01 graphic

to y1[n]=y2[n] dla 0x01 graphic

Przykłady:

  1. idealny system opoźniający0x01 graphic
    , gdzie n0 to liczba całkowita, nazywana opóźnieniem systemu. System bezpamięciowy tylko dla n0=0. Liniowy, stacjonarny, stabilny, przyczynowy jeśli 0x01 graphic

dla n0=3

0x08 graphic
0x08 graphic

0x08 graphic

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic

0 1 2 3 0 1 2 3 4 5 6

  1. system wzmacniający średnią bieżącą:

0x01 graphic

n-ta próbka ciągu wyjściowego jest obliczona jako średnia (M1+M2+1) próbek ciągu wejściowego otaczających próbkę n-tą. Liniowy, stacjonarny, stabilny. Przyczynowy gdy M1 0 i M2 0.

Bez pamięci tylko gdy: M1=M2=0

  1. akumulator:

0x01 graphic

Liniowy, przyczynowy, z pamięcią. Nie jest stabilny.

E25. Jakie warunki należy spełnić by reprezentację sygnału ciągłą i dyskretną można było uważać za równoważne?

Tw. Nyquista o próbkowaniu równomiernym.

Niech xl(t) będzie sygnałem ograniczonym o widmie tzn xl(jΩ)=0 dla |Ω| > ΩM [rad/s]

Wówczas xl(t) jest jednoznacznie określone przez ciąg swoich próbek x[n]=xl(nT) n=0,±1,..

jeżeli ΩP= (2π)/T = 2πFP > 2ΩM lub inaczej (częstotliwość w MHz) FP > 2 FM, czyli jeśli częstotliwość próbkowania jest co najmniej 2-krotnie większa od max częstotliwości w widmie sygnału xl(t).

Tak więc minimalna częst. próbkowania, przy której można jeszcze odtworzyć sygnał x(t) bez zniekształceń, czyli częst. FP=2FM nazywa się szybkością Nyquista. Odpowiadający jej max okres próbkowania T= 1/FP = 1/2FM zazywa się przedziałem Nyquista.

E26. Zdefiniuj podstawowe operacje wykonywane na sygnałach dyskretnych (ciągach liczowych)

-- System opóźniający y[n]=T{x[n]}=x[n-n0] n0- dodatkowa całkowita [opóźnienie systemu]

-- Średnia bieżąca 0x01 graphic
, n-ta próbka jest obliczana jako średnia (M1+M2+1) próbek ciągu wejściowego otaczających próbkę n-tą.

-- Akumulator 0x01 graphic

-- Różnica wprzód y[n] = y[x+1] - x[n]

-- Różnica wstecz y[n] = x[n] - x[n-1]

-- Kompensator y[n] = x[M·n] => wyselekcjonowanie M-tej próbki ciągu wejściowego

E27. Impuls jednostkowy czyli delta Kroneckera

0x01 graphic

Skok jednostkowy 0x01 graphic

Ciąg wykładniczy rzeczywisty x[n]=anu[n]

Ciąg wykładniczy zespolony x[n]=A*ejan

Ciąg oscylacyjny - charakteryzuje się przejściami (przejściem) przez zero. np.:

0x08 graphic
0x01 graphic

Ciąg harmoniczny x[n]=Acos(ωt+φ)

Ciąg sinusoidalny x[n]=Asin(ωt).

E29. Sygnały okresowe

Sygnał okresowy x(t) z czasem ciągłym.

x(t)0x08 graphic
0x01 graphic
Sygnał powtarza się w określonych okresach, jest zdefiniowany w każdej chwili czasowej.

Sygnał okresowy z czasem dyskretnym (dyskretny) - ciąg liczbowy {x[n]}, którego elementami są próbki x[n]sygnału ciągłego (analogowego) pobierane w chwilach tn=nTs, gdzie Ts jest okresem próbkowania.

0x08 graphic
0x01 graphic
x[n], n=0,±1,±2,.. 0x01 graphic

Ciąg x[n] dyskretny nieokresowy: spróbkowana wartość napięcia na kondensatorze w obwodzie prądu elektrycznego.

Ciąg x[n] dyskretny cos, okresowy: spróbkowana cosinusoida z generatora.

E30. Jaki zachodzi związek pomiędzy częstotliwością (lub pulsacją) analogową a dyskretną ?

0x01 graphic
gdzie: 0x01 graphic
- pulsacja dyskretna (unormowana względem T), ω - pulsacja analogowa.

E31. Podaj operatorową definicje systemu dyskretnego

Def: System dyskretny to operator który przekształca ciągi WE x[n] w ciąg WY y[n]

0x08 graphic

y[n] = T{x[n]} {x[n]} => {y[n]}

Operator T{ 0x01 graphic
} określa właściwości systemu dyskretnego, czyli sposób przekształcenia { x[n] } w { y[n] }

W ogólności n - ta próbka y[n] może być funkcja wszystkich próbek x[n]

E32. Zdefiniuj następujące systemy dyskretne: bez pamięci, liniowy, stacjonarny (niezmienny względem przesunięcia), przyczynowy, stabilny

Def. System dyskretny to operator, który przekształca ciąg wejściowy x[n] w ciąg wyjściowy y[n]:

y[n] = T { x[n] }

System bez pamięci: Jeżeli odpowiedz y[n] dla każdego n zależy tylko od wymuszenia x[n] dla tej samej wartości n, to mówimy, ze system nie posiada pamięci.

Np. - System opisany zależnością: y[n] = ( x[n] )2 dla wszystkich n [NIE MA PAM.]

- y[n] = ( x[n-1] )2 [Z PAMIĘCIĄ]

System liniowy: spełnia zasadę superpozycji; jeżeli { x1[n] } => { y1[n] } oraz

{ x2[n] } => { y2[n] } to system dyskretny jest liniowy W.I.T.W gdy:

T{ ax1[n] + bx2[n] } = aT { x1[n]} + bT { x2[n]}

System stacjonarny: spełnia następujący warunek:

Jeżeli x[n] => y[n] oraz x1[n] = x[n-n0]

n0 - całkowite to:

y1[n] = y[n-n0]

System przyczynowy: System jest przyczynowy, jeżeli dla każdej liczby całkowitej n0 próbka odpowiedzi dla n = n0 zależy tylko od próbek wymuszenia o indeksach

n ≤ n0. Krotko zapisujemy:

Jeżeli: x1[n] = x2[n] dla n ≤ n0

To: y1[n] = y2[n] dla n ≤ n0

System stabilny: System dyskretny nazywamy stabilnym W.I.T.W, gdy ograniczony ciąg wejściowy daje w odpowiedzi ciąg ograniczony, tzn.

Jeśli: 0x01 graphic
to 0x01 graphic

Gdzie 0x01 graphic
i 0x01 graphic
< 0x01 graphic

E33. Co oznacza skrót: system DSL?

DLS = DYSKRETNY, LINIOWY, STACJONARNY

DYSKRETNY SYSTEM: to operator, przekształca ciąg wejściowy x[n] w ciąg wyjściowy y[n]:

y[n] = T{x[n]} {x[n]} => {y[n]}

0x08 graphic

Operator T{ 0x01 graphic
} określa właściwości systemu dyskretnego, czyli sposób przekształcenia { x[n] } w { y[n] }

LINIOWY SYSTEM: spełnia zasadę superpozycji, tzn. jeżeli spełnia zasadę superpozycji; jeżeli { x1[n] } => { y1[n] } oraz { x2[n] } => { y2[n] } to system dyskretny jest liniowy W.I.T.W gdy:

T{ ax1[n] + bx2[n] } = aT { x1[n]} + bT { x2[n]}

STACJONARNY SYSTEM: niezależny względem przesunięcia

Jeżeli x[n] => y[n] oraz x1[n] = x[n-n0]

n0 - całkowite to:

y1[n] = y[n-n0]

(opóźnione w czasie pobudzenie => opóźniona odpowiedz)

E34. Podaj definicje idealnego systemu opóźniającego i zilustruj jego działanie na dowolne wylosowanym ciągu. Jakie są właściwości tego systemu?

Idealny system opóźniający jest to system opisany zależnością:

y[n] = T{ x[n] } = x[n-n0]

dla -0x01 graphic
<n<0x01 graphic

n0 - liczba całkowita dodatnia nazwana opóźnieniem systemu

System przesuwa ciąg wejściowy w prawo o n0 próbek. Jeżeli n0 jest liczbą całkowitą ujemna - przesuwamy w lewo oznacza to wyprzedzenie w czasie np. dla ciągu

0x01 graphic

dla opóźnienia n0 = 3 mamy

0x01 graphic

0x08 graphic
0x08 graphic

Właściowści:

E35. Zdefiniuj system wyznaczający średnią bieżącą. Jakie są jego właściwości? Podaj przykład zastosowania tego systemu.

System wyznaczający średnia bieżącą zdefiniowany jest równaniem:

0x01 graphic

n - ta próbka ciągu wejściowego jest obliczana jako średnia (M1+M2+1) próbek ciągu wejściowego otaczających próbkę n - tą.

Właściwości systemu:

E36. Zdefiniuj: akumulator, kompresor oraz system różnicy w przód I wstecz. Jakie są ich zastosowania? Podaj przykład zastosowania systemu.

AKUMULATOR

Def. 0x01 graphic

Właściwości:

- liniowy, przyczynowy, z pamięcią

Przykład odpowiedzi impulsowej (systemów DLS):

0x01 graphic

nieskończona odpowiedz impulsowa

Transmitancja h[n] = u[n]

H(z) = Z{ u[n] } = 0x01 graphic

Zastosowanie jako układ całkujący

KOMPRESOR

Def. 0x01 graphic

M - całkowita liczba dodatnia

Właściwości: ciąg liczbowy wyjściowy jest przez wyselekcjonowanie każdej M-tej próbki ciągu wejściowego. System nie jest stacjonarny, widać to z przykładu.

Niech x1[n] => y1[n] oraz

x1[n] => x[n - n0]

wówczas y1[n] = x1[Mn] = x[Mn - n0]

Dla systemu stacjonarnego musi zachodzić:

y[n - n0] = x[M(n - n0)] ≠ y1[n]

a tu zachodzi jedynie dla n=1.

RÓŻNICA W PRZÓD I WSTECZ

Def. W przód y[n] = x[n + 1] - x[n]

system nie jest przyczynowy, bo występuje x[n+1]

Wstecz y[n] = x[n] - x[n-1]

przyczynowy

Przykład odpowiedzi impulsowej (system DLS)

różnica w przód: h[n] = δ[n + 1] - δ[n]

różnica wstecz: h[n] = δ[n] - δ[n - 1]

Transmitancje:

różnica w przód: h[n] = δ[n + 1] - δ[n] => H(z) = z - 1

różnica wstecz: h[n] = δ[n] - δ[n - 1] => H(z) = 1 - z -1

Zastosowanie różnicy wstecz jako różniczkowanie.

E37. Podaj definicję odpowiedzi impulsowej systemu DLS.

Dla systemu DLS definiuje się odpowiedź impulsową h[n], dla n=±1,...

Df. h[n] to odpowiedź impulsowa systemu DLS na impuls jednostkowy (δ[n] - delta Kroneckera) przy zerowych warunkach początkowych. h[n] całkowicie charakteryzuje system DLS.

Odpowiedź y[n] systemu DLS na dowolne pobudzenie x[n] jest splotem dyskretnym ciągów: pobudzenie x[n] i odpowiedzi impulsowej h[n]. Splot dyskretny jest przemienny.

0x01 graphic

E38. Podaj warunki stabilności i przyczynowości systemu DLS, wyrażone za pośrednictwem odpowiedzi impulsowej.

Układ stabilny definiujemy jako układ, w którym każde ograniczone pobudzenie powoduje ograniczoną odpowiedź. Układy liniowe, niezmienne względem przesunięcia* są stabilne gddy:

0x01 graphic

Układ przyczynowy jest to układ, w którym zmiany na wyjściu nie poprzedzają zmian na wejściu. Oznacza to, że dla układów przyczynowych, jeśli x1[n]=x2[n], n<n0, to y1[n]=y2[n], n<n0.

Układ liniowy niezmienny względem przesunięcia jest przyczynowy gddy jego odp. Impulsowa dla n<0 ma wyrazy równe zeru. Z tego względu wygodnie nazywać ciąg zerowy dla n<0 ciągiem przyczynowym rozumiejąc, że mógłby on być odp. Impulsową układu przyczynowego.

Jako przykład stabilności i przyczynowości można rozpatrzyć układ liniowy niezmienny względem przesunięcia o odp. Impulsowej h[n]=anu[n]. Ponieważ odp. Impulsowa jest zerowa dla n<0, układ jest przyczynowy. Aby określić stabilność musimy obliczyć sumę: 0x01 graphic
0x01 graphic
Jeśli |a|<1, suma postępu geometrycznego istnieje i jest równa 0x01 graphic
. Jednak jeśli |a|>=1, szereg jest rozbieżny, a zatem układ jest stabilny tylko dla |a|<1.

*Układ niezmienny względem przesunięcia charakteryzuje się przez właściwość, że jeśli y[n] jest odp. na pobudzenie x[n] to y[n-k] jest odp. na x[n-k], gdzie k to dodatnia lub ujemna liczba całkowita.

E39. Jaką odpowiedź impulsową ma połączenie kaskadowe, a jaką połączenie równoległe systemów DLS? Narysuj te połączenia.

(na rysunkach na wejściu jest x[n] na wyjściu y[n] )

  1. kaskadowo: odp. impulsowa takiego systemu to: h[n]=h1[n]*h2[n]

0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic

  1. równolegle: odp. impulsowa takiego systemu to h[n]=h1[n]+h2[n]

0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

0x08 graphic
0x08 graphic

E40. W jakim celu stosuje się splot odpowiedzi impulsowej z sygnałem wejściowym systemu dyskretnego? Przytocz odpowiednie wzory.

Celem splotu odp. impulsowej h[n] z sygnałem wej. x[n] jest obliczenie odpowiedzi y[n] systemu.

Odpowiedź y[n] systemu DLS na dowolne pobudzenie x[n] jest splotem dyskretnym ciągów: pobudzenie x[n] i odpowiedzi impulsowej h[n]. Splot dyskretny jest przemienny.

0x01 graphic

E41. Oblicz wynik {y[n]} splotu ciągów: 0x01 graphic
za pomocą

dowolnej metody.

0x08 graphic
Najszybciej można to zrobić metodą mnożenia, ale należy uważać na kropkę (kropkę stawiamy przed wartością współczynnika dla n=0)

0x08 graphic

0x08 graphic

{y[n]}n=-1={3,8,14,8,3}

E42. Zdefiniuj dyskretno-czasowe przekształcenie Fouriera (DTFT) i przekształcenie do niego odwrotne. Objaśnij symbole występujące w tych wzorach.

Dyskretno-czasowe przekształcenie Fouriera (DTFT) 0x01 graphic
, odwrotne: 0x01 graphic
(IDTFT DTFT-1)

0x01 graphic
- widmo sygnału

x[n] - ciąg próbek

E43. Na podstawie DTFT narysuj widma :amplitudowe i fazowe sygnału 0x01 graphic
dla wybranej wartości „a”

0x01 graphic
- transformata ciagu x[n]

Dla a = 0,5 0x01 graphic

0x01 graphic

Widmo amplitudowe:

0x01 graphic

0x08 graphic
0x01 graphic
=0 0x01 graphic
0x01 graphic
=0x01 graphic
=2

0x01 graphic
=0x01 graphic
0x01 graphic
=0x01 graphic
=0,66

0x01 graphic
=0x01 graphic
0x01 graphic
=0x01 graphic

0x08 graphic
Widmo fazowe:

0x01 graphic

0x01 graphic
=0 arg=0

0x01 graphic
=0x01 graphic
arg=0

0x01 graphic
=0x01 graphic
arg=-arctg(0,5)

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic

E44. Zdefiniuj charakterystyki: amplitudowo-fazową, amplitudową i fazową systemu DLS o odpowiedzi impulsowej h[n]=a nu[n]. Narysuj te charakterystyki dla wybranej wartości a takiej, która zapewnia stabilność tego systemu.

Charakterystyka amplitudowo fazowa: 0x01 graphic
dla |d|<1

Amplitudowa: 0x01 graphic
; fazowa: 0x01 graphic

0x01 graphic

E45. Własności DTFT

Jeżeli x[n] ↔ X(e j ω) i y[n] ↔ Y (e j ω) to słuszne są następujące twierdzenia:

1. O liniowości: 0x01 graphic

2. O przesunięciu w dziedzinie czasu 0x01 graphic
; no - całkowite

3. O przesunięciu w dziedzinie częstotliwości: 0x01 graphic

4. O różniczkowaniu względem częstotliwości: 0x01 graphic

5. O energii (Parsevala): 0x01 graphic

6. O splocie: 0x01 graphic

7. O iloczynie: 0x01 graphic

E46. Wiadomo, że DTFT ciągu x[n] wynosi X(ejω). Jakie są transformaty ciągów:

  1. x1[n]=x[n-5], b) x2[n]=x[n]ej5n, c) x3[n]=anx[n], d) x4[n]=nx[n]

  1. x1[n]=x[n-5], z twierdzenia o przesunięciu w dziedzinie czasu

x[n-n0] ⇔ e-jωn0 X(ejω) , n0 - całkowite

x[n-5] ⇔ e-jω5 X(ejω) , n0 - całkowite

  1. x2[n]= x[n]ej5n, z twierdzenia o przesunięciu w dziedzinie częstotliwości

x[n] ejωn0 ⇔ X(ej(ω-ω0))

x[n] ejω5 ⇔ X(ej(ω-5))

  1. x3[n]= anu[n]

X(ejω)=0x01 graphic
, gdzie |a|<1 - w-nek istnienia X(ejω)

  1. x4[n]=nx[n], z twierdzenia o różniczkowaniu względem częstotliwości

nx[n]⇔ j0x01 graphic

Energię każdego z sygnałów można obliczyć korzystając z twierdzenia Parsevala

E=0x01 graphic

E47. Oblicz charakterystykę amplitudowo-fazową i odpowiedź impulsową systemu opisanego równaniem różnicowym y[n]=0.9y[n-1]+x[n]-5x[n-1].

0x01 graphic
; 0x01 graphic

0x01 graphic
; 0x01 graphic

0x01 graphic

0x01 graphic

charakterystyka fazowa 0x01 graphic

charakterystyka fazowa 0x01 graphic

E48. Zdefiniuj przekształcenie Z i oblicz Z-transformaty ciągów: a) x1[n] = δ[n], b) x2[n] = x[n-n0], n i n0 całkowite, c) x3[n] = u[n], d) x4[n] = anu[n]. Podaj odpowiednie obszary zbieżności.

-dwustronna 0x01 graphic
-jednostronna 0x01 graphic

  1. 0x01 graphic
    , 0x01 graphic
    ; 0x01 graphic
    , ROC - wszystkie z

  2. 0x01 graphic
    ; 0x01 graphic
    ,

ROC - gdy no>0 to takie jak dla X(z) z wyjątkiem 0, gdy no<0 to takie jak dla X(z) z wyjątkiem nieskończoności.

  1. 0x01 graphic
    ; 0x01 graphic

ROC: |q|<1 → (1/|z|<1) → |z|>1

  1. 0x01 graphic
    ; 0x01 graphic

ROC: |q|<1 → (|a| / |z|) < 1 → |z| > |a|

E49.Przytocz zapamiętane właściwości przekształcenia Z.

ax1[n]+bx2[n]↔ aX1(z)+bX2(z)

x[n-1]↔ z-1·X(z)+x[-1]

x[n-k]↔ z-1·X(z)+ n-1n x[-1]·z-k+n

x[n-k]·u[n-k] ↔ z-k·X(z)

x[n+1]↔ z·X(z) +x[0]

an·x[n]↔ X(z/a) dla | z | > | a |

n'·x[n]↔ -z ·dx(τ)/dt

x1[n]·x2[n]↔ X1(z)·X2(z)

E50. Oblicz oryginał x[n], gdy dana jest transformata:

0x01 graphic

E51.Jak definiujemy transmitancje H(z) systemu DLS?

0x08 graphic
Def. H(z)=Y(z)/X(z) przy z.w.p. gdzie:x[n] →X (z); y[n] →Y(z);

h[n] →H(z) stąd H(z)=Z{h[n]}= h=-∞ h[n] ·z -n
Warto wiedzieć że dla zadanej H(z) odpowiedź y[n] na x[n] dla systemu DLS przy z.w.p. obliczamy: y[n]=Z-1 {H(z)·X(z)}

E52. Zdefiniuj warunki przyczynowości i stabilności systemu DLS o transmitancji H(z):

Transmitancja H(z):

0x01 graphic
(po faktoryzacji)

gdzie: ci -zera transmitancji ; di - bieguny transmitancji

System nazywamy przyczynowym, jeśli każda próbka jego odpowiedzi zależy tyko i wyłącznie od poprzedzających ją próbek wymuszenia. Jeśli system o transmitancji H(z) jest przyczynowy to:

H[n]<0, dla n<0

i obszar zbieżności ROC jego transmitancji leży na zewnątrz bieguna di o największym promieniu.

System nazywamy stabilnym wtedy i tylko wtedy gdy ograniczony ciąg wejściowy daje w odpowiedzi także ciąg ograniczony. Jeśli odpowiedź impulsowa h[n] systemu jest bezwzględnie sumowalna to jest on systemem stabilnym:

0x01 graphic

Obszar zbieżności ROC stabilnego i przyczynowego systemu H(z) zawiera okrąg jednostkowy. Wszystkie bieguny di takiego systemu leżą wewnątrz okręgu jednostkowego, tzn: 0x01 graphic
, i=1,...,m

Np. 0x01 graphic

0x08 graphic

System FIR jest zawsze stabilny (w odróżnieniu od IIR).

E53. Objaśnić jak zbadać, czy przyczynowy system DLS jest stabilny. Pokaż, czy stabilny jest system o transmitancji H(z) = 1/[(z - 1/8)*(z - ¼)*(z - ½)*(z - 1)]

Jeżeli nie jest stabilny to co należałoby zmienić żeby był stabilny? Odp. Uzasadnij.

ROC stabilnego, przyczynowego systemu H(z) zawiera okrąg jednostkowy. Wszystkie bieguny di takiego systemu leżą wewnątrz okręgu jednostkowego tzn. | di | < 1, i = 1, … ,n.

H(z) = 1/[(z - 1/8)*(z - ¼)*(z - ½)*(z - 1)]

0x08 graphic

Bieguny: d1 = 1/8, d2 = ¼, d3 = ½ d4 = 1

System nie jest stabilny, gdyż biegun d4 nie leży wewnątrz okręgu jednostkowego.

Żeby był stabilny trzeba biegun d4 przesunąć w lewo aby znalazł się wewnątrz

okręgu jednostkowego.

E54. Czy stabilne są następujące, przyczynowe systemy DLS:

a)0x01 graphic

ponieważ jest to transmitancja FIR'u od razu można stwierdzić iż jest to układ stabilny

!!! ”system FIR jest zawsze stabilny ” !!!

0x01 graphic

0x01 graphic

0x01 graphic

zera: ... nieważne [położenie zer nie stanowi o stabilności układu]

bieguny: 0 (podwójny)

0x01 graphic

Ponieważ 0 < 1 to układ jest stabilny

b) 0x01 graphic

0x01 graphic

0x01 graphic

zera: ...

bieguny: 0x01 graphic

0x01 graphic
ponieważ 0x01 graphic
i 0x01 graphic
(bieguny znajdują się wewnątrz okręgu jednostkowego) to układ JEST STABILNY.

E55.Zapisz ogólne transmisje systemów DLS: IIR i FIR, i objaśnij znaczenie występujących w nich symboli.

-system IIR H(z)=Y(z)/X(z) {przy z.w.p.}=( b0+b1·z-1+...+bz-M )/(1+a1·z-1+...+aM·z-M) gdzie: M-rząd systemu

-system FIR H(z)=b0+b1z-1+.....bMz-M = N=0M h[n]·z-N

gdzie: Y(z)- transformata z odpowiedzi; X(z)-transformata z pobudzenia, ai bi -współ. wielomianu( współ. w równaniu różnicowym układu, jako odpowiednie wzmocnienia)

E56. Jakimi sposobami można obliczyć odpowiedź systemu DLS na dowolne pobudzenie, przy zerowych warunkach początkowych?

-na dwa sposoby(a,b): 1)odp. y[n] systemu DLS na dowolne pobudzenie x[n] jest splotem dyskretnym ciągów: pobudzenia x[n] i odp. impulsowej h[n]. Splot dyskretny jest przemienny- a)y[n]=x[n]*h[n]=h[n]*x[n]=Σ(od k=-∞ do +∞) x[k] h[n-k]= Σ(od k=-∞ do +∞) h[k] x[n-k]; 2)transmitancja systemu DLS: H(z)=Y(z)/X(z) przy z.w.p., gdzie x[n]↔X(z), y[n]↔Y(z); ponieważ: h[n]↔H(z) tzn. H(z)=Z{h[n]}= Σ(od n=-∞ do +∞) h[n]z^(-n); czyli dla zadanej H(z): b)y[n]=Z^(-1){X(Z) H(Z)}

E57. Zdefiniuj elementy schematu blokowego systemu DLS. Narysuj schematy blokowe: a) Systemu opóźniającego, b) systemu pierwszej różnicy wstecz, c) akumulatora.

DLS - dyskretny, liniowy, stacjonarny

0x08 graphic

0x08 graphic
0x08 graphic
x[n] y[n] y[n] = T{x[n]}

a) system opozniajacy y[n] = T{x[n]} = x[n-n0]

0x08 graphic
0x08 graphic
0x08 graphic
x[n] x[n-n0]

H(z) = z-n0

b) system pierwszej różnicy wstecz y[n] = x[n] - x[n-1]

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
Dyferator

x[n] y[n]

x[n-1]

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
c) akumulator y[n] =0x01 graphic
= y[n-1] + x[n]

x[n] y[n]

y[n-1]

E58. Narysuj schematy blokowe systemu IIR pierwszego rzędu i drugiego rzędu: a) strukturę bezpośrednią 1D i b) strukturę bezpośrednią 2D i zapisz ich algorytm.

System IIR :

y[n] = 0x01 graphic
bi [ n-i ] -0x01 graphic
ai y [n - i ]

m - rząd systemu

a) pierwszego rzędu IIR => m=1

H(z) = Y(z) / X(z) | zwp = (b0 + b1z-1) / ( 1 + a1 z-1)

Algorytm:

Y[n] = b0 x[n] + b1 x[n - 1] - a1 y [n-1]

0x08 graphic
Schemat blokowy:

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic

b) drugiego rzędu IIR => m=2

H(z) = Y(z) / X(z) | zwp = (b0 + b1z^(-1) + b2 z^(-2)) / ( 1 + a1 z^(-1) + a2 z^(-2))

Algorytm:

Y[n] = b0 x[n] + b1 x[n - 1] + b2 x[n - 2] - a1 y [n-1] - a2 y [n-2]

Schemat blkowy

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x01 graphic

E59. Zdefiniuj DFT i przekształcenie do niego odwrotne.

0x08 graphic
DFT- dyskretne przekształcenie Fouriera; definiowane tylko dla ciągów o skończonej długości:{x[n]}(od n=0 do n=N-1), gdzie N- skończona długość; definicja DFT: X[k]= Σ(od n=0 do N-1) x[n]e^(-j2Πkn/N), dla k=0,1,…N-1; definicja przekształcenia odwrotnego do DFT (IDFT): x[n]= 1/N Σ(od n=0 do N-1) X[k]e^(j2Πkn/N), dla n=0,1,…,N-1;

E60. Jakie zachodza zwiazki pomiedzy DFT,DTFT i transformata Z przyczynowego ciagu x[n] przy zalozeniu, ze obszar zbieznosci X(z) zawiera okrag jednostkowy

0x08 graphic

E61.Wypisz i nazwij zapamiętane właściwości DFT.

-liniowość (oba ciągi o tej samej dł., w przeciwnym razie konieczność uzupełnienia zerami) ax[n]+by[n]→ aX[k]+bY[k]

-przesunięcie kołowe w dziedzinie n (n0<0 opóźnienie całkowite) x[n+n0] N→ X[k]W N -kno

-modulacja x[n] W N kno→ X[[k+k0 ]]N

-splot kołowy w dziedzinie n n=0N-1 x[m] ·y[[n-m]]N → X[k] ·Y[k]

-mnożenie w dziedzinie czasu x[n] ·y[n] → 1/N l=0N-1 X[l] ·Y[k-l]N

-sprzężenie x*[n] → x*[[-k]]N -sprzężenie i obrót x*[[-n]]N→ X*[k]

-symetria x[n] → N·x[[-k]]N

-związek Parsevala l=0N-1 ׀x[n]׀2 = 1/N k=0N-1 ׀x[k]׀2

E62. Czym roznia sie widma ciag x[n] obliczone za pomoca DTFT i DFT? Pokaz na dowolnym przykladzie

0x08 graphic

0x08 graphic

E63 Objasnij zastosoawnie notacji macierzowej w obliczaniu DFT. Mozesz posluzyc sie przykladem.

0x08 graphic

0x08 graphic

0x01 graphic

E65. Na czym polega szybki splot i kiedy warto go stosować?

Splot kołowy (szybki) za pośrednictwem DFT można wykorzystać do obliczeń splotu liniowego. Stosuje się go gdy liczba próbek jest większa niż 32 (próbki zespolone) lub 64 (próbki rzeczywiste).

R1. Podać znane Ci metody zobrazowania analogowego sygnału zespolonego

a) interpretacja geometryczna, graficzna, czyli zobrazowanie sygnału zespolonego na:

- dwa wykresy (w jednakowej skali) obu składowych kartezjańskich x(t)=re u(t) i y(t)=im u(t)

- wykres dwuwymiarowy (2D) na płaszczyźnie zespolonej (płaszczyźnie Arganda) = trajektoria zespolona; wymaga dodatkowego skalowania w czasie

- wykres trójwymiarowy (3D) - spirala wije się wokół osi t

b) sygnał zespolony u(t) zdefiniowany w postaci algebraicznej można wyrażać w postaci trygonometrycznej/biegunowej (we współrzędnych biegunowych):

u(t)=a(t)exp(jφ(t)); a≥0 oraz 0x01 graphic

R2. Na czym polega faktoryzacja biegunowa sygnału zespolonego w dziedzinie czasu.

Faktoryzacja biegunowa sygnału zespolonego 0x01 graphic
; γ(t) to faktor-wskaz (o jednostkowej amplitudzie).

Po stronie częstotliwości: 0x01 graphic
, * to symbol operacji splotu.

Fazor sygnału zespolonego: 0x01 graphic
, reprezentuje cechy/atrybuty strukturalne sygnału u(t), takie jak jego „szybkość kątowa”/pulsację chwilową ω(t) i orientację/fazę początkową φ0.

R3. Zdefiniować amplitudę chwilową, fazę chwilową i pulsację chwilową danego sygnału zespolonego.

Amplituda chwilowa - jest wartością bezwzględną (modułem) zespolonej wartości chwilowej: 0x01 graphic
. Przebieg amplitudy chwilowej często nazywa się obwiednią przebiegów rzeczywistych, zarówno x(t) jak i y(t), gdyż -a(t)≤x(t),y(t)≤a(t)

Faza chwilowa - argument (część urojona) zespolonego logarytmu funkcji zespolonej u(t) zmiennej rzeczywistej t: 0x01 graphic
.

W powyższej definicji fazy chwilowej zarówno arg jak i ln są operatorami (nie funkcjami „bezinercyjnymi” jak np. |( )|, czy exp( )) zdefiniowanymi przez Cauchy'ego jako arg≡im ln oraz 0x01 graphic

Pulsacja chwilowa sygnału u(t) to pochodna jego fazy chwilowej φ(t) względem czasu t: 0x01 graphic
. Korzystając z definicji fazy chwilowej otrzymujemy: 0x01 graphic

R4. Na czym polega rozwijanie fazy zespolonej sygnału analogowego? Kiedy lub po co jest ona stosowana?

Fazę chwilową obliczoną trygonometrycznie/cyklometrycznie ze wzoru: 0x01 graphic
, nazywamy główną lub zwiniętą w odróżnieniu od fazy rozwiniętej obliczanej według wzoru: 0x01 graphic

Przedział 0x01 graphic
. Przedział 0x01 graphic
nazywa się zwojem (lub okresem) głównym fazy funkcji określonej ex(jφ). Faza rozwinięta φ(t) i faza główna φP(t) są powiązane wzorem: 0x01 graphic

Obliczanie φ(t) w oparciu φP(t) nazywa się rozwijaniem fazy, a odwrotnie: φP(t) w oparciu o φ(t) - zwijanie fazy.

R5. Moc chwilowa i twierdzenie Parsevala dla sygnałów analogowych zespolonych. Energia sygnału zespolonego.

Moc chwilowa sygnału zespolonego u(t): 0x01 graphic

Energia (całkowita) sygnału zespolonego: 0x01 graphic
może być skończona albo nieskończona. Z tego wywodzą się dwie klasy sygnałów analogowych:

- o skończonej energii: 0x01 graphic

- o skończonej mocy (chwilowej): 0x01 graphic

Równość Parsevala dla sygnałów zespolonych:

0x01 graphic
, ale też:

0x01 graphic

R6. Zinterpretować sygnał zespolony na płaszczyźnie zespolonej Arganda
DIAGRAM ARGANDA

0x01 graphic

Diagram Arganda jest przedstawia liczby zespolone jako punkty

z=x+jy

oś iksów jako oś części rzeczywistej i oś igreków jako oś części zespolonej

na wykresie powyżej przerywany okrąg reprezentuje zespolony moduł |z| liczby `z'

a kąt 0x01 graphic
reprezentuje jej zespolony argument

Diagram Argand'a tak naprawdę został wymyślony wcześniej przez C.Wessel'a

R7. Omówić znane Ci reprezentacje biegunowe sygnału x(t)+jy(t); x(t), y(t) 0x01 graphic
R 0x01 graphic

Kartezjańska reprezentacja sygnału x(t): u(t)=x(t)+jy(t) 0x01 graphic
równoważnik (sygnału) analityczny

Biegunowe reprezentacje sygnału x(t):

0x01 graphic
(AM - modulacja amplitudy, FM - modulacja częstotliwości)

0x01 graphic
zespolona reprezentacja dynamiczna

0x01 graphic
log obwiednia, poziom chwilowy [dB]

0x01 graphic
pulsacja chwilowa, 0x01 graphic
- częstotliwość chwilowa

0x01 graphic
zespolona pulsacja chwilowa

0x01 graphic
(PM - modulacja fazy)

0x01 graphic
zespolona faza chwilowa

0x01 graphic
”zapis nutowy”; 0x01 graphic
wysokość chwilowa (instantaneous pitch), gdy p=2 mianowana w oktawach, gdy p=21/12 - w półtonach

R8. Omówić główne (kluczowe) właściwości reprezentacji biegunowych

  1. mnożenie (wartości chwilowych) przebiegów zespolonych skutkuje mnożeniem ich amplitud chwilowych oraz dodawaniem ich faz chwilowych i pulsacji chwilowych

  2. niech a(t), ϕ(t) i ω(t) są odpowiednio amplitudą chwilową, fazą chwilową i pulsacją chwilową przebiegu zespolonego, to:

R9.Amplituda chwilowa, faza chwilowa i pulsacja chwilowa analogowego sygnału rzeczywistego - postulaty Valkmana

Amplituda chwilowa (obwiednia), faza i pulsacja chwilowa sygnału rzeczywistych muszą spełniać trzy postulaty Valkmana, które spełniają amplituda chwilowa, faza i pulsacja sygnałów zespolonych.

  1. Małe zmiany w przebiegu wartości chwilowych sygnałów rzeczywistego x(t), co zapisujemy jako: x(t) + δx(t), powinny powodować również małe zmiany w przebiegu jego amplitudy chwilowej (obwiedni) a(t), co zapisujemy jako

Obwiednia {x(t) + δx(t)} ⇒ obwiednia {x(t)}, gdy {δx(t)} 

  1. Faza i częstotliwość chwilowa powinny być niezależne od wielkości skali/sygnału

c0 R faza chwilowa {cx(t)} = faza chwilowa{x(t)}

  1. Amplituda chwilowa, faza chwilowa i pulsacja chwilowa ”czystej” sinusoidy rzeczywistej (a0≠0, ω0>0)

a0cos(ω0t + ϕ0) powinny być równe parametrom tej sinusoidy

amplituda chwilowa {a0cos(ω0t + ϕ0)} = a0

faza chwilowa {a0cos(ω0t + ϕ0)} = ω0t + ϕ0

pulsacja chwilowa {a0cos(ω0t + ϕ0)} = ω0

R10. Amplituda chwilowa, faza chwilowa i pulsacja chwilowa analogowego sygnału rzeczywistego - koncepcja Gabora - Ville (klasyczna)

Koncepcja utworzenia sygnału analitycznego u(t) z sygnału rzeczywistego x(t) poprzez usunięcie „połowy widma” X(jΩ) zajmującego ujemna półoś jΩ wywodzi się z możliwości (czy chęci)m dwukrotnej redukcji „nadmiarowości widmowej” sygnału rzeczywistego x(t), którego widmo jest Hermitowsko symetryczne/parzyste: X(-jΩ)= X*(jΩ). Wzmocnienie dwukrotne widma pozostawionego na dodatniej półosi jΩ ma zapewnić równość mocy i energii danego sygnału rzeczywistego x(t) i jego równoważnika zespolonego u(t) oraz to, by x(t) pozostał częscią rzeczywistą u(t), czyli

x(t)=Reu(t) F X(jΩ) = (U(jΩ) + U*(-jΩ))/2, F - oznacza parę transformat Fouriera

Jedynym sygnałem zespolonym, który w taki (liniowy) sposób jest związany z danym sygnałem rzeczywistym x(t) jest sygnał u(t) którego widmo

U(jΩ)=0x01 graphic

i który Gabor nazwał sygnałem analitycznym a Ville - preobwiednią (preenvelope). Część urojona sygnału x(t) spełnia oczywistą zależność

y(t)=Imu(t) F Y(jΩ) = (U(jΩ) - U*(-jΩ))/(2j)

R11.Zdefiniować (w dziedzinie czasu i w dziedzinie częstotliwości) przekształcenie Hilberta i operator wykonujący to przekształcenie.

y(t) = CPV 0x01 graphic
, 0<a<A<∞

CPV- wartość główna Cauchy'ego

W wyniku całkowania otrzymamy przebieg y(t) nazwany transformatą Hilberta przebiegu x(t). Transformata Hilberta danego przebiegu x(t) powstaje w wyniku filtracji lilniowej splotu liniowego danego przebiegu x(t) z odpowiedzią impulsową

hT(t)=0x01 graphic

Nieprzyczynowy o zerowym opóźnieniu filtr HT o podwyższonej odpowiedzi impulsowej charakterystyce częstotliwościowej

HT(j2πF) = -jsign(F); F=0x01 graphic
nazywamy transformatorem Hilberta.

R12. Omówić własności par transformat Hilberta i sygnałów analitycznych

Własności pary transformat Hilberta = właściwości składowych sygnału analitycznego = właściwości zespolonego równoważnika x(t)+jy(t) sygnału rzeczywistego x(t) wynikają m. in. z faktu, że HT-1=-HT=HT*.

  1. transformator Hilberta nie przenosi składowej stałej

  2. Jeśli y(t)=HT{x(t)}, to HT-1{y(t)}=0x01 graphic

Gdzie 0x01 graphic
<=składowa zmienna sygnału y(t)

0x01 graphic
<= składowa stała, centroid (wartość średnia) y(t)

  1. dwukrotne przejście przebiegu x(t) przez transformator Hilberta HT (=przejście przez kaskadę dwóch transformatorów Hilberta HTHT) zmienia znak sygnału na przeciwny (i usuwa z niego część stałą), czyli HTHT-1{x(t)}=0x01 graphic

  2. Widma: amplitudowe, mocy/energii i opóźnienia grupowego przebiegów x(t) i y(t) są identyczne za wyjątkiem częstotliwości F=0

  3. Widma fazowe przebiegów x(t) i y(t) różnią się o π/2 (90o) => przebiegi x(t) i y(t) są ortogonalne

  4. Niech x(t)=xe(t)+xo(t), gdzie:

xe(t):=(x(t)+x(-t))/2 <= składowa parzysta,

xo(t):=(x(t)-x(-t))/2 <= składowa nieparzysta prebiegu x(t),

Niech y(t)=HT{x(t)}=ye(t)+yo(t). Wtedy HT{xe(t)}=yo(t) i HT{xe(t)}=ye(t)

  1. Twierdzenie Bedrosiana: Niech x(t)=a(t)*b(t) przy czym widma czynników a(t) i b(t) nie zachodzą na siebie (co można zapisać jako A(jΩ)B(jΩ)=0 ∀Ω∈(-∞;+∞)), wtedy y(t)=HT{a(t)b(t)}=a(t)HT{b(t)}=b(t)HT{a(t)}

  2. Trajektoria sygnału analitycznego x(t)+jy(t) „jest dodatnioskrętna”, czyli zarówno lokalnie (wokół chwilowego środka krzywizny) jak i globalnie (wokół „środka ciężkości” = centroid ) wije się w kierunku dodatnim.

R13. Zdefiniować podstawowe parametry analogowego sygnału pasmowego (rzeczywistego i zespolonego): pasmo, szerokość pasma, częstotliwość środkowa.

To F=0 jest częstotliwością środkową sygnału podstawowo pasmowego: FC=0, a szerokość pasma jest równa podwójnej górnej częstotliwości granicznej B=2Fu

0x08 graphic

Zdefiniowany jako w(t)= exp(-2jπFC) * HA {x(t)} =2HLP {exp (-2jπFC) * x(t)}

Gdzie :

HA - zespolony filtr Hilberta

H LP - filtr dolnoprzepustowy o paśmie przepustowym (-B/2, B/2) i paśmie zaporowym powyżej 2FC-B/2

0x08 graphic

R14. Zdefiniować obwiednię zespoloną rzeczywistego sygnału pasmowego.

Definicja obwiedni:

w(t)=exp(-2jπFC) * HA {x(t)} =2HLP {exp (-2jπFC) * x(t)}

R15. Omówić własności obwiedni zespolonej.

Zdefiniowany wyżej sygnał w(t) nazywa się obwiednią zespoloną sygnału rzeczywistego x(t) pasmowego / wąskopasmowego.

Obwiednia zespolona:

w(t)=i(t)+jq(t) ; i(t), q(t)∈R

jest podstawową kartezjańską reprezentacją zespoloną sygnału x(t) o najniższej częstotliwości granicznej FU=B/2, FL(w)=-B/2.

Składowe kartezjańskie obwiedni zespolonej nazywają się odpowiednio:

38

T{ 0x01 graphic
}

x[n]

y[n]

T{ 0x01 graphic
}

x[n]

y[n]

y[n]

n

0

1

2

3

4

5

6

7

x[n]

n

0

1

2

3

4

h2[n]

h1[n]

h1[n]

h2[n]

.3 2 1

x 1 .2 3

9 6 3

6 4 2

+3 2 1

3 .8 14 8 3

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

Arctgx(0x01 graphic
)

Arctg(0,5)

T{•}

z-n0

+

z-1

z-1

+

y[n]

+

0x01 graphic

+

0x01 graphic

z-1

0x01 graphic

x[n]

b0

0x01 graphic

-a1

0x01 graphic

b1

0x01 graphic

0x01 graphic

y[n]

x[n]

+

0x01 graphic

+

0x01 graphic

+

0x01 graphic

+

0x01 graphic

z-2

0x01 graphic

z-1

0x01 graphic

-a1

0x01 graphic

-a2

0x01 graphic

b2

0x01 graphic

b1

0x01 graphic

b0

0x01 graphic



Wyszukiwarka