Inversion of the Laplace Transform of Point Masses

Michael McKenna,Hrushikesh N. Mhaskar,Richard G. Spencer
2024-02-07
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.
Numerical Analysis
What problem does this paper attempt to address?