least square curve fitting on a nonlinear equation , set of data available
9 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Mohammad
el 17 de Feb. de 2023
Editada: Star Strider
el 17 de Feb. de 2023
Hi all, I want to use least square curve fitting on a nonlinear equation , set of data available , I did the analytics on paper and want to be sure using pure code.
to summrize the method;
we first partial derivative with respect to a and equate it to 0
then partial derivative with respect to b and equate it to 0 as well.
we end up with 2 equations ( nonlinear).
how to solve and find a and b using code? I am new to matlab, any help would be appriciated.
the equation is where y and x are the data set and it goes 27 iterations.
Thanks.
0 comentarios
Respuesta aceptada
Star Strider
el 17 de Feb. de 2023
Editada: Star Strider
el 17 de Feb. de 2023
What you appear to be describing is the derivation of a linear least square regression. In that context, ‘linear’ implies ‘linear in the parameters’, such the partial derivatives of the objective function with respect to each parameter are not functions of that parameter or of any other parameters.
Your data are best determined by fminsearch. There are other options, however those require the Optimization Toolbox, the Statistics and Machine Learning Toolbox or Curve Fitting Toolbox.
Example —
xv = 1:10;
yv = sqrt(xv);
objfcn = @(b,x) b(1).*(1-exp(-b(2).*x));
B0 = rand(2,1);
[B,fval] = fminsearch(@(b)norm(yv-objfcn(b,xv)), B0)
figure
plot(xv,yv,'x', 'DisplayName','Data')
hold on
plot(xv, objfcn(B,xv), '-r', 'DisplayName','Regression Fit')
hold off
grid
xlabel('x')
ylabel('y')
legend('Location','best')
text(1,3, sprintf('$y = %.2f\\ (1-e^{-%.2f \\ x})$',B), 'Interpreter','latex')
EDIT — (17 Feb 2023 at 16:58)
Corrected typographical errors.
.
0 comentarios
Ver también
Categorías
Más información sobre Interpolation en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!