3D density inversion of gravity gradient data using the extrapolated Tikhonov regularization

We use the extrapolated Tikhonov regularization to deal with the ill-posed problem of 3D density inversion of gravity gradient data. The use of regularization parameters in the proposed method reduces the deviations between calculated and observed data. We also use the depth weighting function based on the eigenvector of gravity gradient tensor to eliminate undesired effects owing to the fast attenuation of the position function. Model data suggest that the extrapolated Tikhonov regularization in conjunction with the depth weighting function can effectively recover the 3D distribution of density anomalies. We conduct density inversion of gravity gradient data from the Australia Kauring test site and compare the inversion results with the published research results. The proposed inversion method can be used to obtain the 3D density distribution of underground anomalies.

Application of Tikhonov regularization method to wind retrieval from scatterometer data Ⅱ: cyclone wind retrieval with consideration of rain

According to the conclusion of the simulation experiments in paper I, the Tikhonov regularization method is applied to cyclone wind retrieval with a rain-effect-considering geophysical model function （called CMF＋Rain）. The CMF＋Rain model which is based on the NASA scatterometer-2 （NSCAT2） GMF is presented to compensate for the effects of rain on cyclone wind retrieval. With the multiple solution scheme （MSS）, the noise of wind retrieval is effectively suppressed, but the influence of the background increases. It will cause a large wind direction error in ambiguity removal when the background error is large. However, this can be mitigated by the new ambiguity removal method of Tikhonov regularization as proved in the simulation experiments. A case study on an extratropical cyclone of hurricane observed with SeaWinds at 25-km resolution shows that the retrieved wind speed for areas with rain is in better agreement with that derived from the best track analysis for the GMF＋Rain model, but the wind direction obtained with the two-dimensional variational （2DVAR） ambiguity removal is incorrect. The new method of Tikhonov regularization effectively improves the performance of wind direction ambiguity removal through choosing appropriate regularization parameters and the retrieved wind speed is almost the same as that obtained from the 2DVAR.

Some studies on the Tikhonov regularization method with additional assumptions for noise data

In this paper, the Tikhonov regularization method was used to solve the nondegenerate compact hnear operator equation, which is a well-known ill-posed problem. Apart from the usual error level, the noise data were supposed to satisfy some additional monotonic condition. Moreover, with the assumption that the singular values of operator have power form, the improved convergence rates of the regularized solution were worked out.

VIBRATING VELOCITY RECONSTRUC-TION USING IBEM AND TIKHONOV REGULARIZATION

The inverse problem to determine the vibrating velocity from known exteriorfield measurement pressure, involves the solution of a discrete ill-posed problem. To facilitate thecomputation of a meaningful approximate solution possible, the indirect boundary element method(IBEM) code for investigating vibration velocity reconstruction and Tikhonov regularization methodby means of singular value decomposition (SVD) are used. The amount of regularization is determinedby a regularization parameter. Its optimal value is given by the L-curve approach. Numerical resultsindicate the reconstructed normal surface velocity is a good approximation to the real source.

MULTI-PARAMETER TIKHONOV REGULARIZATION FOR LINEAR ILL-POSED OPERATOR EQUATIONS

We consider solving linear ill-posed operator equations. Based on a multi-scale decomposition for the solution space, we propose a multi-parameter regularization for solving the equations. We establish weak and strong convergence theorems for the multi-parameter regularization solution. In particular, based on the eigenfunction decomposition, we develop a posteriori choice strategy for multi-parameters which gives a regularization solution with the optimal error bound. Several practical choices of multi-parameters are proposed. We also present numerical experiments to demonstrate the outperformance of the multiparameter regularization over the single parameter regularization.

Investigation of the Tikhonov Regularization Method in Regional Gravity Field Modeling by Poisson Wavelets Radial Basis Functions

