A Multi-Baseline PolInSAR Forest Height Inversion Method Taking into Account the Model Ill-posed Problem

Affected by the insufficient information of single baseline observation data,the three-stage method assumes the Ground-to-Volume Ratio(GVR)to be zero so as to invert the vegetation height.However,this assumption introduces much biases into the parameter estimates which greatly limits the accuracy of the vegetation height inversion.Multi-baseline observation can provide redundant information and is helpful for the inversion of GVR.Nevertheless,the similar model parameter values in a multi-baseline model often lead to ill-posed problems and reduce the inversion accuracy of conventional algorithm.To this end,we propose a new step-by-step inversion method applied to the multi-baseline observations.Firstly,an adjustment inversion model is constructed by using multi-baseline volume scattering dominant polarization data,and the regularized estimates of model parameters are obtained by regularization method.Then,the reliable estimates of GVR are determined by the MSE(mean square error)analysis of each regularized parameter estimation.Secondly,the estimated GVR is used to extracts the pure volume coherence,and then the vegetation height parameter is inverted from the pure volume coherence by least squares estimation.The experimental results show that the new method can improve the vegetation height inversion result effectively.The inversion accuracy is improved by 26%with respect to the three-stage method and the conventional solution of multi-baseline.All of these have demonstrated the feasibility and effectiveness of the new method.

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.

A mixed Newton-Tikhonov method for nonlinear ill-posed problems

Newton type methods are one kind of the efficient methods to solve nonlinear ill-posed problems, which have attracted extensive attention. However, computational cost of Newton type methods is high because practical problems are complicated. We propose a mixed Newton-Tikhonov method, i.e., one step Newton-Tikhonov method with several other steps of simplified Newton-Tikhonov method. Convergence and stability of this method are proved under some conditions. Numerical experiments show that the proposed method has obvious advantages over the classical Newton method in terms of computational costs.

Arnoldi Projection Fractional Tikhonov for Large Scale Ill-Posed Problems

It is well known that Tikhonov regularization in standard form may determine approximate solutions that are too smooth for ill-posed problems,so fractional Tikhonov methods have been introduced to remedy this shortcoming.And Tikhonov regularization for large-scale linear ill-posed problems is commonly implemented by determining apartial Arnoldi decomposition of the given matrix.In this paper,we propose a new method to compute an approximate solution of large scale linear discrete ill-posed problems which applies projection fractional Tikhonov regularization in Krylov subspace via Arnoldi process.The projection fractional Tikhonov regularization combines the fractional matrices and orthogonal projection operators.A suitable value of the regularization parameter is determined by the discrepancy principle.Numerical examples with application to image restoration are carried out to examine that the performance of the method.

Nonlinear implicit iterative method for solving nonlinear ill-posed problems

In the paper, we extend the implicit iterative method for linear ill-posed operator equations to solve nonlinear ill-posed problems. We show that under some conditions the error sequence of solutions of the nonlinear implicit iterative method is monotonically decreasing and, with this monotonicity, prove convergence of the new method for both the exact and perturbed equations.

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.

A Regularized Randomized Kaczmarz Algorithm for Large Discrete Ill-Posed Problems

Tikhonov regularization is a powerful tool for solving linear discrete ill-posed problems.However,effective methods for dealing with large-scale ill-posed problems are still lacking.The Kaczmarz method is an effective iterative projection algorithm for solving large linear equations due to its simplicity.We propose a regularized randomized extended Kaczmarz(RREK)algorithm for solving large discrete ill-posed problems via combining the Tikhonov regularization and the randomized Kaczmarz method.The convergence of the algorithm is proved.Numerical experiments illustrate that the proposed algorithm has higher accuracy and better image restoration quality compared with the existing randomized extended Kaczmarz(REK)method.

Fast Multilevel Methods for Solving Ill-posed Problems

1 Introduetion Many industrial and engineering applieations require numerieally solving ill-posed Problems. Regularization methods are employed to find approximate solutions of these problems.The choice of regularization

On the Regularization Method for Solving Ill-Posed Problems with Unbounded Operators

Let be a linear, closed, and densely defined unbounded operator, where X and Y are Hilbert spaces. Assume that A is not boundedly invertible. Suppose the equation Au=f is solvable, and instead of knowing exactly f only know its approximation satisfies the condition: In this paper, we are interested a regularization method to solve the approximation solution of this equation. This approximation is a unique global minimizer of the functional , for any , defined as follows: . We also study the stability of this method when the regularization parameter is selected a priori and a posteriori. At the same time, we give an application of this method to the weak derivative operator equation in Hilbert space.

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.

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.

A New Regularized Solution to Ill-Posed Problem in Coordinate Transformation

Coordinates transformation is generally required in GPS applications. If the transformation parameters are solved with the known coordinates in a small area using the Bursa model, the precision of transformed coordinates is generally very poor. Since the translation parameters and rotation parameters are highly correlated in this case, a very large condition number of the coefficient matrix A exists in the linear system of equations. Regularization is required to reduce the effects caused by the intrinsic ill-conditioning of the problem and noises in the data, and to stabilize the solution. Based on advanced regularized methods, we propose a new regularized solution to the ill-posed coordinate transformation problem. Simulation numerical experiments of coordinate transformation are given to shed light on the relationship among different regularization approaches. The results indicate that the proposed new method can obtain a more reasonable resolution with higher precision and/or robustness.

MULTILEVEL ITERATION METHODS FOR SOLVING LINEAR ILL-POSED PROBLEMS

In this paper we develop multilevel iteration methods for solving linear systems resulting from the Galerkin method and Tikhonov regularization for ill-posed problems. The algorithm and its convergence analysis are presented in an abstract framework.

