Convergent Data-Driven Regularizations for CT Reconstruction

The reconstruction of images from their corresponding noisy Radon transform is a typical example of an ill-posed linear inverse problem as arising in the application of computerized tomography(CT).As the(naive)solution does not depend on the measured data continuously,regularization is needed to reestablish a continuous dependence.In this work,we investigate simple,but yet still provably convergent approaches to learning linear regularization methods from data.More specifically,we analyze two approaches:one generic linear regularization that learns how to manipulate the singular values of the linear operator in an extension of our previous work,and one tailored approach in the Fourier domain that is specific to CT-reconstruction.We prove that such approaches become convergent regularization methods as well as the fact that the reconstructions they provide are typically much smoother than the training data they were trained on.Finally,we compare the spectral as well as the Fourier-based approaches for CT-reconstruction numerically,discuss their advantages and disadvantages and investigate the effect of discretization errors at differentresolutions.

The Regularity of Solutions to Mixed Boundary Value Problems of Second-Order Elliptic Equations with Small Angles

This paper considers the regularity of solutions to mixed boundary value problems in small-angle regions for elliptic equations. By constructing a specific barrier function, we proved that under the assumption of sufficient regularity of boundary conditions and coefficients, as long as the angle is sufficiently small, the regularity of the solution to the mixed boundary value problem of the second-order elliptic equation can reach any order.

The Criteria for Reducing Centrally Restricted Three-Body Problem to Two-Body Problem

Our Solar System contains eight planets and their respective natural satellites excepting the inner two planets Mercury and Venus. A satellite hosted by a given Planet is well protected by the gravitational pertubation of much heavier planets such as Jupiter and Saturn if the natural satellite lies deep inside the respective host Planet Hill sphere. Each planet has a Hill radius aH and planet mean radius RP and the ratio R1=RP/aH. Under very low R1 (less than 0.006) the approximation of CRTBP (centrally restricted three-body problem) to two-body problem is valid and planet has spacious Hill lobe to capture a satellite and retain it. This ensures a high probability of capture of natural satellite by the given planet and Sun’s perturbation on Planet-Satellite binary can be neglected. This is the case with Earth, Mars, Jupiter, Saturn, Neptune and Uranus. But Mercury and Venus has R1=RP/aH =0.01 and 5.9862 × 10-3 respectively hence they have no satellites. There is a limit to the dimension of the captured body. It must be a much smaller body both dimensionally as well masswise. The qantitative limit is a subject of an independent study.

Variational regularization method of solving the Cauchy problem for Laplace's equation: Innovation of the Grad–Shafranov(GS) reconstruction

The simplified linear model of Grad-Shafranov （GS） reconstruction can be reformulated into an inverse boundary value problem of Laplace＇s equation. Therefore, in this paper we focus on the method of solving the inverse boundary value problem of Laplace＇s equation. In the first place, the variational regularization method is used to deal with the ill- posedness of the Cauchy problem for Laplace＇s equation. Then, the ＇L-Curve＇ principle is suggested to be adopted in choosing the optimal regularization parameter. Finally, a numerical experiment is implemented with a section of Neumann and Dirichlet boundary conditions with observation errors. The results well converge to the exact solution of the problem, which proves the efficiency and robustness of the proposed method. When the order of observation error δ is 10-1, the order of the approximate result error can reach 10-3.

REGULARIZATION METHODS FOR NONLINEAR ILL-POSED PROBLEMS WITH ACCRETIVE OPERATORS

This article is devoted to the regularization of nonlinear ill-posed problems with accretive operators in Banach spaces. The data involved are assumed to be known approximately. The authors concentrate their discussion on the convergence rates of regular solutions.

THE REGULARIZED SOLUTION APPROXIMATION OF FORWARD/BACKWARD PROBLEMS FOR A FRACTIONAL PSEUDO-PARABOLIC EQUATION WITH RANDOM NOISE

This paper deals with the forward and backward problems for the nonlinear fractional pseudo-parabolic equation ut+(-Δ)^(s1)ut+β(-Δ)^(s2)u=F(u,x,t)subject o random Gaussian white noise for initial and final data.Under the suitable assumptions s1,s2andβ,we first show the ill-posedness of mild solutions for forward and backward problems in the sense of Hadamard,which are mainly driven by random noise.Moreover,we propose the Fourier truncation method for stabilizing the above ill-posed problems.We derive an error estimate between the exact solution and its regularized solution in an E‖·‖Hs22norm,and give some numerical examples illustrating the effect of above method.

