Source author record

Bernd Hofmann

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

14 published item(s)

preprint2022arXiv

Open questions asked to analysis and numerics concerning the Hausdorff moment problem

We address facts and open questions concerning the degree of ill-posedness of the composite Hausdorff moment problem aimed at the recovery of a function $x \in L^2(0,1)$ from elements of the infinite dimensional sequence space $\ell^2$ that characterize moments applied to the antiderivative of $x$. This degree, unknown by now, results from the decay rate of the singular values of the associated compact forward operator $A$, which is the composition of the compact simple integration operator mapping in $L^2(0,1)$ and the non-compact Hausdorff moment operator $B^{(H)}$ mapping from $L^2(0,1)$ to $\ell^2$. There is a seeming contradiction between (a) numerical computations, which show (even for large $n$) an exponential decay of the singular values for $n$-dimensional matrices obtained by discretizing the operator $A$, and \linebreak (b) a strongly limited smoothness of the well-known kernel $k$ of the Hilbert-Schmidt operator $A^*A$. Fact (a) suggests severe ill-posedness of the infinite dimensional Hausdorff moment problem, whereas fact (b) lets us expect the opposite, because exponential ill-posedness occurs in common just for $C^\infty$-kernels $k$. We recall arguments for the possible occurrence of a polynomial decay of the singular values of $A$, even if the numerics seems to be against it, and discuss some issues in the numerical approximation of non-compact operators.

preprint2020arXiv

Convergence results and low order rates for nonlinear Tikhonov regularization with oversmoothing penalty term

For the Tikhonov regularization of ill-posed nonlinear operator equations, convergence is studied in a Hilbert scale setting. We include the case of oversmoothing penalty terms, which means that the exact solution does not belong to the domain of definition of the considered penalty functional. In this case, we try to close a gap in the present theory, where Hölder-type convergence rates results have been proven under corresponding source conditions, but assertions on norm convergence of regularized solutions without source conditions are completely missing. A result of the present work is to provide sufficient conditions for convergence under a priori and a posteriori regularization parameter choice strategies, without any additional smoothness assumption on the solution. The obtained error estimates moreover allow us to prove low order convergence rates under associated (for example logarithmic) source conditions. Some numerical illustrations are also given.

preprint2020arXiv

Two new non-negativity preserving iterative regularization methods for ill-posed inverse problems

Many inverse problems are concerned with the estimation of non-negative parameter functions. In this paper, in order to obtain non-negative stable approximate solutions to ill-posed linear operator equations in a Hilbert space setting, we develop two novel non-negativity preserving iterative regularization methods. They are based on fixed point iterations in combination with preconditioning ideas. In contrast to the projected Landweber iteration, for which only weak convergence can be shown for the regularized solution when the noise level tends to zero, the introduced regularization methods exhibit strong convergence. There are presented convergence results, even for a combination of noisy right-hand side and imperfect forward operators, and for one of the approaches there are also convergence rates results. Specifically adapted discrepancy principles are used as a posteriori stopping rules of the established iterative regularization algorithms. For an application of the suggested new approaches, we consider a biosensor problem, which is modelled as a two dimensional linear Fredholm integral equation of the first kind. Several numerical examples, as well as a comparison with the projected Landweber method, are presented to show the accuracy and the acceleration effect of the novel methods. Case studies of a real data problem indicate that the developed methods can produce meaningful featured regularized solutions.

preprint2019arXiv

Convergence Analysis of (Statistical) Inverse Problems under Conditional Stability Estimates

Conditional stability estimates require additional regularization for obtaining stable approximate solutions if the validity area of such estimates is not completely known. In this context, we consider ill-posed nonlinear inverse problems in Hilbert scales satisfying conditional stability estimates characterized by general concave index functions. For that case, we exploit Tikhonov regularization and provide convergence and convergence rates of regularized solutions for both deterministic and stochastic noise. We further discuss a priori and a posteriori parameter choice rules and illustrate the validity of our assumptions in different model and real world situations.

preprint2018arXiv

On the singular value decomposition of n-fold integration operators

