Fit multiple data with one mathmatical model/function

Question

AHMED FAKHRI 2020-2-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/506852-fit-multiple-data-with-one-mathmatical-model-function

评论： the cyclist 2020-2-24

采纳的回答： the cyclist

mydata.xlsx

Hi

I have the following Excel file which contain my data.

The input data is the number of cycles (N) while the output data is the capacity at different temepratures (278.15, 298.15,308.15, 318.15 ).

I want to fit the input and output data with a single function or mathmatical formula, preferably the following equation:

C=a1*(a2-a3*N^0.5) *exp(((-a4/8.314)*(1/T-1/298.15))-((a5/8.314)*(1/T-1/298.15)^2))

Therfore, what I want to do is to get C that correspond to (T and N) everytime. a1-a5 are free fitting parameters.

For T is is clear that I want it to chage only from 278.15 to 318.15.

Can I solve this using lsqcurvefit? or any other better appraoch?

Thanks

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

请先登录，再进行评论。

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

Answer 1

the cyclist 2020-2-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/506852-fit-multiple-data-with-one-mathmatical-model-function#answer_416721

在 MATLAB Online 中打开

mydata.xlsx

I would restructure my data (see attached), and use the nlinfit function from the Machine Learning and Statistics Toolbox.

This way of structuring the data is known as tidy data.

tbl = readtable('mydata.xlsx');
C = @(a,X) a(1).*(a(2)-a(3).*X(:,1).^0.5).*exp(((-a(4)/8.314)*(1./X(:,2)-1/298.15))-((a(5)/8.314).*(1./X(:,2)-1/298.15).^2));
% Define initial fit coefficients
beta0 = [1 100 1 -50 100];
% Fit the model
mdl = fitnlm(tbl,C,beta0);

There are some subtleties here regarding whether there are correlated errors among these measurements. The method here assumes that the errors in each row are independent from each other.

4 个评论
显示 2更早的评论隐藏 2更早的评论

AHMED FAKHRI 2020-2-21

Thanks. What kind of warning you are referring to that the model produces?

Also, the RMSE now is 2 , do you see this acceptable ?

Many thanks

the cyclist 2020-2-21

在 MATLAB Online 中打开

This is the warning:

Warning: The Jacobian at the solution is ill-conditioned, and some model parameters may not be
estimated well (they are not identifiable).  Use caution in making predictions. 
> In nlinfit (line 384)
  In NonLinearModel/fitter (line 1127)
  In classreg.regr/FitObject/doFit (line 94)
  In NonLinearModel.fit (line 1434)
  In fitnlm (line 99)

I don't think I can really comment on the error question. It's too complex a decision about whether a model is "good enough" for a particular application.

请先登录，再进行评论。

Answer 2

Alex Sha 2020-2-22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/506852-fit-multiple-data-with-one-mathmatical-model-function#answer_416823

在 MATLAB Online 中打开

Hello, AHMED, your fit function is overdeterminted which will lead to multi-solutions, the parameter of "a1" is redundant:

C=a1*(a2-a3*N^0.5)*exp(((-a4/8.314)*(1/T-1/298.15))-((a5/8.314)*(1/T-1/298.15)^2));

the effect of above is exact same as the function of below

C=(a2-a3*N^0.5)*exp(((-a4/8.314)*(1/T-1/298.15))-((a5/8.314)*(1/T-1/298.15)^2));

Root of Mean Square Error (RMSE): 2.01307596538931
Sum of Squared Residual: 405.247484242812
Correlation Coef. (R): 0.962621756201578
R-Square: 0.926640645512611
Adjusted R-Square: 0.920394018937605
Determination Coef. (DC): 0.762037250139209
F-Statistic: 19.6453235739073
Parameter          	Best Estimate
--------------------	-------------
a2	102.235245663151
a3	1.13255810875212
a4	844.097603208569
a5	13460420.8019471

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

请先登录，再进行评论。

Answer 3

AHMED FAKHRI 2020-2-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/506852-fit-multiple-data-with-one-mathmatical-model-function#answer_416853

Thank you both Alex Sha and the cyclist

I have removed a(1) parameter and the error/warning written in the cyclist comment has been mitigated.

I will try to enhance the model further and reduce the RMSE by manipulating the function. After that I will test it.

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

请先登录，再进行评论。

Answer 4

AHMED FAKHRI 2020-2-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/506852-fit-multiple-data-with-one-mathmatical-model-function#answer_417108

Thanks guys, it worked. I only need to reduce the RMSE by modifying the function.

Is there anyway that we give the data and a 'tool' or 'algorithim' suggests a mathmatical model ?

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

AHMED FAKHRI 2020-2-24

Thanks. Any suggestions to improve the RMSE? such as making the nlinfit roboust ? currently the RMSE is 2.18. I am not conveinced exactly for the fit of the 5C ( 278.15) data.

the cyclist 2020-2-24

Usually one chooses the fitting function based on understanding of the underlying process. So, for example, is there some other factor besides number of cycles and temperature, that you are not accounting for? Is there something else that happens at lower temperature, that you can capture with another term?

Of course, it is possible that you just have noise that you cannot really account for. Then your RMSE is simply the best you can do.

请先登录，再进行评论。

Fit multiple data with one mathmatical model/function

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

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（3 个）

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

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

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

另请参阅

类别

标签

Community Treasure Hunt

Fit multiple data with one mathmatical model/function

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

采纳的回答

4 个评论 显示 2更早的评论隐藏 2更早的评论

更多回答（3 个）

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

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

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

另请参阅

类别

标签

Community Treasure Hunt

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

4 个评论
显示 2更早的评论隐藏 2更早的评论

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

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

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