Articulo de referencia

Invariancia de impulsos

La invariancia de impulso es una técnica para diseñar filtros de respuesta de impulso infinita (IIR) de tiempo discreto a partir de filtros de tiempo continuo, en la que se mues...

La invariancia de impulso es una técnica para diseñar filtros de respuesta de impulso infinita (IIR) de tiempo discreto a partir de filtros de tiempo continuo, en la que se muestrea la respuesta de impulso del sistema de tiempo continuo para producir la respuesta de impulso del sistema de tiempo discreto. La respuesta en frecuencia del sistema de tiempo discreto será una suma de copias desplazadas de la respuesta en frecuencia del sistema de tiempo continuo; si el sistema de tiempo continuo tiene un ancho de banda aproximadamente limitado a una frecuencia menor que la frecuencia de Nyquist del muestreo, entonces la respuesta en frecuencia del sistema de tiempo discreto será aproximadamente igual a ella para frecuencias inferiores a la frecuencia de Nyquist.

Discusión

La respuesta impulsional del sistema de tiempo continuo,hdo(t){\displaystyle h_{c}(t)}, se muestrea con período de muestreoT{\displaystyle T}para producir la respuesta impulsional del sistema de tiempo discreto,h[norte]{\displaystyle h[n]}.

h[norte]=Thdo(norteT){\displaystyle h[n]=Th_{c}(nT)\,}

Por lo tanto, las respuestas en frecuencia de los dos sistemas están relacionadas por

H(mijω)=1Tk=THdo(jωT+j2πTk){\displaystyle H(e^{j\omega })={\frac {1}{T}}\sum _{k=-\infty }^{\infty }{TH_{c}\left(j{\frac {\omega }{T}}+j{\frac {2{\pi }}{T}}k\right)}\,}

Si el filtro de tiempo continuo tiene un ancho de banda aproximadamente limitado (es decir,Hdo(jΩ)<δ{\displaystyle H_{c}(j\Omega )<\delta }cuando|Ω|π/T{\displaystyle |\Omega |\geq \pi /T}), entonces la respuesta en frecuencia del sistema de tiempo discreto será aproximadamente la respuesta en frecuencia del sistema de tiempo continuo para frecuencias inferiores a π radianes por muestra (inferiores a la frecuencia de Nyquist 1/(2 T ) Hz):

H(mijω)=Hdo(jω/T){\displaystyle H(e^{j\omega })=H_{c}(j\omega /T)\,}para|ω|π{\displaystyle |\omega |\leq \pi \,}

Comparación con la transformada bilineal

Cabe destacar que se producirá aliasing, incluso por debajo de la frecuencia de Nyquist, siempre que la respuesta del filtro de tiempo continuo sea distinta de cero por encima de dicha frecuencia. La transformación bilineal es una alternativa a la invariancia de impulso que utiliza una función de mapeo diferente: la respuesta en frecuencia del sistema de tiempo continuo, hasta frecuencia infinita, se proyecta en el rango de frecuencias hasta la frecuencia de Nyquist en el caso de tiempo discreto, a diferencia de la invariancia de impulso, que mapea las frecuencias linealmente con solapamiento circular.

Efecto sobre los polos en la función del sistema

Si los polos continuos ens=sk{\displaystyle s=s_{k}}, la función del sistema se puede escribir en una expansión en fracciones parciales como

Hdo(s)=k=1norteAkssk{\displaystyle H_{c}(s)=\sum _{k=1}^{N}{\frac {A_{k}}{s-s_{k}}}\,}

Así, utilizando la transformada inversa de Laplace, la respuesta al impulso es

