prony

The Prony method for designing filters.

Library

EngeeDSP

Syntax

Function call

b,a = prony(h,bord,aord) — returns the coefficients of the numerator and denominator for a causal rational transfer function with an impulse response h, by the order of the numerator bord and the order of the denominator aord.

Arguments

Input arguments

# h — impulse response

+ vector

Details

The impulse response, given as a vector.

Типы данных	`Float32`, `Float64`
Support for complex numbers	Yes

# bord,aord — orders of numerator and denominator

+ positive integer scalar

Details

The orders of the numerator and denominator, given as positive integer scalars. If the length h less max(bord,aord), then the function fills the impulse response with zeros.

If you need a transfer function with poles only, specify for bord value 0.
If you need a transfer function with only zeros, specify for aord value 0.

Типы данных

Float32, Float64

Output arguments

# b,a are the coefficients of the transfer function

+ vector

Details

The coefficients of the transfer function returned as vectors. Vector b has a length of bord+1 and a has a length of aord+1.

Additional information

Transfer function

Details

The transfer function is the Z—transform of the impulse response :

rational transfer function is the ratio of polynomials in negative powers . This equation describes a causal rational transfer function with numerator order and the order of the denominator :

Where .

Literature

Parks, Thomas W., and C. Sidney Burrus. Digital Filter Design. New York, NY, USA: Wiley-Interscience, 1987.