How to fix one point with no error while curve fitting

Question

aneps 2014-6-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/133385-how-to-fix-one-point-with-no-error-while-curve-fitting

评论： curoi 2015-1-14

I have a data to fit linear plot. But I want to fit the first point with no error. I mean, the fitted curve should fix on the first point and then fit to the rest of the point in 'least square fit' method. How can I do that?

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

请先登录，再进行评论。

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

Answer 1

Shashank Prasanna 2014-6-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/133385-how-to-fix-one-point-with-no-error-while-curve-fitting#answer_140081

编辑：Shashank Prasanna 2014-6-10

在 MATLAB Online 中打开

[Edited] removed least squares comment

You could give more importance to the first observation using weights, I do this using the fitlm function in the statistics toolbox:

x = 1:10; y = randn(10,1);
lm1 = fitlm(x,y);
lm2 = fitlm(x,y,'weights',[100, ones(1,length(y)-1)]);
scatter(x,y), hold on
plot(x, lm1.Fitted,'r')
plot(x, lm2.Fitted,'b')

9 个评论
显示 7更早的评论隐藏 7更早的评论

John D'Errico 2014-6-10

编辑：John D'Errico 2014-6-10

在 MATLAB Online 中打开

Sorry, but YES, you can do a least squares fit with a model with no constant term. Effectively, this is what was asked for. I think you misunderstand the idea of least squares. Least squares does NOT require a constant term to be estimated! You might as well say it is not least squares if there is no quadratic term in the model!

Anyway use of a high weight is a poor way to implement the solution. It causes numerical issues if the weight is too large. And even then, you won't get an exact solution that fits the requirements!

To forestall the inevitable argument that the problem as requested asked for a model that passes through a given point is not equivalent to a model with no constant term, a model that passes through the point (x0,y0) is best estimated as

y - y0 = m*(x - x0)

Where x0 and y0 are known constants. That is, a model with NO constant term.

You can estimate that model as:

m = (x(:) - x0))\(y(:) - y0);

This is the OLS estimate, in the MATLAB way.

Of course, lsqlin could also have been used.

John D'Errico 2015-1-14

If there are intrinsically nonlinear parameters to estimate, then no, lsqlin will not suffice. It solves only problems that are linear in the parameters. However, you CAN construct a model which will pass through a point exactly for a sinusoidal model, or have a variety of fixed properties. I assume your question relates to the one you posed recently on answers. I'm sorry, but that was a confusing question that I chose not to answer, partly because it was not at all clear what form your data was in, what form the model took on, etc. It looked like you had only a blurry image, and needed to extract a model of some form from that.

curoi 2015-1-14

Yeah, sorry about that. It does have to do with the question asked but I didn't realize the images were so blurry. I do just want to 'construct a model which will pass through a point exactly for a sinusoidal model' although it may be several points and multiple sine terms if that's possible.

The data is just in a form according to the indices of the image.

请先登录，再进行评论。

Answer 2

Matt J 2014-6-9

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/133385-how-to-fix-one-point-with-no-error-while-curve-fitting#answer_140093

在 MATLAB Online 中打开

It depends on the model function for your curve. Ideally, you would choose the formula for the curve to be such that it passes through the desired point. For example, the following linear curve function has unknown slope parameter m, but is written so that it always passes through the known point (x0,y0)

y=m*(x-x0)+y0

If the formula for your model is not tractable enough for this, you may need to use fmincon to do the curve-fitting and specify your constraint using the nonlcon argument.