Wykresy dwóch funkcji rzeczywistych oraz ich transformaty Fouriera. Drugi sygnał, przesunięty w czasie względem pierwszego, pomimo identycznej postaci, ma inną transformatę Fouriera
Transformacja Fouriera – transformacja całkowa, która przyjmuje funkcję rzeczywistą
f
(
t
)
,
t
∈
R
{\displaystyle f(t),t\in R}
jako daną wejściową i wyprowadza funkcję o wartościach zespolonych, która opisuje stopień, w jakim w funkcji wejściowej występują funkcje zespolone
f
ω
=
e
−
i
ω
t
.
{\displaystyle f_{\omega }=e^{-i\omega t}.}
Transformatą Fouriera nazywa się funkcję, która powstaje w wyniku działania transformacji Fouriera.
Zastosowanie transformacji Fouriera do określenia składowych akordu C dur fortepianu. Pierwsze trzy piki od lewej odpowiadają podstawowym częstotliwościom C, E, G. Pozostałe są wyższymi harmonicznymi.
Przykładowo, obliczanie transformaty Fouriera trwającego jakiś czas sygnału akustycznego oznacza obliczanie amplitud składowych zespolonych tego dźwięku (por. wykresy transformaty obok dla funkcji schodkowych).
Transformacja Fouriera jest operatorem liniowym działającym na przestrzeni funkcyjnej funkcji zmiennej rzeczywistej i dokonującym rozkładu danej funkcji
f
{\displaystyle f}
w bazie ortonormalnej zespolonych funkcji eksponencjalnych
f
ω
=
e
−
i
ω
t
,
ω
∈
R
{\displaystyle f_{\omega }=e^{-i\omega t},\omega \in R}
poprzez liczenie iloczynów skalarnych danej funkcji z funkcjami bazy.
Transformacja Fouriera została nazwana na cześć Jeana Baptiste’a Josepha Fouriera .
Definicja transformacji Fouriera w postaci zespolonej dla funkcji jednej zmiennej [ edytuj | edytuj kod ]
Podana tu zostanie definicja transformacji Fouriera w postaci zespolonej. Obok tej postaci występują także postacie całkowe transformaty Fouriera, wyrażone przez funkcje rzeczywiste sinus i cosinus, bez użycia jednostki urojonej; postacie te można otrzymać z transformaty Fouriera w postaci zespolonej; nie omawia jednak tego zagadnienia ten artykuł.
Zakładamy, że funkcja
f
{\displaystyle f}
spełnia pierwszy i drugi warunek Dirichleta , tj. jest to funkcja rzeczywista
f
:
R
⟶
R
,
{\displaystyle f\colon \mathbb {R} \longrightarrow \mathbb {R} ,}
określona na skończonym przedziale
(
a
,
b
)
{\displaystyle (a,b)}
oraz[1]
f
(
t
)
{\displaystyle f(t)}
jest przedziałami monotoniczna w przedziale
(
a
,
b
)
{\displaystyle (a,b)}
– tzn. jest w nim ograniczona i można podzielić ten przedział na skończoną liczbę podprzedziałów, wewnątrz których funkcja jest monotoniczna (czyli funkcja ta ma skończoną liczbę maksimów lokalnych i minimów lokalnych ),
f
(
t
)
{\displaystyle f(t)}
jest ciągła w przedziale
(
a
,
b
)
{\displaystyle (a,b)}
z wyjątkiem co najwyżej skończonej liczby punktów nieciągłości pierwszego rodzaju , przy czym w każdym punkcie nieciągłości spełniony jest warunek
f
(
t
0
)
=
1
2
[
f
(
t
0
+
)
+
f
(
t
0
−
)
.
{\displaystyle f(t_{0})={\tfrac {1}{2}}[f(t_{0+})+f(t_{0-}).}
Uwagi:
Całkowalność funkcji, gwarantowana przez spełnienie warunków Dirichleta, jest istotnym wymogiem w obliczaniu transformaty Fouriera, która opiera się na obliczaniu całek. (Ten sam wymóg dotyczy szeregów Fouriera , bowiem współczynniki rozwinięcia w szereg Fouriera są wyrażone za pomocą całek.)
Dla funkcji okresowej warunki Dirichleta wystarczy zbadać na dowolnym przedziale
(
a
,
a
+
T
)
{\displaystyle (a,a+T)}
o długości równej okresowi funkcji
T
.
{\displaystyle T.}
Warunki Dirichleta są to warunki wystarczające do tego, by dla funkcji, która je spełnia, istniała całka Fouriera i szereg Fouriera. Jednak nie są to warunki konieczne – istnieje klasa funkcji, które nie spełniają warunków Dirichleta, a mimo to mają całkę i szereg Fouriera. Także dla funkcji uogólnionych, tzw. dystrybucji (np. delta Diraca ), które mają istotne znaczenie w obliczeniach fizyki cząstek elementarnych , definiuje się transformaty Fouriera.
Transformacja z dziedziny czasu
t
{\displaystyle t}
w dziedzinę pulsacji (częstości kołowej )
ω
{\displaystyle \omega }
ma postać[2]
f
^
(
ω
)
=
∫
−
∞
∞
f
(
t
)
e
−
i
ω
t
d
t
,
{\displaystyle {\hat {f}}(\omega )=\int \limits _{-\infty }^{\infty }f(t)e^{-i\omega t}dt,}
gdzie:
f
(
t
)
{\displaystyle f(t)}
– funkcja (oryginał ) w dziedzinie czasu,
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
– transformata Fouriera (widmo Fouriera ) w dziedzinie pulsacji,
i
{\displaystyle i}
– jednostka urojona .
Transformacja odwrotna zadana jest całką[2]
f
(
t
)
=
1
2
π
∫
−
∞
∞
f
^
(
ω
)
e
i
ω
t
d
ω
{\displaystyle f(t)={\frac {1}{2\pi }}\int \limits _{-\infty }^{\infty }{\hat {f}}(\omega )e^{i\omega t}d\omega }
– pozwala odtworzyć funkcję pierwotną z jej widma.
Transformatę
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
funkcji rzeczywistej
f
(
t
)
{\displaystyle f(t)}
można zapisać następująco[3] :
f
^
(
ω
)
=
|
f
^
(
ω
)
|
⋅
e
i
θ
(
ω
)
,
|
θ
(
ω
)
|
⩽
π
.
{\displaystyle {\hat {f}}(\omega )=|{\hat {f}}(\omega )|\cdot e^{i\theta (\omega )}\,,\qquad |\theta (\omega )|\leqslant \pi .}
Używa się następujących nazw:
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
– charakterystyka widmową (gęstością widma, krótko: widmo ) funkcji
f
(
t
)
,
{\displaystyle f(t),}
|
f
^
(
ω
)
|
{\displaystyle |{\hat {f}}(\omega )|}
– charakterystyka amplitudowa (widmo amplitudowe ) funkcji
f
(
t
)
,
{\displaystyle f(t),}
θ
(
ω
)
{\displaystyle \theta (\omega )}
– charakterystyka fazowa (widmo fazowe ) funkcji
f
(
t
)
.
{\displaystyle f(t).}
Widmo amplitudowe jest funkcją rzeczywistą; można je obliczyć ze wzoru:
|
f
^
(
ω
)
|
=
R
e
(
f
^
(
ω
)
)
2
+
I
m
(
f
^
(
ω
)
)
2
.
{\displaystyle |{\hat {f}}(\omega )|={\sqrt {Re({\hat {f}}(\omega ))^{2}+Im({\hat {f}}(\omega ))^{2}}}.}
Widmo fazowe jest funkcją rzeczywistą; można je obliczyć ze wzorów
cos
θ
(
ω
)
=
R
e
(
f
^
(
ω
)
R
e
(
f
^
(
ω
)
)
2
+
I
m
(
f
^
(
ω
)
)
2
,
{\displaystyle \cos \theta (\omega )={\frac {Re({\hat {f}}(\omega )}{\sqrt {Re({\hat {f}}(\omega ))^{2}+Im({\hat {f}}(\omega ))^{2}}}},}
sin
θ
(
ω
)
=
−
I
m
(
f
^
(
ω
)
R
e
(
f
^
(
ω
)
)
2
+
I
m
(
f
^
(
ω
)
)
2
.
{\displaystyle \sin \theta (\omega )=-\,{\frac {Im({\hat {f}}(\omega )}{\sqrt {Re({\hat {f}}(\omega ))^{2}+Im({\hat {f}}(\omega ))^{2}}}}.}
Stąd:
θ
(
ω
)
=
arctg
[
sin
θ
(
ω
)
cos
θ
(
ω
)
]
.
{\displaystyle \theta (\omega )=\operatorname {arctg} {\Big [}{\tfrac {\sin \theta (\omega )}{\cos \theta (\omega )}}{\Big ]}.}
Funkcję
e
i
ω
t
=
cos
ω
t
+
i
sin
ω
t
{\displaystyle e^{i\omega t}=\cos \omega t+i\sin \omega t}
w płaszczyźnie zespolonej reprezentuje wektor o końcu w punkcie
(
cos
ω
t
,
sin
ω
t
)
,
{\displaystyle (\cos \omega t,\sin \omega t),}
który obraca się wokół początku układu współrzędnych przeciwnie do ruchu wskazówek zegara, gdy rośnie wartość zmiennej
t
.
{\displaystyle t.}
Pokazano też rzut tego wektora na oś rzeczywistą, kreślącego cosinusoidę
y
(
t
)
=
cos
ω
t
.
{\displaystyle y(t)=\cos \omega t.}
Sens transformaty Fouriera jest następujący:
Transformata
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
jest liczbą zespoloną , której moduł i faza określają amplitudę i fazę składowej o pulsacji
ω
,
{\displaystyle \omega ,}
wchodzącej w skład sygnału
f
(
t
)
.
{\displaystyle f(t).}
Obliczenie transformaty Fouriera pozwala więc znaleźć amplitudy i fazy wszystkich takich składowych w sygnale.
Przykład obliczeń wraz w wykresami funkcji, opisujących charakterystyki widmową, amplitudową i fazową podano dalej.
Obracający się przeciwnie do ruchu wskazówek zegara wektor
(
cos
t
,
sin
t
)
{\displaystyle (\cos t,\sin t)}
ma dodatnią częstotliwość
ω
=
+
1
{\displaystyle \omega =+1}
rad /s. Nie pokazano wektora
(
cos
(
−
t
)
,
sin
(
−
t
)
)
{\displaystyle (\cos(-t),\sin(-t))}
obracającego się zgodnie z ruchem wskazówek zegara, który ma ujemną częstotliwość
ω
=
−
1
{\displaystyle \omega =-1}
rad/s. Końce obu wektorów obracają się po okręgu jednostkowym, ale w przeciwnych kierunkach.
Unitarną transformacją Fouriera z dziedziny czasu
t
{\displaystyle t}
w dziedzinę pulsacji (częstości kołowej)
ω
{\displaystyle \omega }
nazywa się transformacją zadaną wzorem, który różni się od wzoru wyżej podanego jedynie stałą, stojącą przed znakiem całki, tj.[4]
f
^
(
ω
)
=
1
2
π
∫
−
∞
∞
f
(
t
)
e
−
i
ω
t
d
t
.
{\displaystyle {\hat {f}}(\omega )={\frac {1}{\sqrt {2\pi }}}\int \limits _{-\infty }^{\infty }f(t)e^{-i\omega t}dt.}
Transformacja odwrotna unitarna
f
(
t
)
=
1
2
π
∫
−
∞
∞
f
^
(
ω
)
e
i
ω
t
d
ω
.
{\displaystyle f(t)={\frac {1}{\sqrt {2\pi }}}\int \limits _{-\infty }^{\infty }{\hat {f}}(\omega )e^{i\omega t}d\omega .}
Czynnik
1
2
π
{\displaystyle {\tfrac {1}{\sqrt {2\pi }}}}
przed transformatą i transformatą odwrotną występuje umownie; zamiast takiej postaci może występować czynnik
1
2
π
{\displaystyle {\tfrac {1}{2\pi }}}
przed transformatą prostą, albo (częściej) przed transformacją odwrotną.
Jeżeli jednak czynnik wynosi
1
2
π
,
{\displaystyle {\tfrac {1}{\sqrt {2\pi }}},}
wtedy transformacja i transformacja odwrotna są izometriami przestrzeni
L
2
(
R
)
.
{\displaystyle L^{2}(\mathbb {R} ).}
Pierwsza z dwu powyższych definicji jest popularniejsza, nie posiada jednak własności unitarności .
W definicji transformaty Fouriera w postaci zespolonej mamy częstotliwości
ω
{\displaystyle \omega }
w zakresie od
−
∞
{\displaystyle -\infty }
do
+
∞
,
{\displaystyle +\infty ,}
występujące w funkcji
f
ω
=
e
−
i
ω
t
.
{\displaystyle f_{\omega }=e^{-i\omega t}.}
W klasycznym rozumieniu teorii dotyczącej drgań i fal częstości kołowe
ω
{\displaystyle \omega }
są wielkościami dodatnimi. Co oznaczają więc ujemne wartości częstotliwości? Sytuacja jest tu nieco inna niż w opisie zjawisk drgań i fal za pomocą liczb rzeczywistych, gdzie definiuje się pojęcie częstotliwości. Tu mamy do czynienia z funkcjami
f
ω
{\displaystyle f_{\omega }}
zespolonymi o częstotliwościach
ω
{\displaystyle \omega }
z zakresu liczb rzeczywistych. Każda funkcja może być przedstawiona jako obracający się w płaszczyźnie zespolonej wektor (wielkość obrotu zależy od czasu
t
{\displaystyle t}
występującego w każdej takiej funkcji), którego koniec porusza się po okręgu o promieniu równym 1 (tzw. okrąg jednostkowy).
Dla dodatnich wartości
ω
{\displaystyle \omega }
wektory
f
ω
{\displaystyle f_{\omega }}
obracają się przeciwnie do wskazówek zegara, dla ujemnych wartości
ω
{\displaystyle \omega }
wektory
f
ω
{\displaystyle f_{\omega }}
obracają się zgodnie ze wskazówkami zegara.
Piękna ilustrację transformacji Fouriera, w tym sens ujemnych częstotliwości, przedstawia animowane omówienie pt. But what is the Fourier Transform? A visual introduction
Transformatę Fouriera oznacza się także symbolem[5]
F
(
j
ω
)
≡
f
^
(
ω
)
=
∫
−
∞
∞
f
(
t
)
e
−
j
ω
t
d
t
,
{\displaystyle F(j\omega )\equiv {\hat {f}}(\omega )=\int \limits _{-\infty }^{\infty }f(t)e^{-j\omega t}dt,}
gdzie
j
≡
i
{\displaystyle j\equiv i}
– jednostka urojona (oznaczenie stosowane często w elektrotechnice, teorii przetwarzania sygnałów).
Oznaczenie to podkreśla, iż mamy do czynienia z transformatą Fouriera w postaci zespolonej, by odróżnić ją od transformat Fouriera kosinusowej i sinusowej. Używa się także symboliki
F
−
1
[
F
(
j
ω
)
]
{\displaystyle {\mathcal {F}}^{-1}[F(j\omega )]}
na oznaczenie transformacji odwrotnej; wtedy mamy:
f
(
t
)
=
F
−
1
[
F
(
j
ω
)
]
.
{\displaystyle f(t)={\mathcal {F}}^{-1}[F(j\omega )].}
Transformata Fouriera funkcji
f
(
t
)
=
t
,
|
t
|
<
1
:
{\displaystyle f(t)=t,|t|<1{:}}
jej część rzeczywista, urojona, moduł
Mamy daną funkcję
f
(
t
)
=
{
t
dla
|
t
|
<
1
1
2
⋅
sgn
(
t
)
dla
|
t
|
=
1
0
dla
|
t
|
>
1
{\displaystyle f(t)=\left\{{\begin{aligned}&t&&{\text{dla}}~~|t|<1\\&{\tfrac {1}{2}}\cdot {\text{ sgn}}(t)&&{\text{dla}}~~|t|=1\\&0&&{\text{dla}}~~|t|>1\end{aligned}}\right.}
Obliczamy
f
^
(
ω
)
=
∫
−
∞
∞
f
(
t
)
e
−
i
ω
t
d
t
,
{\displaystyle {\hat {f}}(\omega )=\int \limits _{-\infty }^{\infty }f(t)e^{-i\omega t}dt,}
czyli
f
^
(
ω
)
=
∫
−
1
1
t
e
−
i
ω
t
d
t
=
∫
−
1
1
t
cos
ω
t
d
t
+
i
∫
−
1
1
t
sin
ω
t
d
t
.
{\displaystyle {\hat {f}}(\omega )=\int \limits _{-1}^{1}te^{-i\omega t}dt=\int \limits _{-1}^{1}t\cos \omega t\,dt+i\int \limits _{-1}^{1}t\sin \omega t\,dt.}
Część rzeczywista obliczanej transformaty jest równa zeru z uwagi na nieparzystość funkcji
t
cos
ω
t
{\displaystyle t\cos \omega t}
, natomiast część zespolona ma wartość
∫
−
1
1
t
sin
ω
t
d
t
=
2
∫
0
1
t
sin
ω
t
d
t
=
−
2
t
cos
ω
t
ω
|
0
1
+
2
ω
∫
0
1
cos
ω
t
d
t
=
−
2
cos
ω
ω
+
2
sin
ω
ω
2
.
{\displaystyle \int \limits _{-1}^{1}t\sin \omega t\,dt=2\int \limits _{0}^{1}t\sin \omega t\,dt=-2t{\frac {\cos \omega t}{\omega }}{\Bigg |}_{0}^{1}+{\frac {2}{\omega }}\int \limits _{0}^{1}\cos \omega t\,dt=-2{\frac {\cos \omega }{\omega }}+2{\frac {\sin \omega }{\omega ^{2}}}.}
Stąd
f
^
(
ω
)
=
0
+
i
2
ω
2
(
ω
cos
ω
−
sin
ω
)
.
{\displaystyle {\hat {f}}(\omega )=0+i{\frac {2}{\omega ^{2}}}(\omega \cos \omega -\sin \omega ).}
Widmo amplitudowe czyli moduł transformaty
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
:
|
f
^
(
ω
)
|
=
R
e
(
f
^
(
ω
)
)
2
+
I
m
(
f
^
(
ω
)
)
2
=
|
2
ω
2
(
ω
cos
ω
−
sin
ω
)
|
.
{\displaystyle |{\hat {f}}(\omega )|={\sqrt {Re({\hat {f}}(\omega ))^{2}+Im({\hat {f}}(\omega ))^{2}}}={\Bigg |}{\frac {2}{\omega ^{2}}}(\omega \cos \omega -\sin \omega ){\Bigg |}.}
Wyniki obliczeń pokazano na wykresach.
Obliczanie transformaty Fouriera w ogólnym przypadku nie jest proste i wymaga na ogół całkowania numerycznego. Poniżej podano kod programu w języku Python , który liczy zespoloną transformatę Fouriera wg wzoru całkowego i rysuje wykresy jej części rzeczywistej, urojonej oraz modułu.
Użytkownik może ustalić:
Definicję funkcji
f
(
x
)
{\displaystyle f(x)}
(linia nr 15).
Liczbę
N
{\displaystyle N}
określającą zakres częstotliwości
ω
∈<
−
N
,
N
>
{\displaystyle \omega \in <-N,N>}
(linia nr 32).
Program można testować, korzystając np. z darmowego notatnika colab google online.
import numpy as np
from scipy.integrate import quad
import matplotlib.pyplot as plt
''' Definicja funkcji f(t)=t
def f(t):
if np.abs(t) < 1:
return t
elif np.abs(t) == 1:
return 0.5 * np.sign(t)
else:
return 0
'''
# Definicja funkcji f(t) - impuls prostokątny
def f ( t ):
if abs ( t - a ) < 1 :
return 1
elif abs ( t - a ) == 1 :
return 0.5 * np . sign ( t )
else :
return 0
# Definicja transformacji Fouriera
def fourier_transform ( omega ):
integrand = lambda t : f ( t ) * np . exp ( - 1 j * omega * t )
real_part = quad ( lambda t : np . real ( integrand ( t )), - np . inf , np . inf )[ 0 ]
imag_part = quad ( lambda t : np . imag ( integrand ( t )), - np . inf , np . inf )[ 0 ]
return real_part + 1 j * imag_part
# Zakres częstotliwości omega
N = 10
omega_values = np . linspace ( - N , N , 500 )
fourier_values = np . array ([ fourier_transform ( omega ) for omega in omega_values ])
# Wykresy: Re, Im oraz |F(ω)|
plt . figure ( figsize = ( 14 , 20 ))
# Wykres części rzeczywistej
plt . subplot ( 3 , 1 , 1 )
plt . plot ( omega_values , fourier_values . real , label = r "Re$\{\hat {f} (\omega)\}$" , color = 'blue' )
plt . title ( "Część rzeczywista transformaty Fouriera" , fontsize = 30 )
plt . xlabel ( r "$\omega$" , fontsize = 20 )
plt . ylabel ( r "Re$\{\hat {f} (\omega)\}$" , fontsize = 14 )
plt . grid ( True )
plt . legend ( fontsize = 20 )
# Wykres części urojonej
plt . subplot ( 3 , 1 , 2 )
plt . plot ( omega_values , fourier_values . imag , label = r "Im$\{\hat {f} (\omega)\}$" , color = 'red' )
plt . title ( "Część urojona transformaty Fouriera" , fontsize = 30 )
plt . xlabel ( r "$\omega$" , fontsize = 20 )
plt . ylabel ( r "Im$\{\hat {f} (\omega)\}$" , fontsize = 14 )
plt . grid ( True )
plt . legend ( fontsize = 20 )
# Wykres modułu |F(ω)|
plt . subplot ( 3 , 1 , 3 )
plt . plot ( omega_values , np . abs ( fourier_values ), label = r "$|\hat {f} (\omega)|$" , color = 'green' )
plt . title ( "Moduł transformaty Fouriera" , fontsize = 30 )
plt . xlabel ( r "$\omega$" , fontsize = 20 )
plt . ylabel ( r "$|\hat {f} (\omega)|$" , fontsize = 14 )
plt . grid ( True )
plt . legend ( fontsize = 20 )
plt . tight_layout ()
plt . show ()
Obliczanie transformaty Fouriera zadane jest de facto poprzez całki niewłaściwe , gdyż zmienna
t
{\displaystyle t}
jest określona w granicach od
−
∞
{\displaystyle -\infty }
do
+
∞
.
{\displaystyle +\infty .}
Obliczanie takich całek sprowadza się do liczenia tzw. wartości głównej całki (wg definicji podanej przez Cauchy’ego)[2] , tj. do liczenia jej jako granicy z całek obliczonych na skończonym przedziale całkowania:
f
^
(
ω
)
=
lim
T
→
+
∞
∫
−
T
T
f
(
t
)
e
−
i
t
ω
d
t
.
{\displaystyle {\hat {f}}(\omega )=\lim _{T\to +\infty }\int \limits _{-T}^{T}f(t)\ e^{-it\omega }\,dt.}
Rzeczywiste sygnały zawsze mają skończony czas trwania. Natomiast można je teoretycznie przedłużać do nieskończoności – wtedy ma zastosowanie powyższy wzór.
Także funkcja pierwotna rekonstruowana z jej transformaty
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
poprzez transformację odwrotną jest zadana przez całkę niewłaściwą, gdyż zmienna
ω
{\displaystyle \omega }
jest określona w granicach od
−
∞
{\displaystyle -\infty }
do
+
∞
.
{\displaystyle +\infty .}
I także obliczanie tej całki sprowadza się do liczenia wartości głównej całki[2] , tj. do liczenia jej jako granicy z całek obliczonych na skończonym przedziale całkowania wg wzoru:
f
(
t
)
=
1
2
π
lim
ω
→
+
∞
∫
−
ω
ω
f
^
(
u
)
e
i
t
u
d
u
.
{\displaystyle f(t)={\frac {1}{2\pi }}\lim _{\omega \to +\infty }\int \limits _{-\omega }^{\omega }{\hat {f}}(u)\ e^{itu}\,du.}
W całce Fouriera funkcje harmoniczne
e
−
i
ω
t
=
cos
ω
t
−
i
sin
ω
t
,
ω
∈
R
{\displaystyle e^{-i\omega t}=\cos \omega t-i\sin \omega t,\omega \in R}
mnożone są przez sygnał
f
(
t
)
;
{\displaystyle f(t);}
wynikowa całka dostarcza informacji nt. zawartości poszczególnych harmonicznych, wchodzących w skład sygnału (dokonuje rozkładu sygnału na jego widmo ).
Transformacja Laplace’a wykonuje podobne działanie, ale o bardziej ogólnym charakterze: Za pomocą całki Laplace’a funkcja rzeczywista
f
(
t
)
,
{\displaystyle f(t),}
reprezentująca zmieniający się w czasie
t
{\displaystyle t}
sygnał (np. pole elektryczne fali, przychodzącej do odbiornika), jest transformowana na płaszczyznę S (płaszczyznę zespoloną); dokonuje się to poprzez scałkowanie iloczynu funkcji
f
(
t
)
{\displaystyle f(t)}
z wyrażeniami typu
e
−
s
t
{\displaystyle e^{-st}}
dla czasu
t
{\displaystyle t}
od
−
∞
{\displaystyle -\infty }
do
+
∞
:
{\displaystyle +\infty {:}}
∫
−
∞
∞
f
(
t
)
e
−
s
t
d
t
;
{\displaystyle \int _{-\infty }^{\infty }f(t)e^{-st}\,dt;}
przy tym
s
=
x
+
i
ω
{\displaystyle s=x+i\omega }
jest liczbą zespoloną , stałą w procesie obliczania całki.
Funkcje
e
−
s
t
{\displaystyle e^{-st}}
w całce Laplace’a de facto mają postać
e
−
s
t
=
e
−
x
t
⋅
e
−
i
ω
t
.
{\displaystyle e^{-st}=e^{-xt}\cdot e^{-i\omega t}.}
Dzięki temu pozwalają dokonać nie tylko analizy zawartości harmonicznych
e
−
i
ω
t
{\displaystyle e^{-i\omega t}}
w sygnale
f
(
t
)
,
{\displaystyle f(t),}
ale również efektów zaniku sygnału lub jego wzrastania w czasie, poprzez funkcję
e
−
x
t
.
{\displaystyle e^{-xt}.}
Na przykład krzywa sinusoidalna tłumiona może być modelowana za pomocą transformacji Laplace’a. Transformacja Fouriera stanowi więc szczególny przypadek przekształcenia Laplace’a dla
s
=
i
ω
.
{\displaystyle s=i\omega .}
Podobnie uogólnieniem dyskretnej transformaty Fouriera stanowi transformata Z , z którą powiązana jest transformata Laplace’a (zob. metoda Tustina ).
W tym rozdziale wprowadzono trzeci sposób definiowania transformaty Fouriera za pomocą częstotliwości
ν
,
{\displaystyle \nu ,}
przy czym zachodzi związek
ω
=
2
π
ν
.
{\displaystyle \omega ={2\pi }\nu .}
W przypadku jednowymiarowym funkcja
f
{\displaystyle f}
jest klasy
L
1
,
{\displaystyle L^{1},}
czyli jest całkowalna w przedziale
(
−
∞
,
∞
)
.
{\displaystyle (-\infty ,\infty ).}
f
^
{\displaystyle {\hat {f}}}
jest funkcją ciągłą . Nie musi natomiast być całkowalna w
R
n
.
{\displaystyle \mathbb {R} ^{n}.}
Jeśli
g
(
x
)
=
f
(
x
−
α
)
,
{\displaystyle g(x)=f(x-\alpha ),}
to
g
^
(
ν
)
=
∫
−
∞
∞
g
(
x
)
e
−
2
π
i
x
ν
d
x
=
∫
−
∞
∞
f
(
x
−
α
)
e
−
2
π
i
(
x
−
α
)
ν
e
−
2
π
i
α
ν
d
x
=
e
−
2
π
i
α
ν
∫
−
∞
∞
f
(
t
)
e
−
2
π
i
t
ν
d
t
=
e
−
2
π
i
α
ν
f
^
(
ν
)
.
{\displaystyle {\hat {g}}(\nu )=\int \limits _{-\infty }^{\infty }g(x)\ e^{-2\pi ix\nu }\,dx=\int \limits _{-\infty }^{\infty }f(x-\alpha )\ e^{-2\pi i(x-\alpha )\nu }e^{-2\pi i\alpha \nu }\,dx=e^{-2\pi i\alpha \nu }\int \limits _{-\infty }^{\infty }f(t)\ e^{-2\pi it\nu }\,dt=e^{-2\pi i\alpha \nu }{\hat {f}}(\nu ).}
Jeśli
α
≠
0
{\displaystyle \alpha \neq 0}
i
g
(
t
)
=
f
(
t
α
)
,
{\displaystyle g(t)=f\left({\frac {t}{\alpha }}\right),}
to
g
^
(
ν
)
=
∫
−
∞
∞
g
(
x
)
e
−
2
π
i
x
ν
d
x
=
∫
−
∞
∞
f
(
x
α
)
e
−
2
π
i
x
α
(
α
ν
)
α
d
(
x
α
)
=
α
∫
−
∞
∞
f
(
t
)
e
−
2
π
i
t
(
α
ν
)
d
t
=
α
f
^
(
α
ν
)
.
{\displaystyle {\hat {g}}(\nu )=\int \limits _{-\infty }^{\infty }g(x)\ e^{-2\pi ix\nu }\,dx=\int \limits _{-\infty }^{\infty }f\left({\frac {x}{\alpha }}\right)\ e^{-2\pi i{\frac {x}{\alpha }}(\alpha \nu )}\alpha \,d\left({\frac {x}{\alpha }}\right)=\alpha \int \limits _{-\infty }^{\infty }f(t)\ e^{-2\pi it(\alpha \nu )}\,dt=\alpha {\hat {f}}(\alpha \nu ).}
f
∗
g
^
=
2
π
f
^
g
^
,
{\displaystyle {\widehat {f*g}}={\sqrt {2\pi }}{\hat {f}}{\hat {g}},}
gdzie operacja
∗
{\displaystyle *}
oznacza splot funkcji f i g
Jeśli pochodna funkcji
f
{\displaystyle f}
należy do
L
1
{\displaystyle L^{1}}
i funkcja zeruje się poza pewnym przedziałem skończonym, to z całkowania przez części wynika, że:
f
′
^
(
ν
)
=
∫
−
∞
∞
f
′
(
x
)
e
−
2
π
i
x
ν
d
x
=
f
(
x
)
e
−
2
π
i
x
ν
|
−
∞
∞
+
2
π
i
ν
∫
−
∞
∞
f
(
x
)
e
−
2
π
i
x
ν
d
x
=
2
π
i
ν
f
^
(
ν
)
.
{\displaystyle {\hat {f'}}(\nu )=\int \limits _{-\infty }^{\infty }f'(x)e^{-2\pi ix\nu }\,dx=f(x)e^{-2\pi ix\nu }{\Big |}_{-\infty }^{\infty }+2\pi i\nu \int \limits _{-\infty }^{\infty }f(x)e^{-2\pi ix\nu }\,dx=2\pi i\nu {\hat {f}}(\nu ).}
Funkcja
Transformata Fouriera unitarna, częstotliwość
Transformata Fouriera unitarna, pulsacja (częstość kołowa)
Transformata Fouriera pulsacja (częstość kołowa)
Uwagi
f
(
x
)
{\displaystyle f(x)}
f
^
(
ν
)
=
{\displaystyle {\hat {f}}(\nu )=}
∫
−
∞
∞
f
(
x
)
e
−
2
π
i
x
ν
d
x
{\displaystyle {}\,\int \limits _{-\infty }^{\infty }f(x)e^{-2\pi ix\nu }\,dx}
f
^
(
ω
)
=
{\displaystyle {\hat {f}}(\omega )=}
1
2
π
∫
−
∞
∞
f
(
x
)
e
−
i
ω
x
d
x
{\displaystyle {\frac {1}{\sqrt {2\pi }}}\int \limits _{-\infty }^{\infty }f(x)e^{-i\omega x}\,dx}
f
^
(
ω
)
=
{\displaystyle {\hat {f}}(\omega )=}
∫
−
∞
∞
f
(
x
)
e
−
i
ω
x
d
x
{\displaystyle {}\,\int \limits _{-\infty }^{\infty }f(x)e^{-i\omega x}\,dx}
101
a
⋅
f
(
x
)
+
b
⋅
g
(
x
)
{\displaystyle a\cdot f(x)+b\cdot g(x)}
a
⋅
f
^
(
ν
)
+
b
⋅
g
^
(
ν
)
{\displaystyle a\cdot {\hat {f}}(\nu )+b\cdot {\hat {g}}(\nu )}
a
⋅
f
^
(
ω
)
+
b
⋅
g
^
(
ω
)
{\displaystyle a\cdot {\hat {f}}(\omega )+b\cdot {\hat {g}}(\omega )}
a
⋅
f
^
(
ω
)
+
b
⋅
g
^
(
ω
)
{\displaystyle a\cdot {\hat {f}}(\omega )+b\cdot {\hat {g}}(\omega )}
Liniowość
102
f
(
x
−
a
)
{\displaystyle f(x-a)}
e
−
2
π
i
a
ν
f
^
(
ν
)
{\displaystyle e^{-2\pi ia\nu }{\hat {f}}(\nu )}
e
−
i
a
ω
f
^
(
ω
)
{\displaystyle e^{-ia\omega }{\hat {f}}(\omega )}
e
−
i
a
ω
f
^
(
ω
)
{\displaystyle e^{-ia\omega }{\hat {f}}(\omega )}
Przesunięcie oryginału w dziedzinie „czasu”
103
e
2
π
i
a
x
f
(
x
)
{\displaystyle e^{2\pi iax}f(x)}
f
^
(
ν
−
a
)
{\displaystyle {\hat {f}}\left(\nu -a\right)}
f
^
(
ω
−
2
π
a
)
{\displaystyle {\hat {f}}(\omega -2\pi a)}
f
^
(
ω
−
2
π
a
)
{\displaystyle {\hat {f}}(\omega -2\pi a)}
Przesunięcie widma w dziedzinie częstotliwości, dualne względem 102
104
f
(
a
x
)
{\displaystyle f(ax)}
1
|
a
|
f
^
(
ν
a
)
{\displaystyle {\frac {1}{|a|}}{\hat {f}}\left({\frac {\nu }{a}}\right)}
1
|
a
|
f
^
(
ω
a
)
{\displaystyle {\frac {1}{|a|}}{\hat {f}}\left({\frac {\omega }{a}}\right)}
1
|
a
|
f
^
(
ω
a
)
{\displaystyle {\frac {1}{|a|}}{\hat {f}}\left({\frac {\omega }{a}}\right)}
Dla dużych wartości
|
a
|
,
{\displaystyle |a|,}
f
(
a
x
)
{\displaystyle f(ax)}
zawęża się wokół zera, a
1
|
a
|
f
^
(
ω
a
)
{\displaystyle {\frac {1}{|a|}}{\hat {f}}\left({\frac {\omega }{a}}\right)}
poszerza się i spłaszcza.
105
106
d
n
f
(
x
)
d
x
n
{\displaystyle {\frac {d^{n}f(x)}{dx^{n}}}}
(
2
π
i
ν
)
n
f
^
(
ν
)
{\displaystyle (2\pi i\nu )^{n}{\hat {f}}(\nu )}
(
i
ω
)
n
f
^
(
ω
)
{\displaystyle (i\omega )^{n}{\hat {f}}(\omega )}
(
i
ω
)
n
f
^
(
ω
)
{\displaystyle (i\omega )^{n}{\hat {f}}(\omega )}
Transformata pochodnej
107
x
n
f
(
x
)
{\displaystyle x^{n}f(x)}
(
i
2
π
)
n
d
n
f
^
(
ν
)
d
ν
n
{\displaystyle \left({\frac {i}{2\pi }}\right)^{n}{\frac {d^{n}{\hat {f}}(\nu )}{d\nu ^{n}}}}
i
n
d
n
f
^
(
ω
)
d
ω
n
{\displaystyle i^{n}{\frac {d^{n}{\hat {f}}(\omega )}{d\omega ^{n}}}}
i
n
d
n
f
^
(
ω
)
d
ω
n
{\displaystyle i^{n}{\frac {d^{n}{\hat {f}}(\omega )}{d\omega ^{n}}}}
Ta właściwość jest dualna względem 106
108
(
f
∗
g
)
(
x
)
{\displaystyle (f*g)(x)}
f
^
(
ν
)
g
^
(
ν
)
{\displaystyle {\hat {f}}(\nu ){\hat {g}}(\nu )}
2
π
f
^
(
ω
)
g
^
(
ω
)
{\displaystyle {\sqrt {2\pi }}{\hat {f}}(\omega ){\hat {g}}(\omega )}
f
^
(
ω
)
g
^
(
ω
)
{\displaystyle {\hat {f}}(\omega ){\hat {g}}(\omega )}
Notacja
f
∗
g
{\displaystyle f*g}
oznacza splot funkcji
f
{\displaystyle f}
i
g
{\displaystyle g}
– tę właściwość określamy jako twierdzenie o splocie
109
f
(
x
)
g
(
x
)
{\displaystyle f(x)g(x)}
(
f
^
∗
g
^
)
(
ν
)
{\displaystyle ({\hat {f}}*{\hat {g}})(\nu )}
(
f
^
∗
g
^
)
(
ω
)
2
π
{\displaystyle {\frac {({\hat {f}}*{\hat {g}})(\omega )}{\sqrt {2\pi }}}}
1
2
π
(
f
^
∗
g
^
)
(
ω
)
{\displaystyle {\frac {1}{2\pi }}({\hat {f}}*{\hat {g}})(\omega )}
Właściwość dualna względem 108
110
Dla funkcji
f
(
x
)
{\displaystyle f(x)}
rzeczywistej i parzystej
f
^
(
ν
)
,
{\displaystyle {\hat {f}}(\nu ),}
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
w postaci unitarnej oraz
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
są funkcjami rzeczywistymi i parzystymi .
111
Dla funkcji
f
(
x
)
{\displaystyle f(x)}
rzeczywistej i nieparzystej
f
^
(
ν
)
,
{\displaystyle {\hat {f}}(\nu ),}
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
w postaci unitarnej oraz
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
są funkcjami urojonymi i nieparzystymi .
W przestrzeniach funkcji dobrze określonych, takich jak
L
2
(
R
)
{\displaystyle L^{2}(\mathbb {R} )}
lub przestrzeń funkcji szybko malejących w nieskończoności
S
(
R
)
{\displaystyle {\mathcal {S}}(\mathbb {R} )}
transformacja Fouriera przyporządkowuje wzajemnie jednoznacznie funkcji (oryginałowi) jej transformatę. Oryginał i jego transformata określana są wtedy jako para fourierowska.
W tabeli zestawiony jest oryginał oraz jego transformaty w dziedzinie częstotliwości i pulsacji[6] .
Uwaga: Tu
x
{\displaystyle x}
oznacza zmienną funkcji pierwotnej
f
,
{\displaystyle f,}
zamiast oznaczenia
t
.
{\displaystyle t.}
Funkcja
Transformata Fouriera unitarna, częstotliwość
Transformata Fouriera unitarna, pulsacja (częstość kołowa)
Transformata Fouriera pulsacja (częstość kołowa)
Uwagi
f
(
x
)
{\displaystyle f(x)}
f
^
(
ν
)
=
{\displaystyle {\hat {f}}(\nu )=}
∫
−
∞
∞
f
(
x
)
e
−
2
π
i
x
ν
d
x
{\displaystyle {}\,\int \limits _{-\infty }^{\infty }f(x)e^{-2\pi ix\nu }\,dx}
f
^
(
ω
)
=
{\displaystyle {\hat {f}}(\omega )=}
1
2
π
∫
−
∞
∞
f
(
x
)
e
−
i
ω
x
d
x
{\displaystyle {\frac {1}{\sqrt {2\pi }}}\int \limits _{-\infty }^{\infty }f(x)e^{-i\omega x}\,dx}
f
^
(
ω
)
=
{\displaystyle {\hat {f}}(\omega )=}
∫
−
∞
∞
f
(
x
)
e
−
i
ω
x
d
x
{\displaystyle {}\,\int \limits _{-\infty }^{\infty }f(x)e^{-i\omega x}\,dx}
201
rect
(
a
x
)
{\displaystyle \operatorname {rect} (ax)}
1
|
a
|
⋅
sinc
(
ν
a
)
{\displaystyle {\frac {1}{|a|}}\cdot \operatorname {sinc} \left({\frac {\nu }{a}}\right)}
1
2
π
|
a
|
⋅
sinc
(
ω
2
π
a
)
{\displaystyle {\frac {1}{{\sqrt {2\pi }}|a|}}\cdot \operatorname {sinc} \left({\frac {\omega }{2\pi a}}\right)}
1
|
a
|
⋅
sinc
(
ω
2
π
a
)
{\displaystyle {\frac {1}{|a|}}\cdot \operatorname {sinc} \left({\frac {\omega }{2\pi a}}\right)}
Funkcja prostokątna i normalizowana funkcja sinc , definiowana jako
sinc
(
x
)
=
s
i
n
(
π
x
)
/
(
π
x
)
{\displaystyle \operatorname {sinc} (x)=sin(\pi x)/(\pi x)}
202
sinc
(
a
x
)
{\displaystyle \operatorname {sinc} (ax)}
1
|
a
|
⋅
rect
(
ν
a
)
{\displaystyle {\frac {1}{|a|}}\cdot \operatorname {rect} \left({\frac {\nu }{a}}\right)}
1
2
π
|
a
|
⋅
rect
(
ω
2
π
a
)
{\displaystyle {\frac {1}{{\sqrt {2\pi }}|a|}}\cdot \operatorname {rect} \left({\frac {\omega }{2\pi a}}\right)}
1
|
a
|
⋅
rect
(
ω
2
π
a
)
{\displaystyle {\frac {1}{|a|}}\cdot \operatorname {rect} \left({\frac {\omega }{2\pi a}}\right)}
Relacja dualna do 201. funkcja prostokątna jest filtrem dolnoprzepustowym , a funkcja sinc jest odpowiedzią impulsową dla takiego filtra.
203
sinc
2
(
a
x
)
{\displaystyle \operatorname {sinc} ^{2}(ax)}
1
|
a
|
⋅
tri
(
ν
a
)
{\displaystyle {\frac {1}{|a|}}\cdot \operatorname {tri} \left({\frac {\nu }{a}}\right)}
1
2
π
|
a
|
⋅
tri
(
ω
2
π
a
)
{\displaystyle {\frac {1}{{\sqrt {2\pi }}|a|}}\cdot \operatorname {tri} \left({\frac {\omega }{2\pi a}}\right)}
1
|
a
|
⋅
tri
(
ω
2
π
a
)
{\displaystyle {\frac {1}{|a|}}\cdot \operatorname {tri} \left({\frac {\omega }{2\pi a}}\right)}
Funkcja
tri
(
x
)
=
λ
(
x
)
{\displaystyle \operatorname {tri} (x)=\lambda (x)}
jest funkcją trójkątną
204
tri
(
a
x
)
{\displaystyle \operatorname {tri} (ax)}
1
|
a
|
⋅
sinc
2
(
ν
a
)
{\displaystyle {\frac {1}{|a|}}\cdot \operatorname {sinc} ^{2}\left({\frac {\nu }{a}}\right)}
1
2
π
|
a
|
⋅
sinc
2
(
ω
2
π
a
)
{\displaystyle {\frac {1}{{\sqrt {2\pi }}|a|}}\cdot \operatorname {sinc} ^{2}\left({\frac {\omega }{2\pi a}}\right)}
1
|
a
|
⋅
sinc
2
(
ω
2
π
a
)
{\displaystyle {\frac {1}{|a|}}\cdot \operatorname {sinc} ^{2}\left({\frac {\omega }{2\pi a}}\right)}
Związek dualny względem 203.
205
e
−
a
x
H
(
x
)
{\displaystyle e^{-ax}H(x)}
1
a
+
2
π
i
ν
{\displaystyle {\frac {1}{a+2\pi i\nu }}}
1
2
π
(
a
+
i
ω
)
{\displaystyle {\frac {1}{{\sqrt {2\pi }}(a+i\omega )}}}
1
a
+
i
ω
{\displaystyle {\frac {1}{a+i\omega }}}
H
(
x
)
{\displaystyle H(x)}
jest funkcją skoku Heaviside’a ,
a
>
0.
{\displaystyle a>0.}
206
e
−
α
x
2
{\displaystyle e^{-\alpha x^{2}}}
π
α
⋅
e
−
(
π
ν
)
2
α
{\displaystyle {\sqrt {\frac {\pi }{\alpha }}}\cdot e^{-{\frac {(\pi \nu )^{2}}{\alpha }}}}
1
2
α
⋅
e
−
ω
2
4
α
{\displaystyle {\frac {1}{\sqrt {2\alpha }}}\cdot e^{-{\frac {\omega ^{2}}{4\alpha }}}}
π
α
⋅
e
−
ω
2
4
α
{\displaystyle {\sqrt {\frac {\pi }{\alpha }}}\cdot e^{-{\frac {\omega ^{2}}{4\alpha }}}}
Funkcja Gaussa
exp
(
−
α
x
2
)
{\displaystyle \exp(-\alpha x^{2})}
jest funkcją własną przekształcenia Fouriera dla odpowiednio dobranego
α
.
{\displaystyle \alpha .}
Funkcja jest całkowalna dla
R
e
(
α
)
>
0.
{\displaystyle \mathrm {Re} (\alpha )>0.}
207
e
−
a
|
x
|
{\displaystyle \operatorname {e} ^{-a|x|}}
2
a
a
2
+
4
π
2
ν
2
{\displaystyle {\frac {2a}{a^{2}+4\pi ^{2}\nu ^{2}}}}
2
π
⋅
a
a
2
+
ω
2
{\displaystyle {\sqrt {\frac {2}{\pi }}}\cdot {\frac {a}{a^{2}+\omega ^{2}}}}
2
a
a
2
+
ω
2
{\displaystyle {\frac {2a}{a^{2}+\omega ^{2}}}}
Dla
a
>
0.
{\displaystyle a>0.}
208
J
n
(
x
)
x
{\displaystyle {\frac {J_{n}(x)}{x}}}
2
i
n
(
−
i
)
n
⋅
U
n
−
1
(
2
π
ν
)
{\displaystyle {\frac {2i}{n}}(-i)^{n}\cdot U_{n-1}(2\pi \nu )}
⋅
1
−
4
π
2
ν
2
rect
(
π
ν
)
{\displaystyle \cdot \ {\sqrt {1-4\pi ^{2}\nu ^{2}}}\operatorname {rect} (\pi \nu )}
2
π
i
n
(
−
i
)
n
⋅
U
n
−
1
(
ω
)
{\displaystyle {\sqrt {\frac {2}{\pi }}}{\frac {i}{n}}(-i)^{n}\cdot U_{n-1}(\omega )}
⋅
1
−
ω
2
rect
(
ω
2
)
{\displaystyle \cdot \ {\sqrt {1-\omega ^{2}}}\operatorname {rect} \left({\frac {\omega }{2}}\right)}
2
i
n
(
−
i
)
n
⋅
U
n
−
1
(
ω
)
{\displaystyle {\frac {2i}{n}}(-i)^{n}\cdot U_{n-1}(\omega )}
⋅
1
−
ω
2
rect
(
ω
2
)
{\displaystyle \cdot \ {\sqrt {1-\omega ^{2}}}\operatorname {rect} \left({\frac {\omega }{2}}\right)}
J
n
(
x
)
{\displaystyle J_{n}(x)}
oznacza funkcję Bessela
n
{\displaystyle n}
-tego rzędu, pierwszego rodzaju.
U
n
(
x
)
{\displaystyle U_{n}(x)}
to wielomiany Czebyszewa drugiego rodzaju . (patrz punkty 315 i 316 poniżej).
209
sech
(
a
x
)
{\displaystyle \operatorname {sech} (ax)}
π
a
sech
(
π
2
a
ν
)
{\displaystyle {\frac {\pi }{a}}\operatorname {sech} \left({\frac {\pi ^{2}}{a}}\nu \right)}
1
a
π
2
sech
(
π
2
a
ω
)
{\displaystyle {\frac {1}{a}}{\sqrt {\frac {\pi }{2}}}\operatorname {sech} \left({\frac {\pi }{2a}}\omega \right)}
π
a
sech
(
π
2
a
ω
)
{\displaystyle {\frac {\pi }{a}}\operatorname {sech} \left({\frac {\pi }{2a}}\omega \right)}
Secans hiperboliczny jest funkcją własną transformcji Fouriera
Dystrybucje, określane też jako funkcje uogólnione nie posiadają transformat w sensie określonym przez powyższe definicje. Możliwe jest jednak uogólnienie pojęcia transformaty i przyjęcie, że uzyskujemy w wyniku przejścia do granicy ciągu transformat lub wychodząc od szeregu Fouriera funkcji okresowej .
Funkcja
Transformata Fouriera unitarna, częstotliwość
Transformata Fouriera unitarna, pulsacja (częstość kołowa)
Transformata Fouriera pulsacja (częstość kołowa)
Uwagi
f
(
x
)
{\displaystyle f(x)}
f
^
(
ν
)
=
{\displaystyle {\hat {f}}(\nu )=}
∫
−
∞
∞
f
(
x
)
e
−
2
π
i
x
ν
d
x
{\displaystyle {}\,\int \limits _{-\infty }^{\infty }f(x)e^{-2\pi ix\nu }\,dx}
f
^
(
ω
)
=
{\displaystyle {\hat {f}}(\omega )=}
1
2
π
∫
−
∞
∞
f
(
x
)
e
−
i
ω
x
d
x
{\displaystyle {\frac {1}{\sqrt {2\pi }}}\int \limits _{-\infty }^{\infty }f(x)e^{-i\omega x}\,dx}
f
^
(
ω
)
=
{\displaystyle {\hat {f}}(\omega )=}
∫
−
∞
∞
f
(
x
)
e
−
i
ω
x
d
x
{\displaystyle {}\,\int \limits _{-\infty }^{\infty }f(x)e^{-i\omega x}\,dx}
301
1
{\displaystyle 1}
δ
(
ν
)
{\displaystyle \delta (\nu )}
2
π
⋅
δ
(
ω
)
{\displaystyle {\sqrt {2\pi }}\cdot \delta (\omega )}
2
π
δ
(
ω
)
{\displaystyle 2\pi \delta (\omega )}
δ
(
ω
)
{\displaystyle \delta (\omega )}
oznacza deltę Diraca .
302
δ
(
x
)
{\displaystyle \delta (x)}
1
{\displaystyle 1}
1
2
π
{\displaystyle {\frac {1}{\sqrt {2\pi }}}}
1
{\displaystyle 1}
Co wynika z zasady 301.
303
e
i
a
x
{\displaystyle e^{iax}}
δ
(
ν
−
a
2
π
)
{\displaystyle \delta \left(\nu -{\frac {a}{2\pi }}\right)}
2
π
⋅
δ
(
ω
−
a
)
{\displaystyle {\sqrt {2\pi }}\cdot \delta (\omega -a)}
2
π
δ
(
ω
−
a
)
{\displaystyle 2\pi \delta (\omega -a)}
Co wynika z własności 103 i 301.
304
cos
(
a
x
)
{\displaystyle \cos(ax)}
δ
(
ν
−
a
2
π
)
+
δ
(
ν
+
a
2
π
)
2
{\displaystyle {\frac {\displaystyle \delta \left(\nu -{\frac {a}{2\pi }}\right)+\delta \left(\nu +{\frac {a}{2\pi }}\right)}{2}}}
2
π
⋅
δ
(
ω
−
a
)
+
δ
(
ω
+
a
)
2
{\displaystyle {\sqrt {2\pi }}\cdot {\frac {\delta (\omega -a)+\delta (\omega +a)}{2}}}
π
(
δ
(
ω
−
a
)
+
δ
(
ω
+
a
)
)
{\displaystyle \pi \left(\delta (\omega -a)+\delta (\omega +a)\right)}
Co wynika ze 101 i 303 przy zastosowaniu wzoru Eulera :
cos
(
a
x
)
=
(
e
i
a
x
+
e
−
i
a
x
)
/
2.
{\displaystyle \cos(ax)=(e^{iax}+e^{-iax})/2.}
305
sin
(
a
x
)
{\displaystyle \sin(ax)}
δ
(
ν
+
a
2
π
)
−
δ
(
ν
−
a
2
π
)
2
i
{\displaystyle {\frac {\displaystyle \delta \left(\nu +{\frac {a}{2\pi }}\right)-\delta \left(\nu -{\frac {a}{2\pi }}\right)}{2i}}}
2
π
⋅
δ
(
ω
+
a
)
−
δ
(
ω
−
a
)
2
i
{\displaystyle {\sqrt {2\pi }}\cdot {\frac {\delta (\omega +a)-\delta (\omega -a)}{2i}}}
−
i
π
(
δ
(
ω
+
a
)
−
δ
(
ω
−
a
)
)
{\displaystyle -i\pi \left(\delta (\omega +a)-\delta (\omega -a)\right)}
Co wynika ze 101 i 303, przy zastosowaniu
sin
(
a
x
)
=
(
e
i
a
x
−
e
−
i
a
x
)
/
(
2
i
)
.
{\displaystyle \sin(ax)=(e^{iax}-e^{-iax})/(2i).}
306
cos
(
a
x
2
)
{\displaystyle \cos(ax^{2})}
π
a
cos
(
π
2
ν
2
a
−
π
4
)
{\displaystyle {\sqrt {\frac {\pi }{a}}}\cos \left({\frac {\pi ^{2}\nu ^{2}}{a}}-{\frac {\pi }{4}}\right)}
1
2
a
cos
(
ω
2
4
a
−
π
4
)
{\displaystyle {\frac {1}{\sqrt {2a}}}\cos \left({\frac {\omega ^{2}}{4a}}-{\frac {\pi }{4}}\right)}
π
a
cos
(
ω
2
4
a
−
π
4
)
{\displaystyle {\sqrt {\frac {\pi }{a}}}\cos \left({\frac {\omega ^{2}}{4a}}-{\frac {\pi }{4}}\right)}
307
sin
(
a
x
2
)
{\displaystyle \sin(ax^{2})}
−
π
a
sin
(
π
2
ν
2
a
−
π
4
)
{\displaystyle -{\sqrt {\frac {\pi }{a}}}\sin \left({\frac {\pi ^{2}\nu ^{2}}{a}}-{\frac {\pi }{4}}\right)}
−
1
2
a
sin
(
ω
2
4
a
−
π
4
)
{\displaystyle {\frac {-1}{\sqrt {2a}}}\sin \left({\frac {\omega ^{2}}{4a}}-{\frac {\pi }{4}}\right)}
−
π
a
sin
(
ω
2
4
a
−
π
4
)
{\displaystyle -{\sqrt {\frac {\pi }{a}}}\sin \left({\frac {\omega ^{2}}{4a}}-{\frac {\pi }{4}}\right)}
308
x
n
{\displaystyle x^{n}}
(
i
2
π
)
n
δ
(
n
)
(
ν
)
{\displaystyle \left({\frac {i}{2\pi }}\right)^{n}\delta ^{(n)}(\nu )}
i
n
2
π
δ
(
n
)
(
ω
)
{\displaystyle i^{n}{\sqrt {2\pi }}\delta ^{(n)}(\omega )}
2
π
i
n
δ
(
n
)
(
ω
)
{\displaystyle 2\pi i^{n}\delta ^{(n)}(\omega )}
Gdzie
n
{\displaystyle n}
jest liczbą naturalną a
δ
(
n
)
(
ξ
)
{\displaystyle \delta ^{(n)}(\xi )}
jest
n
{\displaystyle n}
-tą pochodną delty Diraca. Wynika to ze 107 i 301. Stosując tę formułę z 101, można transformować dowolne wielomiany .
309
1
x
{\displaystyle {\frac {1}{x}}}
−
i
π
sgn
(
ν
)
{\displaystyle -i\pi \operatorname {sgn}(\nu )}
−
i
π
2
sgn
(
ω
)
{\displaystyle -i{\sqrt {\frac {\pi }{2}}}\operatorname {sgn}(\omega )}
−
i
π
sgn
(
ω
)
{\displaystyle -i\pi \operatorname {sgn}(\omega )}
Gdzie
sgn
(
ξ
)
{\displaystyle \operatorname {sgn}(\xi )}
to funkcja znaku . Zauważmy, że
1
/
x
{\displaystyle 1/x}
nie jest dystrybucją. Własność przydatna w odniesieniu do transformaty Hilberta.
310
1
x
n
{\displaystyle {\frac {1}{x^{n}}}}
−
i
π
(
−
2
π
i
ν
)
n
−
1
(
n
−
1
)
!
sgn
(
ν
)
{\displaystyle -i\pi {\frac {(-2\pi i\nu )^{n-1}}{(n-1)!}}\operatorname {sgn}(\nu )}
−
i
π
2
⋅
(
−
i
ω
)
n
−
1
(
n
−
1
)
!
sgn
(
ω
)
{\displaystyle -i{\sqrt {\frac {\pi }{2}}}\cdot {\frac {(-i\omega )^{n-1}}{(n-1)!}}\operatorname {sgn}(\omega )}
−
i
π
(
−
i
ω
)
n
−
1
(
n
−
1
)
!
sgn
(
ω
)
{\displaystyle -i\pi {\frac {(-i\omega )^{n-1}}{(n-1)!}}\operatorname {sgn}(\omega )}
Uogólnienie 309.
311
1
|
x
|
{\displaystyle {\frac {1}{\sqrt {|x|}}}}
1
|
ν
|
{\displaystyle {\frac {1}{\sqrt {|\nu |}}}}
1
|
ω
|
{\displaystyle {\frac {1}{\sqrt {|\omega |}}}}
2
π
|
ω
|
{\displaystyle {\frac {\sqrt {2\pi }}{\sqrt {|\omega |}}}}
312
sgn
(
x
)
{\displaystyle \operatorname {sgn}(x)}
1
i
π
ν
{\displaystyle {\frac {1}{i\pi \nu }}}
2
π
⋅
1
i
ω
{\displaystyle {\sqrt {\frac {2}{\pi }}}\cdot {\frac {1}{i\omega }}}
2
i
ω
{\displaystyle {\frac {2}{i\omega }}}
Dualne do 309.
313
H
(
x
)
{\displaystyle H(x)}
1
2
(
1
i
π
ν
+
δ
(
ν
)
)
{\displaystyle {\frac {1}{2}}\left({\frac {1}{i\pi \nu }}+\delta (\nu )\right)}
π
2
(
1
i
π
ω
+
δ
(
ω
)
)
{\displaystyle {\sqrt {\frac {\pi }{2}}}\left({\frac {1}{i\pi \omega }}+\delta (\omega )\right)}
π
(
1
i
π
ω
+
δ
(
ω
)
)
{\displaystyle \pi \left({\frac {1}{i\pi \omega }}+\delta (\omega )\right)}
Funkcja
H
(
x
)
{\displaystyle H(x)}
jest funkcją Heaviside’a; to wynika ze 101, 301 i 312.
314
∑
n
=
−
∞
∞
δ
(
t
−
n
T
)
{\displaystyle \sum _{n=-\infty }^{\infty }\delta (t-nT)}
1
T
∑
k
=
−
∞
∞
δ
(
ν
−
k
T
)
{\displaystyle {\frac {1}{T}}\sum _{k=-\infty }^{\infty }\delta \left(\nu -{\frac {k}{T}}\right)}
2
π
T
∑
k
=
−
∞
∞
δ
(
ω
−
2
π
k
T
)
{\displaystyle {\frac {\sqrt {2\pi }}{T}}\sum _{k=-\infty }^{\infty }\delta \left(\omega -{\frac {2\pi k}{T}}\right)}
2
π
T
∑
k
=
−
∞
∞
δ
(
ω
−
2
π
k
T
)
{\displaystyle {\frac {2\pi }{T}}\sum _{k=-\infty }^{\infty }\delta \left(\omega -{\frac {2\pi k}{T}}\right)}
Funkcja grzebieniowa . Do wyznaczenia z 302 i 102, wraz z faktem, że
∑
n
=
−
∞
∞
e
i
n
x
=
∑
k
=
−
∞
∞
δ
(
x
+
2
π
k
)
{\displaystyle \sum _{n=-\infty }^{\infty }e^{inx}=\sum _{k=-\infty }^{\infty }\delta (x+2\pi k)}
jako dystrybucje.
315
J
0
(
x
)
{\displaystyle J_{0}(x)}
2
rect
(
π
ν
)
1
−
4
π
2
ν
2
{\displaystyle {\frac {2\,\operatorname {rect} (\pi \nu )}{\sqrt {1-4\pi ^{2}\nu ^{2}}}}}
2
π
⋅
rect
(
ω
2
)
1
−
ω
2
{\displaystyle {\sqrt {\frac {2}{\pi }}}\cdot {\frac {\operatorname {rect} \left(\displaystyle {\frac {\omega }{2}}\right)}{\sqrt {1-\omega ^{2}}}}}
2
rect
(
ω
2
)
1
−
ω
2
{\displaystyle {\frac {2\,\operatorname {rect} \left(\displaystyle {\frac {\omega }{2}}\right)}{\sqrt {1-\omega ^{2}}}}}
J
0
(
x
)
{\displaystyle J_{0}(x)}
funkcją Bessela pierwszego rodzaju, rzędu zerowego.
316
J
n
(
x
)
{\displaystyle J_{n}(x)}
2
(
−
i
)
n
T
n
(
2
π
ν
)
rect
(
π
ν
)
1
−
4
π
2
ν
2
{\displaystyle {\frac {2(-i)^{n}T_{n}(2\pi \nu )\operatorname {rect} (\pi \nu )}{\sqrt {1-4\pi ^{2}\nu ^{2}}}}}
2
π
(
−
i
)
n
T
n
(
ω
)
rect
(
ω
2
)
1
−
ω
2
{\displaystyle {\sqrt {\frac {2}{\pi }}}{\frac {(-i)^{n}T_{n}(\omega )\operatorname {rect} \left(\displaystyle {\frac {\omega }{2}}\right)}{\sqrt {1-\omega ^{2}}}}}
2
(
−
i
)
n
T
n
(
ω
)
rect
(
ω
2
)
1
−
ω
2
{\displaystyle {\frac {2(-i)^{n}T_{n}(\omega )\operatorname {rect} \left(\displaystyle {\frac {\omega }{2}}\right)}{\sqrt {1-\omega ^{2}}}}}
Uogólnienie 315. Funkcja
J
n
(
x
)
{\displaystyle J_{n}(x)}
jest funkcją Bessela
n
{\displaystyle n}
-tego rzędu, pierwszego rodzaju. Funkcja
T
n
(
x
)
{\displaystyle T_{n}(x)}
jest wielomianem Czebyszewa pierwszego rodzaju.
Zależność określającą transmitancję widmową
H
(
j
ω
)
{\displaystyle H(j\omega )}
można wyznaczyć:
Pojedyncza zespolona składowa harmoniczna o amplitudzie
A
w
e
,
{\displaystyle A_{we},}
pulsacji
ω
{\displaystyle \omega }
i fazie
p
w
e
{\displaystyle p_{we}}
x
(
t
)
=
A
w
e
e
j
(
ω
t
+
p
w
e
)
,
{\displaystyle x(t)=A_{we}e^{j(\omega t+p_{we})},}
(gdzie
j
{\displaystyle j}
oznacza jednostkę urojoną ), generuje odpowiedź na wyjściu układu w postaci sygnału sinusoidalnego o amplitudzie
A
w
y
{\displaystyle A_{wy}}
i fazie
p
w
y
:
{\displaystyle p_{wy}{:}}
y
(
t
)
=
A
w
y
e
j
(
ω
t
+
p
w
y
)
.
{\displaystyle y(t)=A_{wy}e^{j(\omega t+p_{wy})}.}
Warto zwrócić uwagę na fakt, że częstotliwość
ω
{\displaystyle \omega }
pozostała taka sama, jedynie amplituda i faza sygnału uległy zmianie w układzie. Transmitancja
H
(
j
ω
)
{\displaystyle H(j\omega )}
opisuje charakter tych zmian w całym spektrum częstotliwości (tj. dla każdej częstotliwości
ω
{\displaystyle \omega }
). Moduł transmitancji opisuje wzmocnienie układu:
A
w
y
A
w
e
=
|
H
(
j
ω
)
|
,
{\displaystyle {\frac {A_{wy}}{A_{we}}}=|H(j\omega )|,}
a argument tej transmitancji opisuje przesunięcie fazowe wprowadzane przez układ:
p
w
y
−
p
w
e
=
arg
(
H
(
j
ω
)
)
.
{\displaystyle p_{wy}-p_{we}=\arg(H(j\omega )).}
Transmitancja
H
(
j
ω
)
{\displaystyle H(j\omega )}
Dla przypadku układów dyskretnych wyrażenie
y
(
i
)
u
(
i
)
=
z
−
k
B
(
z
−
1
)
A
(
z
−
1
)
|
z
=
e
j
ω
T
p
=
K
(
e
−
j
ω
T
p
)
{\displaystyle {\frac {y(i)}{u(i)}}=z^{-k}{\frac {B(z^{-1})}{A(z^{-1})}}{\Bigg |}_{z=e^{j\omega T_{p}}}=K(e^{-j\omega T_{p}})}
definiuje dyskretną transmitancję widmową .
Transformatę Fouriera można określić dla funkcji
f
∈
L
1
(
R
n
)
,
{\displaystyle f\in L^{1}(\mathbb {R} ^{n}),}
gdzie
L
1
(
R
n
)
{\displaystyle L^{1}(\mathbb {R} ^{n})}
jest przestrzenią wektorową funkcji całkowalnych na
R
n
,
{\displaystyle \mathbb {R} ^{n},}
za pomocą wzoru[7] :
f
^
(
ω
)
=
∫
R
n
f
(
x
)
e
−
i
(
x
,
ω
)
d
x
,
{\displaystyle {\hat {f}}(\omega )=\int \limits _{\mathbb {R} ^{n}}f(x)\ e^{-i(x,\omega )}\,dx,}
gdzie:
i
{\displaystyle i}
– jednostka urojona
(
i
2
=
−
1
)
,
{\displaystyle (i^{2}=-1),}
(
x
,
ω
)
{\displaystyle (x,\omega )}
– iloczyn skalarny wektorów
x
,
ω
∈
R
n
.
{\displaystyle x,\omega \in \mathbb {R} ^{n}.}
Tw. 1. Transformata
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
jest funkcją istotnie ograniczoną :
f
^
∈
L
∞
(
R
n
)
{\displaystyle {\hat {f}}\in L^{\infty }(\mathbb {R} ^{n})}
(wynika to z twierdzenia Riemanna-Lebesgue’a ).
Tw. 2. W przypadku gdy funkcja
f
{\displaystyle f}
jest ponadto całkowalna z kwadratem (czyli
f
∈
L
1
(
R
n
)
∩
L
2
(
R
n
)
{\displaystyle f\in L^{1}(\mathbb {R} ^{n})\cap L^{2}(\mathbb {R} ^{n})}
), transformata
f
^
(
ω
)
{\displaystyle {\hat {f}}(\omega )}
jest również całkowalna z kwadratem. Innymi słowy, transformacja Fouriera jest odwzorowaniem przestrzeni:
F
:
L
1
(
R
n
)
∩
L
2
(
R
n
)
→
L
2
(
R
n
)
.
{\displaystyle {\mathcal {F}}\colon L^{1}(\mathbb {R} ^{n})\cap L^{2}(\mathbb {R} ^{n})\to L^{2}(\mathbb {R} ^{n}).}
Tw. 3. Twierdzenie Plancherela wówczas mówi, że odwzorowanie to przedłuża się do izometrii przestrzeni
L
2
(
R
n
)
{\displaystyle L^{2}(\mathbb {R} ^{n})}
na siebie.
↑ Żakowski i Leksiński 2017 ↓ , s. 351.
↑ a b c d Żakowski i Leksiński 2017 ↓ , s. 360.
↑ Żakowski i Leksiński 2017 ↓ , s. 361–362.
↑ Żakowski i Leksiński 2017 ↓ , s. 365.
↑ Żakowski i Leksiński 2017 ↓ , s. 360, 365.
↑ David W. D.W. Kammler David W. D.W. , A First Course in Fourier Analysis , Upper Saddle River, NJ: Prentice Hall, 2000, ISBN 0-13-578782-3 , OCLC 43118245 . Brak numerów stron w książce
↑ Hanna Marcinkowska: Dystrybucje, przestrzenie Sobolewa, równania różniczkowe . PWN, 1993, s. 56–60. ISBN 83-01-10864-9 .
Polskojęzyczne
Anglojęzyczne
Eric W. E.W. Weisstein Eric W. E.W. , Fourier Transform , [w:] MathWorld , Wolfram Research [dostęp 2020-12-12] (ang. ) .
Interaktywne wprowadzenie do transformacji Fouriera [Angielski z animacjami]
Grant Sanderson, But what is the Fourier Transform? A visual introduction , kanał 3blue1brown , YouTube , 26 stycznia 2018 [dostęp 2021-03-15].
Aled Walker, F is for Fourier Transform (ang. ) , Oxford University Mathematical Institute , maths.ox.ac.uk, 19 sierpnia 2022 [dostęp 2023-05-29].
transformacje całkowe
inne transformacje
w rachunku prawdopodobieństwa