Z-trasformata - Z-transform

In matematica e nell'elaborazione del segnale , la trasformata Z converte un segnale a tempo discreto , che è una sequenza di numeri reali o complessi , in una rappresentazione complessa nel dominio della frequenza .

Può essere considerato come un equivalente a tempo discreto della trasformata di Laplace . Questa somiglianza è esplorata nella teoria del calcolo su scala temporale .

Storia

L'idea di base ora nota come trasformazione Z era nota a Laplace e fu reintrodotta nel 1947 da W. Hurewicz e altri come un modo per trattare i sistemi di controllo dei dati campionati utilizzati con il radar. Fornisce un modo trattabile per risolvere equazioni alle differenze lineari a coefficiente costante . E 'stato poi soprannominato "la z-trasformata" da Ragazzini e Zadeh nel gruppo di controllo dei dati campionati presso la Columbia University nel 1952.

La trasformazione Z modificata o avanzata è stata successivamente sviluppata e resa popolare da EI Jury .

L'idea contenuta all'interno della Z-trasformata è anche conosciuta nella letteratura matematica come il metodo di generazione di funzioni che può essere fatto risalire già al 1730 quando fu introdotto da de Moivre in concomitanza con la teoria della probabilità. Da un punto di vista matematico, la trasformata Z può anche essere vista come una serie di Laurent in cui si vede la sequenza di numeri in esame come l'espansione (di Laurent) di una funzione analitica.

Definizione

La trasformata Z può essere definita come una trasformazione unilaterale o bilaterale .

Trasformata Z bilaterale

Il bilaterale o due lati Z-trasformata di un segnale discreto è le serie formali definiti

 

 

 

 

( Eq.1 )

dove è un intero ed è, in generale, un numero complesso :

dove è la grandezza di , è l' unità immaginaria ed è l' argomento complesso (indicato anche come angolo o fase ) in radianti .

Trasformata Z unilaterale

In alternativa, nei casi in cui è definita solo per il solo lato o unilaterale trasformata z è definito come

 

 

 

 

( Eq.2 )

Nel trattamento del segnale , questa definizione può essere usata per valutare la Z-trasformata di risposta impulso unitario di un tempo discreto sistema causale .

Un esempio importante della trasformata Z unilaterale è la funzione generatrice di probabilità , dove la componente è la probabilità che una variabile casuale discreta assuma il valore , e la funzione è solitamente scritta come in termini di . Le proprietà delle trasformate Z (sotto) hanno interpretazioni utili nel contesto della teoria della probabilità.

Trasformata Z inversa

L' inversa Z-trasformata è

 

 

 

 

( Eq.3 )

dove C è un percorso chiuso in senso antiorario che circonda l'origine e interamente nella regione di convergenza (ROC). Nel caso in cui il ROC sia causale (vedi Esempio 2 ), questo significa che il cammino C deve circondare tutti i poli di .

Un caso speciale di questo integrale di contorno si verifica quando C è il cerchio unitario. Questo contorno può essere utilizzato quando il ROC include il cerchio unitario, che è sempre garantito quando è stabile, cioè quando tutti i poli sono all'interno del cerchio unitario. Con questo contorno, la trasformata Z inversa si semplifica nella trasformata di Fourier a tempo discreto inversa , o serie di Fourier , dei valori periodici della trasformata Z attorno al cerchio unitario:

La trasformata Z con un intervallo finito di n e un numero finito di valori z uniformemente spaziati può essere calcolata in modo efficiente tramite l'algoritmo FFT di Bluestein . La trasformata discreta di Fourier (DTFT) - da non confondere con la trasformata discreta di Fourier (DFT) - è un caso speciale di tale trasformata Z ottenuta limitando z a giacere sul cerchio unitario.

Regione di convergenza

La regione di convergenza (ROC) è l'insieme dei punti nel piano complesso per cui converge la sommatoria in Z-trasformata.

Esempio 1 (nessun ROC)

