[SciPy-User] deconvolution of 1-D signals

[Lou Pecora]

I'm surprised that no one has brought it up, yet, but deconvolution has to be *very carefully* done.  Look up "ill-posed problems".  Yes, it can be a filter process, but just dividing by the filter FFT coefficients is dangerous since they approach zero (usually) as the frequency increases.  That's the ill-posed part and it has to be controlled.  

"Regularization" is the what people call the controlling process.  It's an extra assumption that one invokes on the deconvolution that limits the problems... if it's done right.  There are surely lots of books and articles on Regularization and Deconvolution.  Start with Google and Wikipedia and Google Scholar.  It's not hard, but it is not a one-step problem, but more like an optimization problem.  
 
-- Lou Pecora,   my views are my own.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20110731/320bcf02/attachment.html>