Transformada Z

A Transformada Z é um método operacional útil no tratamento de sistemas (de tempo) discretos. É também de grande importância na análise de sinais digitais e no projeto de sistemas de controle digital. Sabe-se que a transformada de Laplace (1749-1827) tem sido usada desde longa data na solução de equações diferenciais contínuas e invariantes no tempo.

História da transformada Z

Em engenharia, a ideia por trás do nome “Transformada” consiste, basicamente, em uma operação matemática que tem por finalidade promove algum tipo de simplificação. Dessa forma, o logaritmo consiste, provavelmente, na ferramenta mais antiga de que se tem notı́cia cujo conceito se aproxima da ideia de transformada, uma vez que transforma multiplicações e divisões em somas e subtrações, além de ser útil na resolução de equações cujos expoentes são desconhecidos (Boyer, 1974). Em verdade, o conceito das transformadas vai muito além dos logaritmos no contexto da engenharia, em que desempenham papel importante. Entre as mais conhecidas (e com certeza mais utilizadas) figura a Transformada Z.

Entretanto, métodos para o tratamento de problemas de tempo discreto são relativamente recentes. Um método para a resolução de equações de diferenças lineares e invariantes no tempo foi apresentado por Gardner e Barnes aos seus alunos de engenharia no início da década de 1940. Eles aplicaram tal procedimento, que era baseado principalmente em jump functions (funções usadas para representar uma sequência de dados amostrados), na resolução de linhas de transmissão e aplicações envolvendo funções de Bessel. Tal abordagem era bastante complexa, e, na tentativa de "dividir para simplificar", uma transformação de um sinal amostrado foi proposta em 1947 por Witold Hurewicz (1904-1956). Tal transformação era escrita como função da sequência amostrada f (no domı́nio do tempo) ao invés do número complexo z da notação moderna:

Definição 1 (transformada z bilateral)

Seja $x[n]$ definida para $n\in 𝐙$. A Transformada Z bilateral da função $x[n]$ é dada por:

$X(z)=𝒵\{x[n]\}=\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}x[n]z^{-n}=\cdots +x[-2]z^2+x[-1]z+x[0]+x[1]z^{-1}+x[2]z^{-2}+\cdots$

Em 1952, cinco anos após a tentativa de Hurewicz, a transformação foi batizada de Transformada Z pelo Sampled-data control group, liderada por John Ralph Raggazini (1912-1988), com Eliahu Ibrahim Jury (que, na época, era aluno de doutorado de Raggazini, mas que acabou sendo um dos principais desenvolvedores da teoria), Lotfi Zadeh (famoso pela criação da lógica Fuzzy) e colaboradores da Columbia University, com o artigo “The Analysis of sampled-Data Systems” (1952), considerado um dos pioneiros trabalhos sobre a transformada Z. Ao que tudo indica, o termo “Z" foi provavelmente utilizado porque ser relativamente incomum contexto da Engenharia Elétrica (na década de 1950), ainda que remonte ao nome do próprio Zadeh.

Transformada Z Inversa

$x[n]={𝒵}^{-1}\{X(z)\}=\frac{1}{2\pi j}{\oint }_{C}X(z)z^{n-1}dz$

onde

$C$

é qualquer curva fechada contendo a origem de forma que a integral indicada converge.

Região de convergência da tranformada Z

Região de convergência (RDC)

A região de convergência é a parte do plano complexo onde a Transformada converge:

$RDC=\{z:|\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}x[n]z^{-n}|<\mathrm{\infty }\}$

No caso em que $x[n]=0$, para $n<0$, a série converge para valores de $z$ em módulo, maiores que o raio de convergência $R$:

$\left|z\right|>$$\underset{n\to +\mathrm{\infty }}{\mathrm{lim\; sup}}\sqrt[n]{|x[n]|}=$$R$

Portanto, a série converge absolutamente para todos os pontos do plano $z$ que se encontram fora do círculo de raio $R$, centrado na origem. Esta região é denominada região de convergência (RDC ou ROC, da sigla em inglês Region of Convergence).

Propriedades da Transformada Z bilateral

$\begin{array}{cc}& 𝒵\{g[n]\}=G(z)\\ & 𝒵\{h[n]\}=H(z)\\ \end{array}$

Linearidade

