I am using a model of several gaussians to fit infrared absorption data. For this purpose, I define the resonances I expect at certain wavenumbers (i.e. the wavelength) using "NormalDistribution" and then use the sum of them as my model.
To find starting values, I use an interactive "Manipulate" environment which is working nicely. Using these, Mathematica finds the fit requested rather quick. Some of the parameters are running to values which are not physical. Therefore, I wanted to restrict these parameters using the constraints as described in the Documentation. As soon as I have a constraint on a single parameter, Mathematica takes much longer to compute a fit and the output is much worse or even absolutely stupid (meaning there is no "fit" at all!).
This behaviour I do not understand and would like to know how I can circumvent this situation (either using a different function than "NonlinearModelFit" or some other options.
Any opinion on this is welcome!
PS: I am working on a short example for you guys. I will provide it as soon as possible.