PS Moduł 5: Różnice pomiędzy wersjami

Wersja z 13:27, 5 wrz 2006

Rozdzielczość DTF jest definiowana jako liczba próbek widmowych przypadająca na jeden okres widma (równy $ω_{s}$ w skali pulsacji lub $f_{s}$ w skali częstotliwości). Jest to zarazem odwrotność odległości między kolejnymi próbkami widma.
Rozdzielczość DTF jest cechą zależną tylko od czasu obserwacji $T$ sygnału, bowiem $r_{f} = N / f_{s} = N T_{s} = T$ .
Metoda uzupełniania zerami polega na przyjęciu zerowych wartości dodatkowych próbek sygnału w punktach $n = N, . . ., N_{1} - 1$ i obliczeniu dla tak zdefiniowanego sygnału $N_{1}$ -punktowej DTF.

Dodanie do ośmiu właściwych próbek sygnału ośmiu dodatkowych próbek zerowych zwiększa rozdzielczość DTF dwukrotnie.
Próbki widmowe leżą zawsze na krzywej ciągłej $A (e^{j θ})$ .
Wraz ze wzrostem dodatkowych zer wzrasta nakład obliczeniowy.

Zgodnie z właściwością symetrii wystarczy obliczyć tylko próbki widmowe $X (0), X (1), X (2), X (3)$ . Jak widzimy tylko próbki $X (1)$ i $X (5)$ (lub równoważnie $X (- 1)$ ) są niezerowe.
Gdy analogowy sygnał harmoniczny jest obserwowany dokładnie tylko przez jeden jego okres, wówczas DTF wyznaczone na podstawie jego próbek daje prawidłowy obraz jego widma. Niezerowe próbki widmowe otrzymujemy tylko dla pulsacji unormowanych $θ \pm π / 3$ odpowiadających częstotliwościom nieunormowanym $f = \pm 1 k H z$ .

Punkty $θ_{k} = π k / 4$ odpowiadają częstotliwością nieunormowanym $f_{k} = k \times 0, 75 k H z$ , nie pokrywającymi się z właściwą częstotliwością sygnału $1 k H z$ .
Można mówić o przecieku mocy sygnału skoncentrowanej w prążku widmowym o częstotliwości $1 k H z$ do prążków widmowych występujących w fałszywych punktach osi częstotliwości.
Obserwacja sygnału $x [n]$ w skończonym przedziale czasu $n = N, . . ., N - 1$ odpowiada mnożeniu tego sygnału przez prostokątne okno czasowe. Mówimy wówczas o okienkowaniu sygnału.

Obserwacja sygnału $x [n]$ w skończonym prostokątnym oknie czasowym $g [n]$ odpowiada splataniu widma tego sygnału przez widmo okna opisane wzorem (5.2).
Na zniekształcenia widma amplitudowego nie ma wpływu czynnik wykładniczy widma (5.2), a decyduje o nich drugi czynnik, którego struktura jest listkowa.
W przypadku okna prostokątnego stosunek wysokości pierwszego listka bocznego do wysokości listka głównego nie może być mniejszy od $0, 212$ .

Przykład 5.2 dobrze ilustruje zjawisko przecieku w przypadku sygnału harmonicznego. Dokładne widmo tego sygnału składa się z (okresowo powtarzanej z okresem $2 π$ ) pary dystrybucji w punktach $\pm θ_{0}$ (rys. a). Skończony czas obserwacji powoduje, że dystrybucje te rozmywają się na dwa widma listkowe rozmieszczone wokół punktów $\pm θ_{0}$ .

Na rys. b i c pokazano wykresy widma amplitudowego sygnału harmonicznego z przykładu 5.2 o pulsacji unormowanej $θ_{0} = π / 3$ okienkowanego oknem prostokątnym w dwóch przypadkach: $N = 6$ oraz $N = 8$ .
Wykres wyjaśnia dlaczego w przypadku wyznaczenia 6-punktowej DTF (rys. b) nie występuje zjawisko przecieku, tzn. otrzymujemy prawidłowe próbki widma sygnału harmonicznego w punktach $\pm π / 3$ (dla $k \pm 1$ ). Pozostałe próbki widmowe, składające się na DFT, wypadają bowiem dokładnie w punktach zerowych rozmytego widma $X_{0} (e^{j θ})$ . Natomiast w przypadku wyznaczania 8-punktowej DTF (rys. c) próbki widmowe wypadają między zerami widma i przeciek występuje.
Zjawisko przecieku objawia się w podobny sposób w przypadku, gdy w skończonym oknie czasowym są obserwowane sygnały nieokresowe. Energia (lub moc) sygnału przecieka wówczas nie do obcych prążków, ale w inne ciągłe przedziały osi $θ$ .

