Fitting data in a non-linear function with constraints

Question

Vogel 2018-7-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/411103-fitting-data-in-a-non-linear-function-with-constraints

编辑： Vogel 2018-7-25

Hello,

I need to fit some data (y_data = f(x1_data,x2_data)) in a certain function. I used 'objfcn', but I have to apply some boundary conditions (i.e. the curve must intersect a certain point, the second derivative must be positive, etc.).

Does anybody know if there is a way to do it with Matlab? (Well, I could artificially add the desired points to my data, but I guess this is not cientifically correct?)

I also checked lsqcurvefit, but, if I understood it well, you can only specify the ranges of the coefficients, but not the data.

In the following you find the case I need to solve, just in case my question was not clear.

Thank you in advance!

===================================================================================

I would like to fit some data in the following equation:

y_data = f(x1_data, x2_data)

Where f(x1, x2) = e^((a*x2^3 + b*x2^2 + c*x2 + d )*(-1/x1)) - z

and a,b,c,d and z are the coefficients to be determined.

There are many constraints:

1) 0 =< x2 =< 1

2) 0 =< y =< 1

(My data ranges from x2 = [0.3, 0.95] and y = [0.02,0.95])

3) x2 = 0 --> f = 0

(This is way I added the coefficient z. It is a kind of correction, since e^0 = 1, and there is otherwise no way to make it 0).

4) x2 = 1 --> f = 1

5) x2 = 1 --> f' = 1

6) x2 = 0 --> f'' > 0

To accomplish constraints 3, 4 and 5 I resolved the equations, so that dependent coefficients are obtained:

z = -e^(-d/x1)

c = -x1/(1+e^(-d/x1)) - 3*a -2*b

I tried to solve it as follows:

objfcn = @(b,x1,x2) exp(-1./x1.*( b(1).*x2.^3 + b(2).*x2.^2 + (-x1./(1+exp(b(3).*(-1./(x1))))-3*b(1)-2*b(2)).*x2 + b(3))) - exp(b(3).*(-1./(x1)));
[B,resnorm] = fminsearch(@(b) norm(y - objfcn(b,x1,x2)), [1;1;1]);

However, I don't know how to implement 1, 2 and 6.

*Here my questions:

- Is it possible to introduce constraints 3,4 and 5 with a Matlab function, instead of solving equations?

- Do you know if it is possible to apply constraints 1,2 and 6?*

Thanks!!!

13 个评论
显示 11更早的评论隐藏 11更早的评论

Vogel 2018-7-23

编辑：Vogel 2018-7-23

Mmmm my mathematical knowledge is not so good. I will try to explain shortly the physical relationship:

x1 = temperature

x2 = mass fraction

y = thermodynamic property, which is a function of temperature and mass fraction

x2 and y can only get values between 0 and 1. It is impossible that they have other values. This should occur theoretically for any temperature, so for any x1. If mathematically I need to constrain x1, then I could give my temperature work range: 273.15 - 400 K.

Attached you see the resulting curves for two temperatures and the whole mass fraction range evaluated in the fit function that I managed to generate by using the procedure described in my first post (and adding artificially as measurements points [x2=0, y=0] and [x2=1,y=1] for temperatures 273:10:400).

It is quiet good, the problem is that for low values of x2 I get negative values of y (very small though). Moreover, the fit underestimates slightly the value of 'y'.

This is why I am looking for another way to fit the data.

Thanks again

Vogel 2018-7-24

编辑：Vogel 2018-7-24

It is a mixture of a physical and mathematical model.

The physics say:

1) ln(y)=(-1/x1).*a.*B(x2), where 'a' is a constant and B(x) introduces the influence of mass fraction x2. This is approximately true for some temperature ranges, which can differ for every substance.

2) y = f(x1, x2=0)=0 and y = f(x1, x2=1)=1

3) 0<=y(x1,x2)<=1

4) y'= f'(x1, x2=1) = 1

5) At x2=0, the slope of the curve must be a stright line with a positive slope.

I need to find B(x2). However, there are not physical models to determine it, so here I just looked for mathematical functions that permits implementing all the boundary conditions. I tried many: a*x2^b, polynoms of second, third and fourth order, ....etc. The best one I found was this: B(x2) = b.*x2.^3 + c.*x2.^2 + d.*x2 + f

So that ln(y)=(-a/x1).*(b.*x2.^3 + c.*x2.^2 + d.*x2 + f).

Therefore: y = exp((-a/x1).*(b.*x2.^3 + c.*x2.^2 + d.*x2 + f))

It is observed that, due to the exponential, 'y' can never be eaxct zero. This is why I had to add a correction factor z, which needes to be a function of temperature and mass fraction. I just applied the same form as the principal function.

z = exp((-p/x1).*(q.*x2.^3 + r.*x2.^2 + s.*x2 + t))

Do you know any other function that could pass better? And, in any case, do you know how to apply constraints 2 to 5 to a fit function?

Thanks

Matt J 2018-7-24

在 MATLAB Online 中打开

Well (1) and (4) look impossible to satisfy simultaneously. I assume the f'() is the derivative with respect to x2, not x1. If so, then using (1) we can compute the form of f'() explicitly

y'(x1,x2=1) = (a/x1)^2*B(1)*B'(1)

There is no way the right hand side can equal 1 independently of x1, as demanded by constraint (4).

Vogel 2018-7-25

编辑：Vogel 2018-7-25

Yes, f' is the derivative respect x2 and I calculated it my self and applied the constraint. However, for complicated functions, it is an arduous task. I was wondering if there was a way for Matlab to do it during the fitting process.

请先登录，再进行评论。

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

Answer 1

Matt J 2018-7-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/411103-fitting-data-in-a-non-linear-function-with-constraints#answer_330164

编辑：Matt J 2018-7-24

I need to find B(x2). However, there are not physical models to determine it, so here I just looked for mathematical functions that permits implementing all the boundary conditions.

If you aren't committed to a particular model, it might be useful for you to look at SLM. It allows you to do spline fitting subject to various boundary conditions and montononicity conditions both on the function and on slope.