hdo(t)={k=1norteAkmiskt,t00,de lo contrario{\displaystyle h_{c}(t)={\begin{cases}\sum _{k=1}^{N}{A_{k}e^{s_{k}t}},&t\geq 0\\0,&{\mbox{en otro caso}}\end{cases}}}

La respuesta impulsional del sistema de tiempo discreto correspondiente se define entonces de la siguiente manera:

h[norte]=Thdo(norteT){\displaystyle h[n]=Th_{c}(nT)\,}
h[norte]=Tk=1norteAkmisknorteT[norte]{\displaystyle h[n]=T\sum _{k=1}^{N}{A_{k}e^{s_{k}nT}u[n]}\,}

Al aplicar una transformada Z a la respuesta impulsional en tiempo discreto, se obtiene la siguiente función del sistema en tiempo discreto.

H(z)=Tk=1norteAk1miskTz1{\displaystyle H(z)=T\sum _{k=1}^{N}{\frac {A_{k}}{1-e^{s_{k}T}z^{-1}}}\,}

De este modo, los polos de la función del sistema de tiempo continuo se traducen en polos en z = e s k T . Los ceros, si los hay, no se transforman de forma tan sencilla.

Polos y ceros

Si la función del sistema tiene tanto polos como ceros, se pueden mapear de la misma manera, pero el resultado ya no es una invariancia de impulso: la respuesta al impulso en tiempo discreto no es simplemente igual a muestras de la respuesta al impulso en tiempo continuo. Este método se conoce como el método de la transformada Z adaptada o mapeo de polos y ceros.

Estabilidad y causalidad

Dado que los polos en el sistema de tiempo continuo en s = s k se transforman en polos en el sistema de tiempo discreto en z = exp( s k T ), los polos en la mitad izquierda del plano s se mapean al interior del círculo unitario en el plano z ; por lo tanto, si el filtro de tiempo continuo es causal y estable, entonces el filtro de tiempo discreto también será causal y estable.

Fórmula corregida

Cuando una respuesta de impulso causal en tiempo continuo tiene una discontinuidad ent=0{\displaystyle t=0}, las expresiones anteriores no son consistentes. [ 1 ] Esto se debe a quehdo(0){\displaystyle h_{c}(0)}tiene límites derechos e izquierdos diferentes, y en realidad solo debería aportar su promedio, la mitad de su valor derecho.hdo(0+){\displaystyle h_{c}(0_{+})}, ah[0]{\displaystyle h[0]}.

Hacer esta corrección da

h[norte]=T(hdo(norteT)12hdo(0+)δ[norte]){\displaystyle h[n]=T\left(h_{c}(nT)-{\frac {1}{2}}h_{c}(0_{+})\delta [n]\right)\,}
h[norte]=Tk=1norteAkmisknorteT([norte]12δ[norte]){\displaystyle h[n]=T\sum _{k=1}^{N}{A_{k}e^{s_{k}nT}}\left(u[n]-{\frac {1}{2}}\delta [n]\right)\,}

Al aplicar una transformada Z a la respuesta impulsional en tiempo discreto, se obtiene la siguiente función del sistema en tiempo discreto.

H(z)=Tk=1norteAk1miskTz1T2k=1norteAk.{\displaystyle H(z)=T\sum _{k=1}^{N}{{\frac {A_{k}}{1-e^{s_{k}T}z^{-1}}}-{\frac {T}{2}}\sum _{k=1}^{N}A_{k}}.}

La segunda suma es cero para filtros sin discontinuidad, por lo que ignorarla suele ser seguro.

Véase también

Referencias

  1. Jackson, LB (1 de octubre de 2000). "Una corrección a la invariancia de impulso". IEEE Signal Processing Letters . 7 (10): 273– 275. Bibcode : 2000ISPL....7..273J . doi : 10.1109/97.870677 . ISSN 1070-9908 . 

Otras fuentes

  • Oppenheim, Alan V. y Schafer, Ronald W. con Buck, John R. Procesamiento de señales en tiempo discreto. Segunda edición. Upper Saddle River, Nueva Jersey: Prentice-Hall, 1999.
  • Sahai, Anant. Clase magistral. Ingeniería Eléctrica 123: Procesamiento Digital de Señales. Universidad de California, Berkeley. 5 de abril de 2007.
  • Eitelberg, Ed. «Invariancia de convolución e invariancia de impulso corregida». Procesamiento de señales, vol. 86, número 5, págs.  1116-1120. 2006
  • Transformación invariante de impulso en CircuitDesign.info. Breve explicación, un ejemplo y aplicación a convertidores analógico-digitales sigma-delta de tiempo continuo .