$\begin{array}{cc}𝒵\{\alpha g[n]+\beta h[n]\}& =\sum _{n=0}^{\mathrm{\infty }}(\alpha g[n]+\beta h[n])z^{-n}\\ & =\alpha \sum _{n=0}^{\mathrm{\infty }}g[n]z^{-n}+\beta \sum _{n=0}^{\mathrm{\infty }}h[n]z^{-n}\\ & =\alpha 𝒵\{g[n]\}+\beta 𝒵\{h[n]\}\\ & =\alpha G(z)+\beta H(z)\end{array}$

Teorema do valor inicial

$g[0]=\underset{z\to \mathrm{\infty }}{\lim }G(z)$

Teorema do valor final

$\underset{n\to \mathrm{\infty }}{\lim }g[n]=\underset{z\to 1}{\lim }(z-1)G(z)$

Deslocamento temporal

Atraso

Se

$g[n]$

é um sinal discreto, então

$\begin{array}{cc}𝒵\{g[n-n_0]\}& =z^{-n_0}G(z),n_0\ge 0\\ & =\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}g[n-n_0]z^{-n}\end{array}$

Definindo

$m=n-n_0$

$=\sum _{m=-\mathrm{\infty }}^{\mathrm{\infty }}g[m]z^{-m}{z}^{-n_0}=z^{-n_0}G(z)$

Mudança de Escala

$\begin{array}{cc}𝒵\{{\alpha }^{n}g[n]\}& =G\left(\frac{z}{\alpha }\right)\\ \end{array}$

Derivada da Transformada Z

$\begin{array}{cc}𝒵\{ng[n]\}& =\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}ng[n]z^{-n}\\ & =z\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}ng[n]z^{-n-1}\\ & =-z\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}g[n](-n{z}^{-n-1})\\ & =-z\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}g[n]\frac{d}{dz}(z^{-n})\\ & =-z\frac{dG(z)}{dz}\end{array}$

Transformadas das sucessões de senos e co-senos

Consideremos uma função discreta

$y[n]=(a+\mathrm{i}b{)}^{n}u[n]$

onde $a$ e $b$ são constantes reais. Usando o resultado da (Equação), o qual é também válido para números reais já que a série geométrica também converge no plano complexo, obtemos

$\mathbb{Z}\{(a+\mathrm{i}b{)}^{n}u[n]\}=\frac{z}{z-a-\mathrm{i}b}$

multiplicando o numerador e o denominador pelo complexo conjugado do denominador, podemos separar as partes real e imaginária

$\mathbb{Z}\{(a+\mathrm{i}b{)}^{n}u[n]\}=\frac{z(z-a)}{(z-a{)}^2+b^2}+\mathrm{i}\frac{bz}{(z-a{)}^2+b^2}$

Por outro lado, se usarmos a representação polar do número complexo $a+\mathrm{i}b$ e a linearidade da transformada Z, podemos escrever

$\mathbb{Z}\{(a+\mathrm{i}b{)}^{n}u[n]\}=\mathbb{Z}\{r^n{e}^{\mathrm{i}n\theta }u[n]\}=\mathbb{Z}\{r^n\cos (n\theta )u[n]\}+\mathrm{i}\mathbb{Z}\{r^n\sin (n\theta )u[n]\}$

onde $r$ e $\theta$ são o módulo e ângulo polar do número complexo. Comparando as partes reais e imaginárias das duas últimas equações, obtemos as transformadas de duas sucessões com senos e co-senos

$\begin{array}{cc}& \mathbb{Z}\{r^n\sin (n\theta )u[n]\}=\frac{bz}{(z-a{)}^2+b^2}\\ & \mathbb{Z}\{r^n\cos (n\theta )u[n]\}=\frac{z(z-a)}{(z-a{)}^2+b^2}\end{array}$

onde as constantes $a$ e $b$ são definidas por

$\begin{array}{cc}& a\equiv r\cos \theta \\ & b\equiv r\sin \theta \end{array}$

Multiplicação por exponencial

A multiplicação da sequência $y[n]$ por uma sequência exponencial da forma $x[n]=a^n$ corresponde a uma dilatação no domínio de $z$:

$\mathbb{Z}\{a^{n}h[n]\}=H\left(\frac{z}{a}\right)$

Reversão temporal

