36
COQ 790 – ANÁLISE DE SISTEMAS DA ENGENHARIA QUÍMICA AULA 10: Domínio Discreto; Transformada Z. Universidade Federal do Rio de Janeiro COPPE – Programa de Engenharia Química 2014/1

Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

  • Upload
    others

  • View
    3

  • Download
    0

Embed Size (px)

Citation preview

Page 1: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

COQ 790 – ANÁLISE DE SISTEMAS DA ENGENHARIA QUÍMICA

AULA 10:

Domínio Discreto; Transformada Z.

Universidade Federal do Rio de JaneiroCOPPE – Programa de Engenharia Química

2014/1

Page 2: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Introdução ao Domínio Discreto

Em aplicações digitais, o tempo contínuo é observado em intervalos discretos...

y(t)

t

y(t) y(tk)

Amostrador

y(tk)

t0 t1 t2 t3 ... tk

t

Page 3: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Desejamos aproximar f(t) por uma função discreta (amostrada) f*(t), de modo que:

*k kf (t) f(t ), com t k t, k 0,1,2,

A literatura sugere interpretar a amostragem da função f(t) como sendo realizada a partir de um trem de perturbações do tipo impulso aplicadas nos tempos discretos tk, de modo que:

*

k kk 0

f (t) f(t ) (t t )

onde (t-tk) é a perturbação delta de Dirac.

Aplicando a transformação de Laplace à equação acima, obtemos:

Page 4: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

*k k

k 0

stk k

k 00

stk k

k 0 0

stk k

k 0 0

f (t) f(t ) (t t )

e f(t ) (t t ) dt

e f(t ) (t t )dt

f(t ) e (t t )dt

L L

kst

kk 0

f(t )e

Substituindo tk por kt, temos:

* * sk t

kk 0

f (s) f (t) f(t )eL

Page 5: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Vamos agora comparar esse resultado com a transformada de Laplace de f(t):

* sk t

kk 0

f (s) f(t )e

st

0f(s) f(t)e dt Caso contínuo

Caso discreto

Agora, por conveniência, introduzimos a seguinte notação de modo a simplificar a equação acima para o caso discreto:

s tz e

e nos referirmos ao resultado L[f*(t)] como “transformada z” de f*(t), escrevemos:

* * k

kk 0

f (z) f (t) f(t )z

sendo essa a definição formal da transformada z de um sinal amostrado f*(t).

Page 6: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Propriedades e teoremas da transformada z

Linearidade:

11 2 2 1 1 2 2c f (k) c f (k) c f (k) c f (k)

Convolução:

Se f(k) e g(k) são duas sequências discretas em k, que são representações das funções discretas f*(t) e g*(t), e se

f(k) f(z)

g(k) g(z)

são suas correspondentes transformadas z, então:

i 0 i 0

f(i)g(k i) g(i)f(k i) f(z)g(z)Útil para obter

respostas impulsionais em sistemas amostrados.

Page 7: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

k

i 0 k 0 i 0f(i)g(k i) f(i)g(k i)z

Fazendo n=k-i:

(n i)

i 0 i 0 n if(i)g(k i) f(i)g(n)z

Como g(tn) é zero para n<0, o limite inferior de n pode ser mudado para zero:

(n i) i n

i 0 i 0 n 0 i 0 n 0f(i)g(k i) f(i)g(n)z f(i)z g(n)z f(z)g(z)

Page 8: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Transformação inversa:

1 kC

1 dzf(k) f(z) f(z)z2 j z

A fórmula é dificilmente usada, dando lugar a outras técnicas e tabelas. Devemos, entretanto, lembrar que:

A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua original, f(t), a partir da qual a sequência foi obtida;

A inversa não fornece informação sobre o intervalo de amostragem, t, para a sequência discreta f(k).

Page 9: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Deslocamento à esquerda:

Se f(k) f(z) então:

f(k 1) zf(z) zf(0)

Se f(0) 0 então:

f(k 1) zf(z)

E ainda:

2 2f(k 2) z f(z) z f(0) zf(1)

Generalizando:

m m m 1

m 1m k

k 0

f(k m) z f(z) z f(0) z f(1) zf(m 1)

z f(z) f(k)z

Page 10: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Deslocamento à direita:

mf(k m) z f(z)

Translação complexa:

ate f(t) f(s a)L

Lembrando que:

A forma equivalente em transformada z é:

*at a t a te f(t) f(e z) f(rz), com r e