In theory and practice of inverse problems, linear operator equations $Tx=y$ with compact linear forward operators $T$ having a non-closed range $\mathcal{R}(T)$ and mapping between infinite dimensional Hilbert spaces plays some prominent role. As a consequence of the ill-posedness of such problems, regularization approaches are required, and due to its unlimited qualification spectral cut-off is an appropriate method for the stable approximate solution of corresponding inverse problems. For this method, however, the singular system $\{σ_i(T),u_i(T),v_i(T)\}_{i=1}^\infty$ of the compact operator $T$ is needed, at least for $i=1,2,...,N$, up to some stopping index $N$. In this note we consider $n$-fold integration operators $T=J^n\;(n=1,2,...)$ in $L^2([0,1])$ occurring in numerous applications, where the solution of the associated operator equation is characterized by the $n$-th generalized derivative $x=y^{(n)}$ of the Sobolev space function $y \in H^n([0,1])$. Almost all textbooks on linear inverse problems present the whole singular system $\{σ_i(J^1),u_i(J^1),v_i(J^1)\}_{i=1}^\infty$ in an explicit manner. However, they do not discuss the singular systems for $J^n,\;n \ge 2$. We will emphasize that this seems to be a consequence of the fact that for higher $n$ the eigenvalues $σ^2_i(J^n)$ of the associated ODE boundary value problems obey transcendental equations, the complexity of which is growing with $n$. We present the transcendental equations for $n=2,3,...$ and discuss and illustrate the associated eigenfunctions and some of their properties.

preprint2016arXiv

Conditional stability versus ill-posedness for operator equations with monotone operators in Hilbert space

In the literature on singular perturbation (Lavrentiev regularization) for the stable approximate solution of operator equations with monotone operators in the Hilbert space the phenomena of conditional stability and local well-posedness and ill-posedness are rarely investigated. Our goal is to present some studies which try to bridge this gap. So we discuss the impact of conditional stability on error estimates and convergence rates for the Lavrentiev regularization and distinguish for linear problems well-posedness and ill-posedness in a specific manner motivated by a saturation result. The role of the regularization error in the noise-free case, called bias, is a crucial point in the paper for nonlinear and linear problems. In particular, for linear operator equations general convergence rates, including logarithmic rates, are derived by means of the method of approximate source conditions. This allows us to extend well-known convergence rates results for the Lavrentiev regularization that were based on general source conditions to the case of non-selfadjoint linear monotone forward operators for which general source conditions fail. Examples presenting the self-adjoint multiplication operator as well as the non-selfadjoint fractional integral operator and Cesàro operator illustrate the theoretical results. Extensions to the nonlinear case under specific conditions on the nonlinearity structure complete the paper.

preprint2016arXiv

Elastic-net regularization versus $\ell^1$-regularization for linear inverse problems with quasi-sparse solutions

We consider the ill-posed operator equation $Ax=y$ with an injective and bounded linear operator $A$ mapping between $\ell^2$ and a Hilbert space $Y$, possessing the unique solution \linebreak $x^†=\{x^†_k\}_{k=1}^\infty$. For the cases that sparsity $x^†\in \ell^0$ is expected but often slightly violated in practice, we investigate in comparison with the $\ell^1$-regularization the elastic-net regularization, where the penalty is a weighted superposition of the $\ell^1$-norm and the $\ell^2$-norm square, under the assumption that $x^†\in \ell^1$. There occur two positive parameters in this approach, the weight parameter $η$ and the regularization parameter as the multiplier of the whole penalty in the Tikhonov functional, whereas only one regularization parameter arises in $\ell^1$-regularization. Based on the variational inequality approach for the description of the solution smoothness with respect to the forward operator $A$ and exploiting the method of approximate source conditions, we present some results to estimate the rate of convergence for the elastic-net regularization. The occurring rate function contains the rate of the decay $x^†_k \to 0$ for $k \to \infty$ and the classical smoothness properties of $x^†$ as an element in $\ell^2$.

preprint2016arXiv

Lavrentiev's regularization method in Hilbert spaces revisited

In this paper, we deal with nonlinear ill-posed problems involving monotone operators and consider Lavrentiev's regularization method. This approach, in contrast to Tikhonov's regularization method, does not make use of the adjoint of the derivative. There are plenty of qualitative and quantitative convergence results in the literature, both in Hilbert and Banach spaces. Our aim here is mainly to contribute to convergence rates results in Hilbert spaces based on some types of error estimates derived under various source conditions and to interpret them in some settings. In particular, we propose and investigate new variational source conditions adapted to these Lavrentiev-type techniques. Another focus of this paper is to exploit the concept of approximate source conditions.

preprint2015arXiv

On $\ell^1$-regularization in light of Nashed's ill-posedness concept

Based on the powerful tool of variational inequalities, in recent papers convergence rates results on $\ell^1$-regularization for ill-posed inverse problems have been formulated in infinite dimensional spaces under the condition that the sparsity assumption slightly fails, but the solution is still in $\ell^1$. In the present paper we improve those convergence rates results and apply them to the Cesáro operator equation in $\ell^2$ and to specific denoising problems. Moreover, we formulate in this context relationships between Nashed's types of ill-posedness and mapping properties like compactness and strict singularity.

preprint2015arXiv

Variational regularization of complex deautoconvolution and phase retrieval in ultrashort laser pulse characterization

