image/svg+xml

Discretizaci贸n de sistemas discretos

Para utilizar t茅cnicas de procesamiento digital de se帽ales primero es necesario discretizar los sistemas, pues los modelos que obtenemos son continuos. Para este fin existen dos acercamientos.

El primero ser铆a discretizar el modelo en tiempo continuo y luego por medio de la transformada $\mathcal{Z}$, pasar al dominio de la frecuencia y obtener la funci贸n de transferencia.

\[\large\boxed{h(t)}\overset{\text{discretizar}}{\longrightarrow}\boxed{h(n)}\overset{\normalsize\mathcal{Z}}{\longrightarrow}\boxed{H(z)}\]

La otra alternativa es obtener la funci贸n de transferencia del modelo en tiempo continuo con la transformada de Laplace y luego por medio de una aproximaci贸n, obtener su sistema digital (discreto) equivalente aproximado.

\[\large\boxed{h(t)}\overset{\normalsize \mathcal{L}}{\longrightarrow}\boxed{H(s)}\overset{\approx}{\longrightarrow}\boxed{H(z)}\]

La equivalencia exacta entre el plano $s$ y el plano $z$ es:

\[s = \frac{\ln(z)}{T}\]

Donde $T$ es el tiempo de muestreo.

Implementar la equivalencia exacta para pasar de la funci贸n de transferencia continua a la funci贸n de transferencia discreta, ser铆a demasiado tortuoso. Por eso se usan las aproximaciones siguientes.

Sin embargo, la transformaci贸n que se usar谩 en el curso es la bilineal, pues la aproximaci贸n hacia adelante de Euler puede hacer que sistemas estables de tiempo continuo, pasen a ser inestables en tiempo discreto. Para la aproximaci贸n hacia atr谩s, esto no ocurre as铆, sin embargo puede llegar a deformar significativamente la funci贸n de transferencia. La de Tustin, hace el trabajo muy bien.

Procedimiento

  1. Obtener la funci贸n de transferencia del sistema.
  2. Calcular $T$ a partir de $f_s > 2B$
  3. Sustituir:

    \[s = \dfrac{2(z-1)}{T(z-1)}\]
  4. Desarrollar

En matlab se puede utilizar el comando c2d(H,T,'Tustin') para obtener $H(z)$, donde H, es la funci贸n de transferencia en tiempo continuo $H(s)$ y $T$ el tiempo de muestreo.