MR image reconstruction from undersampled data by using the iterative refinement procedure

Magnetic resonance imaging (MRI) reconstruction from sparsely sampled data has been a difficult problem in medical imaging field. We approach this problem by formulating a cost functional that includes a constraint term that is imposed by the raw measurement data in k-space and the L¹ norm of a sparse representation of the reconstructed image. The sparse representation is usually realized by total variational regularization and/or wavelet transform. We have applied the Bregman iteration to minimize this functional to recover finer scales in our recent work. Here we propose nonlinear inverse scale space methods in addition to the iterative refinement procedure. Numerical results from the two methods are presented and it shows that the nonlinear inverse scale space method is a more efficient algorithm than the iterated refinement method. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the order optimality of the regularization via inexact Newton iterations

Inexact Newton regularization methods have been proposed by Hanke and Rieder for solving nonlinear ill-posed inverse problems. Every such a method consists of two components: an outer Newton iteration and an inner scheme providing increments by regularizing local linearized equations. The method is terminated by a discrepancy principle. In this paper we consider the inexact Newton regularization methods with the inner scheme defined by Landweber iteration, the implicit iteration, the asymptotic regularization and Tikhonov regularization. Under certain conditions we obtain the order optimal convergence rate result which improves the suboptimal one of Rieder. We in fact obtain a more general order optimality result by considering these inexact Newton methods in Hilbert scales.     

基于遗传算法重建多个散射体的组合Newton法

本文研究由单个入射声波或电磁波及其远场数据反演多个柔性散射体边界的逆散射问题.通过建立边界到边界总场的非线性算子及其n6chet导数,本文首先给出了基于单层位势的组合Newton法.将组合Newton法转化为泛响优化问题,从而获得了该方法重建单个散射体的收敛性分析.然后,基于遗传算法和正则化参数选取的模型函数方法,给出...     

A multi-scale image reconstruction algorithm for electrical capacitance tomography

Electrical capacitance tomography (ECT) is considered as a promising process tomography (PT) technology, and its successful applications depend mainly on the precision and speed of the image reconstruction algorithms. In this paper, based on the wavelet multi-scale analysis method, an efficient image reconstruction algorithm is presented. The original inverse problem is decomposed into a sequence of inverse problems, which are solved successively from the largest scale to the smallest scale. At different scales, the inverse problem is solved by a generalized regularized total least squares (TLS) method, which is developed using a combinational minimax estimation method and an extended stabilizing functional, until the solution of the original inverse problem is found. The homotopy algorithm is employed to solve the objective functional. The proposed algorithm is tested by the noise-free capacitance data and the noise-contaminated capacitance data, and excellent numerical performances and satisfactory results are observed. In the cases considered in this paper, the reconstruction results show remarkable improvement in the accuracy. The spatial resolution of the reconstructed images by the proposed algorithm is enhanced and the artifacts in the reconstructed images can be eliminated effectively. As a result, a promising algorithm is introduced for ECT image reconstruction.     

Cubically convergent methods for selecting the regularization parameters in linear inverse problems

We present three cubically convergent methods for choosing the regularization parameters in linear inverse problems. The detailed algorithms are given and the convergence rates are estimated. Our basic tools are Tikhonov regularization and Morozov's discrepancy principle. We prove that, in comparison with the standard Newton method, the computational costs for our cubically convergent methods are nearly the same, but the number of iteration steps is even less. Numerical experiments for an elliptic boundary value problem illustrate the efficiency of the proposed algorithms.     

Semi-smooth Newton methods for the Signorini problem

Semi-smooth Newton methods are analyzed for the Signorini problem. A proper regularization is introduced which guarantees that the semi-smooth Newton method is superlinearly convergent for each regularized problem. Utilizing a shift motivated by an augmented Lagrangian framework, to the regularization term, the solution to each regularized problem is feasible. Convergence of the regularized problems is shown and a report on numerical experiments is given.     

Inexact Newton regularization methods in Hilbert scales

We consider a class of inexact Newton regularization methods for solving nonlinear inverse problems in Hilbert scales. Under certain conditions we obtain the order optimal convergence rate result.     

On decomposition-based block preconditioned iterative methods for half-quadratic image restoration

Image restoration is a fundamental problem in image processing. Except for many different filters applied to obtain a restored image in image restoration, a degraded image can often be recovered efficiently by minimizing a cost function which consists of a data-fidelity term and a regularization term. In specific, half-quadratic regularization can effectively preserve image edges in the recovered images and a fixed-point iteration method is usually employed to solve the minimization problem. In this paper, the Newton method is applied to solve the half-quadratic regularization image restoration problem. And at each step of the Newton method, a structured linear system of a symmetric positive definite coefficient matrix arises. We design two different decomposition-based block preconditioning matrices by considering the special structure of the coefficient matrix and apply the preconditioned conjugate gradient method to solve this linear system. Theoretical analysis shows the eigenvector properties and the spectral bounds for the preconditioned matrices. The method used to analyze the spectral distribution of the preconditioned matrix and the correspondingly obtained spectral bounds are different from those in the literature. The experimental results also demonstrate that the decomposition-based block preconditioned conjugate gradient method is efficient for solving the half-quadratic regularization image restoration in terms of the numerical performance and image recovering quality.     

A modified quasi-boundary value method for an inverse source problem of the time-fractional diffusion equation

In this paper, we consider an inverse source problem for a time-fractional diffusion equation with variable coefficients in a general bounded domain. That is to determine a space-dependent source term in the time-fractional diffusion equation from a noisy final data. Based on a series expression of the solution, we can transform the original inverse problem into a first kind integral equation. The uniqueness and a conditional stability for the space-dependent source term can be obtained. Further, we propose a modified quasi-boundary value regularization method to deal with the inverse source problem and obtain two kinds of convergence rates by using an a priori and an a posteriori regularization parameter choice rule, respectively. Numerical examples in one-dimensional and two-dimensional cases are provided to show the effectiveness of the proposed method.     

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.     

