Substract fitted values from original data

Question

Serhii Tetora 2020-9-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/601858-substract-fitted-values-from-original-data

评论： Ameer Hamza 2020-9-29

Mat.mat

I want to substract fit values from original data. Here is my code. What is wrong? Why Z and fitresult values are not equal? (see plot)

[xData, yData] = prepareCurveData(xdata,zz);
ft = fittype( 'sin4' );
opts = fitoptions( 'Method', 'NonlinearLeastSquares' );
opts.Normalize = 'on';
[fitresult, ~] = fit(xData, yData, ft, opts );
warning('off','curvefit:fit:equationBadlyConditioned')
x = xData;
y = yData;
varnames = coeffnames(fitresult)
varvalues = coeffvalues(fitresult)
for i = 1:length(varnames)
    eval([varnames{i},'=',num2str(varvalues(i)),';']);
end
Z = a1*sin(b1*x+c1) +a2*sin(b2*x+c2) +a3*sin(b3*x+c3) + a4*sin(b4*x+c4);
figure
plot(x,y)
hold on
plot(fitresult)
plot(x,Z)
legend('data','fit function','fit exact values')
%%
z_INT = Z_int-Z;

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

请先登录，再进行评论。

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

Answer 1

Ameer Hamza 2020-9-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/601858-substract-fitted-values-from-original-data#answer_502282

编辑：Ameer Hamza 2020-9-29

在 MATLAB Online 中打开

First, eval is evil; avoid it as much as possible: https://www.mathworks.com/matlabcentral/answers/304528-tutorial-why-variables-should-not-be-named-dynamically-eval . Also, in your case, there is no need, you can directly evaluate fitresult

y_predictions = fitresult(xData)

For the issue in your question. This seems to be a bug in the implementation of fit() function (unless I am overlooking something). For some reason, the fitresult does not report the coefficient value correctly. The issue resolves if you don't pass the opts structure to fit()

[fitresult, ~] = fit(xData, yData, ft);

You may consider filing a bug report: https://www.mathworks.com/support/bugreports/report_bug