Trasformata zeta

In analisi funzionale la trasformata zeta è una trasformata integrale che permette di trasformare una funzione discreta in una funzione più semplice, utilizzata principalmente nella teoria dei segnali.

Storia

Il concetto di trasformata zeta era già noto a Laplace, ma fu reintrodotto nel 1947 da W. Hurewicz come mezzo utile a risolvere equazioni alle differenze lineari a coefficienti costanti.[1] Il termine "trasformata zeta" fu coniato successivamente, nel 1952, da Ragazzini e Zadeh, ricercatori della Columbia University.[2][3] Il nome potrebbe esser derivato dall'idea che la lettera "z" sia somigliante a una lettera "s" campionata/digitalizzata, ove "s" è la lettera spesso usata per indicare la variabile indipendente nella trasformata di Laplace. Un'altra possibile origine è la presenza della lettera "Z" in entrambi i nomi Ragazzini e Zadeh. Questa nomenclatura diverge dall'usanza adottata in ambito scientifico, in cui si associa un metodo o un teorema col nome del principale sviluppatore. La terza probabile origine risiede nel dominio dei segnali discreti, che è solito essere Z {\displaystyle \mathbb {Z} } o un suo sottoinsieme.

Definizione

Trasformata unilatera

Sia x [ n ] {\displaystyle x[n]} una successione di numeri complessi, indicizzata con n N {\displaystyle n\in \mathbb {N} } . La sua trasformata unilatera è definita come la serie formale di potenze complesse

X ( z ) = Z { x [ n ] } = n = 0 x [ n ] z n ,  per  z C . {\displaystyle X(z)={\mathcal {Z}}\{x[n]\}=\sum _{n=0}^{\infty }x[n]z^{-n},\qquad {\mbox{ per }}z\in \mathbb {C} .}

In teoria dei segnali questa definizione è utilizzata per valutare la trasformata della risposta all'impulso unitario di un sistema causale tempo-discreto. Solitamente, in tale ambito la successione x [ n ] {\displaystyle x[n]} rappresenta il campionamento regolare di un segnale f : R C {\displaystyle f:\mathbb {R} \to \mathbb {C} } causale (i.e. f {\displaystyle f} è nulla per tempi negativi), in corrispondenza dei tempi della forma t = n τ {\displaystyle t=n\,\tau } . Il passo di campionamento τ > 0 {\displaystyle \tau >0} è fissato. In altre parole

x [ n ] = f ( n τ ) ,  per ogni  n N . {\displaystyle x[n]=f(n\,\tau ),\qquad {\mbox{ per ogni }}n\in \mathbb {N} .}

Regione di convergenza

La regione di convergenza è la parte di piano complesso dove la serie che definisce la trasformata della funzione converge:

R O C = { z C : | n = 0 x [ n ] z n | < } {\displaystyle ROC=\left\{z\in \mathbb {C} \,:\,\left|\sum _{n=0}^{\infty }x[n]\,z^{-n}\right|<\infty \right\}}

La serie converge per valori di z {\displaystyle z} in modulo maggiori del raggio di convergenza R {\displaystyle R} , definito tramite il criterio della radice come:

R = lim sup n | x [ n ] | n {\displaystyle R=\limsup _{n\to \infty }{\sqrt[{n}]{|x[n]|}}}

Di applicazione meno generale è il criterio del rapporto, poiché esso richiede che i termini siano diversi da zero a partire da un n {\displaystyle n} arbitrario in poi. Nondimeno, spesso è più agevole calcolare il limite tramite tale criterio piuttosto che utilizzando quello della radice. Nel caso entrambi i limiti esistano, essi coincidono. Non bisogna tuttavia prendere il reciproco del limite superiore, in quanto la trasformata zeta unilatera è una serie di potenze con esponente negativo.

Trasformata bilatera

Talvolta, può essere utile definire la trasformata di una successione x [ n ] {\displaystyle x[n]} indicizzata su n Z {\displaystyle n\in \mathbb {Z} } . In tal caso, la sua trasformata bilatera è definita come la serie formale di potenze

X ( z ) = Z { x [ n ] } = n = x [ n ] z n {\displaystyle X(z)={\mathcal {Z}}\{x[n]\}=\sum _{n=-\infty }^{\infty }x[n]z^{-n}}

dove di nuovo z {\displaystyle z} è complesso.

Formula di inversione

L'espressione della trasformata inversa, che può essere ottenuta utilizzando il teorema integrale di Cauchy, è la seguente:

x [ n ] = Z 1 { X ( z ) } = 1 2 π j C X ( z ) z n 1 d z , n N . {\displaystyle x[n]={\mathcal {Z}}^{-1}\{X(z)\}={\frac {1}{2\pi j}}\oint _{C}X(z)z^{n-1}dz,\qquad n\in \mathbb {N} .}

dove C {\displaystyle C} è un percorso antiorario chiuso che è situato nella regione di convergenza di X ( z ) {\displaystyle X(z)} e circonda l'origine del piano. La formula precedente diventa particolarmente utile quando X ( z ) {\displaystyle X(z)} ammette un'estensione a tutto il piano complesso, tranne al più un numero finito di singolarità isolate z 1 , , z {\displaystyle z_{1},\dots ,z_{\ell }} . Infatti, in tal caso si può fare appello al Teorema dei Residui ed ottenere

