Iterated Tikhonov regularization with a general penalty term

Tikhonov regularization is one of the most popular approaches to solving linear discrete ill‐posed problems. The choice of the regularization matrix may significantly affect the quality of the computed solution. When the regularization matrix is the identity, iterated Tikhonov regularization can yield computed approximate solutions of higher quality than (standard) Tikhonov regularization. This paper provides an analysis of iterated Tikhonov regularization with a regularization matrix different from the identity. Computed examples illustrate the performance of this method.     

反向热传导问题的一种新的正则化方法

讨论一个高维反向热传导问题,这是一个经典的严重不适定问题.关于这一问题我们给出一种新的正则化方法-改进的Tikhonov正则化方法,以恢复解对数据的连续依赖性.通过构造一个重要的不等式和提高先验光滑条件,获得正则解在0相似文献   

求解椭圆方程柯西问题的修正吉洪诺夫正则化方法

研究了一类变系数椭圆方程的柯西问题,这类问题出现在很多实际问题领域.由于问题的不适定性,不可能通过经典的数值方法来求解上述问题,必须引入正则化手段.采用了一种修正吉洪诺夫正则化方法来求解上述问题.在一种先验和一种后验参数选取准则下,分别获得了问题的误差估计.数值例子进一步显示方法是稳定有效的.     

约束Fractional Tikhonov正则化的模迭代方法

反问题是现在数学物理研究中的一个热点问题,而反问题求解面临的一个本质性困难是不适定性。求解不适定问题的普遍方法是:用与原不适定问题相“邻近”的适定问题的解去逼近原问题的解,这种方法称为正则化方法.如何建立有效的正则化方法是反问题领域中不适定问题研究的重要内容.当前,最为流行的正则化方法有基于变分原理的Tikhonov正则化及其改进方法,此类方法是求解不适定问题的较为有效的方法,在各类反问题的研究中被广泛采用,并得到深入研究.     

Identification of a time-dependent source term for a time fractional diffusion problem

In this paper, we study an inverse problem of identifying a time-dependent term of an unknown source for a time fractional diffusion equation using nonlocal measurement data. Firstly, we establish the conditional stability for this inverse problem. Then two regularization methods are proposed to for reconstructing the time-dependent source term from noisy measurements. The first method is an integral equation method which formulates the inverse source problem into an integral equation of the second kind; and a prior convergence rate of regularized solutions is derived with a suitable choice strategy of regularization parameters. The second method is a standard Tikhonov regularization method and formulates the inverse source problem as a minimizing problem of the Tikhonov functional. Based on the superposition principle and the technique of finite-element interpolation, a numerical scheme is proposed to implement the second regularization method. One- and two-dimensional examples are carried out to verify efficiency and stability of the second regularization method.     

The Levenberg-Marquardt method applied to a parameter estimation problem arising from electrical resistivity tomography

An efficient and robust electrical resistivity tomographic inversion algorithm based on the Levenberg-Marquardt method is considered to obtain quantities like grain size, spatial scale and particle size distribution of mineralized rocks. The corresponding model in two dimensions is based on the Maxwell equations and leads to a partial differential equation with mixed Dirichlet-Neumann boundary conditions. The forward problem is solved numerically with the finite-difference method. However, the inverse problem at hand is a classical nonlinear and ill-posed parameter estimation problem. Linearizing and applying the Tikhonov regularization method yields an iterative scheme, the Levenberg-Marquardt method. Several large systems of equations have to be solved efficiently in each iteration step which is accomplished by the conjugate gradient method without setting up the corresponding matrix. Instead fast matrix-vector multiplications are performed directly. Therefore, the derivative and its adjoint for the parameter-to-solution map are needed. Numerical results demonstrate the performance of our method as well as the possibility to reconstruct some of the desired parameters.     

Nonnegative ill-conditioned linear systems and GMRES method

In this paper we consider solving ill-conditioned linear systems under nonnegativity constraints with noisy right hand sides. The classical approaches to solve such systems are constrained least square (quadratic programming) and barrier methods. First we present these classical methods. Then a modified version of the GMRES method (NGMRES) is presented. Since we assume that the coefficient matrices are ill-conditioned, then the Tikhonov regularization of the problem is considered. Our computational experiments show that the NGMRES provides us high quality solutions much faster than the other two approaches.     