$\begin{array}{cc}𝒵\{h[-n]\}& =\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}f[-n]z^{-n}\\ & =\sum _{m=-\mathrm{\infty }}^{\mathrm{\infty }}h[m]z^m\\ & =\sum _{m=-\mathrm{\infty }}^{\mathrm{\infty }}h[m]{(z^{-1})}^{-m}\\ & =H(z^{-1})\\ \end{array}$

Convolução em Tempo Discreto

$𝒵\{g[n]*h[n]\}=H(z)G(z)$

Transformada da Derivada

$𝒵\{g[n]-g[n-1]\}=(1-z^{-1})G(z)$

Significado Físico da transformada Z

Seja

$H(z)=z^{-1}$

Função de Transferência

$\begin{array}{cc}U(z)& =z^{-1}E(z)\\ \sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}u[k]z^{-k}& =z^{-1}\sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}e[k]z^{-k}\\ \sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}u[k]z^{-k}& =z^{-1}\sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}e[k-1]z^{-(k-1)}\\ \sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}u[k]z^{-k}& =\sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}e[k-1]z^{-(k-1)}{z}^{-1}\\ \sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}u[k]z^{-k}& =\sum _{k=-\mathrm{\infty }}^{\mathrm{\infty }}e[k-1]z^{-k}\\ u[k]& =e[k-1]\end{array}$

Assim, a função transferência

$z^{-1}$

significa o sinal de entrada atrasado por um período de amostragem, como mostra a figura ao lado

Relação com a Transformada de Laplace

$\begin{array}{cc}x[n]=x(nT)& =\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}x(t)\delta (t-nT)\\ & =\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}x(nT)\delta (t-nT)\\ \end{array}$

onde

$T$

é o tempo de amostragem. A Transformada de Laplace

$X(s)$

do sinal

$x(n)$

é:

$\begin{array}{cc}X(s)& =\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}(\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}x(nT)\delta (t-nT))e^{-st}dt\\ & =\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}x(nT)\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}(\delta (t-nT))e^{-st}dt\\ & =\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}x(nT)e^{-nTs}\\ \end{array}$

Obtemos assim a definição de Transformada Z como a Transformada de Laplace com a mudança de variável

$z=e^{Ts}$

.

$\begin{array}{cc}X(s)& =\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}x(nT)e^{-nTs}\\ & =\sum _{n=-\mathrm{\infty }}^{\mathrm{\infty }}f(nT)z^{-n}\end{array}$

Definição (transformada z unilateral)

Seja $x[n]$ definida para $n\in 𝐙$. A Transformada Z uniateral da função $x[n]$ é dada por:

$X(z)={𝒵}_U\{x[n]\}=\sum _{n=0}^{\mathrm{\infty }}x[n]z^{-n}=x[0]+x[1]z^{-1}+x[2]z^{-2}+\cdots$

Propriedades da Transformada Z unilateral

$\begin{array}{cc}& {𝒵}_U\{g[n]\}=G_U(z)\\ & {𝒵}_U\{h[n]\}=H_U(z)\\ \end{array}$

Deslocamento temporal

$\begin{array}{cc}{𝒵}_U\{g[n-1]\}& =z^{-1}{G}_U(z)+g[-1],\\ {𝒵}_U\{g[n-2]\}& =z^{-2}{G}_U(z)+z^{-1}g[-1]+g[-2],\\ {𝒵}_U\{g[n-3]\}& =z^{-3}{G}_U(z)+z^{-2}g[-1]+z^{-1}g[-2]+g[-3],\\ ⋮\end{array}$

$\begin{array}{cc}{𝒵}_U\{g[n+1]\}& =z{G}_U(z)-zg[0],\\ {𝒵}_U\{g[n+2]\}& =z^2{G}_U(z)-z^{2}g[0]-zg[1],\\ {𝒵}_U\{g[n+3]\}& =z^3{G}_U(z)-z^{3}g[0]-z^{2}g[1]-zg[2],\\ ⋮\end{array}$

Aplicações da Transformada Z

Resolução de equações de diferenças lineares não homogêneas

A propriedade de deslocamento no tempo (atraso ou avanço) da transformada z unilateral é empregada para resolução de equações de diferenças lineares com coeficientes constantes. Converte-se equações de diferenças em equações algébricas e encontram-se as soluções no domínio-z. A transformada inversa determina a solução no domínio do tempo.

$y[n+2]+3y[n+1]+2y[n]=3^n,n\ge 0,\text{com}y[0]=1y[1]=0$