1 kf(rz) r f(k)

Page 11: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Teorema do valor inicial:

*

t 0 k 0 zlim f (t) lim f(k) lim f(z)

Teorema do valor final: |z| < 1

*t k z 1lim f (t) lim f(k) lim z 1 f(z)

Contanto que o valor final da

sequência discreta seja

finito!

t 0 slim f(t) lim s f(s)

t s 0lim f(t) lim s f(s)

Page 12: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Representação gráfica da relação entre as transformadas z e de Laplace

Page 13: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Resolução de Equações Lineares de Diferenças

Como a transformada de Laplace, que é usada na resolução de equações diferenciais ordinárias lineares, a transformada z encontra aplicação na solução de equações lineares de diferenças da forma:

1 2 n 0 1 2 my(k) a y(k 1) a y(k 2) a y(k n) b u(k) b u(k 1) b u(k 2) a u(k m)

Aplicando a transformada z em ambos os lados da equação, obtemos:

1 2 n 1 2 m1 2 n 0 1 2 m1 a z a z a z y(z) b b z b z b z u(z)

De modo que:

1 2 m0 1 2 m

1 2 n1 2 n

b b z b z b zy(z) u(z)

1 a z a z a z

Page 14: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Alguns exemplos

1) Degrau

f(t) a f(k) a

k 1 2 3

k 01

1

f(z) f(k) az a 1 z z z

1 z a a , para z 1z 11 z

2) Rampa

f(t) at f(k) ak t

k 1 1 2

k 01

2 21

f(z) f(k) a t kz a t z 1 2z 3z

z z a t a tz 11 z

Page 15: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

3) Exponencial

btf(t) ae b tkf(k) ae

bk t k

k 02 3b t 1 b t 1 b t 1

b t 1

f(z) f(k) ae z

a 1 e z e z e z

1 a1 e z

Page 16: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Inversão de transformadas z

A transformada z inversa é definida por:

* 1f(k) f (t) f(z) f(z)

A transformada z inversa consiste de valores amostrados f*(t) representados por f(k). A inversão de transformadas z conta com três métodos:

a) Expansão em frações parciaisb) Divisão longac) Integração de linha

a) Expansão em frações parciais

Consideremos que f(z) tem a seguinte forma:

11

1N(z )f(z )D(z )

Page 17: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

onde: N(z) é um polinômio de ordem n em z-1;D(z) é um polinômio de ordem m em z-1.

Admitindo que D(z) possa ser fatorado em m raízes reais distintas (i.e., os pólos de f(z)), denotados por p1, p2, ... pm. Então:

1

11 1 1

1 2 m

N(z )f(z )1 p z 1 p z 1 p z

Expandindo em frações parciais, na forma:

1 1 2 m

1 1 11 2 m

r r rf(z )1 p z 1 p z 1 p z

Cada coeficiente do numerador ri pode ser obtido de maneira similar àquela usada no cálculo de transformadas inversas de Laplace. A transformada inversa é dada, portanto, por:

1 1 11 2 m1 1 1

1 2 m

r r rf(k)1 p z 1 p z 1 p z

Page 18: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Como:

1 kii i1

i

r rp1 p z

Então:

k k k1 1 2 2 m mf(k) r p r p r p

Consideremos, agora, um caso simples:

k1 1f(k) r p

com r1=1 e p1 assumindo diferentes valores.

Page 19: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Real

Imaginário

x(1)

Círculounitário

(1)

x(2)

(2)

x(3)

(3)

(4)

x(4)

x(5)

x(6)

(5)

(6)

Respostas temporais para localizações diferentes do pólo de f(z):

(1)(1)

(2)(2)

(3)

(4)(4)(4)

Page 20: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Respostas temporais para localizações diferentes do pólo de f(z):

Seborg et al. (1989)

Page 21: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Exemplo: Usando expansão em frações parciais, encontre a inversa da função abaixo para um intervalo de amostragem,t=1:

1

11 10.5zf(z )

1 z 1 0.5z

Expandindo em frações parciais:

1

1 1 21 11 1

r r0.5zf(z )1 z 1 0.5z1 z 1 0.5z

Multiplicando por 11 z

e fazendo z=1: 10.5r 10.5

Multiplicando por 11 0.5z

e fazendo z=0.5: 20.5 2r 11 2

11 1

1 1f(z )1 z 1 0.5z

a t ae e 0.5 a 0.693

11 0.693 t 1