Sia x[n] = (0.5) n . Espandendo x[n] sull'intervallo (−∞, ∞) diventa

Guardando la somma

Pertanto, non esistono valori di z che soddisfano questa condizione.

Esempio 2 (ROC causale)

ROC mostrato in blu, il cerchio unitario come un cerchio grigio punteggiato e il cerchio | z | = 0,5 è mostrato come un cerchio nero tratteggiato

Let (dove u è la funzione gradino di Heaviside ). Espandendo x[n] sull'intervallo (−∞, ∞) diventa

Guardando la somma

L'ultima uguaglianza deriva dalla serie geometrica infinita e l'uguaglianza vale solo se |0,5 z −1 | < 1 che può essere riscritto in termini di z come | z | > 0,5. Quindi, il ROC è | z | > 0,5. In questo caso il ROC è il piano complesso con un disco di raggio 0,5 all'origine "punzonato".

Esempio 3 (ROC anti causale)

ROC mostrato in blu, il cerchio unitario come un cerchio grigio punteggiato e il cerchio | z | = 0,5 è mostrato come un cerchio nero tratteggiato

Let (dove u è la funzione gradino di Heaviside ). Espandendo x[n] sull'intervallo (−∞, ∞) diventa

Guardando la somma

Usando ancora la serie geometrica infinita , l'uguaglianza vale solo se |0.5 −1 z | < 1 che può essere riscritto in termini di z come | z | < 0,5. Quindi, il ROC è | z | < 0,5. In questo caso il ROC è un disco centrato nell'origine e di raggio 0,5.

Ciò che differenzia questo esempio dall'esempio precedente è solo il ROC. Questo è intenzionale per dimostrare che il risultato della trasformazione da solo è insufficiente.

Esempi conclusione

Gli esempi 2 e 3 mostrano chiaramente che la trasformata Z X(z) di x[n] è unica quando e solo quando si specifica il ROC. La creazione del grafico polo-zero per il caso causale e anticausale mostra che il ROC per entrambi i casi non include il polo che è a 0,5. Questo si estende ai casi con più poli: il ROC non conterrà mai i poli.

Nell'esempio 2, il sistema causale produce un ROC che include | z | = ∞ mentre il sistema anticausale nell'esempio 3 produce un ROC che include | z | = 0.

ROC mostrato come un anello blu 0,5 < | z | < 0,75

In impianti con più poli è possibile avere un ROC che non includa né | z | = ∞ né | z | = 0. Il ROC crea una banda circolare. Per esempio,

ha poli a 0,5 e 0,75. Il ROC sarà 0,5 < | z | < 0,75, che non include né l'origine né l'infinito. Tale sistema è chiamato sistema di causalità mista poiché contiene un termine causale (0.5) n u [ n ] e un termine anticausale −(0.75) n u [− n −1].

La stabilità di un sistema può essere determinata anche conoscendo solo il ROC. Se il ROC contiene il cerchio unitario (cioè | z | = 1) allora il sistema è stabile. Nei suddetti sistemi il sistema causale (Esempio 2) è stabile perché | z | > 0,5 contiene il cerchio unitario.

Supponiamo che ci venga fornita una trasformata Z di un sistema senza ROC (cioè un x[n] ambiguo ). Possiamo determinare un unico x[n] purché desideriamo quanto segue:

  • Stabilità
  • Causalità

Per la stabilità il ROC deve contenere il cerchio unitario. Se abbiamo bisogno di un sistema causale, allora il ROC deve contenere l'infinito e la funzione del sistema sarà una sequenza di lato destro. Se abbiamo bisogno di un sistema anticausale, il ROC deve contenere l'origine e la funzione del sistema sarà una sequenza di lato sinistro. Se abbiamo bisogno sia di stabilità che di causalità, tutti i poli della funzione del sistema devono essere all'interno del cerchio unitario.

L'unico x[n] può quindi essere trovato.

Proprietà

Proprietà della trasformata z
Dominio del tempo Z-dominio Prova ROC
Notazione
Linearità Contiene ROC 1 ∩ ROC 2
Espansione temporale

