Esta página aún no se ha traducido para esta versión. Puede ver la versión más reciente de esta página en inglés.
Diseño de filtro FIR de fase lineal de menos cuadrados
b = firls(n,f,a)
b = firls(n,f,a,w)
b = firls(n,f,a,'ftype'
)
b = firls(n,f,a,w,'ftype'
)
diseña unfirls
fase lineal Filtro FIR que minimiza el error cuadrado ponderado e integrado entre una función lineal por partes ideal y la respuesta de magnitud del filtro sobre un conjunto de bandas de frecuencia deseadas.
b = firls(n,f,a)
devuelve el vector de fila que contiene los coeficientes del filtro FIR de la orden cuyas características de amplitud de frecuencia coinciden aproximadamente con las dadas por los vectores y .b
n+1
n
f
a
Los coeficientes de filtro de salida, o "toques", en obedecer la relación de simetría.b
Estos son filtros de fase lineal tipo I (impar) y tipo II (par).n
n
Vectores y especifique las características de amplitud de frecuencia del filtro:f
a
es un vector de pares de puntos de frecuencia, especificado en el rango entre 0 y 1, donde 1 corresponde a la frecuencia Nyquist.f
Las frecuencias deben estar en orden creciente. Se permiten puntos de frecuencia duplicados y, de hecho, se pueden utilizar para diseñar un filtro exactamente igual que los devueltos por el y funciones con una ventana rectangular ( ).fir1
fir2
rectwin
es un vector que contiene la amplitud deseada en los puntos especificados en .a
f
La función de amplitud deseada a frecuencias entre pares de puntos ( ( ( ), ( +1)) para impar es el segmento de línea que conecta los puntos ( ( ( ), ) y ( +1), ( +1)).fkfkkfkakfkak
La función de amplitud deseada en frecuencias entre pares de puntos ( ( ), ( +1)) para par no está especificada.fkfkk Estas son regiones de transición o "no le importa".
y tienen la misma longitud.f
a
Esta longitud debe ser un número par.
siempre utiliza un orden de filtro uniforme para configuraciones con una banda de paso en la frecuencia Nyquist.firls
Esto se debe a que para las órdenes impares, la respuesta de frecuencia en la frecuencia Nyquist es necesariamente 0. Si especifica un valor impar , lo incrementa en 1.n
firls
La figura siguiente ilustra la relación entre los vectores y en la definición de una respuesta de amplitud deseada.f
a
b = firls(n,f,a,w)
utiliza el pesos en vector para ponderar el ajuste en cada banda de frecuencia.w
La longitud de es la mitad de la longitud de y , por lo que hay exactamente un peso por banda.w
f
a
b = firls(n,f,a,
Y'ftype'
)
b = firls(n,f,a,w,
especificar un tipo de filtro, donde 'ftype'
)'ftype'
Es:
para filtros de fase lineal con simetría impar (tipo III y tipo IV).'hilbert'
Los coeficientes de producción en obedecer la relaciónb
Esta clase de filtros incluye el Transformador Hilbert, que tiene una amplitud deseada de 1 en toda la banda.
'differentiator'
para filtros tipo III y tipo IV, utilizando una técnica de ponderación especial. Para bandas de amplitud distintas de cero, el error cuadrado integrado tiene un peso de (1/ )f2 de modo que el error a bajas frecuencias sea mucho menor que en frecuencias altas. Para los diferenciadores FIR, que tienen una característica de amplitud proporcional a la frecuencia, los filtros minimizan el error cuadrado integrado relativo (la integral del cuadrado de la relación del error con la amplitud deseada).
Se muestra uno de los siguientes mensajes de diagnóstico cuando se utiliza un argumento incorrecto:
F debe ser de longitud uniforme. F y A deben tener longitudes iguales. Requiere simetría para ser 'hilbert' o 'diferenciador'. Requiere un peso por banda. Las frecuencias en F deben ser no disminuidos. Las frecuencias en F deben estar en el rango [0,1].
Un mensaje de advertencia más serio es
Warning: Matrix is close to singular or badly scaled.
Esto tiende a suceder cuando el producto de la longitud del filtro y el ancho de transición crece. En este caso, es posible que los coeficientes de filtro b no representen el filtro deseado. Puede comprobar el filtro mirando su respuesta de frecuencia.
La referencia describe el enfoque teórico detrás de .[1]firls
La función resuelve un sistema de ecuaciones lineales que implican una matriz de producto interna de tamaño aproximadamente utilizando el operador.n/2
MATLAB®\
Esta función diseña filtros de fase lineal tipo I, II, III y IV. Tipo I y II son los valores predeterminados para n pares e impares respectivamente, mientras que las banderas y producen filtros de tipo III (n par) y IV (n impares).'hilbert'
'differentiator'
Los distintos tipos de filtro tienen diferentes simetrías y restricciones en sus respuestas de frecuencia (consulte para obtener más información).[2]
Tipo de filtro de fase lineal | Orden de filtro | Simetría de coeficientes | Respuesta H(f), f - 0 | Respuesta H(f), f á 1 (Nyquist) |
---|---|---|---|---|
Tipo I | Incluso | Sin restricciones | Sin restricciones | |
Tipo II | Extraño | Sin restricciones | (1) 0H | |
Tipo III | Incluso | (0) á 0H | (1) 0H | |
Tipo IV | Extraño | (0) á 0H | Sin restricciones |
[1] Parks, Thomas W., and C. Sidney Burrus. Digital Filter Design. New York: John Wiley & Sons, 1987, pp. 54–83.
[2] Oppenheim, Alan V., Ronald W. Schafer, and John R. Buck. Discrete-Time Signal Processing. Upper Saddle River, NJ: Prentice Hall, 1999.
fir1
| fir2
| firpm
| rcosdesign