Salta al contenuto
Note per Studenti Serie di Fourier - coefficienti di segnali notevoli

Serie di Fourier - coefficienti di segnali notevoli

In questa pagina 7

In questa nota si calcolano con la formula di analisi (Serie di Fourier - analisi e sintesiUn segnale periodico di periodo $T$ si scrive come somma di esponenziali in relazione armonica, $x(t)=\sum_ka_ke^{jk\omega_0t}$ con $\omega_0=2\pi/T$, e i coefficienti si ottengono per proiezione, $a_k=\frac1T\int_Tx(t)e^{-jk\omega_0t}dt$. L'ortogonalità degli esponenziali dà la formula; la convergenza è in media quadratica (Riesz-Fischer), con il fenomeno di Gibbs nei salti.Serie di Fourier - analisi e sintesi →) i coefficienti dei segnali periodici che compaiono in tutti gli esercizi. Conviene impararli come una piccola tabella, ma soprattutto conoscere i metodi: (1) segnali trigonometrici: si legge con Eulero, senza integrali; (2) segnali a tratti: integrale diretto o derivata generalizzata; (3) segnali ottenuti dai precedenti con le proprietà (Serie di Fourier - proprietà e sistemi LTI in regime periodicoLe proprietà della serie di Fourier (linearità, ribaltamento, coniugio, ritardo $\to$ fattore $e^{-jk\omega_0t_1}$, modulazione $\to$ traslazione dei coefficienti, derivata $\to$ fattore $jk\omega_0$, Parseval $P=\sum|a_k|^2$, convoluzione periodica $\to$ prodotto) permettono di calcolare i coefficienti senza integrare. Un LTI stabile con ingresso periodico produce uscita periodica con coefficienti $b_k=H(k\omega_0)a_k$.Serie di Fourier - proprietà e sistemi LTI in regime periodico →).

1. Costante e sinusoidi: si leggono senza integrare

Se il segnale è già una combinazione di esponenziali, i coefficienti sono i coefficienti degli esponenziali, una volta scelto ω0\omega_0 tale che tutte le pulsazioni siano multiple di ω0\omega_0.

  • Costante x(t)=cx(t)=c: a0=ca_0=c, tutti gli altri nulli.
  • x(t)=Acos⁡(ω0t+φ0)x(t)=A\cos(\omega_0t+\varphi_0): per Eulero =A2ejφ0ejω0t+A2e−jφ0e−jω0t=\frac A2e^{j\varphi_0}e^{j\omega_0t}+\frac A2e^{-j\varphi_0}e^{-j\omega_0t}, quindi a1=A2ejφ0a_1=\frac A2e^{j\varphi_0}, a−1=A2e−jφ0a_{-1}=\frac A2e^{-j\varphi_0}, altri nulli.

Esempio: periodicità e coefficienti di x(t)=Asin⁡(πtΔ)+Acos⁡2(2πt3Δ)x(t)=A\sin\left(\frac{\pi t}\Delta\right)+A\cos^2\left(\frac{2\pi t}{3\Delta}\right). Il periodo comune è T=6ΔT=6\Delta (Segnali periodici, sinusoidi ed esponenziali immaginari puriUn segnale continuo è periodico di periodo T se x(t+T)=x(t); $e^{j2\pi f_0t}$ ha periodo minimo $1/|f_0|$ e una somma di periodici è periodica solo se il rapporto dei periodi è razionale. A tempo discreto $e^{j2\pi\nu n}$ è periodico solo se $\nu$ è razionale, e le pulsazioni che differiscono di $2\pi$ sono lo stesso segnale: ogni sinusoide si riporta alla forma canonica con $\omega\in[0,\pi]$.Segnali periodici, sinusoidi ed esponenziali immaginari puri →), quindi ω0=2π6Δ=π3Δ\omega_0=\frac{2\pi}{6\Delta}=\frac{\pi}{3\Delta}. Il primo termine ha pulsazione πΔ=3ω0\frac{\pi}{\Delta}=3\omega_0; cos⁡2(2πt3Δ)=12+12cos⁡(4πt3Δ)\cos^2\left(\frac{2\pi t}{3\Delta}\right)=\frac12+\frac12\cos\left(\frac{4\pi t}{3\Delta}\right) e 4π3Δ=4ω0\frac{4\pi}{3\Delta}=4\omega_0. Scrivendo con Eulero: x(t)=A2jej3ω0t−A2je−j3ω0t+A2+A4ej4ω0t+A4e−j4ω0t.x(t)=\frac A{2j}e^{j3\omega_0t}-\frac{A}{2j}e^{-j3\omega_0t}+\frac A2+\frac A4e^{j4\omega_0t}+\frac A4e^{-j4\omega_0t}. Quindi a0=A2a_0=\frac A2, a3=A2j=−jA2a_3=\frac{A}{2j}=-\frac{jA}2, a−3=jA2a_{-3}=\frac{jA}{2}, a±4=A4a_{\pm4}=\frac A4, tutti gli altri nulli.

2. Il pettine periodico

x(t)=comb⁡T(t)=∑kδ(t−kT)x(t)=\operatorname{comb}_T(t)=\sum_k\delta(t-kT). Nel periodo [−T2,T2][-\frac T2,\frac T2] c'è un solo impulso in t=0t=0: ak=1T∫−T/2T/2δ(t) e−jkω0t dt=1T e0=1T∀k.a_k=\frac1T\int_{-T/2}^{T/2}\delta(t)\,e^{-jk\omega_0t}\,dt=\frac1T\,e^{0}=\frac1T\quad\forall k. Tutte le armoniche hanno la stessa ampiezza 1T\frac1T: il pettine contiene "tutte le frequenze allo stesso modo" (un impulso contiene tutte le frequenze: si veda anche Trasformata di Fourier di segnali periodici e generalizzatiCon le delta di Dirac in frequenza si trasformano anche segnali non integrabili: $1\leftrightarrow2\pi\delta(\omega)$, $e^{j\omega_0t}\leftrightarrow2\pi\delta(\omega-\omega_0)$, $\cos\omega_0t\leftrightarrow\pi[\delta(\omega-\omega_0)+\delta(\omega+\omega_0)]$, $u(t)\leftrightarrow\pi\delta(\omega)+\frac1{j\omega}$. Un segnale periodico con coefficienti $a_k$ ha spettro a righe $2\pi\sum a_k\delta(\omega-k\omega_0)$: la periodizzazione nel tempo è un campionamento in frequenza.Trasformata di Fourier di segnali periodici e generalizzati →). Il valor medio vale a0=1Ta_0=\frac1T (area 11 per periodo).

3. Onda quadra con duty cycle qualsiasi

x(t)=rep⁡Trect⁡(t2a)x(t)=\operatorname{rep}_T\operatorname{rect}\left(\frac{t}{2a}\right): impulsi rettangolari di ampiezza 11 e durata 2a2a, ripetuti con periodo TT. Il duty cycle è la frazione di periodo in cui il segnale è attivo, d=2aTd=\frac{2a}{T}. Con ω0=2πT\omega_0=\frac{2\pi}T e integrando su [−T2,T2][-\frac T2,\frac T2] (dove il segnale vale 11 solo in ∣t∣<a|t|<a): ak=1T∫−aae−jkω0tdt=1T⋅ejkω0a−e−jkω0ajkω0=1T⋅2sin⁡(kω0a)kω0(k≠0).a_k=\frac1T\int_{-a}^{a}e^{-jk\omega_0t}dt=\frac1T\cdot\frac{e^{jk\omega_0a}-e^{-jk\omega_0a}}{jk\omega_0}=\frac1T\cdot\frac{2\sin(k\omega_0a)}{k\omega_0}\quad(k\ne0). Si porta 2a2a fuori e si ricostruisce un sinc: 2sin⁡(kω0a)T kω0=2aT⋅sin⁡(kω0a)kω0a\frac{2\sin(k\omega_0a)}{T\,k\omega_0}=\frac{2a}{T}\cdot\frac{\sin(k\omega_0a)}{k\omega_0a}. Con kω0a=k2πaT=πkdk\omega_0a=k\frac{2\pi a}T=\pi kd il rapporto è sin⁡πkdπkd=sinc⁡(kd)\frac{\sin\pi kd}{\pi kd}=\operatorname{sinc}(kd):

ak=d sinc⁡(kd)=2aTsinc⁡(2akT)(vale anche per k=0: a0=d).\boxed{a_k=d\,\operatorname{sinc}(kd)=\frac{2a}T\operatorname{sinc}\left(\frac{2ak}{T}\right)}\quad(\text{vale anche per }k=0:\ a_0=d).

Le armoniche si annullano quando kdkd è un intero non nullo.

Esempi numerici.

  • T=2T=2, a=12a=\frac12 (d=12d=\frac12): ak=12sinc⁡(k2)a_k=\frac12\operatorname{sinc}\left(\frac k2\right) (esempio della nota precedente).
  • T=4T=4, x=rep⁡4rect⁡(t2)x=\operatorname{rep}_4\operatorname{rect}\left(\frac t2\right) (a=1a=1, d=12d=\frac12): ak=12sinc⁡(k2)a_k=\frac12\operatorname{sinc}\left(\frac{k}{2}\right); a0=12a_0=\frac12, a±1=1πa_{\pm1}=\frac1\pi, a±2=0a_{\pm2}=0, a±3=−13πa_{\pm3}=-\frac1{3\pi} (tema d'esame febbraio 2026, Esercizio 13 · onda quadra attraverso un filtro passa-basso ideale (tema d'esame febbraio 2026), con ω0=π2\omega_0=\frac\pi2).
  • T=3T=3, x=rep⁡3rect⁡(t2)x=\operatorname{rep}_3\operatorname{rect}\left(\frac t2\right): d=23d=\frac23, ak=23sinc⁡(2k3)a_k=\frac23\operatorname{sinc}\left(\frac{2k}3\right) (Esercizio 14 · segnale periodico determinato da condizioni sui coefficienti (tema d'esame aprile 2024)).

Effetto del duty cycle. Se dd diminuisce (impulsi più stretti) i coefficienti decadono più lentamente: i primi zeri sono in k=1dk=\frac1d, cioè lo spettro si allarga. Al limite (a→0a\to0 con area 2a⋅12a=12a\cdot\frac{1}{2a}=1) si ritrova il pettine, con tutti i coefficienti uguali. È la riprova che un segnale breve nel tempo ha spettro largo (e viceversa).

4. Coseno rettificato

x(t)=∣cos⁡(πt)∣x(t)=|\cos(\pi t)|. Il coseno ha periodo 22, ma il valore assoluto "ribalta" le semionde negative: il periodo diventa T=1T=1 (∣cos⁡π(t+1)∣=∣−cos⁡πt∣|\cos\pi(t+1)|=|-\cos\pi t|) e ω0=2π\omega_0=2\pi. Su [−12,12][-\frac12,\frac12] vale x=cos⁡(πt)≥0x=\cos(\pi t)\ge0:

ak=∫−1/21/2cos⁡(πt) e−j2πktdt=12∫−1/21/2[e−j2πt(k−12)+e−j2πt(k+12)]dt.a_k=\int_{-1/2}^{1/2}\cos(\pi t)\,e^{-j2\pi kt}dt=\frac12\int_{-1/2}^{1/2}\left[e^{-j2\pi t(k-\frac12)}+e^{-j2\pi t(k+\frac12)}\right]dt. Per l'integrale esponenziale (Numeri complessi, formula di Eulero ed esponenziali complessiUn numero complesso si scrive in forma cartesiana a+jb o polare |x|e^{jφ}; il prodotto moltiplica i moduli e somma le fasi. La formula di Eulero e^{jα}=cos α + j sin α lega esponenziali e sinusoidi e permette di trattare tutti i segnali del corso come somme di esponenziali complessi e^{(σ+jω)t}.Numeri complessi, formula di Eulero ed esponenziali complessi →): ∫−1/21/2e−j2πtβdt=ejπβ−e−jπβj2πβ=sin⁡πβπβ=sinc⁡(β)\int_{-1/2}^{1/2}e^{-j2\pi t\beta}dt=\frac{e^{j\pi\beta}-e^{-j\pi\beta}}{j2\pi\beta}=\frac{\sin\pi\beta}{\pi\beta}=\operatorname{sinc}(\beta). Quindi ak=12[sinc⁡(k−12)+sinc⁡(k+12)].a_k=\frac12\left[\operatorname{sinc}\left(k-\tfrac12\right)+\operatorname{sinc}\left(k+\tfrac12\right)\right]. Si può semplificare: sinc⁡(k−12)=sin⁡(πk−π/2)π(k−12)=−cos⁡πkπ(k−12)=−(−1)kπ(k−12)\operatorname{sinc}(k-\frac12)=\frac{\sin(\pi k-\pi/2)}{\pi(k-\frac12)}=\frac{-\cos\pi k}{\pi(k-\frac12)}=-\frac{(-1)^k}{\pi(k-\frac12)} e sinc⁡(k+12)=sin⁡(πk+π/2)π(k+12)=(−1)kπ(k+12)\operatorname{sinc}(k+\frac12)=\frac{\sin(\pi k+\pi/2)}{\pi(k+\frac12)}=\frac{(-1)^k}{\pi(k+\frac12)}. Sommando: ak=(−1)k2π[1k+12−1k−12]=(−1)k2π⋅−1k2−14=2 (−1)k+1π (4k2−1).a_k=\frac{(-1)^k}{2\pi}\left[\frac{1}{k+\frac12}-\frac1{k-\frac12}\right]=\frac{(-1)^k}{2\pi}\cdot\frac{-1}{k^2-\frac14}=\frac{2\,(-1)^{k+1}}{\pi\,(4k^2-1)}. Valori: a0=2π≈0,637a_0=\frac2\pi\approx0{,}637 (è il valor medio di ∣cos⁡∣|\cos|), a±1=23π≈0,212a_{\pm1}=\frac{2}{3\pi}\approx0{,}212, a±2=−215π≈−0,042a_{\pm2}=-\frac{2}{15\pi}\approx-0{,}042, a±3=235π≈0,018a_{\pm3}=\frac2{35\pi}\approx0{,}018. Il decadimento è come 1k2\frac1{k^2}: più rapido dell'onda quadra, perché ∣cos⁡∣|\cos| è continua (ha solo spigoli, non salti). Questi valori sono stati verificati con integrazione numerica.

Grafico interattivo: Coseno rettificato x(t) = |cos(πt)|: periodo fondamentale 1 (la metà del periodo del coseno)

Grafico interattivo: Coefficienti di |cos(πt)|: a_k = 2·(−1)^{k+1}/(π(4k²−1)), con a₀ = 2/π ≈ 0,64 (valor medio)

5. Dente di sega e onda triangolare

Dente di sega x(t)=tx(t)=t per ∣t∣<1|t|<1, periodo T=2T=2 (ω0=π\omega_0=\pi): è dispari. Con l'integrazione per parti, ∫−11te−jkπtdt=[te−jkπt−jkπ]−11+1jkπ∫−11e−jkπtdt\int_{-1}^{1}te^{-jk\pi t}dt=\left[\frac{te^{-jk\pi t}}{-jk\pi}\right]_{-1}^1+\frac{1}{jk\pi}\int_{-1}^{1}e^{-jk\pi t}dt, e l'ultimo integrale è nullo (numero intero di giri); resta e−jkπ+ejkπ−jkπ=2(−1)k−jkπ=2j(−1)kkπ\frac{e^{-jk\pi}+e^{jk\pi}}{-jk\pi}=\frac{2(-1)^k}{-jk\pi}=\frac{2j(-1)^k}{k\pi}. Dividendo per T=2T=2: ak=j(−1)kkπ(k≠0),a0=0.a_k=\frac{j(-1)^k}{k\pi}\quad(k\ne0),\qquad a_0=0. I coefficienti sono immaginari e dispari (segnale reale e dispari). Il decadimento 1k\frac1k è quello tipico di un salto.

Grafico interattivo: Dente di sega x(t) = t per |t| < 1, periodo 2: salto da +1 a −1 ad ogni multiplo dispari di 1

Metodo della derivata (più veloce per i segnali lineari a tratti). Se bkb_k sono i coefficienti di x′x' (derivata generalizzata), la derivata moltiplica per jkω0jk\omega_0 (Serie di Fourier - proprietà e sistemi LTI in regime periodicoLe proprietà della serie di Fourier (linearità, ribaltamento, coniugio, ritardo $\to$ fattore $e^{-jk\omega_0t_1}$, modulazione $\to$ traslazione dei coefficienti, derivata $\to$ fattore $jk\omega_0$, Parseval $P=\sum|a_k|^2$, convoluzione periodica $\to$ prodotto) permettono di calcolare i coefficienti senza integrare. Un LTI stabile con ingresso periodico produce uscita periodica con coefficienti $b_k=H(k\omega_0)a_k$.Serie di Fourier - proprietà e sistemi LTI in regime periodico →): bk=jkω0akb_k=jk\omega_0a_k, quindi ak=bkjkω0a_k=\frac{b_k}{jk\omega_0} per k≠0k\ne0 (il valor medio a0a_0 va calcolato a parte: la derivata lo perde). Per il dente di sega: x′(t)=1−2comb⁡2(t−1)x'(t)=1-2\operatorname{comb}_2(t-1) (pendenza 11 e un salto di −2-2 nei multipli dispari di 1). I coefficienti di comb⁡2(t−1)\operatorname{comb}_2(t-1) sono 12e−jkπ⋅1=12(−1)k\frac12e^{-jk\pi\cdot1}=\frac12(-1)^k, quindi bk=−(−1)kb_k=-(-1)^k per k≠0k\ne0 e ak=−(−1)kjkπ=j(−1)kkπa_k=\frac{-(-1)^k}{jk\pi}=\frac{j(-1)^k}{k\pi} ✓.

Onda triangolare x(t)=rep⁡TΛ(2tT)x(t)=\operatorname{rep}_T\Lambda\left(\frac{2t}T\right) (triangoli di altezza 11 e base TT che si toccano): valor medio 12\frac12 e ak=12sinc⁡2(k2):a0=12, a±1=2π2≈0,203, a±2=0, a±3=29π2≈0,0225.a_k=\frac12\operatorname{sinc}^2\left(\frac k2\right):\quad a_0=\tfrac12,\ a_{\pm1}=\tfrac2{\pi^2}\approx0{,}203,\ a_{\pm2}=0,\ a_{\pm3}=\tfrac{2}{9\pi^2}\approx0{,}0225. (Verificato numericamente.) Derivando una volta si ottiene un'onda quadra ±, derivando due volte due pettini: decadimento 1k2\frac1{k^2}.

Tabella riassuntiva

Segnale (periodo TT, ω0=2π/T\omega_0=2\pi/T) Coefficienti aka_k
costante cc a0=ca_0=c
Acos⁡(ω0t+φ)A\cos(\omega_0t+\varphi) a±1=A2e±jφa_{\pm1}=\frac A2e^{\pm j\varphi}
comb⁡T(t)\operatorname{comb}_T(t) 1T\frac1T per ogni kk
rep⁡Trect⁡(t2a)\operatorname{rep}_T\operatorname{rect}\left(\frac t{2a}\right) 2aTsinc⁡(2akT)\frac{2a}T\operatorname{sinc}\left(\frac{2ak}T\right)
∣cos⁡(πt)∣\lvert\cos(\pi t)\rvert (T=1T=1) 2(−1)k+1π(4k2−1)\frac{2(-1)^{k+1}}{\pi(4k^2-1)}
dente di sega tt su ∣t∣<1\lvert t\rvert<1 (T=2T=2) j(−1)kkπ\frac{j(-1)^k}{k\pi}, a0=0a_0=0
triangolare rep⁡TΛ(2tT)\operatorname{rep}_T\Lambda\left(\frac{2t}T\right) 12sinc⁡2(k2)\frac12\operatorname{sinc}^2\left(\frac k2\right)

Errori comuni

  • Scegliere il periodo sbagliato: per ∣cos⁡(πt)∣|\cos(\pi t)| il periodo è 11, non 22.
  • Perdere il valor medio quando si usa il metodo della derivata.
  • Dimenticare di ricostruire il sinc con il giusto fattore: ak=dsinc⁡(kd)a_k=d\operatorname{sinc}(kd), non sinc⁡(kd)\operatorname{sinc}(kd).
  • Per le sinusoidi, dimenticare il fattore 12\frac12 (cos⁡=12ej+12e−j\cos=\frac12e^{j}+\frac12e^{-j}).

Versione ripasso

  • Metodi: (1) sinusoidi: Eulero, senza integrali, con ω0\omega_0 tale che tutte le pulsazioni siano multiple; (2) a tratti: integrale o derivata generalizzata ak=bkjkω0a_k=\frac{b_k}{jk\omega_0} (a0a_0 a parte); (3) proprietà.
  • Sinusoidi: Acos⁡(ω0t+φ)A\cos(\omega_0t+\varphi): a±1=A2e±jφa_{\pm1}=\frac A2e^{\pm j\varphi}. Es.: Asin⁡πtΔ+Acos⁡22πt3ΔA\sin\frac{\pi t}\Delta+A\cos^2\frac{2\pi t}{3\Delta}, T=6ΔT=6\Delta: a0=A2a_0=\frac A2, a±3=∓jA2a_{\pm3}=\mp\frac{jA}2 (cioè a3=−jA2a_3=-\frac{jA}{2}), a±4=A4a_{\pm4}=\frac A4.
  • Pettine: ak=1Ta_k=\frac1T per ogni kk.
  • Onda quadra rep⁡Trect⁡t2a\operatorname{rep}_T\operatorname{rect}\frac t{2a}, duty cycle d=2aTd=\frac{2a}T: ak=dsinc⁡(kd)a_k=d\operatorname{sinc}(kd). Es.: T=4,a=1T=4,a=1: 12sinc⁡k2\frac12\operatorname{sinc}\frac k2; T=3,a=1T=3,a=1: 23sinc⁡2k3\frac23\operatorname{sinc}\frac{2k}3. Impulsi più stretti ⇒\Rightarrow spettro più largo.
  • ∣cos⁡πt∣|\cos\pi t| (T=1T=1): ak=12[sinc⁡(k−12)+sinc⁡(k+12)]=2(−1)k+1π(4k2−1)a_k=\frac12[\operatorname{sinc}(k-\frac12)+\operatorname{sinc}(k+\frac12)]=\frac{2(-1)^{k+1}}{\pi(4k^2-1)}; a0=2πa_0=\frac2\pi, a±1=23πa_{\pm1}=\frac2{3\pi}, a±2=−215πa_{\pm2}=-\frac2{15\pi}.
  • Dente di sega (tt su ∣t∣<1|t|<1, T=2T=2): ak=j(−1)kkπa_k=\frac{j(-1)^k}{k\pi}, a0=0a_0=0. Triangolare rep⁡TΛ2tT\operatorname{rep}_T\Lambda\frac{2t}T: 12sinc⁡2k2\frac12\operatorname{sinc}^2\frac k2.
  • Decadimento: salti ∼1k\sim\frac1k; solo spigoli ∼1k2\sim\frac1{k^2}.
  • Errori: periodo di ∣cos⁡∣|\cos| (è 1, non 2); media persa con la derivata; dd dimenticato nel sinc. Vedi Serie di Fourier - analisi e sintesiUn segnale periodico di periodo $T$ si scrive come somma di esponenziali in relazione armonica, $x(t)=\sum_ka_ke^{jk\omega_0t}$ con $\omega_0=2\pi/T$, e i coefficienti si ottengono per proiezione, $a_k=\frac1T\int_Tx(t)e^{-jk\omega_0t}dt$. L'ortogonalità degli esponenziali dà la formula; la convergenza è in media quadratica (Riesz-Fischer), con il fenomeno di Gibbs nei salti.Serie di Fourier - analisi e sintesi →.

Esercizi su questo argomento

Teoria collegata