Inversion of the Laplace Transform of Point Masses
<|reference_start|>Inversion of the Laplace Transform of Point Masses: 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.<|reference_end|>
@article{mckenna2024inversion,
    title={Inversion of the Laplace Transform of Point Masses},
    author={Michael McKenna, Hrushikesh N. Mhaskar, and Richard G. Spencer},
    journal={arXiv preprint arXiv:2402.04348},
    year={2024},
    archivePrefix={arXiv},
    eprint={2402.04348},
    primaryClass={math.NA cs.NA}
}