Stability of Least Squares for the Solution of an Ill-posed Inverse Problem of Reconstructing Real Value of the Pcrrnittivity of a Dielectric Layer in a Rectangular Waveguide

Author(s):  
E. A. Sheina ◽  
Y. V. Shestopalov ◽  
A. P. Smitnov
2017 ◽  
Vol 7 (2) ◽  
pp. 325-342
Author(s):  
Xi-Ming Fang ◽  
Fu-Rong Lin ◽  
Chao Wang

AbstractWe consider the nonlinear and ill-posed inverse problem where the Robin coefficient in the Laplace equation is to be estimated using the measured data from the accessible part of the boundary. Two regularisation methods are considered — viz. L2 and H1 regularisation. The regularised problem is transformed to a nonlinear least squares problem; and a suitable regularisation parameter is chosen via the normalised cumulative periodogram (NCP) curve of the residual vector under the assumption of white noise, where information on the noise level is not required. Numerical results show that the proposed method is efficient and competitive.


2020 ◽  
Vol 20 (4) ◽  
pp. 783-798
Author(s):  
Shukai Du ◽  
Nailin Du

AbstractWe give a factorization formula to least-squares projection schemes, from which new convergence conditions together with formulas estimating the rate of convergence can be derived. We prove that the convergence of the method (including the rate of convergence) can be completely determined by the principal angles between {T^{\dagger}T(X_{n})} and {T^{*}T(X_{n})}, and the principal angles between {X_{n}\cap(\mathcal{N}(T)\cap X_{n})^{\perp}} and {(\mathcal{N}(T)+X_{n})\cap\mathcal{N}(T)^{\perp}}. At the end, we consider several specific cases and examples to further illustrate our theorems.


Author(s):  
Karl Kunisch ◽  
Philip Trautmann

AbstractIn this work we discuss the reconstruction of cardiac activation instants based on a viscous Eikonal equation from boundary observations. The problem is formulated as a least squares problem and solved by a projected version of the Levenberg–Marquardt method. Moreover, we analyze the well-posedness of the state equation and derive the gradient of the least squares functional with respect to the activation instants. In the numerical examples we also conduct an experiment in which the location of the activation sites and the activation instants are reconstructed jointly based on an adapted version of the shape gradient method from (J. Math. Biol. 79, 2033–2068, 2019). We are able to reconstruct the activation instants as well as the locations of the activations with high accuracy relative to the noise level.


Proceedings ◽  
2019 ◽  
Vol 33 (1) ◽  
pp. 21
Author(s):  
Fabrizia Guglielmetti ◽  
Eric Villard ◽  
Ed Fomalont

A stable and unique solution to the ill-posed inverse problem in radio synthesis image analysis is sought employing Bayesian probability theory combined with a probabilistic two-component mixture model. The solution of the ill-posed inverse problem is given by inferring the values of model parameters defined to describe completely the physical system arised by the data. The analysed data are calibrated visibilities, Fourier transformed from the ( u , v ) to image planes. Adaptive splines are explored to model the cumbersome background model corrupted by the largely varying dirty beam in the image plane. The de-convolution process of the dirty image from the dirty beam is tackled in probability space. Probability maps in source detection at several resolution values quantify the acquired knowledge on the celestial source distribution from a given state of information. The information available are data constrains, prior knowledge and uncertain information. The novel algorithm has the aim to provide an alternative imaging task for the use of the Atacama Large Millimeter/Submillimeter Array (ALMA) in support of the widely used Common Astronomy Software Applications (CASA) enhancing the capabilities in source detection.


Geophysics ◽  
2021 ◽  
pp. 1-60
Author(s):  
Chuang Li ◽  
Zhaoqi Gao ◽  
Jinghuai Gao ◽  
Feipeng Li ◽  
Tao Yang

Angle-domain common-image gathers (ADCIGs) that can be used for migration velocity analysis and amplitude versus angle analysis are important for seismic exploration. However, because of limited acquisition geometry and seismic frequency band, the ADCIGs extracted by reverse time migration (RTM) suffer from illumination gaps, migration artifacts, and low resolution. We have developed a reflection angle-domain pseudo-extended plane-wave least-squares RTM method for obtaining high-quality ADCIGs. We build the mapping relations between the ADCIGs and the plane-wave sections using an angle-domain pseudo-extended Born modeling operator and an adjoint operator, based on which we formulate the extraction of ADCIGs as an inverse problem. The inverse problem is iteratively solved by a preconditioned stochastic conjugate gradient method, allowing for reduction in computational cost by migrating only a subset instead of the whole dataset and improving image quality thanks to preconditioners. Numerical tests on synthetic and field data verify that the proposed method can compensate for illumination gaps, suppress migration artifacts, and improve resolution of the ADCIGs and the stacked images. Therefore, compared with RTM, the proposed method provides a more reliable input for migration velocity analysis and amplitude versus angle analysis. Moreover, it also provides much better stacked images for seismic interpretation.


Sign in / Sign up

Export Citation Format

Share Document