1 1f(z )1 z 1 e z

0.693k tf(k) 1 e

Page 22: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

f k( ) 1 e 0.693 k

k 0 10

f k( ) 1 e 0.693 k

k 0 10

0 2 4 6 8 100

0.2

0.4

0.6

0.8

1

f k( )

k

Page 23: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

b) Divisão longa

Da definição da transformada z, podemos escrever f(z) como:

1 k 1 2 3k 0 1 2 3

k 0f(z ) f(t )z f(t ) f(t )z f(t )z f(t )z

definindo uma série infinita em z-1. Assim, para qualquer transformada z, f(z), expandida como uma série infinita em z-1, digamos:

1 1 2 30 1 2 3f(z ) z z z

sua inversa, a sequência discreta f(k), pode ser facilmente recuperada através da comparação dos coeficientes das potências em z-1 na equações acima de modo que:

k kf(t ) , k 1, ,n

Assim, para obter f(k), devemos mostrar como uma função arbitrária f(z) pode ser representada como potências de z-1.

Page 24: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Mais uma vez, consideremos f(z) na seguinte forma:

11

1N(z )f(z )D(z )

1 2 n1 0 1 2 n

1 2 m0 1 2 m

b b z b z b zf(z )a a z a z a z

ou

Poderíamos realizar a divisão acima, e obter uma série infinita (longa) em z-1, mas esse procedimento seria muito tedioso. Uma alternativa, menos tediosa, mas totalmente equivalente, consiste de definir:

1 1 2 30 1 2 3(z ) z z z

De modo que:1

1N(z )1 1 1 1D(z )

(z ) (z )D(z ) N(z )

Page 25: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Na forma expandida:

1 2 3 1 2 m0 1 2 3 0 1 2 m

1 2 n0 1 2 n

z z z a a z a z a z

b b z b z b z

Realizando a multiplicação do lado esquerdo da igualdade acima, e comparando os coeficientes das mesmas potências em z-1, observamos que:

00 0 0 0

0

ba ba

E que (potência z-1): 1 11 0 0 1 1 1 0

0 0

b aa a ba a

(potência z-2): 2 1 22 0 1 1 0 2 2 2 1 0

0 0 0

b a aa a a ba a a

Generalizando:k

0k k i k i 0

0 0i 1

b1 b a , k 1,2, com a a

Page 26: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Resumo do procedimento:

1) Represente f(z) como uma razão de polinômios em z-1, obtenha os coeficientes a0, a1, …, an, e b0, b1, …, bn.

2) Com a condição inicial 0=b0/a0, use a expressão geral obtida acima para calcular os demais 1, 2, 3…

3) Observe que f(k)= k.

• Exemplo: Usando a técnica de expansão em uma série infinita (acima), encontre a transformada inversa de:

11

rf(z )1 pz

Observe que: 0

k

0

1

k

b rb 0, k 1,2,a 1a pa 0, k 2,3,

00

0

1 1 1 00

22 2 1 1 2 0

0

kk 1 k 1

b ra1 b a rp

a1 b a a rp

a

a rp

Page 27: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Logo: 1 2 2 3 3f(z) r 1 pz p z p z

De modo que reconhecemos que:

kf(k) rp

Observe que, fazendo: a tp e f(k) representa a função exponencial discreta

Page 28: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Função de transferência no domínio temporal discreto

Lembremos que desenvolvemos uma expressão para a função de transferência impulsional para o caso contínuo:

t

*

0y(t) g(t )u ( )d

Sendo g(t) a resposta impulsional do processo. u*(t) é uma perturbação discreta ao processo, representada na forma de uma série de impulsos:

*

k kk 0

u ( ) u(t ) ( t )

De modo que:

t

k kk 00

t

k kk 0 0

y(t) g(t ) u(t ) ( t ) d

g(t ) ( t )d u(t )

Page 29: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Continuando:

t

k kk 00

t

k kk 0 0

k kk 0

y(t) g(t ) u(t ) ( t ) d

g(t ) ( t )d u(t )

g(t t )u(t )

Em particular, para t=tn:

n n k k

k 0y(t ) g(t t )u(t )

Observe, agora, que uma função contínua y(t) amostrada em intervalos discretos apresenta a seguinte transformada z, por definição:

n

nn 0

y(z) y(t )z

Page 30: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Logo:

n

n k kn 0 k 0

y(z) g(t t )u(t ) z

Fazendo i=n-k:

(i k)i k k k