x [ n ] = j = 1 R e s ( X ( z ) z n 1 , z j ) ,  per ogni  n N {\displaystyle x[n]=\sum _{j=1}^{\ell }\mathrm {Res} (X(z)\,z^{n-1},z_{j}),\qquad {\mbox{ per ogni }}n\in \mathbb {N} }

Inoltre, nel caso in cui le singolarità isolate z 1 , , z {\displaystyle z_{1},\dots ,z_{\ell }} siano dei poli, il calcolo dei residui nella formula precedente risulta particolarmente agevole, usando la formula

R e s ( X ( z ) z n 1 , z j ) = 1 ( m j 1 ) ! lim z z j d m j 1 d z m j 1 ( X ( z ) z n 1 ( z z j ) m j ) {\displaystyle \mathrm {Res} (X(z)\,z^{n-1},z_{j})={\frac {1}{(m_{j}-1)!}}\,\lim _{z\to z_{j}}{\frac {d^{m_{j}-1}}{dz^{m_{j}-1}}}{\Big (}X(z)\,z^{n-1}(z-z_{j})^{m_{j}}{\Big )}}

ove m j {\displaystyle m_{j}} è l'ordine del polo z j {\displaystyle z_{j}} .

Un caso di particolare importanza si presenta quando C {\displaystyle C} è la circonferenza unitaria. In tal caso la trasformata zeta inversa assume la forma della trasformata di Fourier discreta inversa:

x [ n ] = 1 2 π π + π X ( e j ω ) e j ω n d ω .   {\displaystyle x[n]={\frac {1}{2\pi }}\int _{-\pi }^{+\pi }X(e^{j\omega })e^{j\omega n}d\omega .\ }

Proprietà

