Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.math.mathematica

Topic: Some linear constraints seem to be ignored in function NMinimize with
Replies: 0  

Advanced Search

Back to Topic List Back to Topic List  
Daniel Stariolo

Posts: 2
Registered: 2/22/14
Some linear constraints seem to be ignored in function NMinimize with
Posted: Feb 22, 2014 3:01 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

I'm trying to minimize a non-linear function of four variables with some linear constraints. Mathematica 8 is unable to find a good solution giving complex values of the function at some point in the iteration. This implies that one or some contraints are not being enabled in the process. Is this a bug or limitation of the optimization function ?

Function to minimize is

ff[lxw_, lwz_, c_, d_] := - J1 (lxw + lwz) - 2 J2 c +
T (-Log[2] - 1/2 (1 - lxw) Log[(1 - lxw)/4] -
1/2 (1 + lxw) Log[(1 + lxw)/4] -
1/2 (1 - lwz) Log[(1 - lwz)/4] -
1/2 (1 + lwz) Log[(1 + lwz)/4] + 1/2 (1 - d) Log[(1 - d)/16] +
1/8 (1 + 2 c + d - 2 lwz - 2 lxw) Log[
1/16 (1 + 2 c + d - 2 lwz - 2 lxw)])

where

T = 10;
J1 = 1;
J2 = -0.2;

are constant parameters. Then I try

NMinimize[{ff[lxw, lwz, c, d],
2 c + d - 2 lwz - 2 lxw >= -0.999 && -0.999 <= lxw <=
0.999 && -0.999 <= lwz <= 0.999 && -0.999 <= c <= 0.999 &&
d <= 0.9999}, {lxw, lwz, c, d}]

with the result

NMinimize::nrnum: "The function value 5.87777-4.87764 I is not a real number at {c,d,lwz,lxw} = {-0.718817,-1.28595,0.69171,-0.932461}."

I would appreciate if someone can give a hint at what is happening here.




Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.