I like the idea, unfortunately I cannot use it since I do not have (and cant afford) the DWT toolbox. :-/
23 Mar 2012
Peak Fitter
Command-line peak fitter for time-series signals. Version 3.6: February, 2013.
Author: Tom O'Haver
@Tom O'Haver, Nice detailed function, I have yet to study it but look forward to it, thanks for sharing. One thing though, in your function 'ExpBroaden' you list:
function yb = ExpBroaden(y,t)
% ExpBroaden(y,t) convolutes y by an exponential decay of time constant t
% by multiplying Fourier transforms and inverse transforming the result.
fy=fft(y);
a=exp(-[1:length(y)]./t);
fa=fft(a);
fy1=fy.*fa';
yb=real(ifft(fy1))./sum(a);
As it currently is coded, this will return the circular convolution, and not the 'usual' convolution. For a normal convolution, the result must be of length(a) + length(y) - 1, in which case the FFT's need to be zero padded for correct result. Otherwise, result is circular convolution. Was this your intent?
Thanks, Avigdor,
The #9 exppulse function itself (starting on line 1152 in version 3.6) is a unit-height function, as are all the shape functions in peakfit.m. The amplitude (peak height) is supplied by the fitting function ("fitexppulse") beginning on line 1140; it's called PEAKHEIGHT. That's because the amplitude is a linear variable that is calculated by linear regression (in line 1148) rather than by iteration using fminsearch.
Comment only
24 May 2013
Peak Fitter
Command-line peak fitter for time-series signals. Version 3.6: February, 2013.
Great script!
I was wondering if you could give an example of what is needed to modify an existing function. Specifically I am looking to add an amplitude variable to the #9 exppulse function.
Keep up the good work!
Comment only
29 Apr 2013
Peak Fitter
Command-line peak fitter for time-series signals. Version 3.6: February, 2013.
Hi Tom,
Thanks for your reply!
What do you think some basic 2D shape, such as Gaussian, elongated Gaussian (Gabor), high peak and (approximately) pole with half sphere on the top?
I will send you an example data later. The data belong to my old project.
Thanks in advance!
Ben
Comment only
26 Apr 2013
Peak Fitter
Command-line peak fitter for time-series signals. Version 3.6: February, 2013.
Ben, not so far, but it's an interesting possibility. Can you give me an example of a 2D data set with peaks and the kind of information your would like to obtain from a 2D peak fitter program?
Comment only
26 Apr 2013
Peak Fitter
Command-line peak fitter for time-series signals. Version 3.6: February, 2013.
Comment only