Solving Optimal Experimental Design with Sequential Quadratic Programming and Chebyshev Interpolation
Computation
2019-12-30 v2
Abstract
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 . An error analysis guarantees the integrality gap shrinks to zero as , and we apply the algorithm on a two-dimensional advection-diffusion equation, to determine the LIDAR's optimal sensing directions for data collection.
Cite
@article{arxiv.1912.06622,
title = {Solving Optimal Experimental Design with Sequential Quadratic Programming and Chebyshev Interpolation},
author = {Jing Yu and Mihai Anitescu},
journal= {arXiv preprint arXiv:1912.06622},
year = {2019}
}