The SD-SPIDER method for the characterization of ultrashort laser pulses requires the solution of a nonlinear integral equation of autoconvolution type with a device-based kernel function. Taking into account the analytical background of a variational regularization approach for solving the corresponding ill-posed operator equation formulated in complex-valued $L^2$-spaces over finite real intervals, we suggest and evaluate numerical procedures using NURBS and the TIGRA method for calculating the regularized solutions in a stable manner. In this context, besides the complex deautoconvolution problem with noisy but full data, a phase retrieval problem is introduced which adapts to the experimental state of the art in laser optics. For the treatment of this problem facet, which is formulated as a tensor product operator equation, we derive well-posedness of variational regularization methods. Case studies with synthetic and real optical data show the capability of the implemented approach as well as its limitation due to measurement deficits.

preprint2013arXiv

On the interplay of basis smoothness and specific range conditions occurring in sparsity regularization

The convergence rates results in $\ell^1$-regularization when the sparsity assumption is narrowly missed, presented by Burger et al. (2013 Inverse Problems 29 025013), are based on a crucial condition which requires that all basis elements belong to the range of the adjoint of the forward operator. Partly it was conjectured that such a condition is very restrictive. In this context, we study sparsity-promoting varieties of Tikhonov regularization for linear ill-posed problems with respect to an orthonormal basis in a separable Hilbert space using $\ell^1$ and sublinear penalty terms. In particular, we show that the corresponding range condition is always satisfied for all basis elements if the problems are well-posed in a certain weaker topology and the basis elements are chosen appropriately related to an associated Gelfand triple. The Radon transform, Symm's integral equation and linear integral operators of Volterra type are examples for such behaviour, which allows us to apply convergence rates results for non-sparse solutions, and we further extend these results also to the case of non-convex $\ell^q$-regularization with 0<q<1.

preprint2013arXiv

Regularization of an autoconvolution problem in ultrashort laser pulse characterization

An ill-posed inverse problem of autoconvolution type is investigated. This inverse problem occurs in nonlinear optics in the context of ultrashort laser pulse characterization. The novelty of the mathematical model consists in a physically required extension of the deautoconvolution problem beyond the classical case usually discussed in literature: (i) For measurements of ultrashort laser pulses with the self-diffraction SPIDER method, a stable approximate solution of an autocovolution equation with a complex-valued kernel function is needed. (ii) The considered scenario requires complex functions both, in the solution and the rhs of the integral equation. Since, however, noisy data are available not only for amplitude and phase functions of the rhs, but also for the amplitude of the solution, the stable approximate reconstruction of the associated smooth phase function represents the main goal of the paper. An iterative regularization approach is described that is specifically adapted to the physical situation in pulse characterization, using a non-standard stopping rule for the iteration process of computing regularized solutions. Our approach is illustrated by several case studies for synthetic noisy data and physically realistic complex-valued kernel functions. Based on an example with focus on amplitude perturbations, we show that the autoconvolution equation is locally ill-posed everywhere. To date, the analytical treatment of the impact of noisy data on phase perturbations remains an open question. However, we show its influence with the help of numerical experiments. Moreover, we formulate assertions on the non-uniqueness of the complex-valued autoconvolution problem, at least for the simplified case of a constant kernel. The presented results and figures associated with case studies illustrate the ill-posedness phenomena also for the case of non-trivial complex kernel functions.

preprint2012arXiv

Convergence rates in $\mathbf{\ell^1}$-regularization if the sparsity assumption fails

Variational sparsity regularization based on $\ell^1$-norms and other nonlinear functionals has gained enormous attention recently, both with respect to its applications and its mathematical analysis. A focus in regularization theory has been to develop error estimation in terms of regularization parameter and noise strength. For this sake specific error measures such as Bregman distances and specific conditions on the solution such as source conditions or variational inequalities have been developed and used. In this paper we provide, for a certain class of ill-posed linear operator equations, a convergence analysis that works for solutions that are not completely sparse, but have a fast decaying nonzero part. This case is not covered by standard source conditions, but surprisingly can be treated with an appropriate variational inequality. As a consequence the paper also provides the first examples where the variational inequality approach, which was often believed to be equivalent to appropriate source conditions, can indeed go farther than the latter.

preprint2010arXiv

Errors of regularisation under range inclusions using variable Hilbert scales

Based on the variable Hilbert scale interpolation inequality bounds for the error of regularisation methods are derived under range inclusions. In this context, new formulae for the modulus of continuity of the inverse of bounded operators with non-closed range are given. Even if one can show the equivalence of this approach to the version used previously in the literature, the new formulae and corresponding conditions are simpler than the former ones. Several examples from image processing and spectral enhancement illustrate how the new error bounds can be applied.