insieme a

Decimazione ohio-state.edu   o   ee.ic.ac.uk
Ritardo

con e

ROC, eccetto z = 0 se k > 0 ez = ∞ se k < 0
Anticipo del tempo

insieme a

Trasformata Z bilaterale:

Trasformata Z unilaterale:
Prima differenza all'indietro

con x [ n ]=0 per n <0

Contiene l'intersezione di ROC di X 1 (z) e z ≠ 0
Prima differenza in avanti
Inversione del tempo
Ridimensionamento nel dominio z
coniugazione complessa
parte reale
parte immaginaria
Differenziazione ROC, se è razionale;

ROC possibilmente escludendo il confine, se è irrazionale

convoluzione Contiene ROC 1 ∩ ROC 2
Cross-correlazione Contiene l'intersezione di ROC di e
Accumulo
Moltiplicazione -

Teorema di Parseval

Teorema del valore iniziale : se x [ n ] è causale, allora

Teorema del valore finale : se i poli di ( z −1) X ( z ) sono all'interno del cerchio unitario, allora

Tabella delle coppie di trasformate Z comuni

Qui:

è la funzione step unit (o Heaviside) e

è la funzione impulso dell'unità di tempo discreto (cfr. funzione delta di Dirac che è una versione a tempo continuo). Le due funzioni vengono scelte insieme in modo che la funzione a gradino unitario sia l'accumulo (totale parziale) della funzione a impulso unitario.

Segnale, Z-trasformare, ROC
1 1 tutto z
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17 , per intero positivo
18 , per intero positivo
19
20
21
22

Relazione con la serie di Fourier e la trasformata di Fourier

Per valori di nella regione , nota come cerchio unitario , possiamo esprimere la trasformata in funzione di una singola variabile reale, , definendo . E la trasformata bilaterale si riduce a una serie di Fourier :

 

 

 

 

( Eq.4 )

nota anche come trasformata di Fourier a tempo discreto (DTFT) della sequenza. Questo 2 π funzione -periodic è la somma periodica di una trasformata di Fourier , che lo rende uno strumento di analisi ampiamente utilizzato. Per capirlo, sia la trasformata di Fourier di qualsiasi funzione, , i cui campioni ad un certo intervallo, T , sono uguali alla sequenza x[ n ]. Quindi il DTFT della sequenza x [ n ] può essere scritto come segue.

 

 

 

 

( Eq.5 )

Quando T ha unità di secondi, ha unità di hertz . Il confronto delle due serie rivela che     è una frequenza normalizzata con unità di radianti per campione . Il valore ω=2 π corrisponde a Hz. Ed ora, con la sostituzione   Eq.4 può essere espresso in termini della trasformata di Fourier, X(•) : 

 

 

 

 

( Eq.6 )

Al variare del parametro T, i singoli termini dell'Eq.5 si allontanano o si avvicinano lungo l'asse f. In Eq.6 tuttavia, i centri rimangono 2 π a parte, mentre la loro larghezza espandono o contratto. Quando la sequenza x ( nT ) rappresenta la risposta all'impulso di un sistema LTI , queste funzioni sono anche note come risposta in frequenza . Quando la sequenza è periodica, la sua DTFT è divergente a una o più frequenze armoniche e zero a tutte le altre frequenze. Questo è spesso rappresentato dall'uso di funzioni delta di Dirac con variazione di ampiezza alle frequenze armoniche. A causa della periodicità, esiste solo un numero finito di ampiezze uniche, che vengono prontamente calcolate dalla trasformata discreta di Fourier (DFT) molto più semplice . (Vedi DTFT § Dati periodici .)

Relazione con la trasformata di Laplace

Trasformata bilineare

La trasformata bilineare può essere utilizzata per convertire filtri a tempo continuo (rappresentati nel dominio di Laplace) in filtri a tempo discreto (rappresentati nel dominio Z) e viceversa. Viene utilizzata la seguente sostituzione:

