BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Thorsten Hohage

Thorsten Hohage appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath.APComputationmath.FAastro-ph.SRmath-phmath.DGmath.MPmath.OCmath.STMethodologyphysics.comp-phphysics.med-phphysics.opticsq-fin.STStatistics Theory

Catalog footprint

What is connected

22works
17topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Convergence Rates for Oversmoothing Banach Space Regularization

This paper studies Tikhonov regularization for finitely smoothing operators in Banach spaces when the penalization enforces too much smoothness in the sense that the penalty term is not finite at the true solution. In a Hilbert space setting, Natterer (1984) showed with the help of spectral theory that optimal rates can be achieved in this situation. ('Oversmoothing does not harm.') For oversmoothing variational regularization in Banach spaces only very recently progress has been achieved in several papers on different settings, all of which construct families of smooth approximations to the true solution. In this paper we propose to construct such a family of smooth approximations based on $K$-interpolation theory. We demonstrate that this leads to simple, self-contained proofs and to rather general results. In particular, we obtain optimal convergence rates for bounded variation regularization, general Besov penalty terms and $\ell^p$ wavelet penalization with $p<1$ which cannot be treated by previous approaches. We also derive minimax optimal rates for white noise models. Our theoretical results are confirmed in numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A spectrally accurate method for the dielectric obstacle scattering problem and applications to the inverse problem

We analyze the inverse problem to reconstruct the shape of a three dimensional homogeneous dielectric obstacle from the knowledge of noisy far field data. The forward problem is solved by a system of second kind boundary integral equations. For the numerical solution of these coupled integral equations we propose a fast spectral algorithm by transporting these equations onto the unit sphere. We review the differentiability properties of the boundary to far field operator and give a characterization of the adjoint operator of the first Fréchet derivative. Using these results we discuss the implementation of the iteratively regularized Gauss-Newton method for the numerical solution of the inverse problem and give numerical results for star-shaped obstacles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the well-posedness of the damped time-harmonic Galbrun equation and the equations of stellar oscillations

We study the time-harmonic Galbrun equation describing the propagation of sound in the presence of a steady background flow. With additional rotational and gravitational terms these equations are also fundamental in helio- and asteroseismology as a model for stellar oscillations. For a simple damping model we prove well-posedness of these equations, i.e. uniqueness, existence, and stability of solutions under mild conditions on the parameters (essentially subsonic flows). The main tool of our analysis is a generalized Helmholtz decomposition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sweeping preconditioners for stratified media in the presence of reflections

In this paper we consider sweeping preconditioners for time harmonic wave propagation in stratified media, especially in the presence of reflections. In the most famous class of sweeping preconditioners Dirichlet-to-Neumann operators for half-space problems are approximated through absorbing boundary conditions. In the presence of reflections absorbing boundary conditions are not accurate resulting in an unsatisfactory performance of these sweeping preconditioners. We explore the potential of using more accurate Dirichlet-to-Neumann operators within the sweep. To this end, we make use of the separability of the equation for the background model. While this improves the accuracy of the Dirichlet-to-Neumann operator, we find both from numerical tests and analytical arguments that it is very sensitive to perturbations in the presence of reflections. This implies that even if accurate approximations to Dirichlet-to-Neumann operators can be devised for a stratified medium, sweeping preconditioners are limited to very small perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Uniqueness of an inverse source problem in experimental aeroacoustics

This paper is concerned with the mathematical analysis of experimental methods for the estimation of the power of an uncorrelated, extended aeroacoustic source from measurements of correlations of pressure fluctuations. We formulate a continuous, infinite dimensional model describing these experimental techniques based on the convected Helmholtz equation in $\mathbb{R}^3$ or $\mathbb{R}^2$. As a main result we prove that an unknown, compactly supported source power function is uniquely determined by idealized, noise-free correlation measurements. Our framework further allows for a precise characterization of state-of-the-art source reconstruction methods and their interrelations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Density Matrix Reconstructions in Ultrafast Transmission Electron Microscopy: Uniqueness, Stability, and Convergence Rates