Tikhonov regularization method for a backward problem for the time-fractional diffusion equation

This paper is devoted to solve a backward problem for a time-fractional diffusion equation with variable coefficients in a general bounded domain by the Tikhonov regularization method. Based on the eigenfunction expansion of the solution, the backward problem for searching the initial data is changed to solve a Fredholm integral equation of the first kind. The conditional stability for the backward problem is obtained. We use the Tikhonov regularization method to deal with the integral equation and obtain the series expression of solution. Furthermore, the convergence rates for the Tikhonov regularized solution can be proved by using an a priori regularization parameter choice rule and an a posteriori regularization parameter choice rule. Two numerical examples in one-dimensional and two-dimensional cases respectively are investigated. Numerical results show that the proposed method is effective and stable.     

Tikhonov Regularization of Large Linear Problems

Many numerical methods for the solution of linear ill-posed problems apply Tikhonov regularization. This paper presents a new numerical method, based on Lanczos bidiagonalization and Gauss quadrature, for Tikhonov regularization of large-scale problems. An estimate of the norm of the error in the data is assumed to be available. This allows the value of the regularization parameter to be determined by the discrepancy principle.     

The Tikhonov regularization extended to equilibrium problems involving pseudomonotone bifunctions

We extend the Tikhonov regularization method widely used in optimization and monotone variational inequality studies to equilibrium problems. It is shown that the convergence results obtained from the monotone variational inequality remain valid for the monotone equilibrium problem. For pseudomonotone equilibrium problems, the Tikhonov regularized subproblems have a unique solution only in the limit, but any Tikhonov trajectory tends to the solution of the original problem, which is the unique solution of the strongly monotone equilibrium problem defined on the basis of the regularization bifunction.     

LOCAL REGULARIZATION METHODS FOR THE STABILIZATION OF LINEAR ILL-POSED EQUATIONS OF VOLTERRA TYPE

Inverse problems based on first-kind Volterra integral equations appear naturally in the study of many applications, from geophysical problems to the inverse heat conduction problem. The ill-posedness of such problems means that a regularization technique is required, but classical regularization schemes like Tikhonov regularization destroy the causal nature of the underlying Volterra problem and, in general, can produce oversmoothed results. In this paper we investigate a class of local regularization methods in which the original (unstable) problem is approximated by a parameterized family of well-posed, second-kind Volterra equations. Being Volterra, these approximating second-kind equations retain the causality of the original problem and allow for quick sequential solution techniques. In addition, the regularizing method we develop is based on the use of a regularization parameter which is a function (rather than a single constant), allowing for more or less smoothing at localized points in the domain. We take this approach even further by adopting the flexibility of an additional penalty term (with variable penalty function) and illustrate the sequential selection of the penalty function in a numerical example.     

An inverse problem for space‐fractional backward diffusion problem

In this paper, an inverse problem for space‐fractional backward diffusion equation, which is highly ill‐posed, is considered. This problem is obtained from the classical diffusion equation by replacing the second‐order space derivative with a Riesz–Feller derivative of order α ∈ (0,2]. We show that such a problem is severely ill‐posed, and further present a simplified Tikhonov regularization method to deal with this problem. Convergence estimate is presented under a priori choice of regularization parameter. Numerical experiments are given to illustrate the accuracy and efficiency of the proposed method. Copyright © 2013 John Wiley & Sons, Ltd.     

双调和方程柯西问题的修正的Tikhonov正则化方法

本文研究了双调和方程柯西问题,这类是不适定的,即问题的解(如果存在)不连续依赖于测量数据.首先在精确解的先验假设下给出问题的条件稳定性结果.接着利用修正的Tikhonov正则化方法求解此不适定问题.在先验和后验正则化参数选取规则下,给出正则解和精确解之间的误差估计式.最后给出几个数值例子验证此正则化方法求解此类反问题的有效性.     

关于一类广义Tikhonov 正则化方法的饱和效应分析

