multivariable curve fitting with constraints

Question

michael 2016-10-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/308726-multivariable-curve-fitting-with-constraints

编辑： Matt J 2016-10-26

Hello,

I'm confused from all the least squares functions which is supported by Matlab: I have a function which is the basis for curve fitting f(x,y,z), some test input data:

xdata=[x,y,z];

and corresponding output data

ydata=y;

I'd like to perform curve fitting of f(x,y,z) to the test.

Some constraints that I would like to use:

x_a < x < x_b
y_a < y < y_b
z_a < z < z_b
f(x,y,z)>=0 for each x,y,z in the domain

Please assist with function which I should use and the corresponding input parameters for it.

Other option: based on the xdata, ydata and constraints, I'd like to get an analytic function f(x,y,z) which best fits the data and constraints.

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

请先登录，再进行评论。

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

Answer 1

Matt J 2016-10-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/308726-multivariable-curve-fitting-with-constraints#answer_240286

编辑：Matt J 2016-10-23

在 MATLAB Online 中打开

LSQCURVEFIT would be one choice,

x = lsqcurvefit(func,x0,xdata,ydata,lb,ub)

As you can see in the documentation, it readily allows you to specify upper and lower bounds lb and ub on the parameters, x. You have to decide on the parametric model and supply a function handle (func) to it of the form F(x,xdata). You would also have to provide an initial parameter guess x0 for the iterative estimation.

The positivity requirement on the function f(x,y,z)>0 will probably require you to choose a parametric family which is inherently positive, e.g., a Gaussian lobe mixture.

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

michael 2016-10-24

x0 is the initial assumption on unknown values of the function. In case of f(x)=ax+b, the unknowns are a,b and the bounds are on the value of x, not on the xdata!

Matt J 2016-10-26

编辑：Matt J 2016-10-26

在 MATLAB Online 中打开

The bounds that you gave in your posted example are on f(x), namely f(x)>=0. x is an independent variable controlled by you.

Constraining f(x)>=0 over a continuous region of x is not something MATLAB has tools for for multi-dimensional x. For scalar x, you could use fsemininf().

The simplest way to ensure f(x)>=0 is to choose an f(x) that can never be negative for any choice of the unknowns. As I mentioned, an inherently positive parametric family like a mixture of Gaussians would do this.

If you're saying that you're stuck with a more complicated family, then the best solution will likely be specific to that family. For a linear family f(x)=ax+b, you would take advantage that f reaches its extremes at the vertices of the region lb < x < ub. Constraining a*x+b>=0 is then just a simple matter of imposing the linear constraints,