Optimal Shape Factor and Fictitious Radius in the MQ-RBF:Solving Ill-Posed Laplacian Problems

To solve the Laplacian problems,we adopt a meshless method with the multiquadric radial basis function(MQRBF)as a basis whose center is distributed inside a circle with a fictitious radius.A maximal projection technique is developed to identify the optimal shape factor and fictitious radius by minimizing a merit function.A sample function is interpolated by theMQ-RBF to provide a trial coefficient vector to compute the merit function.We can quickly determine the optimal values of the parameters within a preferred rage using the golden section search algorithm.The novel method provides the optimal values of parameters and,hence,an optimal MQ-RBF;the performance of the method is validated in numerical examples.Moreover,nonharmonic problems are transformed to the Poisson equation endowed with a homogeneous boundary condition;this can overcome the problem of these problems being ill-posed.The optimal MQ-RBF is extremely accurate.We further propose a novel optimal polynomial method to solve the nonharmonic problems,which achieves high precision up to an order of 10^(−11).

Frozen Landweber Iteration for Nonlinear Ill-Posed Problems

In this paper we propose a modification of the Landweber iteration termed frozen Landweber iteration for nonlinear ill-posed problems. A convergence analysis for this iteration is presented. The numerical performance of this frozen Landweber iteration for a nonlinear Hammerstein integral equation is compared with that of the Landweber iteration. We obtain a shorter running time of the frozen Landweber iteration based on the same convergence accuracy.

Regularization method with two parameters for nonlinear ill-posed problems

This paper is devoted to the regularization of a class of nonlinear ill-posed problems in Banach spaces. The operators involved are multi-valued and the data are assumed to be known approximately. Under the assumption that the original problem is solvable, a strongly convergent approximation procedure is designed by means of the Tikhonov regularization method with two pa- rameters.

A Restarted Conjugate Gradient Method for Ill-posed Problems

Abstract This paper presents a restarted conjugate gradient iterative algorithm for solving ill-posed problems. The damped Morozov's discrepancy principle is used as a stopping rule. Numerical experiments are given to illustrate the efficiency of the method.

Some results on the regularization of LSQR for large-scale discrete ill-posed problems

LSQR, a Lanczos bidiagonalization based Krylov subspace iterative method, and its mathematically equivalent conjugate gradient for least squares problems(CGLS) applied to normal equations system, are commonly used for large-scale discrete ill-posed problems. It is well known that LSQR and CGLS have regularizing effects, where the number of iterations plays the role of the regularization parameter. However, it has long been unknown whether the regularizing effects are good enough to find best possible regularized solutions. Here a best possible regularized solution means that it is at least as accurate as the best regularized solution obtained by the truncated singular value decomposition(TSVD) method. We establish bounds for the distance between the k-dimensional Krylov subspace and the k-dimensional dominant right singular space. They show that the Krylov subspace captures the dominant right singular space better for severely and moderately ill-posed problems than for mildly ill-posed problems. Our general conclusions are that LSQR has better regularizing effects for the first two kinds of problems than for the third kind, and a hybrid LSQR with additional regularization is generally needed for mildly ill-posed problems. Exploiting the established bounds, we derive an estimate for the accuracy of the rank k approximation generated by Lanczos bidiagonalization. Numerical experiments illustrate that the regularizing effects of LSQR are good enough to compute best possible regularized solutions for severely and moderately ill-posed problems, stronger than our theory, but they are not for mildly ill-posed problems and additional regularization is needed.

An iterative virtual projection method to improve the reconstruction performance for ill-posed emission tomographic problems

In order to improve the reconstruction performance for ill-posed emission tomographic problems with limited projections, a generalized interpolation method is proposed in this paper, in which the virtual lines of projection are fabricated from, but not linearly dependent on, the measured projections. The method is called the virtual projection（VP） method.Also, an iterative correction method for the integral lengths is proposed to reduce the error brought about by the virtual lines of projection. The combination of the two methods is called the iterative virtual projection（IVP） method. Based on a scheme of equilateral triangle plane meshes and a six asymmetrically arranged detection system, numerical simulations and experimental verification are conducted. Simulation results obtained by using a non-negative linear least squares method,without any other constraints or regularization, demonstrate that the VP method can gradually reduce the reconstruction error and converges to the desired one by fabricating additional effective projections. When the mean square deviation of normal error superimposed on the simulated measured projections is smaller than 0.03, i.e., the signal-to-noise ratio（SNR）for the measured projections is higher than 30.4, the IVP method can further reduce the reconstruction error reached by the VP method apparently. In addition, as the regularization matrix in the Tikhonov regularization method is updated by an iterative correction process similar to the IVP method presented in this paper, or the Tikhonov regularization method is used in the IVP method, good improvement is achieved.

STRUCTURED CONDITION NUMBERS FOR THE TIKHONOV REGULARIZATION OF DISCRETE ILL-POSED PROBLEMS

The possibly most popular regularization method for solving the least squares problem rain ‖Ax - b‖2 with a highly ill-conditioned or rank deficient coefficient matrix A is the x Tikhonov regularization method. In this paper we present the explicit expressions of the normwise, mixed and componentwise condition numbers for the Tikhonov regularization when A has linear structures. The structured condition numbers in the special cases of nonlinear structure i.e. Vandermonde and Cauchy matrices are also considered. Some comparisons between structured condition numbers and unstructured condition numbers are made by numerical experiments. In addition, we also derive the normwise, mixed and componentwise condition numbers for the Tikhonov regularization when the coefficient matrix, regularization matrix and right-hand side vector are all perturbed, which generalize the results obtained by Chu et al. [Numer. Linear Algebra Appl., 18 （2011）, 87-103].