The application of Tikhonov regularization method dealing with the ill-conditioned problems in the regional gravity field modeling by Poisson wavelets is studied. In particular, the choices of the regularization matrices as well as the approaches for estimating the regularization parameters are investigated in details. The numerical results show that the regularized solutions derived from the first-order regularization are better than the ones obtained from zero-order regularization. For cross validation, the optimal regularization parameters are estimated from L-curve, variance component estimation（VCE） and minimum standard deviation（MSTD） approach, respectively, and the results show that the derived regularization parameters from different methods are consistent with each other. Together with the firstorder Tikhonov regularization and VCE method, the optimal network of Poisson wavelets is derived, based on which the local gravimetric geoid is computed. The accuracy of the corresponding gravimetric geoid reaches 1.1 cm in Netherlands, which validates the reliability of using Tikhonov regularization method in tackling the ill-conditioned problem for regional gravity field modeling.

Genetic Algorithm Based Tikhonov Regularization Method for Displacement Reconstruction

In this paper,a genetic algorithm based Tikhonov regularization method is proposed for determination of globally optimal regularization factor in displacement reconstruction.Optimization mathematic models are built by using the generalized cross-validation(GCV)criterion,L-curve criterion and Engl’s error minimization(EEM)criterion as the objective functions to prevent the regularization factor sinking into the locally optimal solution.The validity of the proposed algorithm is demonstrated through a numerical study of the frame structure model.Additionally,the influence of the noise level and the number of sampling points on the optimal regularization factor is analyzed.The results show that the proposed algorithm improves the robustness of the algorithm effectively,and reconstructs the displacement accurately.

Performance Improvement of Discrete-time Linear Control Systems Subject to Varying Sampling Rates Using the Tikhonov Regularization Method

Methods to stabilize discrete-time linear control systems subject to variable sampling rates,i.e.,using state feedback controllers,are well known in the literature.Several recent works address the use of the Tikhonov regularization method,originally designed to attenuate the noise effects on ill-posed problems,with the aim of improving performance and stabilizing approximately controllable dynamical systems.Inspired by these works,we propose the use of a feedback controller designed using the Tikhonov method to regularize discrete-time linear systems subject to varying sampling rates.The goal is to minimize an error function,thus improving the performance of the closed loop system and reducing the possibility of instability.Illustrative examples show the effectiveness of the proposed method.

DIRECT IMPLEMENTATION OF TIKHONOV REGULARIZATION FOR THE FIRST KIND INTEGRAL EQUATION

A common way to handle the Tikhonov regularization method for the first kind Fredholm integral equations,is first to discretize and then to work with the final linear system.This unavoidably inflicts discretization errors which may lead to disastrous results,especially when a quadrature rule is used.We propose to regularize directly the integral equation resulting in a continuous Tikhonov problem.The Tikhonov problem is reduced to a simple least squares problem by applying the Golub-Kahan bidiagonalization(GKB)directly to the integral operator.The regularization parameter and the iteration index are determined by the discrepancy principle approach.Moreover,we study the discrete version of the proposed method resulted from numerical evaluating the needed integrals.Focusing on the nodal values of the solution results in a weighted version of GKB-Tikhonov method for linear systems arisen from the Nystr¨om discretization.Finally,we use numerical experiments on a few test problems to illustrate the performance of our algorithms.

Multi-parameter Tikhonov Regularization—An Augmented Approach

We study multi-parameter regularization(multiple penalties) for solving linear inverse problems to promote simultaneously distinct features of the sought-for objects. We revisit a balancing principle for choosing regularization parameters from the viewpoint of augmented Tikhonov regularization, and derive a new parameter choice strategy called the balanced discrepancy principle. A priori and a posteriori error estimates are provided to theoretically justify the principles, and numerical algorithms for efficiently implementing the principles are also provided. Numerical results on deblurring are presented to illustrate the feasibility of the balanced discrepancy principle.

Extrapolated Tikhonov method and inversion of 3D density images of gravity data

Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method（EXTR） is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.

Downward continuation of airborne geomagnetic data based on two iterative regularization methods in the frequency domain

Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed to realize effective continuation. According to the Poisson integral plane approximate relationship between observation and continuation data, the computation formulae combined with the fast Fourier transform（FFT）algorithm are transformed to a frequency domain for accelerating the computational speed. The iterative Tikhonov regularization method and the iterative Landweber regularization method are used in this paper to overcome instability and improve the precision of the results. The availability of these two iterative regularization methods in the frequency domain is validated by simulated geomagnetic data, and the continuation results show good precision.