convertire una funzione nel dominio di Laplace in una funzione nel dominio Z ( trasformazione di Tustin ), o

dal dominio Z al dominio di Laplace. Attraverso la trasformazione bilineare, il piano s complesso (della trasformata di Laplace) viene mappato sul piano z complesso (della trasformata z). Sebbene questa mappatura sia (necessariamente) non lineare, è utile in quanto mappa l'intero asse del piano s sul cerchio unitario nel piano z. Come tale, la trasformata di Fourier (che è la trasformata di Laplace valutata sull'asse) diventa la trasformata di Fourier a tempo discreto. Ciò presuppone che esista la trasformata di Fourier; cioè, che l' asse è nella regione di convergenza della trasformata di Laplace.

Trasformazione speciale

Data una trasformata Z unilaterale, X(z), di una funzione campionata nel tempo, la corrispondente trasformata con stelle produce una trasformata di Laplace e ripristina la dipendenza dal parametro di campionamento, T :

La trasformata inversa di Laplace è un'astrazione matematica nota come funzione campionata a impulsi .

Equazione di differenza lineare costante-coefficiente

L'equazione della differenza di coefficiente costante lineare (LCCD) è una rappresentazione per un sistema lineare basato sull'equazione della media mobile autoregressiva .

Entrambi i lati dell'equazione sopra possono essere divisi per α 0 , se non è zero, normalizzando α 0 = 1 e l'equazione LCCD può essere scritta

Questa forma dell'equazione LCCD è favorevole per rendere più esplicito che l'uscita "corrente" y[n] è una funzione delle uscite passate y[n−p] , dell'ingresso corrente x[n] e degli ingressi precedenti x[n− q] .

Funzione di trasferimento

Prendendo la Z-trasformata dell'equazione di cui sopra (usando linearità e leggi di time-shifting) si ottiene

e riordinando i risultati in

Zeri e poli

Dal teorema fondamentale dell'algebra il numeratore ha M radici (corrispondenti agli zeri di H) e il denominatore ha N radici (corrispondenti ai poli). Riscrivere la funzione di trasferimento in termini di zeri e poli

dove q k è il k -esimo zero e p k è il k -esimo polo. Gli zeri e i poli sono comunemente complessi e quando vengono tracciati sul piano complesso (piano z) viene chiamato grafico polo-zero .

Inoltre, possono esistere anche zeri e poli in z = 0 ez = ∞. Se prendiamo in considerazione questi poli e zeri, nonché zeri e poli di ordine multiplo, il numero di zeri e poli è sempre uguale.

Scomponendo in fattori il denominatore, è possibile utilizzare la scomposizione parziale della frazione , che può quindi essere riconvertita nel dominio del tempo. Ciò comporterebbe la risposta all'impulso e l'equazione di differenza del coefficiente costante lineare del sistema.

Risposta in uscita

Se un tale sistema H(z) è pilotato da un segnale X(z) allora l'uscita è Y(z) = H(z)X(z) . Eseguendo la scomposizione parziale della frazione su Y(z) e quindi prendendo l'inversa Z-trasformata si può trovare l'output y[n] . In pratica, è spesso utile scomporre in modo frazionario prima di moltiplicare quella quantità per z per generare una forma di Y(z) che ha termini con trasformate Z inverse facilmente calcolabili.

Guarda anche

Riferimenti

Ulteriori letture

  • Refaat El Attar, Appunti delle lezioni su Z-Transform , Lulu Press, Morrisville NC, 2005. ISBN  1-4116-1979-X .
  • Ogata, Katsuhiko, Discrete Time Control Systems 2nd Ed , Prentice-Hall Inc, 1995, 1987. ISBN  0-13-034281-5 .
  • Alan V. Oppenheim e Ronald W. Schafer (1999). Elaborazione del segnale a tempo discreto, 2a edizione, serie di elaborazione del segnale Prentice Hall. ISBN  0-13-754920-2 .

link esterno