How can I force Gaussian profile to fit peaks completely?

Hi there,

I have 3 peaks that I fitted with gaussian profile. They are connected and it looks like the gaussian fit is not reaching the dips between the peaks.

I used this line of command to fit my peaks:

f = fit(x, y, 'gauss3');  

Here is the plot:

It would be nice to get the fit down to the dips between peaks.

Thanks

Answer 1

[Image Analyst]

I'm not sure you want to fit 3 Gaussians to this waveform which is not 3 Gaussians. WHY do you think you want to do that? Let's say you did it. What would you do with that information after that? Is it really the dividing line locations, widths, and plateau center locations you want?

For what it's worth, I'm attaching a multiGaussian fitting demo where you can fir some specified number of Gaussians to your data.

Answer 2

[Sam Chak]

Hi @Naif Alsalem

You can try specifying the bounds in fit options to help the algorithm to fit better...

options = fitoptions('gauss3')
options.Lower = [a1 b1 c1 a2 b2 c2 a3 b3 c3];   % specify the values for a1, b1, c1, a2, b2, c2, a3, b3, c3
options.Upper = [a1 b1 c1 a2 b2 c2 a3 b3 c3];

In your case, the amplitudes should be non-negative, and the centers lie around -0.1, -0.05, 0.01.

Manipulate the upper bounds of the non-negative widths with some small values until you get the desired fit to the dips.

However, the wavelength signal has 3 plateaus, while the gaussian function has only peaks.