Fsolve produces comlex numbers

Question

Bismark 2023-5-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1975284-fsolve-produces-comlex-numbers

评论： Walter Roberson 2023-5-30

When i run my code, the result are complex numbers. Is there anyway i can locate the cause or any suggestions on how to obtian non-complex numbers will be helpful.

2 个评论
显示无隐藏无

Rik 2023-5-29

Have a read here and here. It will greatly improve your chances of getting an answer.

Bismark 2023-5-29

Okay Rik. Thank you

请先登录，再进行评论。

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

Answer 1

Torsten 2023-5-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1975284-fsolve-produces-comlex-numbers#answer_1246674

编辑：Torsten 2023-5-29

The usual reasons why solvers like "fsolve" produce complex numbers are:

you take the square-root of a negative number (sqrt(x) with x negative)
you exponentiate a negative number (x^a with x negative)
you take the logarithm of a negative number (log(x) with x negative)
or similar operations with functions that give complex results ...

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

John D'Errico 2023-5-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1975284-fsolve-produces-comlex-numbers#answer_1246684

在 MATLAB Online 中打开

Remember that fsolve does not allow constraints. Not even bound constraints. So when one of the variables tries to go into a bad place, where, perhaps you are taking the sqrt of a negative number, then complex results start getting created. That commonly happens when a log, a sqrt, a fractional power are applied to a negative number. But other operations surely cause problems of this sort too.

I would strongly suggest you learn to use the debugger. Use it to watch your objective function, and see when something strange happens.

You might decide to set breakpoints in your function. Or, you might just decide to dump the current set of parameters at each iteration to the command window, as well as the objective. Watch for when a problem arises, and you will see that what we have said is surely true.

help debug
 DEBUG List MATLAB debugging functions
 
    dbstop     - Set breakpoint.
    dbclear    - Remove breakpoint.
    dbcont     - Resume execution.
    dbdown     - Change local workspace context.
    dbmex      - Enable MEX-file debugging.
    dbstack    - List who called whom.
    dbstatus   - List all breakpoints.
    dbstep     - Execute one or more lines.
    dbtype     - List file with line numbers.
    dbup       - Change local workspace context.
    dbquit     - Quit debug mode.
 
    When a breakpoint is hit, MATLAB goes into debug mode, the debugger
    window becomes active, and the prompt changes to a K>>.  Any MATLAB
    command is allowed at the prompt.  
 
    To resume program execution, use DBCONT or DBSTEP.  
    To exit from the debugger use DBQUIT.

    Documentation for debug
       doc debug

Can you stop fsolve from going to the bad places in your parameter space? Sadly, no. You can use other optimizers. But even there, you need to know where it is trying to go, and what constraints you need to apply. You can also decide to transform the parameters, to avoid the problem from happening. Of course, since we don't see any code, it becomes impossible to offer any useful help.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Walter Roberson 2023-5-30

在 MATLAB Online 中打开

Debugging tip:

On the first executable line after you calculate the value to be returned, put in a conditional breakpoint. The command line form of it would be like

dbstop at LINENUMBER if any(imag(VARIABLE)~=0)

or you can right-click when pointing to the dot before the line in the editor and insert a conditional breakpoint with expression any(imag(VARIABLE)~=0)

This requires that there is at least one executable line in the function before the function is finished. The easiest way to arrange that is to use an end statement matching the function header. (Note that all functions in the same file must be the same style -- either they all have end matching function or else none of them do.)

请先登录，再进行评论。