In the recent paper [17] the first experimental determination of the density matrix of a free electron beam has been reported. The employed method leads to a linear inverse problem with a positive semidefinite operator as unknown. The purpose of this paper is to complement the experimental and algorithmic results in the work mentioned above by a mathematical analysis of the inverse problem concerning uniqueness, stability, and rates of convergence under different types of a-priori information.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Elastic energy regularization for inverse obstacle scattering problems

By introducing a shape manifold as a solution set to solve inverse obstacle scattering problems we allow the reconstruction of general, not necessarily star-shaped curves. The bending energy is used as a stabilizing term in Tikhonov regularization to gain independence of the parametrization. Moreover, we discuss how self-intersections can be avoided by penalization with the Möbius energy and prove the regularizing property of our approach as well as convergence rates under variational source conditions. In the second part of the paper the discrete setting is introduced, and we describe a numerical method for finding the minimizer of the Tikhonov functional on a shape-manifold. Numerical examples demonstrate the feasibility of reconstructing non-star-shaped obstacles.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Regularized Newton Methods for X-ray Phase Contrast and General Imaging Problems

Like many other advanced imaging methods, x-ray phase contrast imaging and tomography require mathematical inversion of the observed data to obtain real-space information. While an accurate forward model describing the generally nonlinear image formation from a given object to the observations is often available, explicit inversion formulas are typically not known. Moreover, the measured data might be insufficient for stable image reconstruction, in which case it has to be complemented by suitable a priori information. In this work, regularized Newton methods are presented as a general framework for the solution of such ill-posed nonlinear imaging problems. For a proof of principle, the approach is applied to x-ray phase contrast imaging in the near-field propagation regime. Simultaneous recovery of the phase- and amplitude from a single near-field diffraction pattern without homogeneity constraints is demonstrated for the first time. The presented methods further permit all-at-once phase contrast tomography, i.e. simultaneous phase retrieval and tomographic inversion. We demonstrate the potential of this approach by three-dimensional imaging of a colloidal crystal at 95 nm isotropic resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convergence Rates for Exponentially Ill-Posed Inverse Problems with Impulsive Noise

This paper is concerned with exponentially ill-posed operator equations with additive impulsive noise on the right hand side, i.e. the noise is large on a small part of the domain and small or zero outside. It is well known that Tikhonov regularization with an $L^1$ data fidelity term outperforms Tikhonov regularization with an $L^2$ fidelity term in this case. This effect has recently been explained and quantified for the case of finitely smoothing operators. Here we extend this analysis to the case of infinitely smoothing forward operators under standard Sobolev smoothness assumptions on the solution, i.e. exponentially ill-posed inverse problems. It turns out that high order polynomial rates of convergence in the size of the support of large noise can be achieved rather than the poor logarithmic convergence rates typical for exponentially ill-posed problems. The main tools of our analysis are Banach spaces of analytic functions and interpolation-type inequalities for such spaces. We discuss two examples, the (periodic) backwards heat equation and an inverse problem in gradiometry.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Verification of a variational source condition for acoustic inverse medium scattering problems

This paper is concerned with the classical inverse scattering problem to recover the refractive index of a medium given near or far field measurements of scattered time-harmonic acoustic waves. It contains the first rigorous proof of (logarithmic) rates of convergence for Tikhonov regularization under Sobolev smoothness assumptions for the refractive index. This is achieved by combining two lines of research, conditional stability estimates via geometrical optics solutions and variational regularization theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Generalization of the Chambolle-Pock Algorithm to Banach Spaces with Applications to Inverse Problems

For a Hilbert space setting Chambolle and Pock introduced an attractive first-order algorithm which solves a convex optimization problem and its Fenchel dual simultaneously. We present a generalization of this algorithm to Banach spaces. Moreover, under certain conditions we prove strong convergence as well as convergence rates. Due to the generalization the method becomes efficiently applicable for a wider class of problems. This fact makes it particularly interesting for solving ill-posed inverse problems on Banach spaces by Tikhonov regularization or the iteratively regularized Newton-type method, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence of infinite element methods for scalar waveguide problems

