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.

prony

Método Prony para el diseño de filtros

Sintaxis

[Num,Den] = prony(impulse_resp,num_ord,denom_ord)

Descripción

[Num,Den] = prony(impulse_resp,num_ord,denom_ord) devuelve los coeficientes numerador y denominador para una función de sistema racional causal con respuesta de impulso.NumDenimpulse_resp La función del sistema tiene orden de numerador y orden de denominador.num_orddenom_ord Las longitudes de y son y .NumDennum_ord+1denom_ord+1 Si la longitud de es menor que el orden más grande ( o ), se rellena con ceros.impulse_respnum_orddenom_ordimpulse_resp Introduzca 0 para una función de sistema de todos los polos.num_ord Para una función de sistema de cero, introduzca un 0 para .denom_ord

Ejemplos

contraer todo

Ajuste un modelo IIR de 4o orden a la respuesta de impulso de un filtro de paso bajo. Trazar las respuestas de impulso originales y diseñados por Prony.

d = designfilt('lowpassiir','NumeratorOrder',4,'DenominatorOrder',4, ...     'HalfPowerFrequency',0.2,'DesignMethod','butter');  impulse_resp = filter(d,[1 zeros(1,31)]); denom_order = 4; num_order = 4; [Num,Den] = prony(impulse_resp,num_order,denom_order);  subplot(2,1,1)  stem(impz(Num,Den,length(impulse_resp))) title 'Impulse Response with Prony Design'  subplot(2,1,2) stem(impulse_resp) title 'Input Impulse Response'

Ajuste un modelo FIR de décimo orden a la respuesta de impulso de un filtro de paso alto. Trazar las respuestas de frecuencia originales y diseñadas por Prony.

d = designfilt('highpassfir', 'FilterOrder', 10, 'CutoffFrequency', .8);  impulse_resp = filter(d,[1 zeros(1,31)]); num_order = 10; denom_order = 0; [Num,Den] = prony(impulse_resp,num_order,denom_order);  fvt = fvtool(Num,Den,d); legend(fvt,'Prony','Original')

Más acerca de

contraer todo

Función del sistema

Es la transformación z de la respuesta de impulsofunción del sistema h[n]:

H(z)=n=h(n)zn.

A es una proporción de polinomios enfunción racional del sistema z–1. Por convención, el polinomio numerador es B(z) y el denominador es A(z). La siguiente ecuación describe una función del sistema racional causal del orden del numerador y el orden del denominador:qp

H(z)=B(z)A(z)=k=0qb(k)zk1+l=1pa(l)zl,

Dónde a[0] = 1.

Referencias

Parks, Thomas W., y C. Sidney Burrus. .Diseño de filtro digital Nueva York: John Wiley & Sons, 1987, págs. 226–228.

Consulte también

| | |

Introducido antes de R2006a