How to apply least square approach in this case?

2 views (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)
end
with fun4 is totally different from fun5, so how to formulate the problem in this case please
  2 Comments
Walter Roberson
Walter Roberson on 2 Feb 2016
In some of those cases your f has 5 parameters, in some it has 4 parameters, in some it has 3 parameters, and in some it has 1 parameter. Those cannot all be the same function.
Zine
Zine on 2 Feb 2016
Edited: Zine on 2 Feb 2016
Thanks Walter for your passing, you are right, I edited the question to be more clear, the idea is one of the parameters is a criteria to choose the objective function

Sign in to comment.

Accepted Answer

Walter Roberson
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.

More Answers (0)

Categories

Find more on Problem-Based Optimization Setup in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!