If you are facing instability in your nonlinear regression using fitnlm, you can try to stabilize the fitting process by setting better initial guesses for your parameters and by constraining the parameter space. Since you have an idea of the expected range of your parameters, you can use this information to set bounds and initial values for the fitting algorithm. fitnlm does not directly support parameter bounds, but you can transform your problem into an optimization problem using lsqcurvefit or fmincon which do support bounds.
Here's an example of how you can use lsqcurvefit to fit your nonlinear model with bounds:
% Define your model function
F = @(b,mu) ((mu.^((2*b(2))-1)).*exp(4*b(1)*((log(mu)).^2) + b(3)));
% Define your data (mu and ydata)
mu = ...; % your mu data
ydata = ...; % your ydata
% Set initial guesses for your parameters
b0 = [0.1, 0.5, 0.1]; % Example initial guesses (b1, b2, b3)
% Set lower and upper bounds for the parameters
lb = [-Inf, 0, 0]; % Lower bounds (b1 can be anything, b2 >= 0, b3 >= 0)
ub = [Inf, 1, Inf]; % Upper bounds (b1 can be anything, b2 <= 1, b3 can be anything)
% Set options for lsqcurvefit
options = optimoptions('lsqcurvefit', 'Algorithm', 'trust-region-reflective', ...
'Display', 'iter', 'MaxIterations', 400, ...
'FunctionTolerance', 1e-9, 'StepTolerance', 1e-9);
% Perform the fitting
[beta,resnorm,residual,exitflag,output] = lsqcurvefit(F, b0, mu, ydata, lb, ub, options);
% Display the results
disp('Fitted parameters:');
disp(beta);
In the above code:
- b0 is the vector of initial guesses for the parameters b1, b2, and b3.
- lb and ub are the vectors of lower and upper bounds for the parameters, respectively.
- options is the set of options for the lsqcurvefit function. You can adjust the MaxIterations, FunctionTolerance, and StepTolerance to improve stability and convergence.
- lsqcurvefit is used to perform the fitting within the specified bounds.
Remember to replace the ... with your actual data for mu and ydata. You can further tune the options to suit your specific problem. If lsqcurvefit does not provide satisfactory results, you can also try using fmincon, which offers even more flexibility in terms of defining constraints and customizing the optimization process.
If this answer helps you, please remember to accept the answer.
