sos2zp: Sos to zero-pole-gain

In gjmvanboxtel/gsignal: Signal Processing

Sos to zero-pole-gain

Description

Convert digital filter second-order section data to zero-pole-gain form.

Usage

sos2zp(sos, g = 1)

Arguments

sos	Second-order section representation, specified as an nrow-by-6 matrix, whose rows contain the numerator and denominator coefficients of the second-order sections: sos <- rbind(cbind(B1, A1), cbind(...), cbind(Bn, An)), where B1 <- c(b0, b1, b2), and A1 <- c(a0, a1, a2) for section 1, etc. The b0 entry must be nonzero for each section.
g	Overall gain factor that effectively scales the output b vector (or any one of the input B_i vectors). Default: 1.

Value

A list of class "Zpg" with the following list elements:

z

complex vector of the zeros of the model (roots of B(z))

p

complex vector of the poles of the model (roots of A(z))

k

overall gain (B(Inf))

Author(s)

Julius O. Smith III [e-mail].
Conversion to R by, Geert van Boxtel [e-mail]

See Also

filter

Examples

sos <- rbind(c(1, 0, 1, 1, 0, -0.81), c(1, 0, 0, 1, 0, 0.49))
zpk <- sos2zp(sos)

gjmvanboxtel/gsignal documentation built on Nov. 22, 2023, 8:19 p.m.