Main Content

La traducción de esta página aún no se ha actualizado a la versión más reciente. Haga clic aquí para ver la última versión en inglés.

ellip

Diseño de filtros elípticos

Descripción

ejemplo

[b,a] = ellip(n,Rp,Rs,Wp) devuelve los coeficientes de la función de transferencia de un filtro digital elíptico paso bajo de n-ésimo orden con una frecuencia de borde de la banda de paso normalizada Wp. El filtro resultante tiene Rp decibelios de ondulación pico a pico de la banda de paso y Rs decibelios de atenuación de la banda de parada desde el valor pico de la banda de paso.

ejemplo

[b,a] = ellip(n,Rp,Rs,Wp,ftype) diseña un filtro elíptico paso bajo, paso alto, paso banda o eliminador de banda, dependiendo del valor de ftype y del número de elementos de Wp. Los diseños de paso banda y de eliminador de banda resultantes son de orden 2n.

Nota: Consulte Limitaciones para obtener información sobre las cuestiones numéricas que afectan a la formación de la función de transferencia.

ejemplo

[z,p,k] = ellip(___) diseña un filtro digital elíptico paso bajo, paso alto, paso banda o eliminador de banda y devuelve sus ceros, polos y ganancia. Esta sintaxis puede incluir cualquiera de los argumentos de entrada de las sintaxis anteriores.

ejemplo

[A,B,C,D] = ellip(___) diseña un filtro digital elíptico paso bajo, paso alto, paso banda o eliminador de banda y devuelve las matrices que especifican su representación en el espacio de estados.

ejemplo

[___] = ellip(___,'s') diseña un filtro analógico elíptico paso bajo, paso alto, paso banda o eliminador de banda con frecuencia angular de borde de la banda de paso Wp, Rp decibelios de ondulación de la banda de paso y Rs decibelios de atenuación de la banda de parada.

Ejemplos

contraer todo

Diseñe un filtro elíptico paso bajo de 6.º orden con 10 dB de ondulación de la banda de paso, 50 dB de atenuación de la banda de parada y una frecuencia de borde de la banda de paso de 300 Hz, que, para datos muestreados a 1000 Hz, corresponde a 0.6π rad/muestra. Represente sus respuestas de magnitud y fase. Utilícelo para filtrar una señal aleatoria de 1000 muestras.

fc = 300;
fs = 1000;

[b,a] = ellip(6,10,50,fc/(fs/2));

freqz(b,a,[],fs)

subplot(2,1,1)
ylim([-100 20])

Figure contains 2 axes objects. Axes object 1 with title Phase, xlabel Frequency (Hz), ylabel Phase (degrees) contains an object of type line. Axes object 2 with title Magnitude, xlabel Frequency (Hz), ylabel Magnitude (dB) contains an object of type line.

dataIn = randn(1000,1);
dataOut = filter(b,a,dataIn);

Diseñe un filtro elíptico eliminador de banda de 6.º orden con frecuencias de borde normalizadas de 0.2π y 0.6π rad/muestra, 5 dB de ondulación de la banda de paso y 50 dB de atenuación de la banda de parada. Represente sus respuestas de magnitud y fase. Utilícelo para filtrar datos aleatorios.

[b,a] = ellip(3,5,50,[0.2 0.6],'stop');
freqz(b,a)

Figure contains 2 axes objects. Axes object 1 with title Phase, xlabel Normalized Frequency (\times\pi rad/sample), ylabel Phase (degrees) contains an object of type line. Axes object 2 with title Magnitude, xlabel Normalized Frequency (\times\pi rad/sample), ylabel Magnitude (dB) contains an object of type line.

dataIn = randn(1000,1);
dataOut = filter(b,a,dataIn);

Diseñe un filtro elíptico paso alto de 6.º orden con una frecuencia de borde de la banda de paso de 300 Hz, que, para datos muestreados a 1000 Hz, corresponde a 0.6π rad/muestra. Especifique 3 dB de ondulación de la banda de paso y 50 dB de atenuación de la banda de parada. Represente las respuestas de magnitud y fase. Convierta los ceros, los polos y la ganancia en secciones de segundo orden para utilizarlos en fvtool.

