Curve fitting toolbox error "NaN computed by model function, fitting cannot continue"

5 次查看（过去 30 天）

显示更早的评论

masoud meskin 2021-8-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/896292-curve-fitting-toolbox-error-nan-computed-by-model-function-fitting-cannot-continue

评论： masoud meskin 2021-8-16

Specimen_mean.csv

在 MATLAB Online 中打开

Hello,

I have a set of experimental data on which i am trying to fit 5 analytical functions (Hyperelastic models) using Levenberg-Marquardt algorithm. The models are as follows:

1) Neo-Hookean: y=2*c1*(x^2-1/x)

2) Yeoh: y=2*(x^2-1/x)*(c1+2*c2*(x^2+2/x-3)+3*c3*(x^2+2/x-3)^2)

3) Mooney-Rivlin: y=2*(x^2-1/x)*(c1+c2/x)

4) Ogden-N3: y=c1*(x^c2-2^(-1+c2)*x^(-c2/2))+c3*(x^c4-2^(-1+c4)*x^(-c4/2))+c5*(x^c6-2^(-1+c6)*x^(-c6/2))

5) Humphrey: y=2*(x^2-1/x)*c1*c2*exp(c2*(x^2+2/x-3))

The data file is attached so you can simply use the custom equation, use the above analytical functions, in the curve fitting toolbox.

clc; clear all; close all;
A = readmatrix('Specimen_mean.csv');
Xdata = A(:,1); % stretch [mm/mm]
Ydata = A(:,2); % Stress [Pa]

The 1,2 and 3 fit very well, though 4 and 5 cannot be fit to the data. The 4 gives the error "NaN computed by model function, fitting cannot continue" and the 5 gives the message "Fit computation did not converge". I have no idea how to deal with the mentioned problems. Any help would be highly appreciated.

Sincerely

Masoud

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

Alan Weiss 2021-8-16

在 MATLAB Online 中打开

I am not sure why you didn't use the Levenberg-Marquardt algorithm in lsqcurvefit when running MultiStart. That might save a step (no curve fitting toolbox needed).

options = optimoptions('lsqcurvefit','Algorithm','levenberg-marquardt');
problem = createOptimProblem('lsqcurvefit',...
'objective',model,'xdata',xData','ydata',yData','x0',[-1e6 -1e-3],'lb',[-1e7 -1],...
'ub',[1e7 1],'options',options);