Main Content

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