Unique reconstruction for discretized inverse problems: a random sketching approach via subsampling

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Jin, Ruhui, Li, Qin, Nair, Anjali, Stechmann, Samuel
Format: Preprint
Veröffentlicht: 2024
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866917002193403904
author Jin, Ruhui
Li, Qin
Nair, Anjali
Stechmann, Samuel
author_facet Jin, Ruhui
Li, Qin
Nair, Anjali
Stechmann, Samuel
contents Theoretical inverse problems are often studied in an ideal infinite-dimensional setting. The well-posedness theory provides a unique reconstruction of the parameter function, when an infinite amount of data is given. Through the lens of PDE-constrained optimization, this means one attains the zero-loss property of the mismatch function in this setting. This is no longer true in computations when we are limited to finite amount of measurements due to experimental or economical reasons. Consequently, one must compromise the goal, from inferring a function, to a discrete approximation. What is the reconstruction power of a fixed number of data observations? How many parameters can one reconstruct? Here we describe a probabilistic approach, and spell out the interplay of the observation size $(r)$ and the number of parameters to be uniquely identified $(m)$. The technical pillar here is the random sketching strategy, in which the matrix concentration inequality and sampling theory are largely employed. By analyzing a randomly subsampled Hessian matrix, we attain a well-conditioned reconstruction problem with high probability. Our main theory is validated in numerical experiments, using an elliptic inverse problem as an example.
format Preprint
id arxiv_https___arxiv_org_abs_2403_05935
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Unique reconstruction for discretized inverse problems: a random sketching approach via subsampling
Jin, Ruhui
Li, Qin
Nair, Anjali
Stechmann, Samuel
Numerical Analysis
Optimization and Control
65M32, 49M41, 65F35
Theoretical inverse problems are often studied in an ideal infinite-dimensional setting. The well-posedness theory provides a unique reconstruction of the parameter function, when an infinite amount of data is given. Through the lens of PDE-constrained optimization, this means one attains the zero-loss property of the mismatch function in this setting. This is no longer true in computations when we are limited to finite amount of measurements due to experimental or economical reasons. Consequently, one must compromise the goal, from inferring a function, to a discrete approximation. What is the reconstruction power of a fixed number of data observations? How many parameters can one reconstruct? Here we describe a probabilistic approach, and spell out the interplay of the observation size $(r)$ and the number of parameters to be uniquely identified $(m)$. The technical pillar here is the random sketching strategy, in which the matrix concentration inequality and sampling theory are largely employed. By analyzing a randomly subsampled Hessian matrix, we attain a well-conditioned reconstruction problem with high probability. Our main theory is validated in numerical experiments, using an elliptic inverse problem as an example.
title Unique reconstruction for discretized inverse problems: a random sketching approach via subsampling
topic Numerical Analysis
Optimization and Control
65M32, 49M41, 65F35
url https://arxiv.org/abs/2403.05935