Tikhonov正则化方法是研究不适定问题最重要的正则化方法之一,但由于这种方法的饱和效应出现的太早,使得无法随着对解的光滑性假设的提高而提高正则逼近解的收敛率,也即对高的光滑性假设,正则解与准确解的误差估计不可能达到阶数最优.Schrroter T 和Tautenhahn U给出了一类广义Tikhonov正则化方法并重点讨论了它的最优误差估计, 但却未能对该方法的饱和效应进行研究.本文对此进行了仔细分析,并发现此方法可以防止饱和效应,而且数值试验结果表明此方法计算效果良好.     

Arnoldi–Tikhonov regularization methods

Tikhonov regularization for large-scale linear ill-posed problems is commonly implemented by determining a partial Lanczos bidiagonalization of the matrix of the given system of equations. This paper explores the possibility of instead computing a partial Arnoldi decomposition of the given matrix. Computed examples illustrate that this approach may require fewer matrix–vector product evaluations and, therefore, less arithmetic work. Moreover, the proposed range-restricted Arnoldi–Tikhonov regularization method does not require the adjoint matrix and, hence, is convenient to use for problems for which the adjoint is difficult to evaluate.     

On nondecreasing sequences of regularization parameters for nonstationary iterated Tikhonov

Nonstationary iterated Tikhonov is an iterative regularization method that requires a strategy for defining the Tikhonov regularization parameter at each iteration and an early termination of the iterative process. A classical choice for the regularization parameters is a decreasing geometric sequence which leads to a linear convergence rate. The early iterations compute quickly a good approximation of the true solution, but the main drawback of this choice is a rapid growth of the error for later iterations. This implies that a stopping criteria, e.g. the discrepancy principle, could fail in computing a good approximation. In this paper we show by a filter factor analysis that a nondecreasing sequence of regularization parameters can provide a rapid and stable convergence. Hence, a reliable stopping criteria is no longer necessary. A geometric nondecreasing sequence of the Tikhonov regularization parameters into a fixed interval is proposed and numerically validated for deblurring problems.     

第一类算子方程的一种新的迭代正则化方法

提出了一种新的解第一类算子方程的迭代正则化方法,与通常的迭代正则化方法相比,提高了j次迭代正则解的渐近阶估计.同时,给出了后验正则化参数的选择.     

The revised generalized Tikhonov method for the backward time-fractional diffusion equation

In this paper, we solve the backward problem for a time-fractional diffusion equation with variable coefficients in a bounded domain by using the revised generalized Tikhonov regularization method. Convergence estimates under an a-priori and a-posteriori regularization parameter choice rules are given. Numerical example shows that the proposed method is effective and stable.     

Error estimates for large-scale ill-posed problems

The computation of an approximate solution of linear discrete ill-posed problems with contaminated data is delicate due to the possibility of severe error propagation. Tikhonov regularization seeks to reduce the sensitivity of the computed solution to errors in the data by replacing the given ill-posed problem by a nearby problem, whose solution is less sensitive to perturbation. This regularization method requires that a suitable value of the regularization parameter be chosen. Recently, Brezinski et al. (Numer Algorithms 49, 2008) described new approaches to estimate the error in approximate solutions of linear systems of equations and applied these estimates to determine a suitable value of the regularization parameter in Tikhonov regularization when the approximate solution is computed with the aid of the singular value decomposition. This paper discusses applications of these and related error estimates to the solution of large-scale ill-posed problems when approximate solutions are computed by Tikhonov regularization based on partial Lanczos bidiagonalization of the matrix. The connection between partial Lanczos bidiagonalization and Gauss quadrature is utilized to determine inexpensive bounds for a family of error estimates. In memory of Gene H. Golub. This work was supported by MIUR under the PRIN grant no. 2006017542-003 and by the University of Cagliari.     

Computation of Control for Linear Approximately Controllable System Using Tikhonov Regularization

Computation of control for a controlled partial differential equation is a di?cult task, especially when the control problem is ill posed. In this paper, we propose a method of computing the regularized control of a diffusion control system using Tikhonov regularization approach when the system is approximately controllable. The method proposed here for choosing regularization parameter guarantees the convergence of the proposed control.