English

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 O(nlog2(n))\mathcal{O}(n\log^2(n)). An error analysis guarantees the integrality gap shrinks to zero as nn\to\infty, and we apply the algorithm on a two-dimensional advection-diffusion equation, to determine the LIDAR's optimal sensing directions for data collection.

Keywords

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}
}
R2 v1 2026-06-23T12:45:28.069Z