[z,p,k] = ellip(6,3,50,300/500,'high');
sos = zp2sos(z,p,k);
fvtool(sos,'Analysis','freq')

Figure Figure 1: Magnitude Response (dB) and Phase Response contains an axes object. The axes object with title Magnitude Response (dB) and Phase Response, xlabel Normalized Frequency ( times pi blank rad/sample), ylabel Magnitude (dB) contains an object of type line.

Diseñe un filtro elíptico paso banda de 20.º orden con una frecuencia de la banda de paso inferior de 500 Hz y una frecuencia de la banda de paso superior de 560 Hz. Especifique una ondulación de la banda de paso de 3 dB, una atenuación de la banda de parada de 40 dB y una tasa de muestreo de 1500 Hz. Utilice la representación del espacio de estados. Diseñe un filtro idéntico utilizando designfilt.

[A,B,C,D] = ellip(10,3,40,[500 560]/750);
d = designfilt('bandpassiir','FilterOrder',20, ...
    'PassbandFrequency1',500,'PassbandFrequency2',560, ...
    'PassbandRipple',3, ...
    'StopbandAttenuation1',40,'StopbandAttenuation2',40, ...
    'SampleRate',1500);

Convierta la representación del espacio de estados en secciones de segundo orden. Visualice las respuestas en frecuencia utilizando fvtool.

sos = ss2sos(A,B,C,D);
fvt = fvtool(sos,d,'Fs',1500);
legend(fvt,'ellip','designfilt')

Figure Figure 1: Magnitude Response (dB) contains an axes object. The axes object with title Magnitude Response (dB), xlabel Frequency (Hz), ylabel Magnitude (dB) contains 2 objects of type line. These objects represent ellip, designfilt.

Diseñe un filtro analógico Butterworth paso bajo de 5.º orden con una frecuencia de corte de 2 GHz. Multiplique por 2π para convertir la frecuencia en radianes por segundo. Calcule la respuesta en frecuencia del filtro en 4096 puntos.

n = 5;
fc = 2e9;

[zb,pb,kb] = butter(n,2*pi*fc,"s");
[bb,ab] = zp2tf(zb,pb,kb);
[hb,wb] = freqs(bb,ab,4096);

Diseñe un filtro Chebyshev Tipo I de 5.º orden con la misma frecuencia de borde y 3 dB de ondulación de banda de paso. Calcule su respuesta en frecuencia.

[z1,p1,k1] = cheby1(n,3,2*pi*fc,"s");
[b1,a1] = zp2tf(z1,p1,k1);
[h1,w1] = freqs(b1,a1,4096);

Diseñe un filtro Chebyshev Tipo II de 5.º orden con la misma frecuencia de borde y 30 dB de atenuación de la banda de parada. Calcule su respuesta en frecuencia.

[z2,p2,k2] = cheby2(n,30,2*pi*fc,"s");
[b2,a2] = zp2tf(z2,p2,k2);
[h2,w2] = freqs(b2,a2,4096);

Diseñe un filtro elíptico de 5.º orden con la misma frecuencia de borde, 3 dB de ondulación de banda de paso y 30 dB de atenuación de banda de parada. Calcule su respuesta en frecuencia.

[ze,pe,ke] = ellip(n,3,30,2*pi*fc,"s");
[be,ae] = zp2tf(ze,pe,ke);
[he,we] = freqs(be,ae,4096);

Diseñe un filtro de Bessel de 5.º orden con la misma frecuencia de borde. Calcule su respuesta en frecuencia.

[zf,pf,kf] = besself(n,2*pi*fc);
[bf,af] = zp2tf(zf,pf,kf);
[hf,wf] = freqs(bf,af,4096);

Represente la atenuación en decibelios. Exprese la frecuencia en gigahercios. Compare los filtros.

plot([wb w1 w2 we wf]/(2e9*pi), ...
    mag2db(abs([hb h1 h2 he hf])))