i k k 0

i ki k

i k k 0

y(z) g(t t )u(t ) z

g(t )z u(t )z

Agora note que g(ti) é zero para i<0; desse modo, o limite inferior de i pode ser mudado para zero e o somatório separado:

i k

i ki 0 k 0

y(z) g(t )z u(t )z y(z) g(z)u(z)

onde g(z) é definida por:

i

i ii 0

g(z) g(t )z , com t i t

, ti-k tk = ti

Page 31: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

i

i ii 0

g(z) g(t )z , com t i t

g(z) é chamada de função de transferência de pulsos (ou pulso) do sistema, erelaciona entrada e saída no domínio discreto da mesma maneira que a funçãode transferência no domínio s relaciona sinais contínuos. Observe que g(z) podeser calculada diretamente a partir de g(t), a resposta impulsional do sistema.

• Relação da função de transferência de pulsos com equações de diferença

1 2 m0 1 2 m

1 2 n0 1 2 n

b b z b z b zy(z) u(z)

a a z a z a z

Mostramos anteriormente que uma equação de diferenças na forma:

0 1 2 n

0 1 2 m

a y(k) a y(k 1) a y(k 2) a y(k n)b u(k) b u(k 1) b u(k 2) a u(k m)

Resulta em:

1 2 m0 1 2 m

1 2 n0 1 2 n

b b z b z b zy(z)g(z)u(z) a a z a z a z

De modo que:

Page 32: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Realização física

Já vimos a noção da possibilidade de realização de funções de transferência nodomínio temporal contínuo. Uma condição análoga pode ser estabelecida para afunção de transferência de pulsos, isto é, que um modelo em tempo discreto nãopode ter um sinal de saída que dependa de entradas futuras. Senão o modelonão é fisicamente realizável.

Desse modelo, o modelo é realizável desde que a00. Para comprovar isso,vamos admitir que a0=0, de modo que a equação de diferenças fique:

1 2 n

0 1 2 m

a y(k 1) a y(k 2) a y(k n)b u(k) b u(k 1) b u(k 2) a u(k m)

Nessa situação, um sinal de entrada futuro, u(k), influenciaria a saída atual,y(k-1), o que não é fisicamente possível!

Page 33: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Conversão de transformadas de Laplace e z

Da definição: s t 1 s tz e z e

Usando uma aproximação de Padé:

s t

1

1

1

2 s te2 s t

2 s tz2 s t

2 1 zst 1 z

Outra aproximação:

2

s t

1 s t

1

(s t)e 1 s t2!

z e 1 s tou

1 zst

Page 34: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

• Retentor de Ordem Zero (ZOH)

A função do retentor é converter um sinal digital em um sinal contínuo. O sinal contínuo de um retentor de ordem zero é uma série de pulsos, ou seja, mantém constante o sinal durante cada intervalo de amostragem.

A função de transferência de um ZOH é:

s t1 eH(s)

s

que é a transformada de Laplace da função h(t) = S(t) – S(t – t), sendo S(t) a função degrau unitário.

u*(t)H(s) G(s)

w(t) y(t) y*(t)t

A transformada z do produto H(s)G(s), representada por HG(z), é:

1 G(s)HG(z) (1 z )

sNote que HG(z) H(z)G(z)

Page 35: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Exemplo: Obter a equação de diferenças de um sistema de primeira ordem com ganho unitário e com um retentor ZOH.

s t1 eH(s)s

1G(s)s 1

1 1HG(z) (1 z )s( s 1)

1 t / 1

1 1 1 1 1s( s 1) s (s 1/ ) 1 z 1 e z

1 t /1

1 t / 1 t / 11 1 z (1 e )HG(z) (1 z )

1 z 1 e z 1 e z

Como: , então:y(z) HG(z)u(z)

yn = a1 yn-1 + (1 – a1) un-1

onde: t /1a e

Page 36: Domínio Discreto; Transformada Z. · 2014. 5. 6. · A transformada z inversa produz apenas a sequência discreta f(k), sendo incapaz de fornecer informação sobre a função contínua

Exemplo: resolva a seguinte equações de diferenças de primeira ordem:

11y(z) u(z)

1 3z

yn + 3 yn-1 = un

y-1 = 0

Aplicando a transformada Z:

y(z) + 3 z-1 y(z) = u(z)

Por exemplo, aplicando a inversa da transformada Z para u(z) = 1 (impulso unitário):

nny ( 3)