Deconvolution of a combination of Gaussians signal

51 visualizaciones (últimos 30 días)
VD
VD el 5 de Mayo de 2021
Comentada: VD el 9 de Jun. de 2021
Hi all,
First of all, thanks for the support in the forum, I`ve learned a lot through the years reading here.
I'm currently trying to deconvolute a signal with three peaks (which are assumed to be gaussian functions with exponential decay). I'm basically trying to do the opposite to the following file Exchange function: https://es.mathworks.com/matlabcentral/fileexchange/74408-fit-multiple-gaussians, by Image Analyst.
The main objective is to find the 'original' signal taking into account the contribution of all peaks. ¿Is there any function in MATLAB, or toolbox, or file-exchange code, dealing with this? I'm trying to do my own code, using the function fit, with gauss3, etc, but I'm still not able to do it, since the fit only retains the first two peaks (maybe due to the start points, etc, that I give to fit.) My data looks something like this:
I'm afraid I cannot provide an exact figure, since it is part of a work still not published.
Any help is appreciated!
Thanks in advance.
  1 comentario
Bjorn Gustavsson
Bjorn Gustavsson el 6 de Mayo de 2021
If you understand your problem well enough you can mock up a similar enough data-set (similar original data convolved with some PSF possibly with an additional noise with similar statistical properties of your secret signal) without revealing your secret data-set.

Iniciar sesión para comentar.

Respuesta aceptada

Bjorn Gustavsson
Bjorn Gustavsson el 6 de Mayo de 2021
If you have a signal that is smeared or convolved with some point-spread-function you can to some extent undo that by using either of the deconvolution-functions in matlab: deconvwnr, deconvreg, deconvlucy. Deconvolution-work is a bit fiddly in the sense that the improvement in resolution is not infinite (noise (either measurement or numerical) will be amplified, and when the PSF/convolving kernel has zeroes in its Fourier-transform those components will be lost.) so both care with respect to estimates of the PSF to deconvolve has to be made and tempered expectations of the deconvolution-result has to be accepted. Have a look at the help and documentation of the above mentioned functions, that should get you started. They are designed for deconvolution on images, but should work also for 1-D signals.
HTH
  1 comentario
VD
VD el 9 de Jun. de 2021
Sorry for the late response, I didn't see the notification of the answer. Thanks for your comments. I managed to try something by myself, but I'll take a look at those functions again.
Thanks again,

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Signal Processing Toolbox en Help Center y File Exchange.

Productos


Versión

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by