Fit data with dependent parameters

2018 7 13
2 个回答
回答已采纳
更新时间：2018 7 25
7 次查看（30 天）

0 个投票

Hi,
There are two rows of data, x and y. I would like to fit y = f(x), where
f(x) = a*x^3 + b*x^2 + (2a+3b)*x,
i.e. parameters are not independent.
I tried to use the function "fittype", but it does not work (Licensing error: -101,147).
I would like to know if there is any other way to solve it.
Thank you!

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

 采纳的回答

Star Strider
Star Strider 2018-7-13

1 个投票

Yours is a linear problem, however the easiest way to estimate the parameters is likely an unconstrained nonlinear solver, such as fminsearch:
x = ...;
y = ...;
objfcn = @(b,x) b(1).*x.^3 + b(2).*x.^2 + (2*b(1) + 3*b(2)).*x;
[B,resnorm] = fminsearch(@(b) norm(y - objfcn(b,x)), [1;1]);
xv = linspace(min(x), max(x));
figure
plot(x, y, 'pb')
hold on
plot(xv, objfcn(B,xv), '-r')
hold off
grid

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

Vogel
Vogel 2018-7-16
Thanks Star Strider. That was a great solution.
Star Strider
Star Strider 2018-7-16
As always, my pleasure!
Vogel
Vogel 2018-7-24
Hi Star Strider,
I have realized that the change of the values of the initial points [1;1] exhibits a considerable influence on the result of the fitted parameters and the function. Do you know how are they usually chosen? By manipulating them, I can get the a more convenient form for the desired curve. But how far is this cheating?
Thanks again
Star Strider
Star Strider 2018-7-24
It is not ‘cheating’. The ‘surface’ (or ‘hypersurface’) of the parameter space can have a number of local minima as well as the global minimum that results in the optimal (lowest error) parameter set in the nonlinear objective function. Nonlinear parameter estimation can be extremely sensitive to the correct choice of initial parameters, because of the presence of local minima.
I usually deal with this by using a genetic algorithm (the Global Optimization Toolbox ga function). With the proper (totally random) and large enough (I usually begin with 500 individuals) initial population, ga has always (in my experience) found the best parameter set. It can take a while for the algorithm to converge, especially if constrained, and for difficult problems. The form of the objective function to use with ga is the same as that used with fminsearch in this example, specifically:
@(b) norm(y - objfcn(b,x))
I structure my code to display the parameter results in the Command Window, so I have a record of them.
Vogel
Vogel 2018-7-25
Thanks a lot for your clear explanation. I am going to check it.
Star Strider
Star Strider 2018-7-25
As always, my pleasure.

请先登录，再进行评论。

更多回答（1 个）

Vogel
Vogel 2018-7-17

0 个投票

Hi again,
I have new questions, which probably have more to do with mathematics. But maybe Matlab can help.
Do you know if it is possible to add some constrains?
- I would like to obtain a fitted curve that intersects the points [x=0,y=0] and [x=1, y=1]. I tried to do it by applying some boundary conditions (this is why the coefficients are interdependent). However, the resulting curve does not intersect exactly the points, it just approximates them. (Actually the equation is much more complex than the posted above. It has exponential functions).
- This is more complicated: Is it possible to impose that the second derivative of the resulting curve must be positive at a certain point x (in my case at x=0)?
Thank you in advance

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

Star Strider
Star Strider 2018-7-17
I have no idea how to do either of those, especially for a nonlinear function.
Vogel
Vogel 2018-7-18
Thank you anyway :)
Matt J
Matt J 2018-7-18
You could probably do it with lsqlin, but better you post this as a new question, detailing your actual model function.
Vogel
Vogel 2018-7-18
Thank you Matt J. I will try it.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Linear Least Squares 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by