Page 1 of 1

The question of optimization function

Posted: Tue Apr 06, 2010 3:15 pm
by heather
hello

i am using vensim to do simulating with my thesis about inventory.
but i can't know how the optimization function calculated even use the trace function
could you tell me the calculating procedure of the Optimization function in Vensim
and what's the different between the modified Powell hill climbing algorithm and the original powell

Thanks for your help

Best Regards

Heather

Posted: Tue Apr 06, 2010 9:18 pm
by bob@vensim.com
Hi Heather,

The Powell search technique is technically not a hill climbing technique but a conjugate gradient technique. The book Numerical Recipes in C that shows up in the bibliography (Appendix B of the Reference Manual) has a reasonable discussion of this. The implementation in Vensim is somewhat modified from this, but largely the same.

Posted: Tue Apr 06, 2010 10:00 pm
by tomfid
A good start is to take a look at the online help; look in the index for "Payoff Computation"

In short, for calibration optimization, the payoff is the sum over time of weight*(model-data)^2. For policy optimization, it's just the sum of weight*modelvar. Kalman filter payoffs are slightly more complicated.

The idea behind the Powell algorithm is described fairly well in Numerical Recipes (see chapter 10 here http://www.fizyka.umk.pl/nrbook/bookcpdf.html ). I don't know what the modifications are in the Vensim code, but they probably have more to do with robustness to floating point errors and boundaries rather than anything fundamental about the approach.

Tom

Posted: Wed Apr 07, 2010 1:31 am
by tomfid
I always thought it was a direction set method. Doesn't conjugate gradient require derivatives (hence the name)?