In order to decrease the sensitivity of the constant scale parameter, adaptively optimize the scale parameter in the iteration regularization model （IRM） and attain a desirable level of applicability for image denoi...In order to decrease the sensitivity of the constant scale parameter, adaptively optimize the scale parameter in the iteration regularization model （IRM） and attain a desirable level of applicability for image denoising, a novel IRM with the adaptive scale parameter is proposed. First, the classic regularization item is modified and the equation of the adaptive scale parameter is deduced. Then, the initial value of the varying scale parameter is obtained by the trend of the number of iterations and the scale parameter sequence vectors. Finally, the novel iterative regularization method is used for image denoising. Numerical experiments show that compared with the IRM with the constant scale parameter, the proposed method with the varying scale parameter can not only reduce the number of iterations when the scale parameter becomes smaller, but also efficiently remove noise when the scale parameter becomes bigger and well preserve the details of images.展开更多
In this paper, a modified Newton type iterative method is considered for ap- proximately solving ill-posed nonlinear operator equations involving m-accretive mappings in Banach space. Convergence rate of the method is...In this paper, a modified Newton type iterative method is considered for ap- proximately solving ill-posed nonlinear operator equations involving m-accretive mappings in Banach space. Convergence rate of the method is obtained based on an a priori choice of the regularization parameter. Our analysis is not based on the sequential continuity of the normalized duality mapping.展开更多
Recently,inverse problems have attracted more and more attention in computational mathematics and become increasingly important in engineering applications.After the discretization,many of inverse problems are reduced...Recently,inverse problems have attracted more and more attention in computational mathematics and become increasingly important in engineering applications.After the discretization,many of inverse problems are reduced to linear systems.Due to the typical ill-posedness of inverse problems,the reduced linear systems are often illposed,especially when their scales are large.This brings great computational difficulty.Particularly,a small perturbation in the right side of an ill-posed linear system may cause a dramatical change in the solution.Therefore,regularization methods should be adopted for stable solutions.In this paper,a new class of accelerated iterative regularization methods is applied to solve this kind of large-scale ill-posed linear systems.An iterative scheme becomes a regularization method only when the iteration is early terminated.And a Morozov’s discrepancy principle is applied for the stop criterion.Compared with the conventional Landweber iteration,the new methods have acceleration effect,and can be compared to the well-known acceleratedν-method and Nesterov method.From the numerical results,it is observed that using appropriate discretization schemes,the proposed methods even have better behavior when comparing withν-method and Nesterov method.展开更多
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 ...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.展开更多
Linear Least Squares（LLS） problems are particularly difficult to solve because they are frequently ill-conditioned, and involve large quantities of data. Ill-conditioned LLS problems are commonly seen in mathematics...Linear Least Squares（LLS） problems are particularly difficult to solve because they are frequently ill-conditioned, and involve large quantities of data. Ill-conditioned LLS problems are commonly seen in mathematics and geosciences, where regularization algorithms are employed to seek optimal solutions. For many problems, even with the use of regularization algorithms it may be impossible to obtain an accurate solution. Riley and Golub suggested an iterative scheme for solving LLS problems. For the early iteration algorithm, it is difficult to improve the well-conditioned perturbed matrix and accelerate the convergence at the same time. Aiming at this problem, self-adaptive iteration algorithm（SAIA） is proposed in this paper for solving severe ill-conditioned LLS problems. The algorithm is different from other popular algorithms proposed in recent references. It avoids matrix inverse by using Cholesky decomposition, and tunes the perturbation parameter according to the rate of residual error decline in the iterative process. Example shows that the algorithm can greatly reduce iteration times, accelerate the convergence,and also greatly enhance the computation accuracy.展开更多
The loping OS-EM iteration is a numerically efficient regularization method for solving ill-posed problems. In this article we investigate the loping OS-EM iterative method in connection with the circular Radon transf...The loping OS-EM iteration is a numerically efficient regularization method for solving ill-posed problems. In this article we investigate the loping OS-EM iterative method in connection with the circular Radon transform. We show that the proposed method converges weakly for the noisy data. Numerical tests are presented for a linear problem related to photoacoustic tomography.展开更多
This article presents a fast convergent method of iterated regularization based on the idea of Landweber iterated regularization, and a method for a-posteriori choice by the Morozov discrepancy principle and the optim...This article presents a fast convergent method of iterated regularization based on the idea of Landweber iterated regularization, and a method for a-posteriori choice by the Morozov discrepancy principle and the optimum asymptotic convergence order of the regularized solution is obtained. Numerical test shows that the method of iterated regularization can quicken the convergence speed and reduce the calculation burden efficiently.展开更多
Cement density monitoring plays a vital role in evaluating the quality of cementing projects,which is of great significance to the development of oil and gas.However,the presence of inhomogeneous cement distribution a...Cement density monitoring plays a vital role in evaluating the quality of cementing projects,which is of great significance to the development of oil and gas.However,the presence of inhomogeneous cement distribution and casing eccentricity in horizontal wells often complicates the accurate evaluation of cement azimuthal density.In this regard,this paper proposes an algorithm to calculate the cement azimuthal density in horizontal wells using a multi-detector gamma-ray detection system.The spatial dynamic response functions are simulated to obtain the influence of cement density on gamma-ray counts by the perturbation theory,and the contribution of cement density in six sectors to the gamma-ray recorded by different detectors is obtained by integrating the spatial dynamic response functions.Combined with the relationship between gamma-ray counts and cement density,a multi-parameter calculation equation system is established,and the regularized Newton iteration method is employed to invert casing eccentricity and cement azimuthal density.This approach ensures the stability of the inversion process while simultaneously achieving an accuracy of 0.05 g/cm^(3) for the cement azimuthal density.This accuracy level is ten times higher compared to density accuracy calculated using calibration equations.Overall,this algorithm enhances the accuracy of cement azimuthal density evaluation,provides valuable technical support for the monitoring of cement azimuthal density in the oil and gas industry.展开更多
In this paper, we analyze the Bregman iterative model using the G-norm. Firstly, we show the convergence of the iterative model. Secondly, using the source condition and the symmetric Bregman distance, we consider the...In this paper, we analyze the Bregman iterative model using the G-norm. Firstly, we show the convergence of the iterative model. Secondly, using the source condition and the symmetric Bregman distance, we consider the error estimations between the iterates and the exact image both in the case of clean and noisy data. The results show that the Bregman iterative model using the G-norm has the similar good properties as the Bregman iterative model using the L2-norm.展开更多
In previous work, we have analyzed the feasibility of the estimation for a source term S（x, y, z） in a transversal section, The present study is concerned with a twodimensional inverse phase change problem. The goal...In previous work, we have analyzed the feasibility of the estimation for a source term S（x, y, z） in a transversal section, The present study is concerned with a twodimensional inverse phase change problem. The goal is the estimation of the dissipated heat flux in the liquid zone （reconstruction of a source term in the energy equation） from experimentally measured temperatures in the solid zone. This work has an application in the electron beam welding of steels of thickness about 8cm. The direct thermometallurgical problem is treated in a quasi steady two-dimensional longitudinal section （x, y）. The beam displacement is normally in the y direction. But in the quasi steady simulation, the beam is steady in the study section. The sample is divided in the axial direction z in few sections. At each section, a source term is defined with a part of the beam and creates a vaporized zone and a fused zone. The goal of this work is the rebuilding of the complete source term with the estimations at each section. In this paper, we analyze the feasibility of the estimation. For this work, we use only the simulated measurements without noise.展开更多
In this paper, by using new analysis techniques, we have studied iterative construc- tion problem for finding zeros of accretive mappings in uniformly smooth Banach spaces, and improved a theorem due to Reich. As its ...In this paper, by using new analysis techniques, we have studied iterative construc- tion problem for finding zeros of accretive mappings in uniformly smooth Banach spaces, and improved a theorem due to Reich. As its application, we have deduced a strong convergence theorem of fixed points for continuous pseudo-contractions.展开更多
基金The National Natural Science Foundation of China(No.60702069)the Research Project of Department of Education of Zhe-jiang Province (No.20060601)+1 种基金the Natural Science Foundation of Zhe-jiang Province (No.Y1080851)Shanghai International Cooperation onRegion of France (No.06SR07109)
文摘In order to decrease the sensitivity of the constant scale parameter, adaptively optimize the scale parameter in the iteration regularization model （IRM） and attain a desirable level of applicability for image denoising, a novel IRM with the adaptive scale parameter is proposed. First, the classic regularization item is modified and the equation of the adaptive scale parameter is deduced. Then, the initial value of the varying scale parameter is obtained by the trend of the number of iterations and the scale parameter sequence vectors. Finally, the novel iterative regularization method is used for image denoising. Numerical experiments show that compared with the IRM with the constant scale parameter, the proposed method with the varying scale parameter can not only reduce the number of iterations when the scale parameter becomes smaller, but also efficiently remove noise when the scale parameter becomes bigger and well preserve the details of images.
文摘In this paper, a modified Newton type iterative method is considered for ap- proximately solving ill-posed nonlinear operator equations involving m-accretive mappings in Banach space. Convergence rate of the method is obtained based on an a priori choice of the regularization parameter. Our analysis is not based on the sequential continuity of the normalized duality mapping.
基金supported by the Natural Science Foundation of China (Nos. 11971230, 12071215)the Fundamental Research Funds for the Central Universities(No. NS2018047)the 2019 Graduate Innovation Base(Laboratory)Open Fund of Jiangsu Province(No. Kfjj20190804)
文摘Recently,inverse problems have attracted more and more attention in computational mathematics and become increasingly important in engineering applications.After the discretization,many of inverse problems are reduced to linear systems.Due to the typical ill-posedness of inverse problems,the reduced linear systems are often illposed,especially when their scales are large.This brings great computational difficulty.Particularly,a small perturbation in the right side of an ill-posed linear system may cause a dramatical change in the solution.Therefore,regularization methods should be adopted for stable solutions.In this paper,a new class of accelerated iterative regularization methods is applied to solve this kind of large-scale ill-posed linear systems.An iterative scheme becomes a regularization method only when the iteration is early terminated.And a Morozov’s discrepancy principle is applied for the stop criterion.Compared with the conventional Landweber iteration,the new methods have acceleration effect,and can be compared to the well-known acceleratedν-method and Nesterov method.From the numerical results,it is observed that using appropriate discretization schemes,the proposed methods even have better behavior when comparing withν-method and Nesterov method.
基金supported by the National Natural Science Foundation of China(41304022,41174026,41104047)the National 973 Foundation(61322201,2013CB733303)+1 种基金the Key laboratory Foundation of Geo-space Environment and Geodesy of the Ministry of Education(13-01-08)the Youth Innovation Foundation of High Resolution Earth Observation(GFZX04060103-5-12)
文摘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.
基金supported by Open Fund of Engineering Laboratory of Spatial Information Technology of Highway Geological Disaster Early Warning in Hunan Province(Changsha University of Science&Technology,kfj150602)Hunan Province Science and Technology Program Funded Projects,China(2015NK3035)+1 种基金the Land and Resources Department Scientific Research Project of Hunan Province,China(2013-27)the Education Department Scientific Research Project of Hunan Province,China(13C1011)
文摘Linear Least Squares（LLS） problems are particularly difficult to solve because they are frequently ill-conditioned, and involve large quantities of data. Ill-conditioned LLS problems are commonly seen in mathematics and geosciences, where regularization algorithms are employed to seek optimal solutions. For many problems, even with the use of regularization algorithms it may be impossible to obtain an accurate solution. Riley and Golub suggested an iterative scheme for solving LLS problems. For the early iteration algorithm, it is difficult to improve the well-conditioned perturbed matrix and accelerate the convergence at the same time. Aiming at this problem, self-adaptive iteration algorithm（SAIA） is proposed in this paper for solving severe ill-conditioned LLS problems. The algorithm is different from other popular algorithms proposed in recent references. It avoids matrix inverse by using Cholesky decomposition, and tunes the perturbation parameter according to the rate of residual error decline in the iterative process. Example shows that the algorithm can greatly reduce iteration times, accelerate the convergence,and also greatly enhance the computation accuracy.
基金supported by the National Natural Science Foundation of China(61271398)K.C.Wong Magna Fund in Ningbo UniversityNatural Science Foundation of Ningbo City(2010A610102)
文摘The loping OS-EM iteration is a numerically efficient regularization method for solving ill-posed problems. In this article we investigate the loping OS-EM iterative method in connection with the circular Radon transform. We show that the proposed method converges weakly for the noisy data. Numerical tests are presented for a linear problem related to photoacoustic tomography.
文摘This article presents a fast convergent method of iterated regularization based on the idea of Landweber iterated regularization, and a method for a-posteriori choice by the Morozov discrepancy principle and the optimum asymptotic convergence order of the regularized solution is obtained. Numerical test shows that the method of iterated regularization can quicken the convergence speed and reduce the calculation burden efficiently.
基金The authors would like to acknowledge the support of the National Natural Science Foundation of China(41974127,42174147).References。
文摘Cement density monitoring plays a vital role in evaluating the quality of cementing projects,which is of great significance to the development of oil and gas.However,the presence of inhomogeneous cement distribution and casing eccentricity in horizontal wells often complicates the accurate evaluation of cement azimuthal density.In this regard,this paper proposes an algorithm to calculate the cement azimuthal density in horizontal wells using a multi-detector gamma-ray detection system.The spatial dynamic response functions are simulated to obtain the influence of cement density on gamma-ray counts by the perturbation theory,and the contribution of cement density in six sectors to the gamma-ray recorded by different detectors is obtained by integrating the spatial dynamic response functions.Combined with the relationship between gamma-ray counts and cement density,a multi-parameter calculation equation system is established,and the regularized Newton iteration method is employed to invert casing eccentricity and cement azimuthal density.This approach ensures the stability of the inversion process while simultaneously achieving an accuracy of 0.05 g/cm^(3) for the cement azimuthal density.This accuracy level is ten times higher compared to density accuracy calculated using calibration equations.Overall,this algorithm enhances the accuracy of cement azimuthal density evaluation,provides valuable technical support for the monitoring of cement azimuthal density in the oil and gas industry.
基金Supported by the National Natural Science Foundation of China(No.10726035,10771065,10801049)Foundation of North China Electric Power University
文摘In this paper, we analyze the Bregman iterative model using the G-norm. Firstly, we show the convergence of the iterative model. Secondly, using the source condition and the symmetric Bregman distance, we consider the error estimations between the iterates and the exact image both in the case of clean and noisy data. The results show that the Bregman iterative model using the G-norm has the similar good properties as the Bregman iterative model using the L2-norm.
文摘In previous work, we have analyzed the feasibility of the estimation for a source term S（x, y, z） in a transversal section, The present study is concerned with a twodimensional inverse phase change problem. The goal is the estimation of the dissipated heat flux in the liquid zone （reconstruction of a source term in the energy equation） from experimentally measured temperatures in the solid zone. This work has an application in the electron beam welding of steels of thickness about 8cm. The direct thermometallurgical problem is treated in a quasi steady two-dimensional longitudinal section （x, y）. The beam displacement is normally in the y direction. But in the quasi steady simulation, the beam is steady in the study section. The sample is divided in the axial direction z in few sections. At each section, a source term is defined with a part of the beam and creates a vaporized zone and a fused zone. The goal of this work is the rebuilding of the complete source term with the estimations at each section. In this paper, we analyze the feasibility of the estimation. For this work, we use only the simulated measurements without noise.
基金the National Natural Science Foundation of China (No.10471033)
文摘In this paper, by using new analysis techniques, we have studied iterative construc- tion problem for finding zeros of accretive mappings in uniformly smooth Banach spaces, and improved a theorem due to Reich. As its application, we have deduced a strong convergence theorem of fixed points for continuous pseudo-contractions.
