Energy resolution calibration is crucial for gamma-ray spectral analysis,as measured using a scintillation detector.A locally constrained regularization method was proposed to determine the resolution calibration para...Energy resolution calibration is crucial for gamma-ray spectral analysis,as measured using a scintillation detector.A locally constrained regularization method was proposed to determine the resolution calibration parameters.First,a Monte Carlo simulation model consistent with an actual measurement system was constructed to obtain the energy deposition distribution in the scintillation crystal.Subsequently,the regularization objective function is established based on weighted least squares and additional constraints.Additional constraints were designed using a special weighting scheme based on the incident gamma-ray energies.Subsequently,an intelligent algorithm was introduced to search for the optimal resolution calibration parameters by minimizing the objective function.The most appropriate regularization parameter was determined through mathematical experiments.When the regularization parameter was 30,the calibrated results exhibited the minimum RMSE.Simulations and test pit experiments were conducted to verify the performance of the proposed method.The simulation results demonstrate that the proposed algorithm can determine resolution calibration parameters more accurately than the traditional weighted least squares,and the test pit experimental results show that the R-squares between the calibrated and measured spectra are larger than 0.99.The accurate resolution calibration parameters determined by the proposed method lay the foundation for gamma-ray spectral processing and simulation benchmarking.展开更多
Seismic data is commonly acquired sparsely and irregularly, which necessitates the regularization of seismic data with anti-aliasing and anti-leakage methods during seismic data processing. We propose a novel method o...Seismic data is commonly acquired sparsely and irregularly, which necessitates the regularization of seismic data with anti-aliasing and anti-leakage methods during seismic data processing. We propose a novel method of 4D anti-aliasing and anti-leakage Fourier transform using a cube-removal strategy to address the combination of irregular sampling and aliasing in high-dimensional seismic data. We compute a weighting function by stacking the spectrum along the radial lines, apply this function to suppress the aliasing energy, and then iteratively pick the dominant amplitude cube to construct the Fourier spectrum. The proposed method is very efficient due to a cube removal strategy for accelerating the convergence of Fourier reconstruction and a well-designed parallel architecture using CPU/GPU collaborative computing. To better fill the acquisition holes from 5D seismic data and meanwhile considering the GPU memory limitation, we developed the anti-aliasing and anti-leakage Fourier transform method in 4D with the remaining spatial dimension looped. The entire workflow is composed of three steps: data splitting, 4D regularization, and data merging. Numerical tests on both synthetic and field data examples demonstrate the high efficiency and effectiveness of our approach.展开更多
Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analy...Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.展开更多
Full waveform inversion(FWI)is a complex data fitting process based on full wavefield modeling,aiming to quantitatively reconstruct unknown model parameters from partial waveform data with high-resolution.However,this...Full waveform inversion(FWI)is a complex data fitting process based on full wavefield modeling,aiming to quantitatively reconstruct unknown model parameters from partial waveform data with high-resolution.However,this process is highly nonlinear and ill-posed,therefore achieving high-resolution imaging of complex biological tissues within a limited number of iterations remains challenging.We propose a multiscale frequency–domain full waveform inversion(FDFWI)framework for ultrasound computed tomography(USCT)imaging of biological tissues,which innovatively incorporates Sobolev space norm regularization for enhancement of prior information.Specifically,we investigate the effect of different types of hyperparameter on the imaging quality,during which the regularization weight is dynamically adapted based on the ratio of the regularization term to the data fidelity term.This strategy reduces reliance on predefined hyperparameters,ensuring robust inversion performance.The inversion results from both numerical and experimental tests(i.e.,numerical breast,thigh,and ex vivo pork-belly tissue)demonstrate the effectiveness of our regularized FWI strategy.These findings will contribute to the application of the FWI technique in quantitative imaging based on USCT and make USCT possible to be another high-resolution imaging method after x-ray computed tomography and magnetic resonance imaging.展开更多
In this paper,C1,1 regularity for solutions to the degenerate dual Orlicz-Minkowski problem is considered.The dual Orlicz-Minkowski problem is a generalization of the Lp dual Minkowski problem in convex geometry.The p...In this paper,C1,1 regularity for solutions to the degenerate dual Orlicz-Minkowski problem is considered.The dual Orlicz-Minkowski problem is a generalization of the Lp dual Minkowski problem in convex geometry.The proof is adapted from Guan-Li[17]and Chen-Tu-Wu-Xiang[11].展开更多
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＋Rai...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.展开更多
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 inv...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.展开更多
In this paper, we continue to construct stationary classical solutions for the incompressible planar flows approximating singular stationary solutions of this problem. This procedure is carried out by constructing sol...In this paper, we continue to construct stationary classical solutions for the incompressible planar flows approximating singular stationary solutions of this problem. This procedure is carried out by constructing solutions for the following elliptic equations{-△u=λ∑1Bδ（x0,j）（u-kj）p＋,in Ω,u=0,onΩ is a bounded simply-connected smooth domain, ki （i = 1,… , k） is prescribed positive constant. The result we prove is that for any given non-degenerate critical pointX0=（x0,1,…,x0,k of the Kirchhoff-Routh function defined on Ωk corresponding to （ k1,……kk ）there exists a stationary classical solution approximating stationary k points vortex solution. Moreover, as λ→＋∞ shrinks to {x05}, and the local vorticity strength near each x0,j approaches kj, j = 1,… , k. This result makes the study of the above problem with p _〉 0 complete since the cases p 〉 1, p = 1, p = 0 have already been studied in [11, 12] and [13] respectively.展开更多
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 discussio...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.展开更多
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.展开更多
Nuclear Magnetic inversion is the basis of NMR Resonance （NMR） T2 logging interpretation. The regularization parameter selection of the penalty term directly influences the NMR T2 inversion result. We implemented b...Nuclear Magnetic inversion is the basis of NMR Resonance （NMR） T2 logging interpretation. The regularization parameter selection of the penalty term directly influences the NMR T2 inversion result. We implemented both norm smoothing and curvature smoothing methods for NMR T2 inversion, and compared the inversion results with respect to the optimal regular- ization parameters （（Xopt） which were selected by the dis- crepancy principle （DP）, generalized cross-validation （GCV）, S-curve, L-curve, and the slope of L-curve methods, respectively. The numerical results indicate that the DP method can lead to an oscillating or oversmoothed solution which is caused by an inaccurately estimated noise level. The （Xopt selected by the L-curve method is occa- sionally small or large which causes an undersmoothed or oversmoothed T2 distribution. The inversion results from GCV, S-curve and the slope of L-curve methods show satisfying inversion results. The slope of the L-curve method with less computation is more suitable for NMR T2 inversion. The inverted T2 distribution from norm smoothing is better than that from curvature smoothing when the noise level is high.展开更多
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 approxima...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.展开更多
Scatterometer is an instrument which provides all-day and large-scale wind field information, and its application especially to wind retrieval always attracts meteorologists. Certain reasons cause large direction erro...Scatterometer is an instrument which provides all-day and large-scale wind field information, and its application especially to wind retrieval always attracts meteorologists. Certain reasons cause large direction error, so it is important to find where the error mainly comes. Does it mainly result from the background field, the normalized radar cross-section （NRCS） or the method of wind retrieval？ It is valuable to research. First, depending on SDP2.0, the simulated ＇true＇ NRCS is calculated from the simulated ＇true＇ wind through the geophysical mode] function NSCAT2. The simulated background field is configured by adding a noise to the simulated ＇true＇ wind with the non-divergence constraint. Also, the simulated ＇measured＇ NRCS is formed by adding a noise to the simulated ＇true＇ NRCS. Then, the sensitivity experiments are taken, and the new method of regularization is used to improve the ambiguity removal with simulation experiments. The results show that the accuracy of wind retrieval is more sensitive to the noise in the background than in the measured NRCS; compared with the two-dimensional variational （2DVAR） ambiguity removal method, the accuracy of wind retrieval can be improved with the new method of Tikhonov regularization through choosing an appropriate regularization parameter, especially for the case of large error in the background. The work will provide important information and a new method for the wind retrieval with real data.展开更多
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.展开更多
In this paper we discuss the edge-preserving regularization method in the reconstruction of physical parameters from geophysical data such as seismic and ground-penetrating radar data. In the regularization method a p...In this paper we discuss the edge-preserving regularization method in the reconstruction of physical parameters from geophysical data such as seismic and ground-penetrating radar data. In the regularization method a potential function of model parameters and its corresponding functions are introduced. This method is stable and able to preserve boundaries, and protect resolution. The effect of regularization depends to a great extent on the suitable choice of regularization parameters. The influence of the edge-preserving parameters on the reconstruction results is investigated and the relationship between the regularization parameters and the error of data is described.展开更多
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 contain...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 new method of constructing a sea level pressure field from satellite microwave scatterometer measurements is presented. It is based on variational assimilation in combination with a regularization method using geost...A new method of constructing a sea level pressure field from satellite microwave scatterometer measurements is presented. It is based on variational assimilation in combination with a regularization method using geostrophic vorticity to construct a sea level pressure field from scatterometer data that are given in this paper, which offers a new idea for the application of scatterometer measurements. Firstly, the geostrophic vorticity from the scatterometer data is computed to construct the observation field, and the vorticity field in an area and the sea level pressure on the borders are assimilated. Secondly, the gradient of sea level pressure （semi-norm） is used as the stable functional to educe the adjoint system, the adjoint boundary condition and the gradient of the cost functional in which a weight parameter is introduced for the harmony of the system and the Tikhonov regularization techniques in inverse problem are used to overcome the ill-posedness of the assimilation. Finally, the iteration method of the sea level pressure field is developed.展开更多
The sea level pressure field can be computed from sea surface winds retrieved from satellite microwave scatterometer measurements, based on variational assimilation in combination with a regularization method given in...The sea level pressure field can be computed from sea surface winds retrieved from satellite microwave scatterometer measurements, based on variational assimilation in combination with a regularization method given in part I of this paper. First, the validity of the new method is proved with a simulation experiment. Then, a new processing procedure for the sea level pressure retrieval is built by combining the geostrophic wind, which is computed from the scatterometer 10-meter wind using the University of Washington planetary boundary layer model using this method. Finally, the feasibility of the method is proved using an actual case study.展开更多
In this paper, a regularization Newton method for mixed complementarity problem(MCP) based on the reformulation of MCP in [1] is proposed. Its global convergence is proved under the assumption that F is a P0-function....In this paper, a regularization Newton method for mixed complementarity problem(MCP) based on the reformulation of MCP in [1] is proposed. Its global convergence is proved under the assumption that F is a P0-function. The main feature of our algorithm is that a priori of the existence of an accumulation point for convergence need not to be assumed.展开更多
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.展开更多
基金supported by the National Natural Science Foundation of China(No.41804141)。
文摘Energy resolution calibration is crucial for gamma-ray spectral analysis,as measured using a scintillation detector.A locally constrained regularization method was proposed to determine the resolution calibration parameters.First,a Monte Carlo simulation model consistent with an actual measurement system was constructed to obtain the energy deposition distribution in the scintillation crystal.Subsequently,the regularization objective function is established based on weighted least squares and additional constraints.Additional constraints were designed using a special weighting scheme based on the incident gamma-ray energies.Subsequently,an intelligent algorithm was introduced to search for the optimal resolution calibration parameters by minimizing the objective function.The most appropriate regularization parameter was determined through mathematical experiments.When the regularization parameter was 30,the calibrated results exhibited the minimum RMSE.Simulations and test pit experiments were conducted to verify the performance of the proposed method.The simulation results demonstrate that the proposed algorithm can determine resolution calibration parameters more accurately than the traditional weighted least squares,and the test pit experimental results show that the R-squares between the calibrated and measured spectra are larger than 0.99.The accurate resolution calibration parameters determined by the proposed method lay the foundation for gamma-ray spectral processing and simulation benchmarking.
文摘Seismic data is commonly acquired sparsely and irregularly, which necessitates the regularization of seismic data with anti-aliasing and anti-leakage methods during seismic data processing. We propose a novel method of 4D anti-aliasing and anti-leakage Fourier transform using a cube-removal strategy to address the combination of irregular sampling and aliasing in high-dimensional seismic data. We compute a weighting function by stacking the spectrum along the radial lines, apply this function to suppress the aliasing energy, and then iteratively pick the dominant amplitude cube to construct the Fourier spectrum. The proposed method is very efficient due to a cube removal strategy for accelerating the convergence of Fourier reconstruction and a well-designed parallel architecture using CPU/GPU collaborative computing. To better fill the acquisition holes from 5D seismic data and meanwhile considering the GPU memory limitation, we developed the anti-aliasing and anti-leakage Fourier transform method in 4D with the remaining spatial dimension looped. The entire workflow is composed of three steps: data splitting, 4D regularization, and data merging. Numerical tests on both synthetic and field data examples demonstrate the high efficiency and effectiveness of our approach.
基金supported by the Na-tional Natural Science Foundation of China(No.52272369).
文摘Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.
基金supported by the National Natural Science Foundation of China(Grant No.12474461)the Basic and Frontier Exploration Project Independently Deployed by Institute of Acoustics,Chinese Academy of Sciences(Grant No.JCQY202402)the Goal-Oriented Project Independently Deployed by Institute of Acoustics,Chinese Academy of Sciences(Grant No.MBDX202113).
文摘Full waveform inversion(FWI)is a complex data fitting process based on full wavefield modeling,aiming to quantitatively reconstruct unknown model parameters from partial waveform data with high-resolution.However,this process is highly nonlinear and ill-posed,therefore achieving high-resolution imaging of complex biological tissues within a limited number of iterations remains challenging.We propose a multiscale frequency–domain full waveform inversion(FDFWI)framework for ultrasound computed tomography(USCT)imaging of biological tissues,which innovatively incorporates Sobolev space norm regularization for enhancement of prior information.Specifically,we investigate the effect of different types of hyperparameter on the imaging quality,during which the regularization weight is dynamically adapted based on the ratio of the regularization term to the data fidelity term.This strategy reduces reliance on predefined hyperparameters,ensuring robust inversion performance.The inversion results from both numerical and experimental tests(i.e.,numerical breast,thigh,and ex vivo pork-belly tissue)demonstrate the effectiveness of our regularized FWI strategy.These findings will contribute to the application of the FWI technique in quantitative imaging based on USCT and make USCT possible to be another high-resolution imaging method after x-ray computed tomography and magnetic resonance imaging.
文摘In this paper,C1,1 regularity for solutions to the degenerate dual Orlicz-Minkowski problem is considered.The dual Orlicz-Minkowski problem is a generalization of the Lp dual Minkowski problem in convex geometry.The proof is adapted from Guan-Li[17]and Chen-Tu-Wu-Xiang[11].
基金Project supported by the National Natural Science Foundation of China (Grant No. 40775023)
文摘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.
基金Project supported by the National Natural Science Foundation of China(Grant No.41175025)
文摘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.
文摘In this paper, we continue to construct stationary classical solutions for the incompressible planar flows approximating singular stationary solutions of this problem. This procedure is carried out by constructing solutions for the following elliptic equations{-△u=λ∑1Bδ（x0,j）（u-kj）p＋,in Ω,u=0,onΩ is a bounded simply-connected smooth domain, ki （i = 1,… , k） is prescribed positive constant. The result we prove is that for any given non-degenerate critical pointX0=（x0,1,…,x0,k of the Kirchhoff-Routh function defined on Ωk corresponding to （ k1,……kk ）there exists a stationary classical solution approximating stationary k points vortex solution. Moreover, as λ→＋∞ shrinks to {x05}, and the local vorticity strength near each x0,j approaches kj, j = 1,… , k. This result makes the study of the above problem with p _〉 0 complete since the cases p 〉 1, p = 1, p = 0 have already been studied in [11, 12] and [13] respectively.
文摘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.
文摘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.
基金funded by Shell International Exploration and Production Inc.(PT45371)the National Natural Science Foundation of China-China National Petroleum Corporation Petrochemical Engineering United Fund(U1262114)the National Natural Science Foundation of China(41272163)
文摘Nuclear Magnetic inversion is the basis of NMR Resonance （NMR） T2 logging interpretation. The regularization parameter selection of the penalty term directly influences the NMR T2 inversion result. We implemented both norm smoothing and curvature smoothing methods for NMR T2 inversion, and compared the inversion results with respect to the optimal regular- ization parameters （（Xopt） which were selected by the dis- crepancy principle （DP）, generalized cross-validation （GCV）, S-curve, L-curve, and the slope of L-curve methods, respectively. The numerical results indicate that the DP method can lead to an oscillating or oversmoothed solution which is caused by an inaccurately estimated noise level. The （Xopt selected by the L-curve method is occa- sionally small or large which causes an undersmoothed or oversmoothed T2 distribution. The inversion results from GCV, S-curve and the slope of L-curve methods show satisfying inversion results. The slope of the L-curve method with less computation is more suitable for NMR T2 inversion. The inverted T2 distribution from norm smoothing is better than that from curvature smoothing when the noise level is high.
基金supported in part by the National Natural Science Foundation of China(No.62073161)the Fundamental Research Funds 2019“Artificial Intelligence+Special Project”of Nanjing University of Aeronautics and Astronautics(No.2019009)
文摘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.
基金supported by the National Natural Science Foundation of China (Grant No. 40775023)
文摘Scatterometer is an instrument which provides all-day and large-scale wind field information, and its application especially to wind retrieval always attracts meteorologists. Certain reasons cause large direction error, so it is important to find where the error mainly comes. Does it mainly result from the background field, the normalized radar cross-section （NRCS） or the method of wind retrieval？ It is valuable to research. First, depending on SDP2.0, the simulated ＇true＇ NRCS is calculated from the simulated ＇true＇ wind through the geophysical mode] function NSCAT2. The simulated background field is configured by adding a noise to the simulated ＇true＇ wind with the non-divergence constraint. Also, the simulated ＇measured＇ NRCS is formed by adding a noise to the simulated ＇true＇ NRCS. Then, the sensitivity experiments are taken, and the new method of regularization is used to improve the ambiguity removal with simulation experiments. The results show that the accuracy of wind retrieval is more sensitive to the noise in the background than in the measured NRCS; compared with the two-dimensional variational （2DVAR） ambiguity removal method, the accuracy of wind retrieval can be improved with the new method of Tikhonov regularization through choosing an appropriate regularization parameter, especially for the case of large error in the background. The work will provide important information and a new method for the wind retrieval with real data.
基金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 in part by the National Natural Science Foundation of China under Grant-in-Aid 40574053the Program for New Century Excellent Talents in University of China (NCET-06-0602)the National 973 Key Basic Research Development Program (No.2007CB209601)
文摘In this paper we discuss the edge-preserving regularization method in the reconstruction of physical parameters from geophysical data such as seismic and ground-penetrating radar data. In the regularization method a potential function of model parameters and its corresponding functions are introduced. This method is stable and able to preserve boundaries, and protect resolution. The effect of regularization depends to a great extent on the suitable choice of regularization parameters. The influence of the edge-preserving parameters on the reconstruction results is investigated and the relationship between the regularization parameters and the error of data is described.
基金This project was supported by TRAPOYT, the Key Project of Chinese Ministry of Education(104126) the NNSF of China(10371046)
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No.41175025)
文摘A new method of constructing a sea level pressure field from satellite microwave scatterometer measurements is presented. It is based on variational assimilation in combination with a regularization method using geostrophic vorticity to construct a sea level pressure field from scatterometer data that are given in this paper, which offers a new idea for the application of scatterometer measurements. Firstly, the geostrophic vorticity from the scatterometer data is computed to construct the observation field, and the vorticity field in an area and the sea level pressure on the borders are assimilated. Secondly, the gradient of sea level pressure （semi-norm） is used as the stable functional to educe the adjoint system, the adjoint boundary condition and the gradient of the cost functional in which a weight parameter is introduced for the harmony of the system and the Tikhonov regularization techniques in inverse problem are used to overcome the ill-posedness of the assimilation. Finally, the iteration method of the sea level pressure field is developed.
基金Project supported by the National Natural Science Foundation of China (Grant No. 41175025)
文摘The sea level pressure field can be computed from sea surface winds retrieved from satellite microwave scatterometer measurements, based on variational assimilation in combination with a regularization method given in part I of this paper. First, the validity of the new method is proved with a simulation experiment. Then, a new processing procedure for the sea level pressure retrieval is built by combining the geostrophic wind, which is computed from the scatterometer 10-meter wind using the University of Washington planetary boundary layer model using this method. Finally, the feasibility of the method is proved using an actual case study.
基金This subject is supported by the NSF of China (10171055,10226022) NSF of Shandong province(Y2003A02)
文摘In this paper, a regularization Newton method for mixed complementarity problem(MCP) based on the reformulation of MCP in [1] is proposed. Its global convergence is proved under the assumption that F is a P0-function. The main feature of our algorithm is that a priori of the existence of an accumulation point for convergence need not to be assumed.
文摘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.
