Uniqueness and nonuniqueness of the solution of boundary-contact acoustics problems

One establishes uniqueness theorems and one gives examples of nonuniqueness for a series of boundary-contact acoustics problems, i.e. diffraction problems under boundary conditions, described by a higher order operator and vanishing at isolated points.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 128, pp. 21–32, 1983.     

Boundary-contact problems for domains with edge singularities

We study boundary-contact problems for elliptic equations (and systems) with interfaces that have edge singularities. Such problems represent continuous operators between weighted edge spaces and subspaces with asymptotics. Ellipticity is formulated in terms of a principal symbolic hierarchy, containing interior, transmission, and edge symbols. We construct parametrices, show regularity with asymptotics of solutions in weighted edge spaces and illustrate the results by boundary-contact problems for the Laplacian with jumping coefficients.     

Some Boundary — Contact Problems of the Elasticity Theory with Mixed Boundary Conditions Outside the Contact Surface

n — Dimensional (n ≥ 2) boundary-contact problems of statics of the elasticity theory for homogeneous anisotropic media are investigated when the contact of two bounded domains occurs from the outside on some part of boundaries with mixed boundary conditions. Theorems on the existence and uniqueness of solutions of boundary-contact problems in Besov and Bessel potential spaces are obtained. The smoothness of solutions is studied in closed domains occupied by elastic media.     

The transmission of a flexural-gravitational wave through several straight obstacles in a floating plate

Periodic wave processes in a thin elastic plate floating on the surface of an incompressible fluid of finite depth are studied. The plate completely covers the fluid surface and executes flexural oscillations under the action of gravitational waves in the fluid. The system of free oscillations in the plate is disrupted along a set of parallel lines. Rigid clamping of the plate, a sliding fastening and an infinitesimally narrow slit are considered as such disruptions. The apparatus used to construct the solution is quite general, and other disruptions in the elastic properties of a plate or its reinforcement, that are realized with linear boundary-contact conditions, can be treated in a similar way. The transmission and reflection of a harmonic flexural-gravitational wave, that is orthogonally incident on the inhomogeneities in the plate, are studied. Exact analytical representations of the wave fields in the plate and in fluid are obtained and the transmission and reflection coefficients for the incident flexural-gravitational wave are determined. The forces developed in the fastenings are found.     

A regularization method for solving the Cauchy problem for the Helmholtz equation

In this paper, we investigate a Cauchy problem associated with Helmholtz-type equation in an infinite “strip”. This problem is well known to be severely ill-posed. The optimal error bound for the problem with only nonhomogeneous Neumann data is deduced, which is independent of the selected regularization methods. A framework of a modified Tikhonov regularization in conjunction with the Morozov’s discrepancy principle is proposed, it may be useful to the other linear ill-posed problems and helpful for the other regularization methods. Some sharp error estimates between the exact solutions and their regularization approximation are given. Numerical tests are also provided to show that the modified Tikhonov method works well.     

Regularization of Discrete Ill-Posed Problems

The paper concerns conditioning aspects of finite-dimensional problems arising when the Tikhonov regularization is applied to discrete ill-posed problems. A relation between the regularization parameter and the sensitivity of the regularized solution is investigated. The main conclusion is that the condition number can be decreased only to the square root of that for the nonregularized problem. The convergence of solutions of regularized discrete problems to the exact generalized solution is analyzed just in the case when the regularization corresponds to the minimal condition number. The convergence theorem is proved under the assumption of the suitable relation between the discretization level and the data error. As an example the method of truncated singular value decomposition with regularization is considered. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the condition number of the antireflective transform

Deconvolution problems with a finite observation window require appropriate models of the unknown signal in order to guarantee uniqueness of the solution. For this purpose it has recently been suggested to impose some kind of antireflectivity of the signal. With this constraint, the deconvolution problem can be solved with an appropriate modification of the fast sine transform, provided that the convolution kernel is symmetric. The corresponding transformation is called the antireflective transform. In this work we determine the condition number of the antireflective transform to first order, and use this to show that the so-called reblurring variant of Tikhonov regularization for deconvolution problems is a regularization method. Moreover, we establish upper bounds for the regularization error of the reblurring strategy that hold uniformly with respect to the size n of the algebraic system, even though the condition number of the antireflective transform grows with n. We briefly sketch how our results extend to higher space dimensions.     

The Fourier method in three-dimensional boundary-contact dynamic problems of elasticity

