求解一维波动方程反问题的消除多次波方法

1．引言一位早期的地球物理学家Noah曾经设想将船和电缆沉入海底采集数据，以避免由海面反射所形成的多次波混响，从而得到更为理想的地震勘探剖面，但是这一方法在实践中是很难实现的．1974年，Ril6y博士［6］从地震波的物理机制出发，提出了一种消除多次波的反卷积方法，可以通过资料处理来实现这个思想．在一继波动方程反问题的研究中，我们重新发现了这一方法［7］；称之为消除多次波方法．方法的基本思想是根据自由表面反射波响应与激发波响应之间的相似关系，通过变换来消除多次波混响，从而近似地得到所需的反射系数信息．与Riley…     

Asymptotic properties of a stochastic predator-prey model with Crowley-Martin functional response

A two-step iterative scheme based on the multiplicative splitting iteration is presented for PageRank computation. The new algorithm is applied to the linear system formulation of the problem. Our method is essentially a two-parameter iteration which can extend the possibility to optimize the iterative process. Theoretical analyses show that the iterative sequence produced by our method is convergent to the unique solution of the linear system, i.e., PageRank vector. An exact parameter region of convergence for the method is strictly proved. In each iteration, the proposed method requires solving two linear sub-systems with the splitting of the coefficient matrix of the problem. We consider using inner iterations to compute approximate solutions of these linear sub-systems. Numerical examples are presented to illustrate the efficiency of the new algorithm.     

波动方程的差分反演模型

为了反演波动方程的系数函数,利用差分离散及扰动假设,推导出一个适合迭代的数值模型．解决了以往方法中正反演模型数值精度不一致问题,以及由此带来的一系列问题．经数值模拟计算说明,该方法是可行的和有效的．     

A procedure for determining a spacewise dependent heat source and the initial temperature

The inverse problem of determining a spacewise dependent heat source, together with the initial temperature for the parabolic heat equation, using the usual conditions of the direct problem and information from two supplementary temperature measurements at different instants of time is studied. These spacewise dependent temperature measurements ensure that this inverse problem has a unique solution, despite the solution being unstable, hence the problem is ill-posed. We propose an iterative algorithm for the stable reconstruction of both the initial data and the source based on a sequence of well-posed direct problems for the parabolic heat equation, which are solved at each iteration step using the boundary element method. The instability is overcome by stopping the iterations at the first iteration for which the discrepancy principle is satisfied. Numerical results are presented for a typical benchmark test example, which has the input measured data perturbed by increasing amounts of random noise. The numerical results show that the proposed procedure gives accurate numerical approximations in relatively few iterations.     

强非线性波动方程孤子行波解

研究了一个强非线性波动方程.利用泛函分析变分迭代方法，首先构造了一个变分, 求出相应的Lagrange乘子；其次构造一个解的变分迭代, 选取初始孤子波；最后利用迭代方法依次求出各次孤子波的近似解.该方法是一个简单可行的近似求解非线性方程的方法     

Asymptotic analysis and accurate approximate solutions for strongly nonlinear conservative symmetric oscillators

A new accurate iterative and asymptotic method is introduced to construct analytical approximate solutions to strongly nonlinear conservative symmetric oscillators. The method is based on applying a second-order expansion with the harmonic balance method and it excludes the requirement of solving a set of coupled nonlinear algebraic equations. Newton's iteration or the linearized model may be readily deduced by considering only the first-order terms in the model. According to this iterative approach, only Fourier series expansions of restoring force function, its first- and second-order derivatives for each iteration are required. It is concluded here that by only one single iteration, very brief and yet accurate analytical approximate solutions can be attained. Three physical examples are solved and accurate solutions are presented to illustrate the physics of the system and the effectiveness of the proposed asymptotic method.     

Fast iteration method in the problem of waves interacting with a set of thin screens

A new iteration method is proposed for the wave equation describing the scattering of a harmonic wave from an arbitrary configuration in the form of an array of thin straight barriers. The problem is reduced to a system of boundary integral equations, which are discretized by applying the Belotserkovskii-Lifanov method. In discrete form, a finite number of systems with Toeplitz matrices (the number of systems is equal to the number of barriers) are solved at each iteration step by applying special fast methods. The algorithm is tested on several geometries, and its convergence in these cases is analyzed.     

Iterative solution of a singular convection-diffusion perturbation problem

