Re: FindMinimum constraints
- To: mathgroup at smc.vnet.net
- Subject: [mg92986] Re: FindMinimum constraints
- From: dantimatter <google at dantimatter.com>
- Date: Wed, 22 Oct 2008 05:37:33 -0400 (EDT)
- References: <gdja4r$idt$1@smc.vnet.net> <gdjq6u$mtu$1@smc.vnet.net>
Thanks for the quick response. Join[] did indeed do the trick. But adding this constraint introduced another problem. By fixing a[100] == 0.5*(a[0] + a[25], I (inadvertently) allowed FindMinimum to choose solutions that seem to have discontinuous derivatives (i.e., there are really sharp points at a[100], a[0], and a[25]). I guess I should also somehow force the solution to be sufficiently 'smooth'. Is there an easy way to do that? Thanks again for all the help! Dan On Oct 21, 1:43 am, Jean-Marc Gulliet <jeanmarc.gull... at gmail.com> wrote: > dantimatter wrote: > > I have a large and complicated model for which I'm using the > > FindMinimum function to find a solution. The solution consists of 10= 1 > > numbers a[i] (i goes from 0 to 100), all greater than zero, and the > > search for each a[i] begins at 0.001: > > > sol = FindMinimum[{model, Table[a[i] > 0,{i, 0, 100}]}, Table[{a[i], > > 0.001},{i, 0, 100}]]; > > > The line of code that I'm showing here works great. My problem now i= s > > that I'd like to add an additional constraint: that a[100] == > > 0.5*(a[0]+a[25]). I can't for the life of me figure out how to > > include this added constraint in the line of code that already works. > > I get errors when I use Append[] to the Table. I would appreciate any > > thoughts on this. > > Using *Join* should do the trick. For instance, > > sol = FindMinimum[{model, > Join[{a[100] == 0.5*(a[0] + a[25])}, > Table[a[i] > 0, {i, 0, 100}]]}, Table[{a[i], 0.001}, {i, 0, 10= 0}]] > > Regards, > -- Jean-Marc