Dominio del tempo Dominio Z Dimostrazione ROC
Notazione x [ n ] = Z 1 { X ( z ) } {\displaystyle x[n]={\mathcal {Z}}^{-1}\{X(z)\}} X ( z ) = Z { x [ n ] } {\displaystyle X(z)={\mathcal {Z}}\{x[n]\}} ROC: r 2 < | z | < r 1   {\displaystyle r_{2}<|z|<r_{1}\ }
Linearità a 1 x 1 [ n ] + a 2 x 2 [ n ]   {\displaystyle a_{1}x_{1}[n]+a_{2}x_{2}[n]\ } a 1 X 1 ( z ) + a 2 X 2 ( z )   {\displaystyle a_{1}X_{1}(z)+a_{2}X_{2}(z)\ } X ( z ) = n = ( a 1 x 1 [ n ] + a 2 x 2 [ n ] ) z n = a 1 n = x 1 [ n ] z n + a 2 n = x 2 [ n ] z n = a 1 X 1 ( z ) + a 2 X 2 ( z ) {\displaystyle {\begin{array}{ll}X(z)&=\sum _{n=-\infty }^{\infty }(a_{1}\,x_{1}[n]+a_{2}\,x_{2}[n])\,z^{-n}\\&=a_{1}\sum _{n=-\infty }^{\infty }x_{1}[n]\,z^{-n}\\&+a_{2}\sum _{n=-\infty }^{\infty }x_{2}[n]\,z^{-n}\\&=a_{1}\,X_{1}(z)+a_{2}\,X_{2}(z)\end{array}}} Almeno la regione di intersezione di ROC1 e ROC2
Espansione temporale x ( k ) [ n ] = { x [ r ] , n = r k 0 , n r k {\displaystyle x_{(k)}[n]={\begin{cases}x[r],&n=rk\\0,&n\not =rk\end{cases}}}

r {\displaystyle r} intero

X ( z k )   {\displaystyle X(z^{k})\ } X k ( z ) = n = x k [ n ] z n = r = x [ r ] z r k = r = x [ r ] ( z k ) r = X ( z k ) {\displaystyle {\begin{array}{ll}X_{k}(z)&=\sum _{n=-\infty }^{\infty }x_{k}[n]\,z^{-n}\\&=\sum _{r=-\infty }^{\infty }x[r]\,z^{-r\,k}\\&=\sum _{r=-\infty }^{\infty }x[r]\,(z^{k})^{-r}\\&=X(z^{k})\end{array}}} r 1 / k {\displaystyle r^{1/k}}
Traslazione temporale x [ n k ] {\displaystyle x[n-k]} z k X ( z ) {\displaystyle z^{-k}X(z)} Z { x [ n k ] } = n = 0 x [ n k ] z n {\displaystyle Z\{x[n-k]\}=\sum _{n=0}^{\infty }x[n-k]z^{-n}}

Posto j = n k {\displaystyle j=n-k} si ha:
n = 0 x [ n k ] z n = j = k x [ j ] z ( j + k ) = j = k x [ j ] z j z k {\displaystyle \sum _{n=0}^{\infty }x[n-k]z^{-n}=\sum _{j=-k}^{\infty }x[j]z^{-(j+k)}=\sum _{j=-k}^{\infty }x[j]z^{-j}z^{-k}} = z k j = k x [ j ] z j {\displaystyle =z^{-k}\sum _{j=-k}^{\infty }x[j]z^{-j}} = z k j = 0 x [ j ] z j {\displaystyle =z^{-k}\sum _{j=0}^{\infty }x[j]z^{-j}}
essendo x [ β ] = 0 {\displaystyle x[\beta ]=0} se β < 0 {\displaystyle \beta <0} . Da cui:
z k j = 0 x [ j ] z j = z k X ( z ) {\displaystyle z^{-k}\sum _{j=0}^{\infty }x[j]z^{-j}=z^{-k}X(z)}

ROC, eccetto z = 0   {\displaystyle z=0\ } se k > 0 {\displaystyle k>0\,} e z = {\displaystyle z=\infty } se k < 0   {\displaystyle k<0\ }
Segnali periodici x [ n + m ] = x [ n ] {\displaystyle x[n+m]=x[n]} X ( z ) = z m z m 1 k = 0 m 1 x [ k ] z k {\displaystyle X(z)={\frac {z^{m}}{z^{m}-1}}\,\sum _{k=0}^{m-1}{\frac {x[k]}{z^{k}}}}
Scalatura nel dominio z a n x [ n ]   {\displaystyle a^{n}x[n]\ } X ( a 1 z )   {\displaystyle X(a^{-1}z)\ } Z { a n x [ n ] } = n = a n x ( n ) z n = n = x ( n ) ( a 1 z ) n = X ( a 1 z ) {\displaystyle {\begin{array}{lcl}Z\{a^{n}x[n]\}=&\\\sum _{n=-\infty }^{\infty }a^{n}x(n)z^{-n}&\\=\sum _{n=-\infty }^{\infty }x(n)(a^{-1}z)^{-n}&\\=X(a^{-1}z)&\\\end{array}}} | a | r 2 < | z | < | a | r 1   {\displaystyle |a|r_{2}<|z|<|a|r_{1}\ }
Inversione temporale x [ n ]   {\displaystyle x[-n]\ } X ( z 1 )   {\displaystyle X(z^{-1})\ } Z { x ( n ) } = n = x ( n ) z n   = m = x ( m ) z m   = m = x ( m ) ( z 1 ) m   = X ( z 1 ) {\displaystyle {\begin{array}{lcl}{\mathcal {Z}}\{x(-n)\}=&\\\sum _{n=-\infty }^{\infty }x(-n)z^{-n}\ &\\=\sum _{m=-\infty }^{\infty }x(m)z^{m}\ &\\=\sum _{m=-\infty }^{\infty }x(m){(z^{-1})}^{-m}\ &\\=X(z^{-1})&\\\end{array}}} 1 r 1 < | z | < 1 r 2   {\displaystyle {\frac {1}{r_{1}}}<|z|<{\frac {1}{r_{2}}}\ }
Coniugazione complessa x [ n ]   {\displaystyle x^{*}[n]\ } X ( z )   {\displaystyle X^{*}(z^{*})\ } Z { x ( n ) } = n = x ( n ) z n   = n = [ x ( n ) ( z ) n ]   = [ n = x ( n ) ( z ) n   ] = X ( z ) {\displaystyle {\begin{array}{lcl}Z\{x^{*}(n)\}=&\\\sum _{n=-\infty }^{\infty }x^{*}(n)z^{-n}\ &\\=\sum _{n=-\infty }^{\infty }[x(n)(z^{*})^{-n}]^{*}\ &\\=[\sum _{n=-\infty }^{\infty }x(n)(z^{*})^{-n}\ ]^{*}&\\=X^{*}(z^{*})&\\\end{array}}} ROC
Parte reale Re { x [ n ] }   {\displaystyle \operatorname {Re} \{x[n]\}\ } 1 2 [ X ( z ) + X ( z ) ] {\displaystyle {\frac {1}{2}}\left[X(z)+X^{*}(z^{*})\right]} ROC
Parte immaginaria Im { x [ n ] }   {\displaystyle \operatorname {Im} \{x[n]\}\ } 1 2 j [ X ( z ) X ( z ) ] {\displaystyle {\frac {1}{2j}}\left[X(z)-X^{*}(z^{*})\right]} ROC
Differenziazione n x [ n ]   {\displaystyle nx[n]\ } z d X ( z ) d z {\displaystyle -z{\frac {dX(z)}{dz}}} Z { n x ( n ) } = n = n x ( n ) z n   = z n = n x ( n ) z n 1   = z n = x ( n ) ( n z n 1 )   = z n = x ( n ) d d z ( z n )   = z d X ( z ) d z {\displaystyle {\begin{array}{lcl}Z\{nx(n)\}=&\\\sum _{n=-\infty }^{\infty }nx(n)z^{-n}\ &\\=z\sum _{n=-\infty }^{\infty }nx(n)z^{-n-1}\ &\\=-z\sum _{n=-\infty }^{\infty }x(n)(-nz^{-n-1})\ &\\=-z\sum _{n=-\infty }^{\infty }x(n){\frac {d}{dz}}(z^{-n})\ &\\=-z{\frac {dX(z)}{dz}}&\\\end{array}}} ROC
Convoluzione x 1 [ n ] x 2 [ n ]   {\displaystyle x_{1}[n]*x_{2}[n]\ } X 1 ( z ) X 2 ( z )   {\displaystyle X_{1}(z)X_{2}(z)\ } Z { x 1 ( n ) x 2 ( n ) } = Z { l = x 1 ( l ) x 2 ( n l ) }   = n = [ l = x 1 ( l ) x 2 ( n l ) ] z n   = l = x 1 ( l ) n = x 2 ( n l ) z n ]   = [ l = x 1 ( l ) z l ] [ n = x 2 ( n ) z n ]   = X 1 ( z ) X 2 ( z ) {\displaystyle {\begin{array}{lcl}{\mathcal {Z}}\{x_{1}(n)*x_{2}(n)\}=&\\{\mathcal {Z}}\{\sum _{l=-\infty }^{\infty }x_{1}(l)x_{2}(n-l)\}\ &\\=\sum _{n=-\infty }^{\infty }[\sum _{l=-\infty }^{\infty }x_{1}(l)x_{2}(n-l)]z^{-n}\ &\\=\sum _{l=-\infty }^{\infty }x_{1}(l)\sum _{n=-\infty }^{\infty }x_{2}(n-l)z^{-n}]\ &\\=[\sum _{l=-\infty }^{\infty }x_{1}(l)z^{-l}][\sum _{n=-\infty }^{\infty }x_{2}(n)z^{-n}]\ &\\=X_{1}(z)X_{2}(z)&\\\end{array}}} Almeno la regione di intersezione di ROC1 e ROC2
Cross-correlazione r x 1 , x 2 = x 1 [ n ] x 2 [ n ]   {\displaystyle r_{x_{1},x_{2}}=x_{1}^{*}[-n]*x_{2}[n]\ } R x 1 , x 2 ( z ) = X 1 ( 1 / z ) X 2 ( z )   {\displaystyle R_{x_{1},x_{2}}(z)=X_{1}^{*}(1/z^{*})X_{2}(z)\ } Almeno la regione di intersezione di ROC of X 1 ( 1 / z ) {\displaystyle X_{1}(1/z^{*})} e X 2 ( z ) {\displaystyle X_{2}(z)}
Prima differenza x [ n ] x [ n 1 ]   {\displaystyle x[n]-x[n-1]\ } ( 1 z 1 ) X ( z )   {\displaystyle (1-z^{-1})X(z)\ } Almeno la regione di intersezione di ROC of X1(z) e | z | > 0 {\displaystyle |z|>0}
Accumulazione k = n x [ k ]   {\displaystyle \sum _{k=-\infty }^{n}x[k]\ } 1 1 z 1 X ( z ) {\displaystyle {\frac {1}{1-z^{-1}}}X(z)} n = k = n x [ k ] z n = n = ( x [ n ] + x [ n 1 ] + x [ n 2 ] x [ ] ) z n = X [ z ] ( 1 + z 1 + z 2 + z 3 ) = X [ z ] j = 0 z j = X [ z ] 1 1 z 1 {\displaystyle {\begin{array}{lcl}\sum _{n=-\infty }^{\infty }\sum _{k=-\infty }^{n}x[k]\cdot z^{-n}\\=\sum _{n=-\infty }^{\infty }(x[n]+x[n-1]+\\x[n-2]\cdots x[-\infty ])z^{-n}\\=X[z](1+z^{-1}+z^{-2}+z^{-3}\cdots )\\=X[z]\sum _{j=0}^{\infty }z^{-j}\\=X[z]{\frac {1}{1-z^{-1}}}\end{array}}}
Moltiplicazione x 1 [ n ] x 2 [ n ]   {\displaystyle x_{1}[n]x_{2}[n]\ } 1 j 2 π C X 1 ( v ) X 2 ( z v ) v 1 d v   {\displaystyle {\frac {1}{j2\pi }}\oint _{C}X_{1}(v)X_{2}({\frac {z}{v}})v^{-1}\mathrm {d} v\ } -
Teorema di Parseval n = x 1 [ n ] x 2 [ n ]   {\displaystyle \sum _{n=-\infty }^{\infty }x_{1}[n]x_{2}^{*}[n]\ } 1 j 2 π C X 1 ( v ) X 2 ( 1 v ) v 1 d v   {\displaystyle {\frac {1}{j2\pi }}\oint _{C}X_{1}(v)X_{2}^{*}({\frac {1}{v^{*}}})v^{-1}\mathrm {d} v\ }