An iterative domain decomposition method is developed to solve a singular perturbation problem. The problem consists of a convection-diffusion equation with a discontinuous (piecewise-constant) diffusion coefficient, and the problem domain is decomposed into two subdomains, on each of which the coefficient is constant. After showing that the boundary value problem is well posed, we indicate a specific numerical implementation of the iterative technique that combines the finite element method on one subdomain with the method of matched asymptotic expansions on the other subdomain. This procedure extends work by Carlenzoli and Quarteroni, which was originally intended for a boundary layer problem with an outer region and an inner region. Our extension carries over to a problem where the domain consists of the outer and inner boundary layer regions plus a region in which the diffusion coefficient is constant and significant in magnitude. An unexpected benefit of our new implementation is its efficiency, which is due to the fact that at each iteration the problem needs to be solved explicitly only on one subdomain. It is only when the final approximation on the entire domain is desired that the matched asymptotic expansions approximation need be computed on the second subdomain. Two-dimensional convergence results and numerical results illustrating the method for a two-dimensional test problem are given.     

Variational iterative method and initial-value problems

The variational iterative method is revisited for initial-value problems in ordinary or partial differential equation. A distributional characterization of the Lagrange multiplier - the keystone of the method - is proposed, that may be interpreted as a retarded Green function. Such a formulation makes possible the simplification of the iteration formula into a Picard iterative scheme, and facilitates the convergence analysis. The approximate analytical solution of a nonlinear Klein-Gordon equation with inhomogeneous initial data is proposed.     

Computational identification of the lowest space-wise dependent coefficient of a parabolic equation

In the present work, we consider a nonlinear inverse problem of identifying the lowest coefficient of a parabolic equation. The desired coefficient depends on spatial variables only. Additional information about the solution is given at the final time moment, i.e., we consider the final redefinition. An iterative process is used to evaluate the lowest coefficient, where at each iteration we solve the standard initial-boundary value problem for the parabolic equation. On the basis of the maximum principle for the solution of the differential problem, the monotonicity of the iterative process is established along with the fact that the coefficient is approached from above. The possibilities of the proposed computational algorithm are illustrated by numerical examples for a model two-dimensional problem.     

Iterative solution of a singular convection-diffusion perturbation problem

An iterative domain decomposition method is developed to solve a singular perturbation problem. The problem consists of a convection-diffusion equation with a discontinuous (piecewise-constant) diffusion coefficient, and the problem domain is decomposed into two subdomains, on each of which the coefficient is constant. After showing that the boundary value problem is well posed, we indicate a specific numerical implementation of the iterative technique that combines the finite element method on one subdomain with the method of matched asymptotic expansions on the other subdomain. This procedure extends work by Carlenzoli and Quarteroni, which was originally intended for a boundary layer problem with an outer region and an inner region. Our extension carries over to a problem where the domain consists of the outer and inner boundary layer regions plus a region in which the diffusion coefficient is constant and significant in magnitude. An unexpected benefit of our new implementation is its efficiency, which is due to the fact that at each iteration the problem needs to be solved explicitly only on one subdomain. It is only when the final approximation on the entire domain is desired that the matched asymptotic expansions approximation need be computed on the second subdomain. Two-dimensional convergence results and numerical results illustrating the method for a two-dimensional test problem are given.Received: February 12, 2004     

Identification of the Exchange Coefficient from Indirect Data for a Coupled Continuum Pipe-Flow Model

Calibration and identification of the exchange effect between the karst aquifers and the underlying conduit network are important issues in order to gain a better understanding of these hydraulic systems. Based on a coupled continuum pipe-flow （CCPF for short） model describing flows in karst aquifers, this paper is devoted to the identification of an exchange rate function, which models the hydraulic interaction between the fissured volume （matrix） and the conduit, from the Neumann boundary data, i.e., matrix/conduit seepage velocity. The authors formulate this parameter identification problem as a nonlinear operator equation and prove the compactness of the forward mapping. The stable approximate solution is obtained by two classic iterative regularization methods, namely, the Landweber iteration and Levenberg-Marquardt method. Numerical examples on noisefree and noisy data shed light on the appropriateness of the proposed approaches.     

Optimal control computation for parabolic systems with boundary conditions involving time delays

In this paper, we consider a class of optimal control problems involving a second-order, linear parabolic partial differential equation with Neumann boundary conditions. The time-delayed arguments are assumed to appear in the boundary conditions. A necessary and sufficient condition for optimality is derived, and an iterative method for solving this optimal control problem is proposed. The convergence property of this iterative method is also investigated.On the basis of a finite-element Galerkin's scheme, we convert the original distributed optimal control problem into a sequence of approximate problems involving only lumped-parameter systems. A computational algorithm is then developed for each of these approximate problems. For illustration, a one-dimensional example is solved.     

New iterative scheme with nonexpansive mappings for equilibrium problems and variational inequality problems in Hilbert spaces