We consider the numerical solution of scalar wave equations in domains which are the union of a bounded domain and a finite number of infinite cylindrical waveguides. The aim of this paper is to provide a new convergence analysis of both the Perfectly Matched Layer (PML) method and the Hardy space infinite element method in a unified framework. We treat both diffraction and resonance problems. The theoretical error bounds are compared with errors in numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Convergence Rates for Inverse Problems with Impulsive Noise

We study inverse problems F(f) = g with perturbed right hand side g^{obs} corrupted by so-called impulsive noise, i.e. noise which is concentrated on a small subset of the domain of definition of g. It is well known that Tikhonov-type regularization with an L^1 data fidelity term yields significantly more accurate results than Tikhonov regularization with classical L^2 data fidelity terms for this type of noise. The purpose of this paper is to provide a convergence analysis explaining this remarkable difference in accuracy. Our error estimates significantly improve previous error estimates for Tikhonov regularization with L^1-fidelity term in the case of impulsive noise. We present numerical results which are in good agreement with the predictions of our analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Generalization of the noise model for time-distance helioseismology

In time-distance helioseismology, information about the solar interior is encoded in measurements of travel times between pairs of points on the solar surface. Travel times are deduced from the cross-covariance of the random wave field. Here we consider travel times and also products of travel times as observables. They contain information about e.g. the statistical properties of convection in the Sun. The basic assumption of the model is that noise is the result of the stochastic excitation of solar waves, a random process which is stationary and Gaussian. We generalize the existing noise model (Gizon and Birch 2004) by dropping the assumption of horizontal spatial homogeneity. Using a recurrence relation, we calculate the noise covariance matrices for the moments of order 4, 6, and 8 of the observed wave field, for the moments of order 2, 3 and 4 of the cross-covariance, and for the moments of order 2, 3 and 4 of the travel times. All noise covariance matrices depend only on the expectation value of the cross-covariance of the observed wave field. For products of travel times, the noise covariance matrix consists of three terms proportional to $1/T$, $1/T^2$, and $1/T^3$, where $T$ is the duration of the observations. For typical observation times of a few hours, the term proportional to $1/T^2$ dominates and $Cov[τ_1 τ_2, τ_3 τ_4] \approx Cov[τ_1, τ_3] Cov[τ_2, τ_4] + Cov[τ_1, τ_4] Cov[τ_2, τ_3]$, where the $τ_i$ are arbitrary travel times. This result is confirmed for $p_1$ travel times by Monte Carlo simulations and comparisons with SDO/HMI observations. General and accurate formulae have been derived to model the noise covariance matrix of helioseismic travel times and products of travel times. These results could easily be generalized to other methods of local helioseismology, such as helioseismic holography and ring diagram analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On parameter identification in stochastic differential equations by penalized maximum likelihood

In this paper we present nonparametric estimators for coefficients in stochastic differential equation if the data are described by independent, identically distributed random variables. The problem is formulated as a nonlinear ill-posed operator equation with a deterministic forward operator described by the Fokker-Planck equation. We derive convergence rates of the risk for penalized maximum likelihood estimators with convex penalty terms and for Newton-type methods. The assumptions of our general convergence results are verified for estimation of the drift coefficient. The advantages of log-likelihood compared to quadratic data fidelity terms are demonstrated in Monte-Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Iterative Estimation of Solutions to Noisy Nonlinear Operator Equations in Nonparametric Instrumental Regression

This paper discusses the solution of nonlinear integral equations with noisy integral kernels as they appear in nonparametric instrumental regression. We propose a regularized Newton-type iteration and establish convergence and convergence rate results. A particular emphasis is on instrumental regression models where the usual conditional mean assumption is replaced by a stronger independence assumption. We demonstrate for the case of a binary instrument that our approach allows the correct estimation of regression functions which are not identifiable with the standard model. This is illustrated in computed examples with simulated data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence rates in expectation for Tikhonov-type regularization of Inverse Problems with Poisson data

