How can I force Gaussian profile to fit peaks completely?

11 次查看（过去 30 天）
Naif Alsalem
Naif Alsalem 2022-6-18
回答： Image Analyst 2022-6-18
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

请先登录，再回答此问题。

采纳的回答

Image Analyst
Image Analyst 2022-6-18
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.

更多回答（1 个）

Sam Chak
Sam Chak 2022-6-18
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.

类别

在 Help Center 和 File Exchange 中查找有关 Get Started with Curve Fitting Toolbox 的更多信息

标签

产品


版本

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by