grpdelay

Retardo promedio del filtro (retardo de grupo)

Sintaxis

[gd,w] = grpdelay(b,a,n)

[gd,w] = grpdelay(sos,n)

[gd,w] = grpdelay(d,n)

[gd,w] = grpdelay(___,'whole')

[gd,f] = grpdelay(___,n,fs)

[gd,f] = grpdelay(___,n,'whole',fs)

gd = grpdelay(___,win)

gd = grpdelay(___,fin,fs)

grpdelay(___)

Descripción

[gd,w] = grpdelay(b,a,n) devuelve la respuesta de retardo de grupo de n puntos, gd, y el correspondiente vector de frecuencia angular, w, para el filtro digital con los coeficientes de la función de transferencia almacenados en b y a.

[gd,w] = grpdelay(sos,n) devuelve la respuesta de retardo de grupo de n puntos correspondiente a la matriz de secciones de segundo orden sos.

[gd,w] = grpdelay(d,n) devuelve la respuesta de retardo de grupo de n puntos del filtro digital d.

[gd,w] = grpdelay(___,'whole') devuelve el retardo de grupo en n puntos de muestra alrededor de todo el círculo de la unidad.

[gd,f] = grpdelay(___,n,fs) devuelve el vector de retardo de grupo gd y el correspondiente vector de frecuencia física f de un filtro digital diseñado para filtrar las señales muestreadas a una tasa fs.

[gd,f] = grpdelay(___,n,'whole',fs) devuelve el vector de frecuencia en n puntos que van de 0 a fs.

gd = grpdelay(___,win) devuelve el vector de retardo de grupo gd evaluado en las frecuencias normalizadas suministradas en win.

gd = grpdelay(___,fin,fs) devuelve el vector de retardo de grupo gd evaluado en las frecuencias físicas suministradas en fin.

grpdelay(___) sin argumentos de salida representa la respuesta de retardo de grupo del filtro.

Ejemplos

Retardo de grupo de un filtro Butterworth

Diseñe un filtro Butterworth de sexto orden con 3 dB normalizados de frecuencia de $$0.2\pi$$ rad/muestra. Utilice grpdelay para mostrar el retardo de grupo.

[z,p,k] = butter(6,0.2);
sos = zp2sos(z,p,k);

grpdelay(sos,128)

Represente tanto el retardo de grupo como el retardo de fase del sistema en la misma figura.

gd = grpdelay(sos,512);

[h,w] = freqz(sos,512);
pd = -unwrap(angle(h))./w;

plot(w/pi,gd,w/pi,pd)
grid
xlabel 'Normalized Frequency (\times\pi rad/sample)'
ylabel 'Group and phase delays'
legend('Group delay','Phase delay')

Respuesta de retardo de grupo de un filtro digital Butterworth

Utilice designfilt para diseñar un filtro Butterworth de sexto orden con 3 dB normalizados de frecuencia de $$0.2\pi$$ rad/muestra. Muestre su respuesta de retardo de grupo.

d = designfilt('lowpassiir','FilterOrder',6, ...
    'HalfPowerFrequency',0.2,'DesignMethod','butter');
grpdelay(d)

Respuesta de retardo de grupo de un filtro FIR de respuesta de magnitud arbitraria

Diseñe un filtro FIR de 88.º orden de respuesta de magnitud arbitraria. El filtro tiene dos bandas de paso y dos bandas de parada. La banda de paso de menor frecuencia tiene el doble de ganancia que la frecuencia de paso de mayor frecuencia. Especifique una tasa de muestreo de 200 Hz. Visualice la respuesta de magnitud y la respuesta de fase del filtro de 10 Hz a 78 Hz.

fs = 200;
d = designfilt('arbmagfir', ...
       'FilterOrder',88, ...
       'NumBands',4, ...
       'BandFrequencies1',[0 20], ...
       'BandFrequencies2',[25 40], ...
       'BandFrequencies3',[45 65], ...
       'BandFrequencies4',[70 100], ...
       'BandAmplitudes1',[2 2], ...
       'BandAmplitudes2',[0 0], ...
       'BandAmplitudes3',[1 1], ...
       'BandAmplitudes4',[0 0], ...
       'SampleRate',fs);
freqz(d,10:1/fs:78,fs)

Calcule y muestre la respuesta de retardo de grupo del filtro a lo largo del mismo rango de frecuencia. Compruebe que sea una mitad del orden del filtro.

filtord(d)

ans = 88

grpdelay(d,10:1/fs:78,fs)

Argumentos de entrada

b, a — Coeficientes de la función de transferencia

vectores

Coeficientes de la función de transferencia, especificados como vectores. Exprese la función de transferencia en términos de b y a como

$$H(e^{j\omega })=\frac{B(e^{j\omega })}{A(e^{j\omega })}=\frac{\text{b(1)}+\text{b(2)}{e}^{-j\omega }+\text{b(3)}{e}^{-j2\omega }+\cdots +\text{b(M)}{e}^{-j(M-1)\omega }}{\text{a(1)}+\text{a(2)}{e}^{-j\omega }+\text{a(3)}{e}^{-j2\omega }+\cdots +\text{a(N)}{e}^{-j(N-1)\omega }}.$$

