Transformação Invariante ao Impulso: mudanças entre as edições

De MediaWiki do Campus São José
Ir para navegação Ir para pesquisar
Criou página com '*Filtros Analógicos: A função de transferência do filtro analógico :<math> H_a(s) = \frac {c_0 + c_1 s + c_2 s^2 + ... + c_m s^m} {d_0 + d_1 s + d_2 s^2 + ... + d_n s^n}</m...'
 
Moecke (discussão | contribs)
 
(5 revisões intermediárias por 2 usuários não estão sendo mostradas)
Linha 1: Linha 1:
*Filtros Analógicos:
=Base teórica=
A função de transferência do filtro analógico  
A função de transferência do filtro analógico  
:<math> H_a(s) = \frac {c_0 + c_1 s + c_2 s^2 + ... + c_m s^m} {d_0 + d_1 s + d_2 s^2 + ... + d_n s^n}</math> onde <math> m \le n </math>
:<math> H_a(s) = \frac {c_0 + c_1 s + c_2 s^2 + ... + c_m s^m} {d_0 + d_1 s + d_2 s^2 + ... + d_n s^n}</math> onde <math> m \le n </math>
Linha 9: Linha 9:
dado a [https://en.wikipedia.org/wiki/List_of_Laplace_transforms transformada de Laplace] da exponencial decrescente
dado a [https://en.wikipedia.org/wiki/List_of_Laplace_transforms transformada de Laplace] da exponencial decrescente


: <math> e^{-\alpha t} \cdot u(t) \longleftrightarrow { 1 \over s+\alpha } </math>
: <math> { 1 \over s - \alpha }  \longleftrightarrow  e^{\alpha t} \cdot u(t) </math>
 
a resposta de cada termo da fração parcial é obtida como:
 
<math> \frac{r_k}{s - p_k}  \longleftrightarrow r_k  e^{p_k t} \cdot u(t) </math>
 
e portanto a resposta ao impulso do filtro analógico é dada por
 
: <math>h_a(t) = \sum_{k=1}^{K} \left( r_k  e^{p_k t} \cdot \right) u(t)</math>
 
Amostrando <math>h_a(t) \ </math> periodicamente em um período <math>T \ </math>, é possível obter a resposta ao impulso do filtro analógico digitalizado.
 
: <math>h_d(n) = h_a(n T)= \sum_{k=1}^{K} \left( r_k  e^{p_k n T} \cdot  \right) u(n T)</math>
 
Considerando o par de [https://en.wikipedia.org/wiki/Z-transform#Table_of_common_Z-transform_pairs transformada Z]
 
:<math>a^n \cdot u[n]  \longleftrightarrow \frac{1}{1-a z^{-1}} </math>
ou
:<math>a^n \cdot u[n]  \longleftrightarrow  \frac{z}{z-a} </math>
 
considerando
:<math>a = e^{p_k T}</math>
 
obtém-se a função de transferência do filtro digital para a substituição exata de :<math>z = e^{sT} \ </math>
:<math>H_d(z) = \sum_{k=1}^{K}  r_k  \frac{1}{1-e^{p_k T}z^{-1}} </math>
ou
:<math>H_d(z) = \sum_{k=1}^{K}  r_k  \frac{z}{z-e^{p_k T}} </math>
 
Portanto cada polo <math>p_k \ </math> do filtro analógico é transformado em um polo <math>a = e^{p_k T} \ </math> no filtro digital
 
=Método de transformação digital=
O método de transformação do filtro analógico em digital consiste básica de:
:1) fazer a expansão em frações parciais da função de transferência <math> H_a(s) \ </math>
:2) calcular os polos <math> e^{p_k T} \ </math> da função de transferência <math> H_d(z) \ </math>
:3) obter a divisão de polinômios para <math> H_d(z) = \frac{num(z)}{den(z)}\ </math>, para obter os coeficientes do filtro digital.
 
 
;Importante:
*Como esse tipo de transformação digital resulta em um enrolamento das frequências do filtro analógico (eixo imaginário no plano ''s'') no circulo unitário do plano ''z'', ocorre a repetição periódica das respostas de frequência e ''aliasing'' de frequência.  Por isso é necessário que o filtro analógico seja limitado em frequência, e que a atenuação na banda de rejeição seja monotonicamente decrescente.  Assim apenas filtros do tipo '''passa-baixas''' e '''passa-faixas''' com aproximação do tipo '''Butterworth''' e '''Chebyshev tipo 1''' podem ser utilizados com esse método.
*A dedução acima considerou apenas polos não múltiplos, e necessita de pequenos ajustes no caso da existência de polos múltiplos (situação raramente encontrada em filtros analógicos).

Edição atual tal como às 01h12min de 23 de abril de 2020

1 Base teórica

A função de transferência do filtro analógico

Ha(s)=c0+c1s+c2s2+...+cmsmd0+d1s+d2s2+...+dnsn onde mn

pode ser expandida em frações parciais, considerando pk os polos não múltiplos de Ha(s):

Ha(s)=k=1K(rkspk)

dado a transformada de Laplace da exponencial decrescente

1sαeαtu(t)

a resposta de cada termo da fração parcial é obtida como:

rkspkrkepktu(t)

e portanto a resposta ao impulso do filtro analógico é dada por

ha(t)=k=1K(rkepkt)u(t)

Amostrando ha(t)  periodicamente em um período T , é possível obter a resposta ao impulso do filtro analógico digitalizado.

hd(n)=ha(nT)=k=1K(rkepknT)u(nT)

Considerando o par de transformada Z

anu[n]11az1

ou

anu[n]zza

considerando

a=epkT

obtém-se a função de transferência do filtro digital para a substituição exata de :z=esT 

Hd(z)=k=1Krk11epkTz1

ou

Hd(z)=k=1KrkzzepkT

Portanto cada polo pk  do filtro analógico é transformado em um polo a=epkT  no filtro digital

2 Método de transformação digital

O método de transformação do filtro analógico em digital consiste básica de:

1) fazer a expansão em frações parciais da função de transferência Ha(s) 
2) calcular os polos epkT  da função de transferência Hd(z) 
3) obter a divisão de polinômios para Hd(z)=num(z)den(z) , para obter os coeficientes do filtro digital.


Importante
  • Como esse tipo de transformação digital resulta em um enrolamento das frequências do filtro analógico (eixo imaginário no plano s) no circulo unitário do plano z, ocorre a repetição periódica das respostas de frequência e aliasing de frequência. Por isso é necessário que o filtro analógico seja limitado em frequência, e que a atenuação na banda de rejeição seja monotonicamente decrescente. Assim apenas filtros do tipo passa-baixas e passa-faixas com aproximação do tipo Butterworth e Chebyshev tipo 1 podem ser utilizados com esse método.
  • A dedução acima considerou apenas polos não múltiplos, e necessita de pequenos ajustes no caso da existência de polos múltiplos (situação raramente encontrada em filtros analógicos).