A modified Tikhonov regularization method for a Cauchy problem of a time fractional diffusion equation

In this paper,we consider a Cauchy problem of the time fractional diffusion equation(TFDE)in x∈[0,L].This problem is ubiquitous in science and engineering applications.The illposedness of the Cauchy problem is explained by its solution in frequency domain.Furthermore,the problem is formulated into a minimization problem with a modified Tikhonov regularization method.The gradient of the regularization functional based on an adjoint problem is deduced and the standard conjugate gradient method is presented for solving the minimization problem.The error estimates for the regularized solutions are obtained under Hp norm priori bound assumptions.Finally,numerical examples illustrate the effectiveness of the proposed method.

New Regularization Algorithms for Solving the Deconvolution Problem in Well Test Data Interpretation

Two new regularization algorithms for solving the first-kind Volterra integral equation, which describes the pressure-rate deconvolution problem in well test data interpretation, are developed in this paper. The main features of the problem are the strong nonuniform scale of the solution and large errors (up to 15%) in the input data. In both algorithms, the solution is represented as decomposition on special basic functions, which satisfy given a priori information on solution, and this idea allow us significantly to improve the quality of approximate solution and simplify solving the minimization problem. The theoretical details of the algorithms, as well as the results of numerical experiments for proving robustness of the algorithms, are presented.

Comparisons between Isotropic and Anisotropic TV Regularizations in Inverse Acoustic Scattering

This article compares the isotropic and anisotropic TV regularizations used in inverse acoustic scattering. It is observed that compared with the traditional Tikhonov regularization, isotropic and anisotropic TV regularizations perform better in the sense of edge preserving. While anisotropic TV regularization will cause distortions along axes. To minimize the energy function with isotropic and anisotropic regularization terms, we use split Bregman scheme. We do several 2D numerical experiments to validate the above arguments.

A Preconditioned Fractional Tikhonov Regularization Method for Large Discrete Ill-posed Problems

The generalized Tikhonov regularization method is one of the most classical methods for the solution of linear systems of equations that arise from the discretization of linear ill-posed problems.However,the approximate solution obtained by the Tikhonov regularization method in general form may lack many details of the exact solution.Combining the fractional Tikhonov method with the preconditioned technique,and using the discrepancy principle for determining the regularization parameter,we present a preconditioned projected fractional Tikhonov regularization method for solving discrete ill-posed problems.Numerical experiments illustrate that the proposed algorithm has higher accuracy compared with the existing classical regularization methods.

Regularization Methods to Approximate Solutions of Variational Inequalities

In this paper, we study the regularization methods to approximate the solutions of the variational inequalities with monotone hemi-continuous operator having perturbed operators arbitrary. Detail, we shall study regularization methods to approximate solutions of following variational inequalities: and with operator A being monotone hemi-continuous form real Banach reflexive X into its dual space X*, but instead of knowing the exact data (y0, A), we only know its approximate data satisfying certain specified conditions and D is a nonempty convex closed subset of X;the real function f defined on X is assumed to be lower semi-continuous, convex and is not identical to infinity. At the same time, we will evaluate the convergence rate of the approximate solution. The regularization methods here are different from the previous ones.

Quasi-Reversibility Regularization Method for Solving a Backward Heat Conduction Problem

Non-standard backward heat conduction problem is ill-posed in the sense that the solution(if it exists) does not depend continuously on the data. In this paper, we propose a regularization strategy-quasi-reversibility method to analysis the stability of the problem. Meanwhile, we investigate the roles of regularization parameter in this method. Numerical result show that our algorithm is effective and stable.

A STRONG POSITIVITY PROPERTY AND A RELATED INVERSE SOURCE PROBLEM FOR MULTI-TERM TIME-FRACTIONAL DIFFUSION EQUATIONS

In this article,we consider the diffusion equation with multi-term time-fractional derivatives.We first derive,by a subordination principle for the solution,that the solution is positive when the initial value is non-negative.As an application,we prove the uniqueness of solution to an inverse problem of determination of the temporally varying source term by integral type information in a subdomain.Finally,several numerical experiments are presented to show the accuracy and efficiency of the algorithm.