The basic three-dimensional boundary-contact dynamic problems are considered for a piecewise-homogeneous isotropic elastic medium bounded by several closed surfaces. Using the Fourier method, the considered problems are proved to be solvable under much weaker restrictions on the initial data of the problems as compared with other methods.     

Asymptotic Distribution of Eigenfunctions and Eigenvalues of the Basic Boundary-Contact Oscillation Problems of the Classical Theory of Elasticity

The basic boundary-contact oscillation problems are considered for a three-dimensional piecewise-homogeneous isotropic elastic medium bounded by several closed surfaces. Using Carleman's method, the asymptotic formulas for the distribution of eigenfunctions and eigenvalues are obtained.     

Solution of the Basic Boundary Value Problems of Stationary Thermoelastic Oscillations for Domains Bounded by Spherical Surfaces

The boundary value problems of stationary thermoelastic oscillations are investigated for the entire space with a spherical cavity, when the limit values of a displacement vector and temperature or of a stress vector and heat flow are given on the boundary. Also, consideration is given to the boundary-contact problems when a nonhomogeneous medium fills up the entire space and consists of several homogeneous parts with spherical interface surfaces. Given on an interface surface are differences of the limit values of displacement and stress vectors, also of temperature and heat flow, while given on a free boundary are the limit values of a displacement vector and temperature or of a stress vector and heat flow. Solutions of the considered problems are represented as absolutely and uniformly convergent series.     

约束Fractional Tikhonov正则化的模迭代方法

反问题是现在数学物理研究中的一个热点问题,而反问题求解面临的一个本质性困难是不适定性。求解不适定问题的普遍方法是:用与原不适定问题相“邻近”的适定问题的解去逼近原问题的解,这种方法称为正则化方法.如何建立有效的正则化方法是反问题领域中不适定问题研究的重要内容.当前,最为流行的正则化方法有基于变分原理的Tikhonov正则化及其改进方法,此类方法是求解不适定问题的较为有效的方法,在各类反问题的研究中被广泛采用,并得到深入研究.     

Necessary and sufficient conditions for linear convergence of ℓ1‐regularization

Motivated by the theoretical and practical results in compressed sensing, efforts have been undertaken by the inverse problems community to derive analogous results, for instance linear convergence rates, for Tikhonov regularization with ℓ¹‐penalty term for the solution of ill‐posed equations. Conceptually, the main difference between these two fields is that regularization in general is an uncon strained optimization problem, while in compressed sensing a constrained one is used. Since the two methods have been developed in two different communities, the theoretical approaches to them appear to be rather different: In compressed sensing, the restricted isometry property seems to be central for proving linear convergence rates, whereas in regularization theory range or source conditions are imposed. The paper gives a common meaning to the seemingly different conditions and puts them into perspective with the conditions from the respective other community. A particularly important observation is that the range condition together with an injectivity condition is weaker than the restricted isometry property. Under the weaker conditions, linear convergence rates can be proven for compressed sensing and for Tikhonov regularization. Thus existing results from the literature can be improved based on a unified analysis. In particular, the range condition is shown to be the weakest possible condition that permits the derivation of linear convergence rates for Tikhonov regularization with a priori parameter choice. © 2010 Wiley Periodicals, Inc.     

On the convergence of a regularization method for variational inequalities

For variational inequalities in a finite-dimensional space, the convergence of a regularization method is examined in the case of a nonmonotone basic mapping. It is shown that a fairly general sufficient condition for the existence of solutions to the original problem also guarantees the convergence and existence of solutions to perturbed problems. Examples of applications to problems on order intervals are presented.     

REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems

The package REGULARIZATION TOOLS consists of 54 Matlab routines for analysis and solution of discrete ill-posed problems, i.e., systems of linear equations whose coefficient matrix has the properties that its condition number is very large, and its singular values decay gradually to zero. Such problems typically arise in connection with discretization of Fredholm integral equations of the first kind, and similar ill-posed problems. Some form of regularization is always required in order to compute a stabilized solution to discrete ill-posed problems. The purpose of REGULARIZATION TOOLS is to provide the user with easy-to-use routines, based on numerical robust and efficient algorithms, for doing experiments with regularization of discrete ill-posed problems. By means of this package, the user can experiment with different regularization strategies, compare them, and draw conclusions from these experiments that would otherwise require a major programming effert. For discrete ill-posed problems, which are indeed difficult to treat numerically, such an approach is certainly superior to a single black-box routine. This paper describes the underlying theory gives an overview of the package; a complete manual is also available.This work was supported by grants from Augustinus Fonden, Knud Højgaards Fond, and Civ. Ing. Frants Allings Legat.     

