deconv(b, a)q and r.deconv(b,a) deconvolves vector a out of vector b.
The quotient is returned in vector q and the remainder in vector
r such that b = conv(a,q)+r. If b and a are vectors of polynomial coefficients,
convolving them is equivalent to multiplying the two polynomials,
and deconvolution is polynomial division.
conv, polymulb <- c(10, 40, 100, 160, 170, 120)
a <- c(1, 2, 3, 4)
p <- deconv(b, a)
p$q #=> 10 20 30
p$r #=> 0 0 0Run the code above in your browser using DataLab