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 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.

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.

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.

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).

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.

Dirac method for nonlinear and non-homogenous boundary value problems of plates

The boundary value problem plays a crucial role in the analytical investigation of continuum dynamics. In this paper, an analytical method based on the Dirac operator to solve the nonlinear and non-homogeneous boundary value problem of rectangular plates is proposed. The key concept behind this method is to transform the nonlinear or non-homogeneous part on the boundary into a lateral force within the governing function by the Dirac operator, which linearizes and homogenizes the original boundary, allowing one to employ the modal superposition method for obtaining solutions to reconstructive governing equations. Once projected into the modal space, the harmonic balance method(HBM) is utilized to solve coupled ordinary differential equations(ODEs)of truncated systems with nonlinearity. To validate the convergence and accuracy of the proposed Dirac method, the results of typical examples, involving nonlinearly restricted boundaries, moment excitation, and displacement excitation, are compared with those of the differential quadrature element method(DQEM). The results demonstrate that when dealing with nonlinear boundaries, the Dirac method exhibits more excellent accuracy and convergence compared with the DQEM. However, when facing displacement excitation, there exist some discrepancies between the proposed approach and simulations;nevertheless, the proposed method still accurately predicts resonant frequencies while being uniquely capable of handling nonuniform displacement excitations. Overall, this methodology offers a convenient way for addressing nonlinear and non-homogenous plate boundaries.

Highly Accurate Golden Section Search Algorithms and Fictitious Time Integration Method for Solving Nonlinear Eigenvalue Problems

This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.

Numerical Exploration of Asymmetrical Impact Dynamics: Unveiling Nonlinearities in Collision Problems and Resilience of Reinforced Concrete Structures

This study provides a comprehensive analysis of collision and impact problems’ numerical solutions, focusing ongeometric, contact, and material nonlinearities, all essential in solving large deformation problems during a collision.The initial discussion revolves around the stress and strain of large deformation during a collision, followedby explanations of the fundamental finite element solution method for addressing such issues. The hourglassmode’s control methods, such as single-point reduced integration and contact-collision algorithms are detailedand implemented within the finite element framework. The paper further investigates the dynamic responseand failure modes of Reinforced Concrete (RC) members under asymmetrical impact using a 3D discrete modelin ABAQUS that treats steel bars and concrete connections as bond slips. The model’s validity was confirmedthrough comparisons with the node-sharing algorithm and system energy relations. Experimental parameterswere varied, including the rigid hammer’s mass and initial velocity, concrete strength, and longitudinal and stirrupreinforcement ratios. Findings indicated that increased hammer mass and velocity escalated RC member damage,while increased reinforcement ratios improved impact resistance. Contrarily, increased concrete strength did notsignificantly reduce lateral displacement when considering strain rate effects. The study also explores materialnonlinearity, examining different materials’ responses to collision-induced forces and stresses, demonstratedthrough an elastic rod impact case study. The paper proposes a damage criterion based on the residual axialload-bearing capacity for assessing damage under the asymmetrical impact, showing a correlation betweendamage degree hammer mass and initial velocity. The results, validated through comparison with theoreticaland analytical solutions, verify the ABAQUS program’s accuracy and reliability in analyzing impact problems,offering valuable insights into collision and impact problems’ nonlinearities and practical strategies for enhancingRC structures’ resilience under dynamic stress.

A FINITE DIMENSIONAL METHOD FOR SOLVINGNONLINEAR ILL-POSED PROBLEMS

We propose a finite dimensional method to compute the solution of nonlinear ill-posed problems approximately and show that under certain conditions, the convergence can be guaranteed. Moreover, we obtain the rate of convergence of our method provided that the true solution satisfies suitable smoothness condition. Finally, we present two examples from the parameter estimation problems of differential equations and illustrate the applicability of our method.

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.

Maximum entropy method for solving nonlinear ill-posed problems

We are specially interested in the case that problem (1) is ill-posed; that is, the solutions of (1) do not depend continuously on the data. Now the regularization techniques are required. The traditional method is Tikhonov regularization. In recent years, the concept of entropy was introduced into the study of ill-posed problems and developed the maximum entropy method. It is found that the maximum entropy method has its

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.

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.

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.

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 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 multi-objective optimization framework for ill-posed inverse problems

Many image inverse problems are ill-posed for no unique solutions. Most of them have incommensurable or mixed-type objectives. In this study, a multi-objective optimization framework is introduced to model such ill-posed inverse problems. The conflicting objectives are designed according to the properties of ill-posedness and certain techniques. Multi-objective evolutionary algorithms have capability to optimize multiple objectives simultaneously and obtain a set of trade-off solutions. For that reason, we use multi-objective evolutionary algorithms to keep the trade-off between these objectives for image ill-posed problems. Two case studies of sparse reconstruction and change detection are imple- mented. In the case study of sparse reconstruction, the measurement error term and the sparsity term are optimized by multi-objective evolutionary algorithms, which aims at balancing the trade-off between enforcing sparsity and reducing measurement error. In the case study of image change detection, two conflicting objectives are constructed to keep the trade-off between robustness to noise and preserving the image details. Experimental results of the two case studies confirm the multi-objective optimization framework for ill-posed inverse problems in image processing is effective.