axis([0 5 -45 5])
grid
xlabel("Frequency (GHz)")
ylabel("Attenuation (dB)")
legend(["butter" "cheby1" "cheby2" "ellip" "besself"])

Figure contains an axes object. The axes object with xlabel Frequency (GHz), ylabel Attenuation (dB) contains 5 objects of type line. These objects represent butter, cheby1, cheby2, ellip, besself.

Los filtros Butterworth y Chebyshev Tipo II tienen bandas de paso planas y bandas de transición anchas. Los filtros Chebyshev de tipo I y elíptico se desinflan más rápido pero tienen ondulaciones en la banda de paso. La entrada de frecuencia a la función de diseño Chebyshev Tipo II establece el comienzo de la banda de parada en lugar del final de la banda de paso. El filtro de Bessel tiene un retardo de grupo aproximadamente constante a lo largo de la banda de paso.

Argumentos de entrada

contraer todo

Orden del filtro, especificado como un escalar entero. En los diseños de paso banda y eliminador de banda, n representa la mitad del orden del filtro.

Tipos de datos: double

Ondulación pico a pico de banda de paso, especificada como un escalar positivo en decibelios.

Si su especificación, ℓ, está en unidades lineales, puede convertirla a decibelios utilizando Rp = 40 log10((1+ℓ)/(1–ℓ)).

Tipos de datos: double

Atenuación de la banda de parada desde el valor pico de la banda de paso, especificada como un escalar positivo en decibelios.

Si su especificación, ℓ, está en unidades lineales, puede convertirla a decibelios utilizando Rs = –20 log10ℓ.

Tipos de datos: double

Frecuencia de borde de banda de paso, especificada como un escalar o un vector de dos elementos. La frecuencia de borde de la banda de paso es la frecuencia a la que la respuesta de magnitud del filtro es Rp decibelios. Con valores menores de ondulación de la banda de paso, Rp, y valores mayores de atenuación de la banda de parada, Rs, se obtienen bandas de transición más anchas.

  • Si Wp es un escalar, ellip diseña un filtro paso bajo o paso alto con una frecuencia de borde Wp.

    Si Wp es el vector de dos elementos [w1 w2], donde w1 < w2, ellip diseña un filtro paso banda o eliminador de banda con una frecuencia de borde inferior w1 y una frecuencia de borde superior w2.

  • Para los filtros digitales, las frecuencias de borde de banda de parada deben estar entre 0 y 1, donde 1 corresponde a la tasa de Nyquist, la mitad de la tasa de muestreo o π rad/muestra.

    Para los filtros analógicos, las frecuencias de borde de banda de paso deben expresarse en radianes por segundo y pueden tomar cualquier valor positivo.

Tipos de datos: double

Tipo de filtro, especificado como uno de los siguientes:

  • 'low' especifica un filtro paso bajo con frecuencia de borde de banda de paso Wp. 'low' es el valor predeterminado para el escalar Wp.

  • 'high' especifica un filtro paso alto con frecuencia de borde de banda de paso Wp.

  • 'bandpass' especifica un filtro paso banda de orden 2n si Wp es un vector de dos elementos. 'bandpass' es el valor por defecto cuando Wp tiene dos elementos.

  • 'stop' especifica un filtro eliminador de banda de orden 2n si Wp es un vector de dos elementos.

Argumentos de salida

contraer todo

Coeficientes de la función de transferencia del filtro, devueltos como vectores fila de longitud n + 1 para los filtros paso bajo y paso alto y 2n + 1 para los filtros paso banda y eliminador de banda.

  • Para los filtros digitales, la función de transferencia se expresa en términos de b y a como

    H(z)=B(z)A(z)=b(1)+b(2)z1++b(n+1)zna(1)+a(2)z1++a(n+1)zn.

  • Para los filtros analógicos, la función de transferencia se expresa en términos de b y a como

    H(s)=B(s)A(s)=b(1)sn+b(2)sn1++b(n+1)a(1)sn+a(2)sn1++a(n+1).

Tipos de datos: double

