How to apply least square approach in this case?
1 view (last 30 days)
I think I have a connected problem, it is given to me the experimental data (xdata,ydata), and I am asked to fit them to a model to find the parameters (a,b,c,d) with a model function as follows:
Y=f(xdata, A,B,C,D) with
A = fun1(xdata,a,b,c,d)
B = fun2(xdata,a,b,c,d)
C = fun3(xdata,a,b,c,d)
if xdata < d
D = fun4(xdata,a,b,c,d)
elseif xdaat> d
D = fun5(xdata,a,b,c,d)
with fun4 is totally different from fun5, so how to formulate the problem in this case please
Walter Roberson on 2 Feb 2016
Edited: Walter Roberson on 2 Feb 2016
The objective function for the minimizers can have if/else clauses. You would minimize over a, b, c, d, pass xdata in as an additional parameter, calculate your A, B, C, D, and use those to calculate f.
However, any time you have an if/else, or a min() or max(), you run the risk that you are introducing a discontinuity in the objective function or its first derivative, and either of those will typically lead to incorrect results for all of the minimizers for fmincon and all of the other minimizers in the Optimization Toolbox. As well, some of the minimizers for fmincon might have problems if the second derivative is not continuous.
If your objective or uts first derivative are not continuous then you need to use a function from the Global Optimization toolbox, or you need to split the inputs into subranges each of which has a continuous first derivative, run those separately, and take the best of the results.