Zero-pole-gain to second-order section format
Convert digital filter zero-pole-gain data to second-order section form.
zp2sos(z, p, g = 1, order = c("down", "up"))z |
complex vector of the zeros of the model (roots of |
p |
complex vector of the poles of the model (roots of |
g |
overall gain ( |
order |
row order, specified as:
The ordering influences round-off noise and the probability of overflow. |
A list with the following list elements:
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.
Overall gain factor that effectively scales the output b
vector (or any one of the input Bi vectors).
Julius O. Smith III, jos@ccrma.stanford.edu.
Conversion to R by Geert van Boxtel, gjmvanboxtel@gmail.com
zpk <- tf2zp (c(1, 0, 0, 0, 0, 1), c(1, 0, 0, 0, 0, .9)) sosg <- zp2sos (zpk$z, zpk$p, zpk$g)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.