Search code examples
rpolynomial-mathdeconvolution

Deconvolution (polynomial division) in R


I need to perform a deconvolution of two vectors. This is equivalent to polynomial division. For example, in Matlab I would use deconv function. Is there a similar/equivalent function in R?


Solution

  • There is deconv in package pracma

    c = c(10, 40, 100, 160, 170, 120)
    u = c(1, 2, 3, 4)
    
    #install.packages("pracma")
    library(pracma)
    deconv(c, u)
    #$q
    #[1] 10 20 30
    
    #$r
    #[1] 0 0 0