Widma wszystkich pokazanych na wykresach okien czasowych mają strukturę listkową. Zniekształcenia widma sygnału spowodowane skończonym czasem jego obserwacji są tym mniejsze im mniejszy jest poziom listków bocznych i jednocześnie węższy jest listek główny. Wymagania te są jednak sprzeczne.
Dobierając inne od prostokątnego okna czasowe, zmniejszamy poziom listków bocznych. Powoduje to jednak jednoczesny wzrost szerokości listka głównego.
Na wykresach pokazano przykłady okien czasowych różnych od prostokątnego dla przypadku $N = 10$ .

Zauważmy, że obie krańcowe wartości okna Hamminga (rys. e) są różne od zeru (równe $0, 08$ ). Ta drobna różnica powoduje, że poziom listków bocznych jest dla tego okna mniejszy w porównaniu np. z oknem Hanna (rys. d) lub oknem Blackmana (rys. f).
W porównaniu z oknem prostokątnym wszystkie okna pokazane na rys. a-f zapewniają mniejszy poziom listków bocznych, ale kosztem wzrostu szerokości listka głównego. Dla okien tych wielkości te przy zadanym $N$ są ustalone i nie można na nie wpływać. Istnieją jednak bardziej złożone okna, w których poziom listków bocznych i szerokość listka głównego można regulować i dobierać ich wartości w drodze kompromisu.
Im poziom listków bocznych widma okna jest mniejszy, tym większa jest tzw. rozróżnialność częstotliwościowa DTF. Natomiast im szerokość listka głównego jest mniejsza, tym większa jest tzw. rozróżnialność amplitudowa DTF.

W szacowaniu złożoności obliczeniowej DTF zakładamy, że sygnał $x [n]$ jest zespolony.
Liczba $W = e^{j 2 π / N}$ ma ciekawe właściwości. Jest ona równa pierwiastkowi $N$ -tego z jedności. Jej $N$ -ta potęga jest równa jedności: $W^{n} = (e^{j 2 π / N})^{N} = e^{j 2 π} = 1$ . Ponadto w wyniku próbkowania $N$ razy w okresie w chwilach $n T_{s}$ zespolonego sygnału harmonicznego $e^{j ω_{0} t}$ otrzymujemy próbki będące kolejnymi potęgami liczby $W : e^{j ω_{0} n T_{s}} = e^{j 2 π n T_{s} / T_{0}} = (e^{j 2 π / N})^{n} = W^{n}$ .
O złożoności obliczeniowej DTF decyduje liczba mnożeń. Dlatego przy jej szacowaniu nie uwzględnia się dodawań. Liczba mnożeń wynosi dokładnie $(N - 1)^{2}$ , gdyż obliczenie próbki widmowej dla $k = 0$ wymaga wykonania tylko dodawań. Dla większych $N$ można dla prostoty szacowania przyjąć, że liczba mnożeń jest równa $N^{2}$ .

@@ Linia 96: / Linia 96: @@
 |width="500px" valign="top"|[[Grafika:PS_M5_Slajd10.png]]
 |valign="top"|
+*W szacowaniu złożoności obliczeniowej DTF zakładamy, że sygnał <math>x[n]\,</math>  jest zespolony.
+*Liczba <math>W=e^{j2\pi/N}</math> ma ciekawe właściwości. Jest ona równa pierwiastkowi <math>N\,</math> -tego z jedności. Jej <math>N\,</math> -ta potęga jest równa jedności: <math>W^n=(e^{j2\pi/N})^N=e^{j2\pi}=1</math> . Ponadto w wyniku próbkowania <math>N\,</math>  razy w okresie w chwilach <math>nT_s\,</math>  zespolonego sygnału harmonicznego <math>e^{j\omega_0 t}\,</math>  otrzymujemy próbki będące kolejnymi potęgami liczby <math>W\, :\, e^{j\omega_0 nT_s}=e^{j2\pi nT_s/T_0}=(e^{j2\pi /N})^n=W^n</math> .
+*O złożoności obliczeniowej DTF decyduje liczba mnożeń. Dlatego przy jej szacowaniu nie uwzględnia się dodawań. Liczba mnożeń wynosi dokładnie <math>(N-1)^2\,</math> , gdyż obliczenie próbki widmowej dla <math>k=0</math>  wymaga wykonania tylko dodawań. Dla większych <math>N\,</math>  można dla prostoty szacowania przyjąć, że liczba mnożeń jest równa <math>N^2\,</math> .
 |}