Impulse invariance method for A/D filter conversion
Convert analog filter with coefficients b and a to digital, conserving impulse response.
impinvar(b, ...) ## S3 method for class 'Arma' impinvar(b, ...) ## Default S3 method: impinvar(b, a, fs = 1, tol = 1e-04, ...)
b
: coefficients of numerator polynomial...
: additional arguments (not used)a
: coefficients of denominator polynomialfs
: sampling frequency (Default: 1 Hz)tol
: tolerance. Default: 0.0001A list of class Arma
containing numerator and denominator polynomial filter coefficients of the A/D converted filter.
Because impinvar
is generic, it can also accept input of class Arma
.
f <- 2 fs <- 10 but <- butter(6, 2 * pi * f, 'low', 's') zbut <- impinvar(but, fs) freqz(zbut, n = 1024, fs = fs)
invimpinvar
Tony Richardson, arichard@stark.cc.oh.us ,
Ben Abbott, bpabbott@mac.com ,
adapted by John W. Eaton.
Conversion to R by Geert van Boxtel, G.J.M.vanBoxtel@gmail.com