Hi Jun Yang,
To the best of my understanding you wish to solve a non-linear optimization problem with “SQP” algorithm offered by “fmincon” function, however in some runs of the function the execution apparently freezes this may happen because of the followings -
- At each iteration, an approximation is made of the Hessian of the Lagrangian function using a quasi-Newton updating method. This is then used to generate a QP subproblem whose solution is used to form a search direction for a line search procedure. During its iterations, the SQP algorithm can attempt to take a step that fails. This means an objective function or nonlinear constraint function you supply returns a value of Inf, NaN, or a complex value.
- If the nonlinear constraints are not satisfied, and an attempted step causes the constraint violation to grow. The SQP algorithm attempts to obtain feasibility using a second-order approximation to the constraints. This second-order technique can lead to a feasible solution. However, this technique can slow the solution by requiring more evaluations of the nonlinear constraint functions, more large the constraints the more slow the iteration.
The SQP algorithm involves solving a sequence of quadratic subproblems, which can be sensitive to numerical issues such as round-off errors, ill-conditioning, or numerical instability In such cases, it is advisable to try alternative optimization algorithms, adjust the algorithm parameters, or reformulate the problem to improve convergence.
Refer to the following documentation for a more detailed description of “fmincon sqp” -https://in.mathworks.com/help/optim/ug/constrained-nonlinear-optimization-algorithms.html#brnox01
Regards,
Vinayak Luha
