În matematică , reprezentarea spectrală a semnalelor este o descriere formală a semnalelor ( funcții în timp) în domeniul frecvenței , adică în ceea ce privește frecvența lor, care este utilizată în multe domenii ale științei, cum ar fi ingineria și fizica . În această descriere, fiecare frecvență din care este compus un semnal se numește armonic și, din punct de vedere matematic, fiecare armonică este făcută să corespundă unui vector al unei baze a unui spațiu vectorial infinit-dimensional cu produs intern (produs scalar) pe câmp complex, care este baza unui spațiu Hilbert . Semnalul este apoi scris ca o combinație liniară în acel spațiu. Analiza de frecvență a comportamentului unui sistem dinamic se numește răspuns în frecvență al sistemului dinamic.
Un spațiu Hilbert este un spațiu vectorial cu un produs scalar pe câmpulreal sau complex care este complet în raport cu distanța indusă de acel produs scalar. Având în vedere un set de vectori {\ displaystyle \ {\ mathbf {s} _ {1}, \ mathbf {s} _ {2}, \ dots \}} a unui spațiu Hilbert complex {\ displaystyle {\ mathcal {H}}} prin urmare, suma și produsul pentru un scalar mențin acești vectori în spațiu:
cu {\ displaystyle \ alpha \ in \ mathbb {R}} (sau {\ displaystyle \ mathbb {C}} ). Mai mult, există doar un invers al sumei{\ displaystyle - \ mathbf {s}} astfel încât {\ displaystyle \ mathbf {s} + (- \ mathbf {s}) = \ mathbf {0}} . În acest context putem defini dependența și independența liniară a vectorilor și conceptul de bază . O bază este un set de vectori liniar independenți care sunt, de asemenea, un sistem de generatori, adică un sistem de vectori {\ displaystyle \ {\ mathbf {u} _ {1}, \ mathbf {u} _ {2}, \ dots \}} este liniar independent și formează un sistem complet astfel încât orice alt vector poate fi reprezentat ca o combinație liniară, posibil infinită, de vectori ai unei baze:
unde este {\ displaystyle c_ {i} \ in \ mathbb {C}} sunt coeficienții combinației liniare. Un spațiu Hilbert este un spațiu normat , adică este definită norma unui vector, este un număr real astfel încât:
Reprezentarea spectrală se bazează pe faptul că orice funcție (semnal) definită într-un interval{\ displaystyle [t_ {1}, t_ {2}]} poate fi dezvoltat în seria Fourier ca o combinație liniară de vectori {\ displaystyle \ mathbf {u} _ {i} (t)} (la rândul său funcții ale timpului) aparținând unei baze ortonormale:
Semnalele periodice sunt de așa natură încât {\ displaystyle s (t) = s (t + T)} , unde este {\ displaystyle T} este perioada : acestea sunt semnalele care se repetă identic după un timp {\ displaystyle T} . Luați în considerare un semnal periodic {\ displaystyle s (t)} continuu, a cărui serie Fourier este:
unde este {\ displaystyle c_ {m} = (s, u_ {m})} sunt coeficienți determinabili cu produsul scalar și {\ displaystyle \ {u_ {m} = e ^ {imt}, m \ in \ mathbb {Z} \}} este baza ortonormală a funcțiilor exponențiale. De sine {\ displaystyle \ omega _ {1} = 2 \ pi / T} este pulsația fundamentală (cea mai mică frecvență a semnalului), sumarea anterioară ia forma: [3]
{\ displaystyle s (t) = {\ frac {a_ {0}} {2}} + \ sum _ {n = 1} ^ {\ infty} \ left [a_ {n} \ cos (n \ omega _ { 1} t) + b_ {n} \ sin (n \ omega _ {1} t) \ right]}
Primul termen este constant și toți ceilalți termeni sunt o combinație liniară de coeficienți adecvați {\ displaystyle a_ {n}} Și {\ displaystyle b_ {n}} a funcțiilor exponențiale. Pentru a determina coeficienții {\ displaystyle a_ {n}} Și {\ displaystyle b_ {n}} produsul dot este în general utilizat.
Constanta {\ displaystyle a_ {0} / 2} este egal cu valoarea medie a semnalului în perioada de definiție, de fapt:
{\ displaystyle \ int _ {- T / 2} ^ {T / 2} s (t) \, dt = \ int _ {- T / 2} ^ {T / 2} {\ frac {a_ {0}} {2}} \, dt + \ int _ {- T / 2} ^ {T / 2} \ sum _ {n = 1} ^ {\ infty} \ left [a_ {n} \ cos (n \ omega _ {1} t) + b_ {n} \ sin (n \ omega _ {1} t) \ right] dt = {\ frac {a_ {0}} {2}} T}
unde a doua integrală a celui de-al doilea termen dispare deoarece integrala pe o perioadă a funcțiilor exponențiale este zero prin simetrie. Prin urmare, avem:
{\ displaystyle {\ frac {a_ {0}} {2}} = {\ frac {1} {T}} \ int _ {- T / 2} ^ {T / 2} s (t) \, dt}
adică valoarea medie a semnalului pe parcursul perioadei {\ displaystyle T} . Pentru a determina coeficienții rămași {\ displaystyle a_ {n}} produsul scalar este executat:
{\ displaystyle b_ {n} = {\ frac {2} {T}} \ int _ {- T / 2} ^ {T / 2} s (t) \ cdot \ sin (n \ omega _ {1} t ) \, dt}
Proprietățile reprezentării seriei Fourier
În reprezentarea semnalului prin seria Fourier, un semnal periodic este descompus într-un set infinit de frecvențe multiple ale celei fundamentale {\ displaystyle \ omega _ {1}} , adică{\ displaystyle \ omega _ {n} = n \ omega _ {1}} , și sunt numite armonice (termen care nu are legătură cu conceptul de funcție armonică ). Fiecare dintre aceste componente spectrale are o amplitudine egală cu:
Dacă semnalul este o funcție uniformă a timpului, adică dacă {\ displaystyle s (t) = s (-t)} , apoi toate armonicile care conțin sinusul (care este o funcție ciudată) se anulează. Deci seria devine:
{\ displaystyle a_ {n} = {\ frac {2} {T}} \ int _ {- T / 2} ^ {T / 2} s (t) \ cdot \ cos (n \ omega _ {1} t ) \, dt}
În mod similar, dacă semnalul este o funcție ciudată a timpului, adică dacă {\ displaystyle s (t) = - s (-t)} , toate armonicile care conțin cosinusul se anulează (deci și valoarea medie) și seria devine:
{\ displaystyle s (t) = \ sum _ {n = 1} ^ {\ infty} b_ {n} \ cdot \ sin (n \ omega _ {1} t)}
cu coeficienți:
{\ displaystyle b_ {n} = {\ frac {2} {T}} \ int _ {- T / 2} ^ {T / 2} s (t) \ cdot \ sin (n \ omega _ {1} t ) \, dt}
Formă complexă a seriei Fourier
Formulele lui Euler pot fi folosite în continuare:
{\ displaystyle \ cos z = {\ frac {e ^ {iz} + e ^ {- iz}} {2}} \, \, \, \, \, \ sin z = {\ frac {e ^ {iz } -e ^ {- iz}} {2i}}}
pentru a obține o formă alternativă la seria Fourier:
Prin aceste transformări matematice putem rescrie seria Fourier ca:
{\ displaystyle s (t) = \ sum _ {n = - \ infty} ^ {\ infty} c_ {n} \ cdot e ^ {în \ omega _ {1} t}}
unde este:
{\ displaystyle c_ {n} = {\ frac {1} {T}} \ int _ {- T / 2} ^ {T / 2} s (t) e ^ {- în \ omega _ {1} t} \, dt}
Rețineți că seria este definită și pentru {\ displaystyle n} negativ.
Reprezentarea semnalelor non-periodice
Reprezentarea semnalelor non-periodice se realizează și folosind baza ortonormală formată din funcțiile armonice, cu condiția ca funcția non-periodică să scadă la infinit cu o regularitate suficientă. Această constrângere se datorează faptului că metoda utilizată pentru reprezentarea în frecvență constă în construirea unui semnal periodic dat de repetarea infinită a unui semnal non-periodic, care trebuie definit într-un interval de timp dincolo de care este nul.
unde este {\ displaystyle f} înseamnă {\ displaystyle s} sau {\ displaystyle \ phi} . Dacă această condiție este validă, transformarea {\ displaystyle s} și anti-transformarea {\ displaystyle \ phi} sunt funcții continue, mărginite și următoarele sunt valabile:
pentru fiecare {\ displaystyle f, g \ în L ^ {1} (\ mathbb {R})} Și {\ displaystyle \ alpha, \ beta \ in \ mathbb {C}} .
Spectrul derivatei și al integralei
Derivata semnalului în timp corespunde, în domeniul frecvenței, multiplicării cu {\ displaystyle i \ omega} a transformării semnalului nederivat. Într-adevăr, fie el {\ displaystyle s (t)} un semnal e {\ displaystyle S (\ omega)} sa transformat. Atunci derivata semnalului este:
{\ displaystyle \ int _ {- \ infty} ^ {+ \ infty} {\ frac {ds (t)} {dt}} e ^ {- i \ omega t} \, dt = \ left [if ^ {- i \ omega t} \ right] _ {- \ infty} ^ {+ \ infty} - \ int _ {- \ infty} ^ {+ \ infty} s (t) \ left [-i \ omega e ^ {- i \ omega t} \ right] dt = 0 + i \ omega S (\ omega)}
De aici și transformarea {\ displaystyle ds (t) / dt} Și {\ displaystyle i \ omega S (\ omega)} și transformarea {\ displaystyle d ^ {n} s (t) / dt ^ {n}} Și{\ displaystyle (i \ omega) ^ {n} S (\ omega)} .
Spectrul integralei unui semnal este dat în schimb de împărțirea prin {\ displaystyle i \ omega} a transformării semnalului (neintegrată). Este {\ displaystyle s (t)} un semnal e {\ displaystyle S (\ omega)} transformarea sa, apoi transformarea integralei semnalului:
{\ displaystyle \ int _ {- \ infty} ^ {t} s (\ tau) d \ tau}
O caracteristică deosebit de utilă a reprezentării spectrale este că convoluția a două funcții în timp este echivalentă cu produsul algebric al transformărilor lor în domeniul frecvenței. Într-adevăr, scriind transformarea produsului {\ displaystyle s (t) = u (t) \ cdot v (t)} a două semnale precum:
{\ displaystyle S (\ omega) = \ int _ {- \ infty} ^ {+ \ infty} u (t) \ cdot v (t) e ^ {- i \ omega t} dt = {\ frac {1} {2 \ pi}} \ int _ {- \ infty} ^ {+ \ infty} u (t) \ left [\ int _ {- \ infty} ^ {+ \ infty} V (\ xi) e ^ {i \ xi t} \, d \ xi \ right] e ^ {- i \ omega t} \, dt = {\ frac {1} {2 \ pi}} \ int _ {- \ infty} ^ {+ \ infty } V (\ xi) \ left [\ int _ {- \ infty} ^ {+ \ infty} u (t) e ^ {- i (\ omega - \ xi) t} \, dt \ right] d \ xi }
unde în primul pas funcția de pornire a fost scrisă ca antitransformă a transformării (între paranteze pătrate), în timp ce în al doilea termenul între paranteze pătrate este transformarea {\ displaystyle U (\ omega - \ xi)} a funcției {\ displaystyle u} tradus prin multiplicare prin exponențial. Prin urmare:
{\ displaystyle S (\ omega) = {\ frac {1} {2 \ pi}} \ int _ {- \ infty} ^ {+ \ infty} V (\ xi) \ cdot U (\ omega - \ xi) d \ xi}
această integrală este un produs al convoluției și este scrisă simbolic ca:
{\ displaystyle s (t) = u (t) \ cdot v (t) \ qquad S (\ omega) = U (\ omega) * V (\ omega)}
Inversul este, de asemenea, adevărat, dacă avem produsul obișnuit din două spectre:
Transformarea unui semnal real {\ displaystyle s (t)} poate fi scris generic ca:
{\ displaystyle S (\ omega) = \ int _ {- \ infty} ^ {+ \ infty} s (t) \, \ cos (\ omega t) \, dt-i \ cdot \ int _ {- \ infty } ^ {+ \ infty} s (t) \, \ sin (\ omega t) \, dt = A (\ omega) -iB (\ omega)}
cu {\ displaystyle A (\ omega) = A (- \ omega)} este partea reală a transformării și este o funcție uniformă, în timp ce {\ displaystyle B (\ omega) = - B (- \ omega)} este partea imaginară a spectrului și este o funcție ciudată . Dacă se efectuează antitransformarea, semnalul real se obține din nou în timp:
{\ displaystyle s (t) = {\ frac {1} {2 \ pi}} \ int _ {- \ infty} ^ {+ \ infty} [A (\ omega) -iB (\ omega)] \ cdot [ \ cos (\ omega t) + i \ cdot \ sin (\ omega t)] \, d \ omega}
Prin dezvoltarea produsului în cadrul integralei avem:
{\ displaystyle s (t) = {\ frac {1} {2 \ pi}} \ int _ {- \ infty} ^ {+ \ infty} [A (\ omega) \ cos (\ omega t) + B ( \ omega) \ sin (\ omega t)] + i \ cdot [A (\ omega) \ sin (\ omega t) -B (\ omega) \ cos (\ omega t)] \, d \ omega}
Pentru ca semnalul să fie real, trebuie să se întâmple în mod necesar ca partea reală și imaginară să fie ambele nule, adică:
{\ displaystyle \ int _ {- \ infty} ^ {+ \ infty} A (\ omega) \ sin (\ omega t) \, d \ omega = 0 \, \, \, \, \, \ int _ { - \ infty} ^ {+ \ infty} B (\ omega) \ cos (\ omega t) \, d \ omega t = 0}
iar această condiție este îndeplinită numai dacă partea reală este pară și partea imaginară este ciudată. Reversul este, de asemenea, adevărat, deci dacă partea reală a transformării unui semnal este pară și partea imaginară este impară, atunci se obține un semnal real.
( EN ) Walter Rudin, Analiză reală și complexă , Mladinska Knjiga, McGraw-Hill, 1970, ISBN0-07-054234-1 .
( EN ) Michael Reed, Barry Simon, Methods of Modern Mathematical Physics, Vol. 1: Functional Analysis , ed. A II-a, San Diego, California, Academic press inc., 1980, ISBN0-12-585050-6 .