Teorema del valore iniziale e del valore finale

Analogamente alla trasformata di Laplace, anche per la trasformata zeta si possono enunciare due teoremi che permettono di conoscere il valore iniziale e il valore finale del campionamento partendo dalla sua trasformata.

Il teorema del valore iniziale afferma che:

x [ 0 ] = lim z X ( z )   {\displaystyle x[0]=\lim _{z\rightarrow \infty }X(z)\ }

se x [ n ] {\displaystyle x[n]} è causale (ovvero nulla per n negativi).

Se la successione x [ n ] {\displaystyle x[n]} ammette limite finito, allora X ( z ) {\displaystyle X(z)} è una funzione analitica all'esterno del disco di raggio 1 {\displaystyle 1} centrato nell'origine e il teorema del valore finale afferma che:

x [ ] = lim R z 1 + ( z 1 ) X ( z )   {\displaystyle x[\infty ]=\lim _{\mathbb {R} \ni z\rightarrow 1^{+}}(z-1)\,X(z)\ }

Il risultato è falso senza l'ipotesi che x [ n ] {\displaystyle x[n]} ammetta limite, come si vede facilmente prendendo la successione x [ n ] = ( 1 ) n {\displaystyle x[n]=(-1)^{n}} , la cui trasformata zeta è data da

