Main Content

poly2rc

Convertir polinomio de filtro de predicción en coeficientes de reflexión

Sintaxis

k = poly2rc(a)
[k,r0] = poly2rc(a,efinal)

Descripción

k = poly2rc(a) convierte el polinomio del filtro de predicción en los coeficientes de reflexión de la estructura de celosía correspondiente. puede ser real o complejo, y no puede ser 0.aaa(1) Si no es igual a , poly2rc normaliza el polinomio del filtro de predicción por . es un vector de fila de tamaño .a(1)1a(1)klength(a)-1

[k,r0] = poly2rc(a,efinal) devuelve la autocorrelación de cero vueltas, , basada en el error de predicción final, .r0efinal

Ejemplos

contraer todo

Dado un polinomio de filtro de predicción, , y un error de predicción final, , determinar los coeficientes de reflexión de la estructura de celosía correspondiente y la autocorrelación de cero vueltas.aefinal

a = [1.0000 0.6149 0.9899 0.0000 0.0031 -0.0082]; efinal = 0.2; [k,r0] = poly2rc(a,efinal)
k = 5×1

    0.3090
    0.9801
    0.0031
    0.0081
   -0.0082

r0 = 5.6032 

Limitaciones

Si es por alguno, encontrar los coeficientes de reflexión es un problema mal condicionado. devuelve algunos s y proporciona un mensaje de advertencia en esos casos.abs(k(i)) == 1ipoly2rcNaN

Sugerencias

Una forma sencilla y rápida de comprobar si tiene todas sus raíces dentro del círculo de la unidad es comprobar si cada uno de los elementos de tiene una magnitud menor que 1.ak

stable = all(abs(poly2rc(a))<1) 

Algoritmos

implementa esta relación recursiva:poly2rc

k(n)=an(n)an1(m)=an(m)k(n)an(nm)1k(n)2,m=1,2,,n1

Esta relación se basa en la recursividad de Levinson.[1] Para implementarlo, recorre en bucle en orden inverso después de descartar su primer elemento.poly2rca Para cada iteración de bucle, la función:i

  1. Establece igual ak(i)a(i)

  2. Aplica la segunda relación anterior a los elementos 1 a través del vector .ia

    a = (a-k(i)*fliplr(a))/(1-k(i)^2); 

Referencias

[1] Kay, Steven M. Modern Spectral Estimation. Englewood Cliffs, NJ: Prentice-Hall, 1988.

Capacidades ampliadas

Consulte también

| | | | |

Introducido antes de R2006a