Title:Inversion of the Laplace Transform of Point Masses

Abstract:Motivated by applications in magnetic resonance relaxometry, we consider the following problem: Given samples of a function $t\mapsto \sum_{k=1}^K A_k\exp(-t\lambda_k)$, where $K\ge 2$ is an integer, $A_k\in\mathbb{R}$, $\lambda_k>0$ for $k=1,\cdots, K$, determine $K$, $A_k$'s and $\lambda_k$'s. Our approach is to transform this function into another function of the same form where $\lambda_k$'s are replaced by $i\lambda_k$. For this purpose, we study the least square approximation using polynomials weighted by the Gaussian weight, and use the fact that Hermite functions are eigenfunctions of the Fourier transform. We provide a detailed analysis of the effect of noise in the data.