Efficient multiuser detector based on box-constrained dichotomous coordinate descent and regularization

The presented iterative multiuser detection technique was based on joint deregularized and box-constrained solution to quadratic optimization with iterations similar to that used in the nonstationary Tikhonov iterated algorithm.The deregularization maximized the energy of the solution,which was opposite to the Tikhonov regularization where the energy was minimized.However,combined with box-constraints,the deregularization forced the solution to be close to the binary set.It further exploited the box-constrained dichotomous coordinate descent algorithm and adapted it to the nonstationary iterative Tikhonov regularization to present an efficient detector.As a result,the worst-case and average complexity are reduced down as K2.8 and K2.5 floating point operation per second,respectively.The development improves the "efficient frontier" in multiuser detection,which is illustrated by simulation results.In addition,most operations in the detector are additions and bit-shifts.This makes the proposed technique attractive for fixed-point hardware implementation.

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.

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.

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.

Thrust estimate method of an on-orbit Hall thruster using Hall drift current

In order to realize the thrust estimation of the Hall thruster during its flight mission,this study establishes an estimation method based on measurement of the Hall drift current.In this method,the Hall drift current is calculated from an inverse magnetostatic problem,which is formulated according to its induced magnetic flux density detected by sensors,and then the thrust is estimated by multiplying the Hall drift current with the characteristic magnetic flux density of the thruster itself.In addition,a three-wire torsion pendulum micro-thrust measurement system is utilized to verify the estimate values obtained from the proposed method.The errors were found to be less than 8%when the discharge voltage ranged from 250 V to 350 V and the anode flow rate ranged from 30 sccm to 50 sccm,indicating the possibility that the proposed thrust estimate method could be practically applied.Moreover,the measurement accuracy of the magnetic flux density is suggested to be lower than 0.015 mT and improvement on the inverse problem solution is required in the future.

Measuring internal residual stress in Al-Cu alloy forgings by crack compliance method with optimized parameters

Measuring the internal stress of Al alloy forgings accurately is critical for controlling the deformation during the subsequent machine process.In this work,the crack compliance method was used to calculate the internal residual stress of Al-Cu high strength alloys,and the effect of various model parameters of crack compliance method on the calculated precision was studied by combining the numerical simulation and experimental method.The results show that the precision first increased and then decreased with increasing the crack range.The decreased precision when using a high crack range was due to the strain fluctuation during the machining process,and the optimized crack range was 71%of the thickness of forgings.Low orders of Legendre polynomial can result in residual stress curve more smooth,while high orders led to the occurrence of distortion.The Tikhonov regularization method effectively suppressed the distortion of residual stress caused by the fluctuation of strain data,which significantly improved the precision.In addition,The crack compliance method with optimized parameters was used to measure the residual stress of Al-Cu alloy with different quenching methods.The calculated results demonstrated that the distribution of residual stress was obtained accurately.

Inversion and III-Posed Problem Solutions in Atmospheric Remote Sensing

With the swift advances in earth observation,satellite remote sensing and application of atmospheric radiation theory have been developed in the past decades,atmospheric sensing inversion with its algorithms is getting more and more importance.It is known that since a remote sensing equation falls into an integral equation of the first kind,thus leading to the fact that it is ill-posed and particularly the solution is unsteady,tremendous difficulties arise from the retrieval.This paper will present a simple review on the inversion techniques with some necessary remarks,before introducing the successful efforts with respect to such equations and the encouraging solutions achieved in recent decades by researchers of the world.

IDENTIFYING AN UNKNOWN SOURCE IN SPACE-FRACTIONAL DIFFUSION EQUATION

In this paper, we identify a space-dependent source for a fractional diffusion equation. This problem is ill-posed, i.e., the solution （if it exists） does not depend continuously on the data. The generalized Tikhonov regularization method is proposed to solve this problem. An a priori error estimate between the exact solution and its regularized approximation is obtained. Moreover, an a posteriori parameter choice rule is proposed and a stable error estimate is also obtained, Numerical examples are presented to illustrate the validity and effectiveness of this method.