A regularized limited memory BFGS method for nonconvex unconstrained minimization

The limited memory BFGS method (L-BFGS) is an adaptation of the BFGS method for large-scale unconstrained optimization. However, The L-BFGS method need not converge for nonconvex objective functions and it is inefficient on highly ill-conditioned problems. In this paper, we proposed a regularization strategy on the L-BFGS method, where the used regularization parameter may play a compensation role in some sense when the condition number of Hessian approximation tends to become ill-conditioned. Then we proposed a regularized L-BFGS method and established its global convergence even when the objective function is nonconvex. Numerical results show that the proposed method is efficient.     

L1 regularization method in electrical impedance tomography by using the L1-curve (Pareto frontier curve)

Electrical impedance tomography (EIT), as an inverse problem, aims to calculate the internal conductivity distribution at the interior of an object from current-voltage measurements on its boundary. Many inverse problems are ill-posed, since the measurement data are limited and imperfect. To overcome ill-posedness in EIT, two main types of regularization techniques are widely used. One is categorized as the projection methods, such as truncated singular value decomposition (SVD or TSVD). The other categorized as penalty methods, such as Tikhonov regularization, and total variation methods. For both of these methods, a good regularization parameter should yield a fair balance between the perturbation error and regularized solution. In this paper a new method combining the least absolute shrinkage and selection operator (LASSO) and the basis pursuit denoising (BPDN) is introduced for EIT. For choosing the optimum regularization we use the L1-curve (Pareto frontier curve) which is similar to the L-curve used in optimising L2-norm problems. In the L1-curve we use the L1-norm of the solution instead of the L2 norm. The results are compared with the TSVD regularization method where the best regularization parameters are selected by observing the Picard condition and minimizing generalized cross validation (GCV) function. We show that this method yields a good regularization parameter corresponding to a regularized solution. Also, in situations where little is known about the noise level σ, it is also useful to visualize the L1-curve in order to understand the trade-offs between the norms of the residual and the solution. This method gives us a means to control the sparsity and filtering of the ill-posed EIT problem. Tracing this curve for the optimum solution can decrease the number of iterations by three times in comparison with using LASSO or BPDN separately.     

Boundary-contact problems for domains with conical singularities

We study boundary-contact problems for elliptic equations (and systems) with interfaces that have conical singularities. Such problems represent continuous operators between weighted Sobolev spaces and subspaces with asymptotics. Ellipticity is formulated in terms of extra transmission conditions along the interfaces with a control of the conormal symbolic structure near conical singularities. We show regularity and asymptotics of solutions in weighted spaces, and we construct parametrices. The result will be illustrated by a number of explicit examples.     

Uniqueness of the solution of boundary-contact problems of acoustics for a plate-liquid system

Conditions guaranteeing the uniqueness of the solution of boundary-contact problems for a system consisting of plates immersed in an acoustic medium are formulated on the basis of the law of conservation of energy.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 99, pp. 43–56, 1980.The author is grateful to B. P. Belinskii for his attention to the work and for useful remarks.     

Operator-difference schemes for a class of systems of evolution equations

For a special system of evolution equations of first order, discrete time approximations for the approximate solution of the Cauchy problem are considered. Such problems arise after the spatial approximation in the Schrödinger equation and the subsequent separation of the imaginary and real parts and in nonstationary problems of acoustics and electrodynamics. Unconditionally stable two time level operator-difference weighted schemes are constructed. The second class of difference schemes is based on the formal passage to explicit operator-difference schemes for evolution equations of second order when explicit-implicit approximation is used for isolated equations of the system. The regularization of such schemes in order to obtain unconditionally stable operator difference schemes is discussed. Splitting schemes involving the solution of simplest problems at each time step are constructed.     

The Nonclassical Boundary-contact Problems of Elasticity for Homogeneous Anisotropic Media

The existence and uniqueness of solutions of the nonclassical boundary-contact problems (i.e., problems with a contact on some part of the boundaries) of elasticity for homogeneous anisotropic media are investigated in Besov and Bessel potential spaces using methods of potential theory and the theory of pseudodifferential equations on manifolds with boundary. The smoothness of the solutions obtained is studied.