On the discrepancy principle for some Newton type methods for solving nonlinear inverse problems

We consider the computation of stable approximations to the exact solution of nonlinear ill-posed inverse problems F(x) = y with nonlinear operators F : X → Y between two Hilbert spaces X and Y by the Newton type methods

On some questions arising in numerical realization of amplitude-phase methods

For computing rapidly oscillating solutions of certain second order differential equations a new version of amplitude-phase methods has recently been proposed in [11]. Error estimates were given to approximate solutions for large arguments in [15]. One of the most important points in these methods is the introduction of Prüfer transformation modified by auxiliary functions. Their appropriate choice makes the methods applicable and efficient. When implementing and applying the methods to practical problems, we face some further questions. In this paper we describe and try to answer them. Efficiency of the methods is confirmed by numerical experiments on concrete problems. This revised version was published online in June 2006 with corrections to the Cover Date.     

Numerical methods for some inverse problems of heart electrophysiology

We consider numerical methods for solving inverse problems that arise in heart electrophysiology. The first inverse problem is the Cauchy problem for the Laplace equation. Its solution algorithm is based on the Tikhonov regularization method and the method of boundary integral equations. The second inverse problem is the problem of finding the discontinuity surface of the coefficient of conductivity of a medium on the basis of the potential and its normal derivative given on the exterior surface. For its numerical solution, we suggest a method based on the method of boundary integral equations and the assumption on a special representation of the unknown surface.     

Wavelet-based numerical techniques for 1D peristatic problems in infinite domain

Peristatics is an important branch of continuum mechanics dealing with nonlocal effects in solid structures. Some 1D peristatic problems with four different micromodulus functions have been studied here. A wavelet-based collocation method has been developed to find multiscale approximate solutions of the peristatic problems that avoid cumbersome evaluation of singular type integrals present in analytical solution available in literature. The obtained results offer important insights into applications and simulations of peristatic models and seem to be useful for investigation in the field of peridynamics in one or higher dimensions.     

On some ranks of infinite groups

Abstract A group G has finite Hirsch-Zaicev rank r_hz(G) = r if G has an ascending series whose factors are either infinite cyclic or periodic and if the number of infinite cyclic factors is exactly r. The authors discuss groups with finite Hirsch-Zaicev rank and the connection between this and groups having finite section p-rank for some prime p, or p=0. Groups all of whose abelian subgroups are of bounded rank are also discussed. Keywords: p-rank, locally generalized radical group, Hirsch-Zaicev rank, torsion-free rank, rank Mathematics Subject Classification (2000): 20F19, 20E25, 20E15     

A numerical method for solving an inverse thermoacoustic problem

In this paper, we consider an inverse problem of determining the initial condition of an initial boundary value problem for the wave equation with some additional information about solving a direct initial boundary value problem. The information is obtained from measurements at the boundary of the solution domain. The purpose of our paper is to construct a numerical algorithm for solving the inverse problem by an iterative method called a method of simple iteration (MSI) and to study the resolution quality of the inverse problem as a function of the number and location of measurement points. Three two-dimensional inverse problem formulations are considered. The results of our numerical calculations are presented. It is shown that the MSI decreases the objective functional at each iteration step. However, due to the ill-posedness of the inverse problem the difference between the exact and approximate solutions decreases up to some fixed number k _min, and then monotonically increases. This shows the regularizing properties of the MSI, and the iteration number can be considered a regularization parameter.     

On the inverse problem for three-dimensional potential scattering

We consider a linear perturbation for the wave equation $\text{□u = 0}\text{in}\text{Ω = E}{}^3$ by “repulsive” smooth potentials q(y) that are small at infinity and suitably small (in magnitude). We use a time-dependent approach to prove that the scattering operator S(q) determines uniquenes uniquely the scatterer q (at least in this class). Energy inequalities will play a central role in our discussion.     

On the inverse problem for a logarithmic potential

Siberian Mathematical Journal -     

Direct numerical method for an inverse problem of hyperbolic equations

A coefficient inverse problem of the one-dimensional hyperbolic equation with overspecified boundary conditions is solved by the finite difference method. The computation is carried out in the x direction instead of the usual t direction. The original boundary condition and the overspecified boundary data are used as the new initial conditions, and the original data at t = 0 are used to compute the coefficient directly. The computation time used by this scheme is almost equal to that for solving the hyperbolic equation in the same region once, even though the inverse problem is essentially nonlinear and hence more difficult to solve. An error estimate is obtained that guarantees the stability of the scheme marching in the x direction. Several numerical experiments are carried out to show the convergence and other properties of the scheme. © 1992 John Wiley & Sons, Inc.     

A comparative study of some reconstruction methods for linear inverse problems

In this note, we compare several reconstruction methods to solve a linear ill-conditioned problem, a finite Haussdorf moment problem. The methods are the usual L₂-regularization method, the linear programming method, and two maxentropic reconstruction methods. The scale seems to lean toward the maxentropic reconstructions methods.     

On some iterative numerical methods for a Volterra functional integral equation of the second kind

In this paper, we study an iterative numerical method for approximating solutions of a certain type of Volterra functional integral equations of the second kind (Volterra integral equations where both limits of integration are variables). The method uses the contraction principle and a suitable quadrature formula. Under certain conditions, we prove the existence and uniqueness of the solution and give error estimates for our approximations. We also included a numerical example which illustrates the fast approximations.     

On some inverse problems for elliptic equations and systems

Under study are the inverse problems of determining the right-hand side of a particular form and the solution for elliptic systems, including a series of elasticity systems. (On the boundary of the domain the solution satisfies either the Dirichlet conditions or mixed Dirichlet-Neumann conditions.) We assume that on a system of planes the normal derivatives of the solution can have discontinuities of the first kind. The conjugating boundary conditions on the discontinuity surface are analogous to the continuity conditions for the fields of displacements and stresses for a horizontally laminated medium. The overdetermination conditions are integral (the average of the solution over some domain is specified) or local (the values of the solution on some lines are specified). We study the solvability conditions for these problems and their Fredholm property.     

On 5- and 10-dissections for some infinite products

The Ramanujan Journal - Quite recently, Xia and Zhao established the 10-dissections for Hirschhorn’s two infinite q-series products by using two MAPLE packages and the theory of modular...     

On some new criteria for infinite differentiability of periodic functions

We study the set of infinitely differentiable periodic functions in terms of generalized -derivatives defined by a pair of sequences ψ ₁ and ψ ₂. It is shown that every function ƒ from the set has at least one derivative whose parameters ψ ₁ and ψ ₂ decrease faster than any power function and, at the same time, for any function ƒ ∈ different from a trigonometric polynomial, there exists a pair ψ whose parameters ψ ₁ and ψ ₂ have the same rate of decrease and for which the -derivative no longer exists. __________ Translated from Ukrains’kyi Matematychnyi Zhurnal, Vol. 59, No. 10, pp. 1399–1409, October, 2007.     

On estimate for numerical radius of some contractions

For the numerical radius of an arbitrary nilpotent operator T on a Hilbert space H, Haagerup and de la Harpe proved the inequality , where n ≥ 2 is the nilpotency order of the operator T. In the present paper, we prove a Haagerup-de la Harpe-type inequality for the numerical radius of contractions from more general classes. Published in Ukrains’kyi Matematychnyi Zhurnal, Vol. 58, No. 10, pp. 1335–1339, October, 2006.