Solving Optimal Experimental Design with Sequential Quadratic Programming and Chebyshev Interpolation

12/13/2019
by   Jing Yu, et al.
0

We propose an optimization algorithm to compute the optimal sensor locations in experimental design in the formulation of Bayesian inverse problems, where the parameter-to-observable mapping is described through an integral equation and its discretization results in a continuously indexed matrix whose size depends on the mesh size n. By approximating the gradient and Hessian of the objective design criterion from Chebyshev interpolation, we solve a sequence of quadratic programs and achieve the complexity O(nlog^2(n)). An error analysis guarantees the integrality gap shrinks to zero as n→∞, and we apply the algorithm on a two-dimensional advection-diffusion equation, to determine the LIDAR's optimal sensing directions for data collection.

READ FULL TEXT

Please sign up or login with your details

Forgot password? Click here to reset