Usando a expressão que obtivemos para a transformada de $y[n]$, podemos escrever

$\begin{array}{cc}& \mathbb{Z}\{y[n+1]\}=z\overline{y}-z\\ & \mathbb{Z}\{y[n+2]\}=z\mathbb{Z}\{y[n+1]\}-zy[1]=z^2\overline{y}-z^2\end{array}$

vemos que

$\mathbb{Z}\{3^{n}u[n]\}=\frac{z}{z-3}$

Assim, a transformada da equação de diferenças será

$(z^2+3z+2)\overline{y}-z^2-3z=\frac{z}{z-3}$

e daí obtemos

$\overline{y}=\frac{z^2+3z}{(z+1)(z+2)}+\frac{z}{(z+1)(z+2)(z-3)}$

O cálculo da transformada inversa é feito em forma análoga as transformadas inversas de Laplace, usando expansão em frações parciais, mas deixando de fora um fator $z$ no numerador, que será necessário manter em todas as frações parciais.

Consequentemente, as frações que deverão ser expandidas são:

$\begin{array}{cc}& \frac{z+3}{(z+1)(z+2)}=\frac{2}{z+1}-\frac{1}{z+2}\\ & \frac{1}{(z+1)(z+2)(z-3)}=\frac{1}{5(z+2)}-\frac{1}{4(z+1)}+\frac{1}{20(z-3)}\end{array}$

multiplicando cada fração parcial pelo fator $z$ que deixamos de fora, obtemos o lado direito da (Equação)

$\overline{y}=\frac{7z}{4(z+1)}-\frac{4z}{5(z+2)}+\frac{z}{20(z-3)}$

assi,, encontramos a solução do problema de valores iniciais

$y[n]=(\frac{7}{4}(-1{)}^n-\frac{4}{5}(-2{)}^n+\frac{3^n}{20})u[n]$

Circuito elétrico

Circuito escada simétrico de (k+1) malhas

Considerando um circuito elétrico de forma escada, conforme figura ao lado, constituído por (k+1) malhas fechadas.

Sendo todas as resistências iguais, $R=\frac{V}{i}$, onde $V$ represente a tensão elétrica medida sobre a resistência e $i$, a intensidade de corrente que passa nessa resistência.

Analisando a primeira malha, tem-se que:

$V=R{i}_0+R(i_0-i_1)$

, ou ainda,

Primeira malha do circuito

$i_1=2i_0-\frac{V}{R}$

Da segunda malha segue:

$R(i_1-i_0)+R{i}_1+R(i_1-i_2)=0$

, ou seja

$3i_1-i_0-i_2=0$

Assim percebe-se que, não é preciso conhecer

$V$

para obter

$i_2$

, pois

$i_2=3i_1-i_0$

k_ésima malha do circuito

Tem-se ainda, uma série:

$\begin{array}{c}{i}_{n+2}=3i_{n+1}-i_n\\ i_{n+2}-3i_{n+1}+i_n=0\end{array}$

, que é uma equação de diferenças de segunda ordem, cuja solução da o valor de

$i_n$

para qualquer elemento

$n$

do circuito.Predefinição:Quote

$z^2𝒵[i_n](z)-z^2{i}_0-z{i}_1-3z𝒵[i_n](z)+3z{i}_0+𝒵[i_n](z)=0$

ou seja,

