Deconvolution of a combination of Gaussians signal
조회 수: 53 (최근 30일)
이전 댓글 표시
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
Bjorn Gustavsson
2021년 5월 6일
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.
채택된 답변
Bjorn Gustavsson
2021년 5월 6일
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
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Signal Processing Toolbox에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!