A UNIFIED TREATMENT OF REGULARIZATION METHODS FOR LINEAR ILL-POSED PROBLEMS

By presenting a general framework, some regularization methods for solving linear ill-posed problems are considered in a unified manner. Applications to some specific approaches are illustrated.

WEAK REGULARIZATION FOR A CLASS OF ILL-POSED CAUCHY PROBLEMS

This article is concerned with the ill-posed Cauchy problem associated with a densely defined linear operator A in a Banach space. A family of weak regularizing operators is introduced. If the spectrum of A is contained in a sector of right-half complex plane and its resolvent is polynomially bounded, the weak regularization for such ill-posed Cauchy problem can be shown by using the quasi-reversibilky method and regularized semigroups. Finally, an example is given.

A Differential Continuation－Regularization Method for Solving Inverse Problems of Acoustic Wave Equations

ＡＤｉｆｆｅｒｅｎｔｉａｌＣｏｎｔｉｎｕａｔｉｏｎ－ＲｅｇｕｌａｒｉｚａｔｉｏｎＭｅｔｈｏｄｆｏｒＳｏｌｖｉｎｇＩｎｖｅｒｓｅＰｒｏｂｌｅｍｓｏｆＡｃｏｕｓｔｉｃＷａｖｅＥｑｕａｔｉｏｎｓ￥（张大力）（韩波）（刘家琦）（姜立功）ＺＨＡＮＧＤａｌｉ；Ｈ...

REGULARIZATION OF NEARLY SINGULAR INTEGRALS IN THE BOUNDARY ELEMENT METHOD OF POTENTIAL PROBLEMS

A general algorithm is applied to the regularization of nearly singular integrals in the boundary element method of planar potential problems. For linear elements, the strongly singular and hypersingular integrals of the interior points very close to boundary were categorized into two forms. The factor leading to the singularity was transformed out of the integral representations with integration by parts, so non-singular regularized formulas were presented for the two forms of integrals. Furthermore, quadratic elements are used in addition to linear ones. The quadratic element very close to the internal point can be divided into two linear ones, so that the algorithm is still valid. Numerical examples demonstrate the effectiveness and accuracy of this algorithm. Especially for problems with curved boundaries, the combination of quadratic elements and linear elements can give more accurate results.

REGULARIZATION OF AN ILL-POSED HYPERBOLIC PROBLEM AND THE UNIQUENESS OF THE SOLUTION BY A WAVELET GALERKIN METHOD

We consider the problem K（x）Uxx = utt , 0 〈 x 〈 1, t 〉 0, with the boundary condition u（O,t） = g（t） E LZ（R） and ux（O,t） = 0, where K（x） is continuous and 0 〈α≤ K （x） 〈 ＋∞. This is an ill-posed problem in the sense that, if the solution exists, it does not depend continuously on g. Considering the existence of a solution u（x, .） E H2（R） and using a wavelet Galerkin method with Meyer multiresolution analysis, we regularize the ill-posedness of the problem. Furthermore we prove the uniqueness of the solution for this problem.

A Newton-type method for nonlinear ill-posed problems with A-smooth regularization

In this paper we present a regularized Newton-type method for ill-posed problems, by using the A-smooth regularization to solve the linearized ill-posed equations. For noisy data a proper a posteriori stopping rule is used that yields convergence of the Newton iteration to a solution, as the noise level goes to zero, under certain smoothness conditions on the nonlinear operator. Some appropriate assumptions on the closedness and smoothness of the starting value and the solution are shown to lead to optimal convergence rates.

TWO REGULARIZATION METHODS FOR IDENTIFYING THE SOURCE TERM PROBLEM ON THE TIME-FRACTIONAL DIFFUSION EQUATION WITH A HYPER-BESSEL OPERATOR

In this paper,we consider the inverse problem for identifying the source term of the time-fractional equation with a hyper-Bessel operator.First,we prove that this inverse problem is ill-posed,and give the conditional stability.Then,we give the optimal error bound for this inverse problem.Next,we use the fractional Tikhonov regularization method and the fractional Landweber iterative regularization method to restore the stability of the ill-posed problem,and give corresponding error estimates under different regularization parameter selection rules.Finally,we verify the effectiveness of the method through numerical examples.