$\begin{array}{cc}(z^2-3z+1)𝒵[i_n](z)& =z^2{i}_0+z{i}_1-3z{i}_0\\ z^2{i}_0+z{i}_1-3z{i}_0& =i_0(z^2-3z)+(2i_0-\frac{V}{R})z\\ & =i_0(z^2-z)-z\frac{V}{R}\\ & =i_0(z^2-z((1+\frac{V}{i_{0}R}))\end{array}$

Assim,

$𝒵[i_n](z)=\frac{z^2-z(1+\frac{V}{i_{0}R})}{z^2-3z+1}{i}_0$

Consequentemente,

$i_n={𝒵}^{-1}\left[\frac{z^2-z(1+\frac{V}{i_{0}R})}{z^2-3z+1}{i}_0\right]$

$\begin{array}{cc}{i}_n& ={𝒵}^{-1}[\frac{z^2-z\cosh (w)}{z^2-2z\cosh (w)+1}{i}_0+\frac{z\cosh (w)-z(1+\frac{V}{i_{0}R})}{z^2-2z\cosh (w)+1}]\\ & =i_0\cosh (nw)+{𝒵}^{-1}\left[\frac{z(\frac{3}{2}-1-\frac{V}{i_{0}R})}{z^2-2z\cosh (w)+1}\right]\\ & =i_0\cosh (nw)+(\frac{1}{2}-\frac{V}{i_{0}R})i_0{𝒵}^{-1}\left[\frac{\frac{2}{\sqrt{5}}z\sinh (w)}{z^2-2z\cosh (w)+1}\right]\\ & =i_0\cosh (nw)+(\frac{1}{2}-\frac{V}{i_{0}R})i_0\frac{2}{\sqrt{5}}\sinh (nw)\end{array}$

$i_n=i_0\cosh (nw)+(\frac{i_0}{\sqrt{5}}-\frac{V}{R\sqrt{5}})\sinh (nw)$

Tabela de Transformadas Z selecionadas

A tabela a seguir provê as transformadas Z para as funções mais comuns de uma variável. Para definições e exemplos, veja a nota explicatória no fim da tabela.

Função	$f[n]={𝒵}^{-1}\{F(z)\}$	$F(z)=𝒵\{f[n]\}$
impulso unitário	$\delta [n]$	$1$
impulso atrasado	$\delta [n-n_0]$	$z^{-n_0}$
degrau unitário	$u[n]$	$\frac{z}{z-1}$
rampa	$n\cdot u[n]$	$\frac{z}{(z-1{)}^2}$
rampa quadrática	$n^2\cdot u[n]$	$\frac{z(z+1)}{(z-1{)}^3}$
rampa cúbica	$n^3\cdot u[n]$	$\frac{z(z^2+4z+1)}{(z-1{)}^4}$
exponencial	$a^n\cdot u[n]$	$\frac{z}{z-a}$
exponencial atrasada	$a^{n-1}\cdot u[n-1]$	$\frac{1}{z-a}$
rampa exponencial	$n\cdot a^n\cdot u[n]$	$\frac{az}{(z-a{)}^2}$
rampa quadrática exponencial	$n^2\cdot a^n\cdot u[n]$	$\frac{az(z+a)}{(z-a{)}^3}$
rampa quadrática exponencial	$\frac{n(n-1)(n-2)...(n-q=1)}{a^q\cdot q!}\cdot a^n\cdot u[n]$	$\frac{z}{(z-a{)}^{q+1}}$
Seno exponencial	$b^n\cdot \sin (\omega n)\cdot u[n]$	$\frac{zb\sin (\omega )}{z^2-2b\cos (\omega )z+b^2}$
Cosseno exponencial (I)	$b^n\cdot \cos (\omega n)\cdot u[n]$	$\frac{z[z-b\cos (\omega )]}{z^2-2b\cos (\omega )z+b^2}$
Cosseno exponencial (II)	$b^n\cdot \cos (\omega n+\theta )\cdot u[n]$	$\frac{z[z\cos (\theta )-b\cos (\omega -\theta )]}{z^2-2b\cos (\omega )z+b^2}$
Cosseno exponencial (III)	$b^n\cdot \cos (\omega n+\theta )\cdot u[n]$	$\frac{e^{j\theta }z}{2(z-b)}-\frac{e^{-j\theta }z}{2(z-b^{*})}$
Cosseno exponencial (IV)	$C{b}^n\cdot \cos (\omega n+\theta )\cdot u[n]$	$\frac{z(Az+B)}{z^2+2az+b^2}$ $C=\sqrt{\frac{A^2{b}^2-2AaB+B^2}{b^2-a^2}}$ $\omega ={\cos }^{-1}(-\frac{a}{b})$ $\theta ={\tan }^{-1}\left(\frac{Aa-B}{A\sqrt{b^2-a^2}}\right)$
Nota explicatória: u[n] representa a Função de Heaviside. $\delta [n]$ representa a Delta de Dirac. ${*}^{}$ denota o conjugado complexo. n, um número inteiro, tipicamente representa tempo discreto, embora possa representar qualquer dimensão independente. Em geral, n > 0. z é a Frequência angular complexa discreta. a, b, ω e θ são números reais, sendo b > 0, 0 < ω < π e 0 < θ < 2π. m e q são números inteiros positivos.