Fit NonLinearLeastSquares to data and constrain curve to pass through point (0,1)

Question

Zachary Nunn 2019-8-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/476324-fit-nonlinearleastsquares-to-data-and-constrain-curve-to-pass-through-point-0-1

编辑： Matt J 2019-8-19

采纳的回答： Jyotsna Talluri

在 MATLAB Online 中打开

I have data I need to fit to an equation, which I can do, but I want the equation to pass through point (0,1).

The equation is an exponential y = a*exp(-x/b)+ c*exp(-x/d)+ e

Currently, this is my code:

f3 = fitoptions('Method','NonlinearLeastSquares','Startpoint',[1,100,1,100,1])
newrelax = fittype('q*exp(-x/r)+ s*exp(-x/t)+u','options',f3);
[h,gof] = fit(timeR,stressR,newrelax)
plot(h,timeR,stressR)

Thank you,

-Zach

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Torsten 2019-8-19

Use

y = a*exp(-x/b) + c*exp(-x/d) + (1-(a+c))

as equation.

请先登录，再进行评论。

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

Answer 1

Jyotsna Talluri 2019-8-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/476324-fit-nonlinearleastsquares-to-data-and-constrain-curve-to-pass-through-point-0-1#answer_388163

在 MATLAB Online 中打开

You can use a “lsqlin” function from optimization toolbox.

t=coeffvalues(h); 
C=[exp(-x/t(2)) exp(-x/t(4)) ones(size(x))]; 
D=y; 
A = []; % No inequality constraint 
B= []; % No inequality constraint 
%Set the specified points in C and D to find equality constraints Aeq, Beq..By substituting (0,1) they turned out to be 
Aeq =[1 1 1];Beq=[1]; 
l=lsqlin(C,D,A,B,Aeq,Beq);     %generates coefficients of the curve fitted passing through a specified point  
Plot(timeR,C*l); 
 

Refer to the below link

https://www.mathworks.com/help/optim/ug/lsqlin.html

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Matt J 2019-8-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/476324-fit-nonlinearleastsquares-to-data-and-constrain-curve-to-pass-through-point-0-1#answer_388191

编辑：Matt J 2019-8-19

在 MATLAB Online 中打开

This answer incorporates Torsten's advice, but I also think you should re-formulate the model to make the exponential terms asymmetric. Otherwise, the solver cannot decide which exponential term belongs to q and which belongs to s.

f3 = fitoptions('Method','NonlinearLeastSquares',...
     'Startpoint',[1,100,1,100,1],'Lower',[-inf,-inf,-inf,0]);
newrelax = fittype('q*exp(-x/r)+ s*exp(-x/(r+d))+(1-(q+s))','options',f3);