Ceros, polos y ganancia del filtro, devueltos como dos vectores columna de longitud n (2n en los diseños de paso banda y eliminador de banda) y un escalar.

  • Para los filtros digitales, la función de transferencia se expresa en términos de z, p y k como

    H(z)=k(1z(1)z1)(1z(2)z1)(1z(n)z1)(1p(1)z1)(1p(2)z1)(1p(n)z1).

  • Para los filtros analógicos, la función de transferencia se expresa en términos de z, p y k como

    H(s)=k(sz(1))(sz(2))(sz(n))(sp(1))(sp(2))(sp(n)).

Tipos de datos: double

Representación de espacios de estados del filtro, devuelta como matrices. Si m = n para los diseños de paso bajo y paso alto, y m = 2n para los filtros paso banda y eliminador de banda, entonces A es m × m, B es m × 1, C es 1 × m y D es 1 × 1.

  • Para los filtros digitales, las matrices del espacio de estado relacionan el vector de estado x, la entrada u y la salida y a través de

    x(k+1)=Ax(k)+Bu(k)y(k)=Cx(k)+Du(k).

  • Para los filtros analógicos, las matrices del espacio de estado relacionan el vector de estado x, la entrada u y la salida y a través de

    x˙=Ax+Buy=Cx+Du.

Tipos de datos: double

Más acerca de

contraer todo

Limitaciones

Inestabilidad numérica de la sintaxis de la función de transferencia

En general, utilice la sintaxis [z,p,k] para diseñar filtros IIR. Para analizar o implementar su filtro, puede entonces utilizar la salida [z,p,k] con zp2sos. Si diseña el filtro utilizando la sintaxis [b,a] puede encontrar problemas numéricos. Estos problemas se deben a errores de redondeo y pueden ocurrir para n tan bajos como 4. El siguiente ejemplo ilustra esta limitación.

n = 6; 
Rp = 0.1;
Rs = 80;
Wn = [2.5e6 29e6]/500e6; 
ftype = 'bandpass'; 

% Transfer Function design 
[b,a] = ellip(n,Rp,Rs,Wn,ftype);            % This filter is unstable 

% Zero-Pole-Gain design 
[z,p,k] = ellip(n,Rp,Rs,Wn,ftype);
sos = zp2sos(z,p,k);

% Plot and compare the results
hfvt = fvtool(b,a,sos,'FrequencyScale','log'); 
legend(hfvt,'TF Design','ZPK Design')

Figure Figure 1: Magnitude Response (dB) contains an axes object. The axes object with title Magnitude Response (dB), xlabel Normalized Frequency ( times pi blank rad/sample), ylabel Magnitude (dB) contains 2 objects of type line. These objects represent TF Design, ZPK Design.

Algoritmos

Los filtros elípticos ofrecen características de descenso más pronunciadas que los filtros Butterworth o Chebyshev, pero son equiripple tanto en la banda de paso como en la banda de parada. En general, los filtros elípticos cumplen determinadas especificaciones de rendimiento con el menor orden de cualquier tipo de filtro.

ellip utiliza un algoritmo de cinco pasos:

  1. Encuentra los polos, ceros y ganancia del prototipo analógico de paso bajo utilizando la función ellipap.

  2. Convierte los polos, los ceros y la ganancia en forma de espacio de estado.

  3. Si es necesario, utiliza una transformación del espacio de estados para convertir el filtro paso bajo en un filtro paso banda, paso alto o eliminador de banda con las restricciones de frecuencia deseadas.

  4. Para el diseño de filtros digitales, utiliza bilinear para convertir el filtro analógico en un filtro digital a través de una transformación bilineal con preconfiguración de frecuencia. Un cuidadoso ajuste de la frecuencia permite que los filtros analógicos y los digitales tengan la misma magnitud de respuesta en frecuencia en Wp o w1 y w2.

  5. Convierte el filtro de espacio de estados de nuevo en una función de transferencia o en un formato cero-polo-ganancia, según sea necesario.

Capacidades ampliadas

Historial de versiones

Introducido antes de R2006a