gsignal (version 0.3-7)

grpdelay: Group delay

Description

Compute the average delay of a filter (group delay).

Usage

grpdelay(filt, ...)
# S3 method for grpdelay
print(x, ...)
# S3 method for grpdelay
plot(
  x,
  xlab = if (x$HzFlag) "Frequency (Hz)" else "Frequency (rad/sample)",
  ylab = "Group delay (samples)",
  type = "l",
  ...
)
# S3 method for default
grpdelay(filt, a = 1, n = 512, whole = FALSE, fs = NULL, ...)
# S3 method for Arma
grpdelay(filt, ...)
# S3 method for Ma
grpdelay(filt, ...)
# S3 method for Sos
grpdelay(filt, ...)
# S3 method for Zpg
grpdelay(filt, ...)

Value

A list of class grpdelay with items:

gd

the group delay, in units of samples. It can be converted to seconds by multiplying by the sampling period (or dividing by the sampling rate fs).

w

frequencies at which the group delay was calculated.

ns

number of points at which the group delay was calculated.

Hzflag

TRUE for frequencies in Hz, FALSE for frequencies in radians.

Arguments

filt

for the default case, the moving-average coefficients of an ARMA model or filter. Generically, filt specifies an arbitrary model or filter operation.

...

for methods of grpdelay, arguments are passed to the default method. For plot.grpdelay, additional arguments are passed through to plot.

x

object to be plotted.

xlab, ylab, type

as in plot, but with more sensible defaults.

a

the autoregressive (recursive) coefficients of an ARMA filter.

n

number of points at which to evaluate the frequency response. If n is a vector with a length greater than 1, then evaluate the frequency response at these points. For fastest computation, n should factor into a small number of small primes. Default: 512.

whole

FALSE (the default) to evaluate around the upper half of the unit circle or TRUE to evaluate around the entire unit circle.

fs

sampling frequency in Hz. If not specified, the frequencies are in radians.

Author

Paul Kienzle, ,
Julius O. Smith III, .
Conversion to R by Tom Short,
adapted by Geert van Boxtel,

Details

If the denominator of the computation becomes too small, the group delay is set to zero. (The group delay approaches infinity when there are poles or zeros very close to the unit circle in the z plane.)

References

https://ccrma.stanford.edu/~jos/filters/Numerical_Computation_Group_Delay.html
https://en.wikipedia.org/wiki/Group_delay

Examples

# Two Zeros and Two Poles
b <- poly(c(1 / 0.9 * exp(1i * pi * 0.2), 0.9 * exp(1i * pi * 0.6)))
a <- poly(c(0.9 * exp(-1i * pi * 0.6), 1 / 0.9 * exp(-1i * pi * 0.2)))
gpd <- grpdelay(b, a, 512, whole = TRUE, fs = 1)
print(gpd)
plot(gpd)