Solving an inverse problem for a generalized time-delayed Burgers-Fisher equation by Haar wavelet method

In this paper, a numerical method consists of combining Haar wavelet method and Tikhonov regularization method to determine unknown boundary condition and unknown nonlinear source term for the generalized time-delayed Burgers-Fisher equation using noisy data is presented. A stable numerical solution is determined for the problem. We also show that the rate of convergence of the method is as exponential $\Bigl(O\left(\frac{1}{2^{J+1}}\right)\Bigr)$, where $J$ is maximal level of resolution of wavelet. Some numerical results are reported to show the efficiency and robustness of the proposed approach for solving the inverse problems.     

Regularized state-constrained boundary optimal control of the Navier-Stokes equations

The numerical solution of the Dirichlet boundary optimal control problem of the Navier-Stokes equations in presence of pointwise state constraints is investigated. Two different regularization techniques are considered. First, a Moreau-Yosida regularization of the problem is studied. Optimality conditions are derived and the convergence of the regularized solutions towards the original one is proved. A source representation of the control combined with a Lavrentiev type regularization strategy is also presented. The analysis concerning optimality conditions and convergence of the regularized solutions is carried out. In the last part of the paper numerical experiments are presented. For the numerical solution of each regularized problem a semi-smooth Newton method is applied.     

An Infinite-Dimensional Semismooth Newton Method for Elasto-Plastic Contact Problems

A Fenchel dualization scheme for the one-step time-discretized elasto-plastic contact problem with kinematic or isotropic hardening is considered. The associated path is induced by a combined Moreau-Yosida / Tichonov regularization of the dual problem. The sequence of solutions to the regularized problems is shown to converge strongly to the solution of the original problem. This property relies on the density of the intersection of certain convex sets. The corresponding conditions are worked out and customary regularization approaches are shown to be valid in this context. It is also argued that without higher regularity assumptions on the data the resulting problems possess Newton differentiable optimality systems in infinite dimensions [2]. Consequently, each regularized subsystem can be solved mesh-independently at a local superlinear rate of convergence [6]. Numerically the problems are solved using conforming finite elements. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

动力系统方法解决无界线性算子方程

针对反问题中出现的第一类算子方程Au=f,其中A是实Hilbert空间H上的一个无界线性算子利用动力系统方法和正则化方法,求解上述问题的正则化问题的解:u'(t)=-A~*(Au(t)-f)利用线性算子半群理论可以得到上述正则化问题的解的半群表示,并证明了当t→∞时,所得的正则化解收敛于原问题的解.     

Globally Convergent BFGS Method for Nonsmooth Convex Optimization1

We propose an implementable BFGS method for solving a nonsmooth convex optimization problem by converting the original objective function into a once continuously differentiable function by way of the Moreau–Yosida regularization. The proposed method makes use of approximate function and gradient values of the Moreau-Yosida regularization instead of the corresponding exact values. We prove the global convergence of the proposed method under the assumption of strong convexity of the objective function.     

The Sinc‐Galerkin method for solving an inverse parabolic problem with unknown source term

The inverse problem of determining an unknown source term depending on space variable in a parabolic equation is considered. A numerical algorithm is presented for recovering the unknown function and obtaining a solution of the problem. As this inverse problem is ill‐posed, Tikhonov regularization is used for finding a stable solution. For solving the direct problem, a Galerkin method with the Sinc basis functions in both the space and time domains is presented. This approximate solution displays an exponential convergence rate and is valid on the infinite time interval. Finally, some examples are presented to illustrate the ability and efficiency of this numerical method. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

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.     

An implicit Landweber method for nonlinear ill-posed operator equations

In this paper, we are interested in the solution of nonlinear inverse problems of the form F(x)=y. We propose an implicit Landweber method, which is similar to the third-order midpoint Newton method in form, and consider the convergence behavior of the implicit Landweber method. Using the discrepancy principle as a stopping criterion, we obtain a regularization method for ill-posed problems. We conclude with numerical examples confirming the theoretical results, including comparisons with the classical Landweber iteration and presented modified Landweber methods.     

A semismooth Newton method for a class of semilinear optimal control problems with box and volume constraints

In this paper we consider optimal control problems subject to a semilinear elliptic state equation together with the control constraints 0≤u≤1 and ∫u=m. Optimality conditions for this problem are derived and reformulated as a nonlinear, nonsmooth equation which is solved using a semismooth Newton method. A regularization of the nonsmooth equation is necessary to obtain the superlinear convergence of the semismooth Newton method. We prove that the solutions of the regularized problems converge to a solution of the original problem and a path-following technique is used to ensure a constant decrease rate of the residual. We show that, in certain situations, the optimal controls take 0–1 values, which amounts to solving a topology optimization problem with volume constraint.     

A one-parametric class of smoothing functions and an improved regularization Newton method for the NCP

In this paper, we introduce a one-parametric class of smoothing functions, which enjoys some favourable properties and includes two famous smoothing functions as special cases. Based on this class of smoothing functions, we propose a regularization Newton method for solving the non-linear complementarity problem. The main feature of the proposed method is that it uses a perturbed Newton equation to obtain the direction. This not only allows our method to have global and local quadratic convergences without strict complementarity conditions, but also makes the regularization parameter converge to zero globally Q-linearly. In addition, we use a new non-monotone line search scheme to obtain the step size. Some numerical results are reported which confirm the good theoretical properties of the proposed method.