Recently, Ceng, Guu and Yao introduced an iterative scheme by viscosity-like approximation method to approximate the fixed point of nonexpansive mappings and solve some variational inequalities in Hilbert space (see Ceng et al. (2009) [9]). Takahashi and Takahashi proposed an iteration scheme to solve an equilibrium problem and approximate the fixed point of nonexpansive mapping by viscosity approximation method in Hilbert space (see Takahashi and Takahashi (2007) [12]). In this paper, we introduce an iterative scheme by viscosity approximation method for finding a common element of the set of a countable family of nonexpansive mappings and the set of an equilibrium problem in a Hilbert space. We prove the strong convergence of the proposed iteration to the unique solution of a variational inequality.     

Determination of a spacewise dependent heat source

This paper investigates the inverse problem of determining a spacewise dependent heat source in the parabolic heat equation using the usual conditions of the direct problem and information from a supplementary temperature measurement at a given single instant of time. The spacewise dependent temperature measurement ensures that the inverse problem has a unique solution, but this solution is unstable, hence the problem is ill-posed. For this inverse problem, we propose an iterative algorithm based on a sequence of well-posed direct problems which are solved at each iteration step using the boundary element method (BEM). The instability is overcome by stopping the iterations at the first iteration for which the discrepancy principle is satisfied. Numerical results are presented for various typical benchmark test examples which have the input measured data perturbed by increasing amounts of random noise.     

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.     

Existence of Solutions to a Class of Fractional Differential Equations

In this paper, the existence of solutions to a class of fractional differential equations $D_{0+}^{\alpha}u(t)=h(t)f(t, u(t), D_{0+}^{\theta}u(t))$ is obtained by an efficient and simple monotone iteration method. At first, the existence of a solution to the problem above is guaranteed by finding a bounded domain $D_M$ on functions $f$ and $g$. Then, sufficient conditions for the existence of monotone solution to the problem are established by applying monotone iteration method. Moreover, two efficient iterative schemes are proposed, and the convergence of the iterative process is proved by using the monotonicity assumption on $f$ and $g$. In particular, a new algorithm which combines Gauss-Kronrod quadrature method with cubic spline interpolation method is adopted to achieve the monotone iteration method in Matlab environment, and the high-precision approximate solution is obtained. Finally, the main results of the paper are illustrated by some numerical simulations, and the approximate solutions graphs are provided by using the iterative method.     

A nonasymptotic method for singular perturbation problems

In this paper, an approximate method for the numerical integration of singularly perturbed two-point boundary-value problems with a boundary layer on the left end of the underlying interval is presented. The method is distinguished by the following fact: the original second-order differential equation is replaced by an approximate first-order differential equation with a small deviating argument and is solved efficiently by employing the Simpson rule, coupled with the discrete invariant imbedding algorithm. The proposed method is iterative on the deviating argument. Several numerical examples have been solved to demonstrate the applicability of the method.     

Regularizing properties of a truncated newton-cg algorithm for nonlinear inverse problems

This paper develops truncated Newton methods as an appropriate tool for nonlinear inverse problems which are ill-posed in the sense of Hadamard. In each Newton step an approximate solution for the linearized problem is computed with the conjugate gradient method as an inner iteration. The conjugate gradient iteration is terminated when the residual has been reduced to a prescribed percentage. Under certain assumptions on the nonlinear operator it is shown that the algorithm converges and is stable if the discrepancy principle is used to terminate the outer iteration. These assumptions are fulfilled, e.g., for the inverse problem of identifying the diffusion coefficient in a parabolic differential equation from distributed data.     

Numerical solution of the inverse electrocardiography problem with the use of the Tikhonov regularization method

The inverse electrocardiography problem related to medical diagnostics is considered in terms of potentials. Within the framework of the quasi-stationary model of the electric field of the heart, the solution of the problem is reduced to the solution of the Cauchy problem for the Laplace equation in R ³. A numerical algorithm based on the Tikhonov regularization method is proposed for the solution of this problem. The Cauchy problem for the Laplace equation is reduced to an operator equation of the first kind, which is solved via minimization of the Tikhonov functional with the regularization parameter chosen according to the discrepancy principle. In addition, an algorithm based on numerical solution of the corresponding Euler equation is proposed for minimization of the Tikhonov functional. The Euler equation is solved using an iteration method that involves solution of mixed boundary value problems for the Laplace equation. An individual mixed problem is solved by means of the method of boundary integral equations of the potential theory. In the study, the inverse electrocardiography problem is solved in region Ω close to the real geometry of the torso and heart.