Is there any way to exclude the complex values ​​of the iterations with lsqnonlin ?

6 次查看(过去 30 天)
I'm using lsqnonlin to solve a non linear system, and my solution must be real, so I wanted to know if I could exclude complex iterations so as to make more efficient the search.
Thanks in advance!

采纳的回答

Matt J
Matt J 2014-9-16
编辑:Matt J 2014-9-16
Just to be clear, the complex numbers are being introduced by the objective function that you have provided, not code elsewhere in lsqnonlin.
Make sure that your objective is not executing expressions like log(z) for negative z and similar things like that. If you have regions like that where the objective function is undefined, you should set the value there to NaN or Inf and you should probably use fmincon to do the optimization instead, with the sqp algorithm. The sqp algorithm has the capability of course-correcting when it strays into regions of NaNs or Infs in the objective or constraints.
  4 个评论
Jurgen
Jurgen 2014-11-11
编辑:Matt J 2014-11-12
In my case I have a system of equations (without objective function), it enough to put zero in the part of the objective function?
fmincon(0,x0,[],[],[],[],lb,ub,nonlcon,options)
Matt J
Matt J 2014-11-12
编辑:Matt J 2014-11-12
I don't know why you would do that. If your equations have an exact solution, then using a least squares objective function should find it. However, you can also replace bad values returned by nonlcon with nans or infs in a similar manner.

请先登录,再进行评论。

更多回答(1 个)

Jurgen
Jurgen 2014-11-15
编辑:Jurgen 2014-11-15
I'm a little confused now. Just to be clear, I need to solve a system containing non-linear equations and inequalities. Considering this I have two questions:
1) What is the best way to solve my system ?. I have tried lsqnonlin, Newton-Raphson and now I'm trying with fmincon.
2) Newton-Raphson and lsqnonlin solve it, but gives me complex numbers, which is useless for me.
Thanks a lot Matt!
  2 个评论
Matt J
Matt J 2014-11-15
编辑:Matt J 2014-11-15
I don't see how you would have tried any of those. lsqnonlin doesn't handle nonlinear inequalities while Newton Raphson doesn't support inequality constraints of any kind.
Any of the fmincon formulations, using the sqp algorithm, that we were discussing should work, so long as a solution does indeed, exist. If an exact solution doesn't exist, you will have to compromise and accept a least squares solution. In that case, you need to decide which of your equations/inequalities need to be exactly satisfied and which you can relax.
Jurgen
Jurgen 2014-11-15
Forget explain that. I use a gap, for example if I wanted X> Y, what I did was put Y-X + H = 0 with H greater than or equal to zero.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Solver-Based Nonlinear Optimization 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by