Transformacja Fouriera - Fourier Transformation

Analiza harmonicznych

Cyfrowa obróbka sygnałów Digital Signal Processing (DSP) wykorzystuje tzw. szybką transformatę Fouriera Fast Fourier Transform (FFT).
Funkcje DSP mogą być używane w symulacji czasowej (Transient) oraz częstotliwościowej (AC).
Jednym z zastosowań szybkiej transformaty Fouriera jest analiza widma częstotliwościowego sygnałów, czyli harmonicznych.
Procedura FFT wymaga dwóch podstawowych parametrów: Parametry TP, N, fmax mogą być też zaprogramowane na karcie DSP Parameters.

Poniższa tabelka przedstawia zakres obliczanych harmonicznych:
Składowa harmonicznaCzęstotliwość harmonicznejUwagi
1F = 0Hzskładowa stała (DC)
2F = FPskładowa podstawowa (pierwsza harmoniczna)
3F = 2×FPdruga harmoniczna
4F = 3×FPtrzecia harmoniczna
.........
NF = (N-1)×FP(N-1) harmoniczna
Składowe harmoniczne począwszy od drugiej nazywany wyższymi harmonicznymi.
Dla celów analizy harmonicznych będziemy używali funkcji HARM(u), gdzie u jest to analizowany przebieg czasowy.

Programowanie analizy harmonicznych (przykład 1)

Przedstaw harmoniczne przebiegu sinusoidalnego o częstotliwości f=1kHz i amplitudzie A=4V.

Schemat układu do symulacji


Programujemy odpowiednie źródło sinusoidalne. Model źródła nazwiemy np. SIN. Kartę źródła wypełnij według rysunku (żółte pola).


Programujemy symulację czasową wraz z analizą harmonicznych. Kartę Transient Analysis Limits wypełnij według rysunku (żółte pola).
Okres przebiegu T=1/f=1/1kHz=1ms. Chcemy wyświetlić np. pięć okresów przebiegu więc w polu Time Range wpisaliśmy wartość 5m (5×1ms=5ms).
Wobec tego dla analizy harmonicznych mamy:
tmax=5ms,
TP=tmax=5ms,
FP=1/TP=1/5ms=200Hz.
Chcemy wyświetlić prążki o częstotliwościach do np. 5kHz.

W wyniku symulacji otrzymujemy wykresy:
Symulator obliczył prążki:
Częstotliwość [Hz]Amplituda [V]
00
200 = FP0
400 = 2×FP0
600 = 3×FP0
800 = 4×FP0
1000 = 5×FP4
1200 = 6×FP0
......
Częstotliwość 1kHz jest więc piątą harmoniczną częstotliwości próbkowania FP=200Hz. Częstotliwość FP jest traktowana jako składowa podstawowa.

Aby częstotliwość próbkowania była równa częstotliwości analizowanego przebiegu (FP=1kHz) należy więc odpowiednio dobrać okres próbkowania TP=T=1ms (tak, aby był równy okresowi analizowanego przebiegu).

W tym celu otwieramy kartę DSP Parameters (patrz rysunek).

Do analizy harmonicznych wybieramy jeden pełny okres przebiegu - np. ostatni. W naszym przypadku otrzymamy:
TP=5ms-4ms=1ms,
FP=1/TP=1/1ms=1kHz
N=1024 (wartość domyślna).

Nie zapomnij włączyć - On Status


Po ponownym uruchomieniu symulacji otrzymujemy teraz wykresy:
Symulator obliczył prążki:
Częstotliwość [kHz]Amplituda [V]
00
1 = FP4
2 = 2×FP0
3 = 3×FP0
4 = 4×FP0
......
Częstotliwość 1kHz jest więc tym razem składową podstawową. Symulacja została prawidłowo zaprogramowana (oczywiście idealny przebieg sinusoidalny nie posiada żadnych harmonicznych tylko samą składową podstawową - tylko jeden prążek).

Programowanie analizy harmonicznych - zniekształcenia skrośne
(przykład 2)

Przedstaw zniekształcenia skrośne we wzmacniaczu mocy klasy B oraz oblicz współczynnik zawartości harmonicznych h uwzględniając pięć pierwszych harmonicznych.

Zniekształcenia skrośne są jednym z rodzajów zniekształceń nieliniowych.
Wielkość zniekształceń nieliniowych podaje się w postaci współczynnika zawartości harmonicznych h, który zdefiniowany jest następująco:

gdzie:
U1 - wartość skuteczna składowej podstawowej
U2, U3, U4,... - wartości skuteczne wyższych harmonicznych
A1 - amplituda składowej podstawowej
A2, A3, A4,... - amplitudy wyższych harmonicznych
- wartość skuteczna całego przebiegu zniekształconego
- wartość skuteczna wyższych harmonicznych
- związek między wartością skuteczną a amplitudą harmonicznej numer i.


Schemat układu do symulacji


Do wejścia wzmacniacza doprowadzamy napięcie sinusoidalne o częstotliwości f=1kHz. Przeprowadzimy dwa razy symulację zmieniając amplitudę modelu źródła gen_1kHz (A=2.5V oraz A=7V).


W wyniku pierwszej symulacji otrzymujemy wykresy:

Przy doprowadzonym niewielkim napięciu wejściowym o amplitudzie 2.5V (napięcie przewodzenia tranzystora wynosi około 0.7V), na wyjściu widzimy wyraźnie zniekształcenia skrośne.

Amplitudy prążków:
f[kHz]A[V]
00
11.39
20
30.29
40
50.11

Obliczamy współczynnik zawartości harmonicznych uwzględniając maksymalnie piątą harmoniczną




W wyniku drugiej symulacji otrzymujemy wykresy:

Przy doprowadzonym większym napięciu wejściowym o amplitudzie 7V (napięcie przewodzenia tranzystora wynosi około 0.7V), na wyjściu widzimy znacznie mniej wyraźnie zniekształcenia skrośne.

Amplitudy prążków:
f[kHz]A[V]
00
15.82
20
30.37
40
50.20

Obliczamy współczynnik zawartości harmonicznych uwzględniając maksymalnie piątą harmoniczną


Opracował: Stanisław Pelczar PZNr10 SME Kęty