Ejemplo: b = [1 3 3 1]/6 y a = [3 0 1 0]/3 especifican un filtro Butterworth de tercer orden con una frecuencia normalizada de 3 dB de 0,5π rad/muestra.

Tipos de datos: double | single
Soporte de números complejos: Sí

n — Número de puntos de evaluación

512 (predeterminado) | escalar entero positivo

Número de puntos de evaluación, especificado como un escalar entero positivo no inferior a 2. Cuando n está ausente, el valor predeterminado es 512. Para obtener los mejores resultados, establezca n en un valor mayor que el orden del filtro.

sos — Coeficientes de la sección de segundo orden

matriz

Coeficientes de la sección de segundo orden, especificados como una matriz. sos es una matriz de K por 6, donde el número de secciones, K, debe ser mayor o igual a 2. Si el número de secciones es inferior a 2, la función trata la entrada como un vector de numerador. Cada fila de sos corresponde a los coeficientes de un filtro de segundo orden (bicuadrado). La fila i de sos corresponde a [bi(1) bi(2) bi(3) ai(1) ai(2) ai(3)].

Ejemplo: s = [2 4 2 6 0 2;3 3 0 6 0 0] especifica un filtro Butterworth de tercer orden con una frecuencia normalizada de 3 dB de 0,5π rad/muestra.

Tipos de datos: double | single
Soporte de números complejos: Sí

d — Filtro digital

objeto digitalFilter

Filtro digital, especificado como un objeto digitalFilter. Utilice designfilt para generar un filtro digital según las especificaciones de frecuencia-respuesta.

Ejemplo: d = designfilt('lowpassiir','FilterOrder',3,'HalfPowerFrequency',0.5) especifica un filtro Butterworth de tercer orden con una frecuencia normalizada de 3 dB de 0,5π rad/muestra.

fs — Tasa de muestreo

escalar positivo

Tasa de muestreo, especificada como un escalar positivo. Cuando la unidad de tiempo es el segundo, fs se expresa en hercios.

Tipos de datos: double

win — Frecuencias angulares

vector

Frecuencias angulares, especificadas como un vector y expresadas en rad/muestra. win debe tener al menos dos elementos porque, de lo contrario, la función lo interpreta como n. win = π corresponde a la frecuencia de Nyquist.

fin — Frecuencias

vector

Frecuencias, especificadas como un vector. fin debe tener al menos dos elementos porque, de lo contrario, la función lo interpreta como n. Cuando la unidad de tiempo es el segundo, fin se expresa en hercios.

Tipos de datos: double

Argumentos de salida

gd — Respuesta de retardo de grupo

vector

Respuesta de retardo de grupo, devuelta como vector. Si especifica n, gd tiene una longitud n. Si no especifica n, o especifica n como el vector vacío, gd tiene una longitud 512.

Si la entrada de grpdelay es de precisión simple, la función calcula la respuesta de retardo de grupo mediante aritmética de precisión simple. La salida h es de una precisión simple.

w — Frecuencias angulares

vector

Frecuencias angulares, devueltas como un vector. w tiene valores que van de 0 a π. Si especifica 'whole' en su entrada, los valores de w van de 0 a 2π. Si especifica n, w tiene una longitud n. Si no especifica n, o especifica n como el vector vacío, w tiene una longitud 512.

f — Frecuencias

vector

Frecuencias, devueltas como un vector expresado en hercios. f tiene valores que van de 0 a fs/2 Hz. Si especifica 'whole' en su entrada, los valores de f van de 0 a fs Hz. Si especifica n, f tiene una longitud n. Si no especifica n, o especifica n como el vector vacío, f tiene una longitud 512.

Más acerca de

Retardo de grupo

La respuesta de retardo de grupo de un filtro es una medida del retardo promedio del filtro como función de frecuencia. Es la primera derivada negativa de la respuesta de fase del filtro. Si la respuesta en frecuencia de un filtro es H(e^jω), el retardo de grupo es

$${\tau }_g(\omega )=-\frac{d\theta (\omega )}{d\omega },$$

donde θ(ω) es la fase o argumento de H(e^jω).

Capacidades ampliadas

Generación de código C/C++

Genere código C y C++ mediante MATLAB® Coder™.

Notas y limitaciones de uso:

• El tamaño de la primera entrada debe ser una constante en tiempo de compilación.

Historial de versiones

Introducido antes de R2006a

Consulte también

cceps | designfilt | digitalFilter | fft | freqz | FVTool | hilbert | icceps | phasedelay | rceps

×

Abrir ejemplo

Tiene una versión modificada de este ejemplo. ¿Desea abrir este ejemplo con sus modificaciones?

No, sobrescribir la versión modificada Sí