TE-072 Processamento Digital de Sinais I - UFPR 3. Transformada Z 3.1. Definição Seja um sistema discreto LTI: x[n] Com: h[n] y[n] x[n] z n , z Complexos A saída y[n] pode ser calculada como: 1 TE-072 Processamento Digital de Sinais I - UFPR y[n] x[n] * h[n] y[n] h[k ] x[n k ] k y[n] h[k ] z nk k y[n] z n k h [ k ] z k 2 y[n] z n h[k ] z TE-072 Processamento Digital de Sinais I - UFPR k k Definindo: H ( z) k h [ k ]. z k Temos que: y[n] H ( z) z n Cte complexa Logo: zn é autofunção do sistema Discreto LTI e H(z) seu autovalor correspondente. 3 TE-072 Processamento Digital de Sinais I - UFPR Logo, definimos Transformada Z do sinal discreto x[n] como: X ( z ) Z{x[n]} n x [ n ]. z n x[n] X ( z) Z Transformada Z Unilateral: X ( z ) x[n].z n n 0 Equivalente à TZ bilateral quando x[n]=0 n<0 ; Usada p/ analisar EDCC com condições iniciais não nulas. 4 TE-072 Processamento Digital de Sinais I - UFPR Escrevendo o número complexo z na sua forma polar: z r.e Temos: j x[n].r X r.e j n .e jn n Logo: Se r=1: X r.e j F x[n].r n X e j Fx[n] Transformada de Fourier 5 TE-072 Processamento Digital de Sinais I - UFPR Logo: Transformada Z pode ser obtida a partir da Transformada de Fourier fazendo: z e j X ( z) Fx[n]e j z O inverso nem sempre é verdade!!! Pois: Z x[n] F x[n].r n Pode fazer com que alguns sinais se tornem convergentes 6 TE-072 Processamento Digital de Sinais I - UFPR Analogia Contínuo Discreto Contínuo: Transformada de Laplace: st x ( t ). e dt X ( s) s j Fazendo: s j Obtemos a Transformada de Fourier: X ( ) x(t ).e jt j dt s Se eixo j ROC X ( s) x(t ).e .e t 0 jt dt Eixo j 7 TE-072 Processamento Digital de Sinais I - UFPR Discreto: Transformada Z: X ( z) x[n].z n z r.e j n ze Fazendo: j Obtemos a Transformada de Fourier: X () x[n].e jn Im{z} z n -1 Se circulounitário ROC X ( z ) x[n]r n .e jn n r0 1 Re{z} Circulo Unitário 8 TE-072 Processamento Digital de Sinais I Ex.: - UFPR x[n] a u[n] n X ( z) n n a u [ n ]. z n X ( z ) a .z n n n 0 PG: =a.z-1 a0=1 n= a.z 1 n n 0 1 n S n a0 1 1 (a.z 1 ) X ( z ) 1 1 a.z 1 Converge se: a. z 1 1 za 9 TE-072 Processamento Digital de Sinais I Neste caso: - UFPR 1 0 X ( z) 1 a.z 1 z X ( z) , z a za Im{z} z -1 1 a Re{z} Se |a|<1 T.Fourier |a|>1 T.Fourier 10 Ex.2: x[n] a u[n 1] TE-072 Processamento Digital de Sinais I n X ( z) - UFPR n n a u [ n 1 ]. z n 1 X ( z ) a .z n n n PG: n 1 =a-1.z a.z n X ( z ) a .z 1 1 1 n a.z 1 n n 1 n a0= a-1.z n= 1 n S n a0 1 1 1 ( a . z ) 1 X ( z ) a z 1 a 1.z Converge se: 1 a .z 1 za 11 TE-072 Processamento Digital de Sinais I Neste caso: - UFPR 1 0 X ( z ) a z 1 a 1.z z X ( z) , za za 1 Im{z} z -1 1 a Re{z} Se |a|>1 T.Fourier |a|<1 T.Fourier 12 TE-072 Processamento Digital de Sinais I - UFPR Conclusão: Sinais diferentes podem ter a mesma expressão algébrica de X(z). Logo uma Transformada Z só é completamente definida se especificarmos: - Expressão algébrica de X(z) - Região de Convergência(ROC) 13 TE-072 Processamento Digital de Sinais I - UFPR Qualquer sinal que pode ser representado como um somatório de exponenciais complexas, poderá ser representado por uma Transformada Z composta de uma razão de dois polinômios. N ( z ) aN z N aN 1 z N 1 aN 2 z N 2 ... a1 z1 a0 z 0 X ( z) D( z ) bM z M bM 1 z M 1 bM 2 z M 2 ... b1 z1 b0 z 0 Raízes de N(z) Zeros da X(z) Raízes de D(z) Pólos da X(z) Nos exemplos: Valores que fazem X(z) igual a ZERO Valores que fazem X(z) igual a INFINITO z X ( z) za Zero: z=0 Pólo: z=a 14 TE-072 Processamento Digital de Sinais I - UFPR Diagrama de pólos e zeros: Representação gráfica no plano z dos pólos e zeros. Im{z} z -1 1 a Re{z} z X ( z) , z a za 15 TE-072 Processamento Digital de Sinais I n Ex.3: - UFPR n 1 1 x[n] u[n] u[n] 2 3 n 1 n n 1 X ( z ) u[n] u[n] z 3 n 2 n n 1 1 n X ( z ) z n 0 2 3 n n 1 1 X ( z ) z n z n 3 n 0 2 n 0 n 1 1 X ( z ) z 1 z 1 n 0 3 n 0 2 n 16 TE-072 Processamento Digital de Sinais I X ( z) 1 1 12 z 1 1 z 2 1 1 13 z 1 - UFPR z z z 12 z 13 1 z 3 2 z 2 16 z 2 zz 121 X ( z) 2 1 1 1 z 2 z 3 z 6 z 16 17 TE-072 Processamento Digital de Sinais I z z 12 z z 13 1 z 2 Im{z} -1 1/2 1 -1 Re{z} z 1 -1/3 Im{z} 2 z z 121 X ( z) z 12 z 13 -1/3 Re{z} z 1/12 -1 UFPR 1 z 3 Im{z} z - 1/2 1 Re{z} 1 z 2 18 TE-072 Processamento Digital de Sinais I n - UFPR n Ex.4: x[n] 1 u[n] 1 u[n 1] 3 2 Usando os resultados das análises anteriores e a Propriedade de Linearidade da Transformada Z. n z 1 Z u[n] z 13 3 n z 1 Z u[n 1] z 12 2 Logo: X ( z) z z z 13 z 12 1 z 3 1 z 2 ROC : 1 1 z 3 2 19 TE-072 Processamento Digital de Sinais I 2 z z 121 z z X ( z) 1 1 z 3 z 2 z 13 z 12 Im{z} 1 1 ROC : z 3 2 -1 UFPR z 1 1/12 -1/3 - 1/2 Re{z} 20 TE-072 Processamento Digital de Sinais I Ex.5: - UFPR n a , 0 n N 1 n x[n] a u[n] u[n N ] 0 , outros X ( z) N 1 n n 0 n n n x [ n ]. z a .z N 1 X ( z ) a.z n 0 PG: a0=1 =a.z-1 n=N 1 n 1 n S n a0 1 1 N 1 a.z X ( z ) 1 1 a.z 1 X(z) converge se Isto é: a.z 1 N a e z 0 21 TE-072 Processamento Digital de Sinais I - UFPR 1 N 1 a.z 1 a N .z N X ( z ) 1 1 1 a.z 1 a.z 1 zN aN N N N 1 z a X ( z) z N 1 za z za z Pólos da X(z): Zeros da X(z): z N 1 0 N 1 pólos em z 0 z a 0 pólo em z a z N a N 0 Zeros raízes do polinôm io zeros em z a.e j 2k / N k 0,1, ..., N 1 Quando k=0: zero em z=a logo cancela com o pólo em z=a 22 TE-072 Processamento Digital de Sinais I - UFPR Logo tem-se: N-1 pólos em z=0 N-1 zeros distribuídos uniformemente sobre um círculo de raio a p/ N=8 Im{z} -1 (7) z 2k/8 1 a Re{z} ROC: Todo plano z com exceção de z=0 23 TE-072 Processamento Digital de Sinais I - UFPR 3.2. Propriedades da ROC Considerando X(z) uma função racional em z e x[n] finito p/n finito 1) A ROC de X(z) é um anel ou disco centrado na origem (z=0) 2) A Transformada de Fourier de x[n] converge absolutamente se e somente se a ROC de X(z) inclui a circunferência unitária. 3) A ROC não contém pólos de X(z) 4) Se x[n] tem duração finita, x[n]0 p/ -<N1nN2<, a ROC é todo plano z com possíveis exceções em z=0 e z= 24 TE-072 Processamento Digital de Sinais I - UFPR 5) Se x[n] é definida à direita, x[n]=0 p/ n<N1<, a ROC estende-se de |z|=r0 (maior pólo) até , incluindo ou não z= 6) Se x[n] é definida à esquerda, x[n]=0 p/ n>N2>, a ROC será |z|<r0 (menor pólo), incluindo ou não z=0 7) Se x[n] é definida à esquerda e à direita, a ROC será um anel compreendido entre 2 pólos. 8) A ROC deve ser uma região conexa. 25 TE-072 Processamento Digital de Sinais I - UFPR 3.3. Transformada Z Inversa Demonstração da fórmula de inversão. X ( z ) F x[n].r n F x[n].r x[n].r F X r.e 1 x[n] r X r.e .e 2 X r.e j n n 1 j j n jn .d 2 1 x[n] 2 X r.e . re .d j j n 2 26 TE-072 Processamento Digital de Sinais I 1 x[n] 2 - UFPR X r.e . re .d j j n 2 z r.e j Mudança de variáveis: dz jr.e j .d 1 1 d j dz j r.e Variando de 0 a 2 z varia sobre uma circunferência de raio r. |z|=r ROC de X(z) x[n] x[n] 1 2 j r zn X ( z ) dz z 1 2 j X ( z ) z n1 dz r Resolve-se utilizando o Teorema dos Resíduos 27 TE-072 Processamento Digital de Sinais I Pares de Transformadas Z x[n] X (z ) [n] 1 Todo Plano Z z z 1 z z 1 z 1 u[n] u[n 1] [n n0 ] n a u[n] a u[n 1] n z n0 z za z za - UFPR ROC z 1 Todo Plano Z exceto: z 0 p / n0 0 z p / n0 0 za za 28 TE-072 Processamento Digital de Sinais I X (z ) x[n] n n.a u[n] 2 n n .a u[n] n 1.a u[n] n a.z z a 2 a.z.z a 3 z a z2 z a 2 - UFPR ROC za za za cos0 n.u[n] zz cos0 z 2 2 cos 0 .z 1 z 1 sin0n.u[n] z sin 0 z 2 2 cos 0 .z 1 z 1 29 TE-072 Processamento Digital de Sinais I x[n] X (z ) r cos0n.u[n] z r r sin0n.u[n] z.r. sin 0 z 2 2r cos0 .z r 2 z r n a , 0 n N 1 0 , outros n zN aN N 1 z z a UFPR ROC zz r. cos0 z 2 2r cos0 .z r 2 n - z 0 30 TE-072 Processamento Digital de Sinais I - UFPR 3.3.1. Inversão por Inspeção Consiste no uso eficiente das tabelas e propriedades da Transformada Z 31 TE-072 Processamento Digital de Sinais I - UFPR 3.3.2. Expansão em Frações Parciais Revisão: Dado G(v) função racional em v com grau N(v) < grau D(v) N (v ) G (v ) D (v ) r Pode ser escrita na forma i Aik G (v) k i 1 k 1 v pi Onde: r = número de pólos i = multiplicidade do pólo i Aik = coeficiente relativo a k-ésima parcela do pólo i 32 TE-072 Processamento Digital de Sinais I - UFPR Onde: i k 1 d i Aik i k v pi .G(v) i k ! dv v pi Ex.: s2 4 H ( s) ( s 2) 2 ( s 2 2s 2) Pólos: duplos em s=-2 e complexo s=-1j 1=2 2=1 3=1 r=4 33 TE-072 Processamento Digital de Sinais I s2 4 A B Cs D H ( s) 2 2 2 2 (s 2) (s 2s 2) (s 2) s 2 s 2s 1 - UFPR 1 d 2 2 2 A 22 s 2 H ( s) 2 2! ds s 2 s2 4 8 A 2 4 s 2s 2 s 2 2 1 d 21 2 B 21 s 2 H ( s ) 2 1! ds s 2 d s2 4 2s 2 4s 8 B 2 4 2 3 2 ds s 2 s 2 s 2 s 4 s 8s 8s 4 s 2 34 TE-072 Processamento Digital de Sinais I Pólo complexo: No caso: - UFPR Cs D A' B' s 2 2s 2 s R jI s R jI Cs D ( A' B' )s A' ( R jI ) B' ( R jI ) 1 d 11 A' 11 ( s 1 j ) H ( s ) (1 1)! ds s 1 j s2 4 42j A' 1 0.5 j 2 ( s 2) ( s 1 j ) s 1 j 4 1 d 11 B' 11 ( s 1 j ) H ( s ) (1 1)! ds s 1 j s2 4 42j B' 1 0 .5 j 2 ( s 2) ( s 1 j ) s 1 j 4 35 TE-072 Processamento Digital de Sinais I - UFPR P/ funções com coeficientes reais: sempre B’=A’* Logo: Assim: C 2. Re{A'} D 2Re{ p}. Re{A'} Im{p}. Im{A'} C 2.(1) 2 D 2(1)(1) (0.5)(1) 3 s2 4 4 2 2s 3 H ( s) 2 2 2 2 (s 2) (s 2s 2) (s 2) s 2 s 2s 1 Matlab: função residue [r,p,k]=residue(n,d) 36 TE-072 Processamento Digital de Sinais I - UFPR No caso específico da Transformada Z Como as funções básicas são na forma: z za A expansão em frações parciais não pode ser aplicada diretamente na X(z). Soluções: 1) Aplicar o método na função: 2) Aplicar o método na função: X ( z) z X z 1 Matlab: função residuez [r,p,k]=residuez(n,d) 37 3z z Ex.: X ( z ) z 14 z 13 2 Método 1: TE-072 Processamento Digital de Sinais I 5 6 - UFPR 3z 56 X ( z) z 14 z 13 z X ( z) A B 1 z z 4 z 13 3z 56 1 A z 1 4 z 14 z 13 1 z 4 3z 56 1 2 B z 3 1 1 z 4 z 3 1 z 3 38 TE-072 Processamento Digital de Sinais I Logo: - UFPR X ( z) 1 2 1 z z 4 z 13 z 2z X ( z) 1 z 4 z 13 Por tabela temos: x[n] u[n] 2 u[n] 1 n 4 1 n 3 39 TE-072 Processamento Digital de Sinais I - UFPR 3z 2 56 z X ( z) z 14 z 13 Método 2: 5 1 3 A B 1 6 z X (z ) 1 1 1 1 1 1 1 1 1 4 z 1 3 z 1 4 z 1 3 z 3 56 z 1 1 1 A 1 4 z 1 1 1 1 1 4 z 1 3 z 3 56 z 1 1 1 B 1 3 z 1 1 1 1 1 4 z 1 3 z 1 z 14 2 z 13 40 TE-072 Processamento Digital de Sinais I 1 1 X (z ) 1 z 1 4 1 - UFPR 2 1 z 1 3 1 Por tabela temos: x[n] u[n] 2 u[n] 1 n 4 1 n 3 41 TE-072 Processamento Digital de Sinais I - UFPR 3.3.2. Expansão em Série de Potência Definição da Transformada Z Série de Laurent X ( z) n x [ n ]. z n X ( z) ... x[2].z 2 x[1].z x[0] x[1].z 1 x[2].z 2 ... 42 TE-072 Processamento Digital de Sinais I Ex.: - UFPR z X ( z) z 0.5 x[n] (0.5)n u[n] Sabemos por tabela: n 0 1 2 3 4 ... Isto é: x[n] 1 0.5 0.25 0.125 0.0625.... Podemos calcular a série de potência de uma razão de polinômios por divisões sucessivas: 43 TE-072 Processamento Digital de Sinais I - UFPR z 0.5 z 1 0.5 z 1 0.25z 2 0.125z 3 .... z 0.5 0.5 0.5 0.25z 1 0.25z 1 0.25z 1 0.125z 2 0.125z 2 0.125z 2 0.0625z 3 .... 44 Ex.2: X ( z) z 2 1 12 z 1 1 z 1 1 z 1 TE-072 Processamento Digital de Sinais I - UFPR Pólos somente em z=0, frações parciais não é apropriado. Multiplicando todos os termos: X ( z) z z 1 z 2 1 2 1 2 1 De tabela temos: x[n] [n 2] 12 [n 1] [n] 12 [n 1] Ex.3: X ( z) log 1 az1 , z a 45 TE-072 Processamento Digital de Sinais I 3.4.Propriedades da Transformada Z - UFPR 46 TE-072 Processamento Digital de Sinais I - UFPR Exercícios: 1) z2 X ( z) , z 14 z 12 2) z 2z 1 X ( z) 2 3 z 2 z 12 z 12 2 3) z3 2z 2 z X ( z) 2 3 z 2 z 12 Verificar as respostas usando Divisões Sucessivas 47