Downloads & Free Reading Options - Results

We consider the problem of finding optimal parameters under simulation optimization setup. For a $p$-dimensional parameter optimization, the classical Kiefer-Wolfowitz Finite Difference Stochastic Approximation (FDSA) scheme uses $p+1$ or $2p$ simulations of the system feedback for one-sided and two-sided gradient estimates respectively. The dependence on the dimension $p$ makes FDSA impractical for high dimensional problems. An alternative approach for gradient estimation in high dimensional problems is the simultaneous perturbation technique that appears in [1],[2]. The main idea in this approach is to estimate the gradient by using only two settings of the $p$-dimensional parameter being optimized. The two settings of the parameter are obtained by simultaneously perturbing all the components of the parameter by adding a random direction. A drawback of using random directions for the gradient estimate is the very large or possibly infinite range of these random directions (for e.g. $\pm 1$ symmetric Bernoulli perturbations typically used in 1SPSA algorithm has a range of cardinality $2^p$ ). In this article we consider deterministic perturbations with a range of cardinality $p+1$ to improve the convergence of these algorithms. A novel construction of deterministic perturbations based on specially chosen circulant matrix is proposed. Convergence analysis of the proposed algorithms is presented along with numerical experiments.