X ( z ) = z z + 1 {\displaystyle X(z)={\frac {z}{z+1}}}

Trasformata di alcune funzioni notevoli

Siano:

  • u [ n ] = { 1 , n 0 0 , n < 0 {\displaystyle u[n]={\begin{cases}1,&n\geq 0\\0,&n<0\end{cases}}}
  • δ [ n ] = { 1 , n = 0 0 , n 0 {\displaystyle \delta [n]={\begin{cases}1,&n=0\\0,&n\neq 0\end{cases}}}
Funzione, x [ n ] {\displaystyle x[n]} Trasformata Z, X ( z ) {\displaystyle X(z)} ROC
δ [ n ] {\displaystyle \delta [n]\,} 1 {\displaystyle 1\,} ogni  z {\displaystyle {\mbox{ogni }}z\,}
δ [ n n 0 ] {\displaystyle \delta [n-n_{0}]\,} z n 0 {\displaystyle z^{-n_{0}}\,} z 0 {\displaystyle z\neq 0\,}
u [ n ] {\displaystyle u[n]\,} 1 1 z 1 {\displaystyle {\frac {1}{1-z^{-1}}}} | z | > 1 {\displaystyle |z|>1\,}
e α n u [ n ] {\displaystyle \,e^{-\alpha n}u[n]} 1 1 e α z 1 {\displaystyle 1 \over 1-e^{-\alpha }z^{-1}} | z | > | e α | {\displaystyle |z|>|e^{-\alpha }|\,}
u [ n 1 ] {\displaystyle -u[-n-1]\,} 1 1 z 1 {\displaystyle {\frac {1}{1-z^{-1}}}} | z | < 1 {\displaystyle |z|<1\,}
n u [ n ] {\displaystyle nu[n]\,} z 1 ( 1 z 1 ) 2 {\displaystyle {\frac {z^{-1}}{(1-z^{-1})^{2}}}} | z | > 1 {\displaystyle |z|>1\,}
n u [ n 1 ] {\displaystyle -nu[-n-1]\,} z 1 ( 1 z 1 ) 2 {\displaystyle {\frac {z^{-1}}{(1-z^{-1})^{2}}}} | z | < 1 {\displaystyle |z|<1\,}
n 2 u [ n ] {\displaystyle n^{2}u[n]\,} z 1 ( 1 + z 1 ) ( 1 z 1 ) 3 {\displaystyle {\frac {z^{-1}(1+z^{-1})}{(1-z^{-1})^{3}}}} | z | > 1 {\displaystyle |z|>1\,}
n 2 u [ n 1 ] {\displaystyle -n^{2}u[-n-1]\,} z 1 ( 1 + z 1 ) ( 1 z 1 ) 3 {\displaystyle {\frac {z^{-1}(1+z^{-1})}{(1-z^{-1})^{3}}}} | z | < 1 {\displaystyle |z|<1\,}
n 3 u [ n ] {\displaystyle n^{3}u[n]\,} z 1 ( 1 + 4 z 1 + z 2 ) ( 1 z 1 ) 4 {\displaystyle {\frac {z^{-1}(1+4z^{-1}+z^{-2})}{(1-z^{-1})^{4}}}} | z | > 1 {\displaystyle |z|>1\,}
n 3 u [ n 1 ] {\displaystyle -n^{3}u[-n-1]\,} z 1 ( 1 + 4 z 1 + z 2 ) ( 1 z 1 ) 4 {\displaystyle {\frac {z^{-1}(1+4z^{-1}+z^{-2})}{(1-z^{-1})^{4}}}} | z | < 1 {\displaystyle |z|<1\,}
a n u [ n ] {\displaystyle a^{n}u[n]\,} 1 1 a z 1 {\displaystyle {\frac {1}{1-az^{-1}}}} | z | > | a | {\displaystyle |z|>|a|\,}
a n u [ n 1 ] {\displaystyle -a^{n}u[-n-1]\,} 1 1 a z 1 {\displaystyle {\frac {1}{1-az^{-1}}}} | z | < | a | {\displaystyle |z|<|a|\,}
n a n u [ n ] {\displaystyle na^{n}u[n]\,} a z 1 ( 1 a z 1 ) 2 {\displaystyle {\frac {az^{-1}}{(1-az^{-1})^{2}}}} | z | > | a | {\displaystyle |z|>|a|\,}
n a n u [ n 1 ] {\displaystyle -na^{n}u[-n-1]\,} a z 1 ( 1 a z 1 ) 2 {\displaystyle {\frac {az^{-1}}{(1-az^{-1})^{2}}}} | z | < | a | {\displaystyle |z|<|a|\,}
n 2 a n u [ n ] {\displaystyle n^{2}a^{n}u[n]\,} a z 1 ( 1 + a z 1 ) ( 1 a z 1 ) 3 {\displaystyle {\frac {az^{-1}(1+az^{-1})}{(1-az^{-1})^{3}}}} | z | > | a | {\displaystyle |z|>|a|\,}
n 2 a n u [ n 1 ] {\displaystyle -n^{2}a^{n}u[-n-1]\,} a z 1 ( 1 + a z 1 ) ( 1 a z 1 ) 3 {\displaystyle {\frac {az^{-1}(1+az^{-1})}{(1-az^{-1})^{3}}}} | z | < | a | {\displaystyle |z|<|a|\,}
cos ( ω 0 n ) u [ n ] {\displaystyle \cos(\omega _{0}n)u[n]\,} z 2 z cos ( ω 0 ) z 2 2 z cos ( ω 0 ) + 1 {\displaystyle {\frac {z^{2}-z\cos(\omega _{0})}{z^{2}-2z\cos(\omega _{0})+1}}} | z | > 1 {\displaystyle |z|>1\,}
sin ( ω 0 n ) u [ n ] {\displaystyle \sin(\omega _{0}n)u[n]\,} z sin ( ω 0 ) z 2 2 z cos ( ω 0 ) + 1 {\displaystyle {\frac {z\sin(\omega _{0})}{z^{2}-2z\cos(\omega _{0})+1}}} | z | > 1 {\displaystyle |z|>1\,}
a n cos ( ω 0 n ) u [ n ] {\displaystyle a^{n}\cos(\omega _{0}n)u[n]\,} z 2 a z cos ( ω 0 ) z 2 2 a z cos ( ω 0 ) + a 2 {\displaystyle {\frac {z^{2}-az\cos(\omega _{0})}{z^{2}-2az\cos(\omega _{0})+a^{2}}}} | z | > | a | {\displaystyle |z|>|a|\,}
a n sin ( ω 0 n ) u [ n ] {\displaystyle a^{n}\sin(\omega _{0}n)u[n]\,} a z sin ( ω 0 ) z 2 2 a z cos ( ω 0 ) + a 2 {\displaystyle {\frac {az\sin(\omega _{0})}{z^{2}-2az\cos(\omega _{0})+a^{2}}}} | z | > | a | {\displaystyle |z|>|a|\,}

Relazione con la trasformata di Laplace

Lo stesso argomento in dettaglio: Trasformata di Laplace.

La trasformata zeta unilatera è la trasformata di Laplace di un segnale campionato in modo ideale con la sostituzione:

z   = d e f   e s T   {\displaystyle z\ {\stackrel {\mathrm {def} }{=}}\ e^{sT}\ }

dove T = 1 / f s   {\displaystyle T=1/f_{s}\ } è il periodo di campionamento, con f s {\displaystyle f_{s}} la frequenza di campionamento (misurata in campioni per secondo o in hertz).

Sia:

Δ T ( t )   = d e f   n = 0 δ ( t n T ) {\displaystyle \Delta _{T}(t)\ {\stackrel {\mathrm {def} }{=}}\ \sum _{n=0}^{\infty }\delta (t-nT)}

un treno di impulsi e sia:

x q ( t ) = d e f   x ( t ) Δ T ( t ) = x ( t ) n = 0 δ ( t n T ) = n = 0 x ( n T ) δ ( t n T ) = n = 0 x [ n ] δ ( t n T ) {\displaystyle {\begin{aligned}x_{q}(t)&{\stackrel {\mathrm {def} }{=}}\ x(t)\Delta _{T}(t)=x(t)\sum _{n=0}^{\infty }\delta (t-nT)\\&=\sum _{n=0}^{\infty }x(nT)\delta (t-nT)=\sum _{n=0}^{\infty }x[n]\delta (t-nT)\end{aligned}}}

la rappresentazione tempo-continua del segnale x [ n ]   = d e f   x ( n T ) {\displaystyle x[n]\ {\stackrel {\mathrm {def} }{=}}\ x(nT)} ottenuto campionando x ( t ) {\displaystyle x(t)} . La trasformata di Laplace di x q ( t ) {\displaystyle x_{q}(t)} è data da:

X q ( s ) = 0 x q ( t ) e s t d t = 0 n = 0 x [ n ] δ ( t n T ) e s t d t = n = 0 x [ n ] 0 δ ( t n T ) e s t d t = n = 0 x [ n ] e n s T {\displaystyle {\begin{aligned}X_{q}(s)&=\int _{0^{-}}^{\infty }x_{q}(t)e^{-st}\,dt\\&=\int _{0^{-}}^{\infty }\sum _{n=0}^{\infty }x[n]\delta (t-nT)e^{-st}\,dt\\&=\sum _{n=0}^{\infty }x[n]\int _{0^{-}}^{\infty }\delta (t-nT)e^{-st}\,dt\\&=\sum _{n=0}^{\infty }x[n]e^{-nsT}\end{aligned}}}

Si tratta della definizione della trasformata zeta unilatera della funzione tempo-discreta x [ n ]   {\displaystyle x[n]\ } , ovvero:

X ( z ) = n = 0 x [ n ] z n {\displaystyle X(z)=\sum _{n=0}^{\infty }x[n]z^{-n}}

con la sostituzione z e s T {\displaystyle z\leftarrow e^{sT}} . Confrontando le ultime due relazioni si ottiene quindi la relazione tra la trasformata zeta unilatera e la trasformata di Laplace del segnale campionato:

X q ( s ) = X ( z ) | z = e s T {\displaystyle X_{q}(s)=X(z){\Big |}_{z=e^{sT}}}

Relazione tra il piano s e il piano z

Per quanto detto la variabile s può essere riscritta utilizzando la rappresentazione rettangolare come:

z = e s T = e T σ e j T ω = e T σ e j T ( ω + 2 k π T ) k R {\displaystyle z=e^{sT}=e^{T\sigma }e^{jT\omega }=e^{T\sigma }e^{jT(\omega +{\frac {2k\pi }{T}})}\qquad k\in \mathbb {R} }

L'ultima identità deriva dal fatto che l'esponenziale complesso è una funzione periodica di periodo i2π.

Da questa relazione si possono fare alcune considerazioni importanti

  • ogni punto sul piano s la cui parte immaginaria differisce di un multiplo intero della pulsazione di campionamento viene trasformato nello stesso punto sul piano z
  • ogni punto sul piano s appartenente al semipiano negativo viene trasformato in un punto interno alla circonferenza di raggio 1 poiché | z | = e T σ {\displaystyle |z|=e^{T\sigma }}
  • ogni punto sul piano s appartenente al semipiano positivo viene trasformato in un punto esterno alla circonferenza di raggio unitario
  • ogni punto appartenente all'asse immaginario viene trasformato in un punto sulla circonferenza di raggio unitario

In virtù di queste considerazioni ha senso definire anche una striscia primaria e più strisce complementari nel piano s. La striscia primaria comprende tutti i numeri complessi con parte immaginaria compresa tra ± j ω s / 2 {\displaystyle \pm j\omega _{s}/2} , le strisce complementari si ottengono, a partire da quella primaria, per traslazione verticale di un multiplo intero della pulsazione di campionamento. Per quanto detto è possibile far corrispondere ogni punto del piano z con un punto della striscia primaria.

Al pari di quanto avviene nel piano s è possibile, anche nel piano z, tracciare dei luoghi a δ {\displaystyle \delta } e ω {\displaystyle \omega } costante.

Campionamento

Si consideri un segnale tempo-continuo x ( t ) {\displaystyle x(t)} , la cui trasformata è:

L { x ( t ) } X ( s ) 0 x ( t ) e s t d t {\displaystyle L\{x(t)\}\equiv X(s)\equiv \int _{0}^{\infty }{x(t)e^{-st}dt}}

Se x ( t ) {\displaystyle x(t)} è campionato uniformemente con un treno di impulsi in modo da ottenere un segnale discreto x [ k ] = x ( k T ) {\displaystyle x^{*}[k]=x(kT)} (supponendo il processo ideale), allora può essere rappresentato come:

x [ k ] = x ( k T ) = k = 0 x ( t ) δ ( t k T ) {\displaystyle x^{*}[k]=x(kT)=\sum _{k=0}^{\infty }{x(t)\delta (t-kT)}}

dove T {\displaystyle T} è l'intervallo di campionamento. In tale contesto la trasformata di Laplace è data da:

L { x ( k T ) } = X ( s ) = 0 k = 0 x ( t ) . δ ( t k T ) e s t d t = k = 0 x ( k ) . z k , z = e s T L { x ( k T ) } | s = ln ( z ) T = X ( s ) | s = ln ( z ) T = Z { x ( k ) } {\displaystyle {\begin{array}{l l l l l l}L\{x(kT)\}&=&X^{*}(s)&=&\int _{0}^{\infty }{\sum _{k=0}^{\infty }{x(t).\delta (t-kT)}e^{-st}dt}\\&=&\sum _{k=0}^{\infty }{x^{*}(k).z^{-k}},z=e^{sT}\\\left.L\{x(kT)\}\right|_{s={\frac {\ln {(z)}}{T}}}&=&\left.X^{*}(s)\right|_{s={\frac {\ln {(z)}}{T}}}&=&Z\{x^{*}(k)\}\end{array}}}

Trasformata di Fourier a tempo discreto

Lo stesso argomento in dettaglio: Trasformata di Fourier a tempo discreto.

La trasformata di Fourier a tempo discreto è un caso particolare della trasformata zeta:

X ( z ) = n = x [ n ] z n {\displaystyle X(z)=\sum _{n=-\infty }^{\infty }x[n]\,z^{-n}}

che si ottiene ponendo z = e i ω {\displaystyle z=e^{i\omega }\,} . Dal momento che | e i ω | = 1 {\displaystyle |e^{i\omega }|=1\,} , la trasformata di Fourier a tempo discreto è la valutazione della trasformata zeta sul cerchio unitario nel piano complesso.

Modello autoregressivo a media mobile

Lo stesso argomento in dettaglio: Modello autoregressivo a media mobile.

Un sistema basato sul modello autoregressivo a media mobile è rappresentato dall'equazione:

p = 0 N y [ n p ] α p = q = 0 M x [ n q ] β q   {\displaystyle \sum _{p=0}^{N}y[n-p]\alpha _{p}=\sum _{q=0}^{M}x[n-q]\beta _{q}\ }

dove entrambi i membri possono essere divisi per α 0   {\displaystyle \alpha _{0}\ } , se è diversa da zero, normalizzando α 0 = 1   {\displaystyle \alpha _{0}=1\ } . In questo modo l'equazione assume la forma:

y [ n ] = q = 0 M x [ n q ] β q p = 1 N y [ n p ] α p {\displaystyle y[n]=\sum _{q=0}^{M}x[n-q]\beta _{q}-\sum _{p=1}^{N}y[n-p]\alpha _{p}}

Tale scrittura consente di visualizzare il fatto che l'uscita al tempo attuale y [ n ] {\displaystyle y[{n}]} è funzione del valore dell'uscita y [ n p ] {\displaystyle y[{n-p}]} a un tempo precedente, dell'ingresso attuale x [ n ] {\displaystyle x[{n}]} e dei precedenti valori x [ n q ]   {\displaystyle x[{n-q}]\ } . Considerando la trasformata zeta della precedente equazione, dalle proprietà di linearità e traslazione temporale si ha:

Y ( z ) p = 0 N z p α p = X ( z ) q = 0 M z q β q   {\displaystyle Y(z)\sum _{p=0}^{N}z^{-p}\alpha _{p}=X(z)\sum _{q=0}^{M}z^{-q}\beta _{q}\ }

che può essere scritta in modo da evidenziare la funzione di trasferimento:

H ( z ) = Y ( z ) X ( z ) = q = 0 M z q β q p = 0 N z p α p = β 0 + z 1 β 1 + z 2 β 2 + + z M β M α 0 + z 1 α 1 + z 2 α 2 + + z N α N {\displaystyle H(z)={\frac {Y(z)}{X(z)}}={\frac {\sum _{q=0}^{M}z^{-q}\beta _{q}}{\sum _{p=0}^{N}z^{-p}\alpha _{p}}}={\frac {\beta _{0}+z^{-1}\beta _{1}+z^{-2}\beta _{2}+\cdots +z^{-M}\beta _{M}}{\alpha _{0}+z^{-1}\alpha _{1}+z^{-2}\alpha _{2}+\cdots +z^{-N}\alpha _{N}}}}

Dal teorema fondamentale dell'algebra il numeratore ha M radici, corrispondenti agli zeri di H {\displaystyle H} , e il denominatore ha N radici, corrispondenti ai poli di H {\displaystyle H} . Riscrivendo la funzione di trasferimento in modo da evidenziare questo fatto si ha:

H ( z ) = ( 1 q 1 z 1 ) ( 1 q 2 z 1 ) ( 1 q M z 1 ) ( 1 p 1 z 1 ) ( 1 p 2 z 1 ) ( 1 p N z 1 )   {\displaystyle H(z)={\frac {(1-q_{1}z^{-1})(1-q_{2}z^{-1})\cdots (1-q_{M}z^{-1})}{(1-p_{1}z^{-1})(1-p_{2}z^{-1})\cdots (1-p_{N}z^{-1})}}\ }

dove q k   {\displaystyle q_{k}\ } è il k-esimo zero e p k   {\displaystyle p_{k}\ } il k-esimo polo. Se il sistema descritto da H ( z )   {\displaystyle H(z)\ } è pilotato dal segnale X ( z )   {\displaystyle X(z)\ } allora l'uscita è data da Y ( z ) = H ( z ) X ( z ) {\displaystyle Y(z)=H(z)X(z)} .

Note

  1. ^ E. R. Kanasewich, Time sequence analysis in geophysics, 3rd, University of Alberta, 1981, pp. 185–186, ISBN 978-0-88864-074-1.
  2. ^ J. R. Ragazzini and L. A. Zadeh, The analysis of sampled-data systems, in Trans. Am. Inst. Elec. Eng., vol. 71, II, 1952, pp. 225–234.
  3. ^ Cornelius T. Leondes, Digital control systems implementation and computational techniques, Academic Press, 1996, p. 123, ISBN 978-0-12-012779-5.

Bibliografia

  • El Jury Theory and Applications of the z-Transform Method (John Wiley & Sons, NY, 1964)
  • Yutaka Yamamoto Digital Control Wiley Encyclopedia of Electrical and Electronics Engineering, 5, 445–457 (1999). PDF

Voci correlate

Collegamenti esterni

  • (EN) Eric W. Weisstein, Trasformata zeta, su MathWorld, Wolfram Research. Modifica su Wikidata
  • (EN) J. H. Matthews e R. W. Howell Introduction to the Z-transform (California State University, Fullerton)
  • (EN) Springer Encyclopedia of Mathematics: Z-Transform, su eom.springer.de.
Controllo di autoritàLCCN (EN) sh85149533 · J9U (ENHE) 987007534267905171
  Portale Controlli automatici
  Portale Ingegneria
  Portale Matematica