In this paper we study a Tikhonov-type method for ill-posed nonlinear operator equations $\gdag = F(\udag)$ where $\gdag$ is an integrable, non-negative function. We assume that data are drawn from a Poisson process with density $t\gdag$ where $t>0$ may be interpreted as an exposure time. Such problems occur in many photonic imaging applications including positron emission tomography, confocal fluorescence microscopy, astronomic observations, and phase retrieval problems in optics. Our approach uses a Kullback-Leibler-type data fidelity functional and allows for general convex penalty terms. We prove convergence rates of the expectation of the reconstruction error under a variational source condition as $t\to\infty$ both for an a priori and for a Lepski{\uı}-type parameter choice rule.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Iteratively regularized Newton-type methods for general data misfit functionals and applications to Poisson data

We study Newton type methods for inverse problems described by nonlinear operator equations $F(u)=g$ in Banach spaces where the Newton equations $F'(u_n;u_{n+1}-u_n) = g-F(u_n)$ are regularized variationally using a general data misfit functional and a convex regularization term. This generalizes the well-known iteratively regularized Gauss-Newton method (IRGNM). We prove convergence and convergence rates as the noise level tends to 0 both for an a priori stopping rule and for a Lepski{\uı}-type a posteriori stopping rule. Our analysis includes previous order optimal convergence rate results for the IRGNM as special cases. The main focus of this paper is on inverse problems with Poisson data where the natural data misfit functional is given by the Kullback-Leibler divergence. Two examples of such problems are discussed in detail: an inverse obstacle scattering problem with amplitude data of the far-field pattern and a phase retrieval problem. The performence of the proposed method for these problems is illustrated in numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Riesz bases and Jordan form of the translation operator in semi-infinite periodic waveguides

We study the propagation of time-harmonic acoustic or transverse magnetic (TM) polarized electromagnetic waves in a periodic waveguide lying in the semi-strip $(0,\infty)\times(0,L)$. It is shown that there exists a Riesz basis of the space of solutions to the time-harmonic wave equation such that the translation operator shifting a function by one periodicity length to the left is represented by an infinite Jordan matrix which contains at most a finite number of Jordan blocks of size $> 1$. Moreover, the Dirichlet-, Neumann- and mixed traces of this Riesz basis on the left boundary also form a Riesz basis. Both the cases of frequencies in a band gap and frequencies in the spectrum and a variety of boundary conditions on the top and bottom are considered.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High order Curl-conforming Hardy space infinite elements for exterior Maxwell problems

A construction of prismatic Hardy space infinite elements to discretize wave equations on unbounded domains $Ω$ in $H^1_{loc}(Ω)$, $H_{loc}(curl;Ω)$ and $H_{loc}(div;Ω)$ is presented. As our motivation is to solve Maxwell's equations we take care that these infinite elements fit into the discrete de Rham diagram, i.e. they span discrete spaces, which together with the exterior derivative form an exact sequence. Resonance as well as scattering problems are considered in the examples. Numerical tests indicate super-algebraic convergence in the number of additional unknowns per degree of freedom on the coupling boundary that are required to realize the Dirichlet to Neumann map.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Iteratively Regularized Gauß-Newton Method with Convex Constraints and Applications in 4Pi-Microscopy

This paper is concerned with the numerical solution of nonlinear ill-posed operator equations involving convex constraints. We study a Newton-type method which consists in applying linear Tikhonov regularization with convex constraints to the Newton equations in each iteration step. Convergence of this iterative regularization method is analyzed if both the operator and the right hand side are given with errors and all error levels tend to zero. Our study has been motivated by the joint estimation of object and phase in 4Pi microscopy, which leads to a semi-blind deconvolution problem with nonnegativity constraints. The performance of the proposed algorithm is illustrated both for simulated and for three-dimensional experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Acceleration techniques for regularized Newton methods applied to electromagnetic inverse medium scattering problems

We study the construction and updating of spectral preconditioners for regularized Newton methods and their application to electromagnetic inverse medium scattering problems. Moreover, we show how a Lepskiĭ-type stopping rule can be implemented efficiently for these methods. In numerical examples, the proposed method compares favorably with other iterative regularization method in terms of work-precision diagrams for exact data. For data perturbed by random noise, the Lepskiĭ-type stopping rule performs considerably better than the commonly used discrepancy principle.

0 citations0 reviewsNo code linkedOpen access