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 de numerador y denominador para una función de sistema causal racional con respuesta de impulso.NumDenimpulse_resp La función del sistema dispone de orden 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 en para una función de sistema All-Pole.num_ord Para una función de sistema All-Zero, introduzca un 0 para.denom_ord

Ejemplos

contraer todo

Ajuste un modelo IIR de 4º orden a la respuesta de impulso de un filtro de paso bajo. Trace las respuestas de impulso originales y diseñadas 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 10 órdenes a la respuesta de impulso de un filtro de paso alto. Trace 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

La transformación z de la respuesta al 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 de sistema causal racional de orden numerador y orden denominador:qp

H(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