求解微分方程初值问题的一种弧长法

对于连续介质力学问题中导出的微分方程初值问题,常常具有解奇异性,如不连续、Ｓｔｉｆ性质或激波间断·本文通过在相应空间,引入一个或数个弧长参数变量,克服解的奇异性·对于常微分方程组引入弧长参数变量后,奇异性得以消除和削弱,应用一般的解常微分方程组的方法（如Ｒｕｎｇｅ＿Ｋｕｔａ法）求解·对于偏微分方程引入弧长参数变量后,在相应的空间离散成常微分方程组,用解奇异性常微分方程组相同的方法即可求解·本文给出了两个算例     

矩阵方程AXA^T＋BYB^T=C的双对称最小二乘解及其最佳逼近

基于共轭梯度法的思想,通过特殊的变形,建立了一类求矩阵方程AXA^T＋BYB^T=C的双对称最小二乘解的迭代算法．对任意的初始双对称矩阵．在没有舍人误差的情况下,经过有限步迭代得到它的双对称最小二乘解;在选取特殊的初始双对称矩阵时,能得到它的的极小范数双对称最小二乘解．另外,给定任意矩阵,利用此方法可得到它的最佳逼近双对称解,数值例子表明,这种方法是有效的．     

MUSIC ALGORITHM FOR LOCATING POINT-LIKE SCATTERERS CONTAINED IN A SAMPLE ON FLAT SUBSTRATE

In this paper, we consider a MUSIC algorithm for locating point-like scatterers contained in a sample on flat substrate. Based on an asymptotic expansion of the scattering amplitude proposed by Ammari et al., the reconstruction problem can be reduced to a calculation of Green function corresponding to the background medium. In addition, we use an explicit formulation of Green function in the MUSIC algorithm to simplify the calculation when the cross-section of sample is a half-disc. Numerical experiments are included to demonstrate the feasibility of this method.     

Fourth‐order compact scheme with local mesh refinement for option pricing in jump‐diffusion model

The value of a contingent claim under a jump‐diffusion process satisfies a partial integro‐differential equation. A fourth‐order compact finite difference scheme is applied to discretize the spatial variable of this equation. It is discretized in time by an implicit‐explicit method. Meanwhile, a local mesh refinement strategy is used for handling the nonsmooth payoff condition. Moreover, the numerical quadrature method is exploited to evaluate the jump integral term. It guarantees a Toeplitz‐like structure of the integral operator such that a fast algorithm is feasible. Numerical results show that this approach gives fourth‐order accuracy in space. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011     

Finite difference approximations and dynamics simulations for the Lévy Fractional Klein‐Kramers equation

The Klein‐Kramers equation describes position and velocity distribution of Langevin dynamics, the diffusion equation and Fokker‐Planck equation are its special cases for characterizing position distribution and velocity distribution, respectively. Incorporating the mechanisms of Lévy flights into the Klein‐Kramers formalism leads to the Lévy fractional Klein‐Kramers equation, which can effectively describe Lévy flights in the presence of an external force field in the phase space. For numerically solving the Lévy fractional Klein‐Kramers equation, this article presents the explicit and implicit finite difference schemes. The discrete maximum principle is generalized, using this result the detailed stability and convergence analyses of the schemes are given. And the extrapolation and some other possible techniques for improving the convergent rate or making the schemes efficient in more general cases are also discussed. The extensive numerical experiments are performed to confirm the effectiveness of the numerical schemes or simulate the superdiffusion processes. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

An asymptotic method for calculating limits of maximal means

For a continuous almost periodic function , we show that the function

Oscillation and Nonoscillation Criteria for a Second Order Linear Equation

New oscillation and nonoscillation criteria are established for the equation

具有未知源的某些非线性伪抛物型方程的反问题

考 虑具有未 知源项的 某些非 线性伪 抛物 型方程 的反演 问题． 首先 将伪抛 物型 方程初 边值问 题化为非线 性发展方 程 Ｃｏｕｃｈ ｙ 问题,然 后,利用半 群理论,论 证发展 方程反问 题解的存 在唯一 性,最后, 利用不 动点方法得到 伪抛物型方程反 问题的可解性     

三维地磁流方程的整体吸引子

在本文中,我们通过引入一类Ｆｒｅｃｈｅｔ空间,得到了三维有界光滑区域上的地磁流方程的弱解的整体吸引子的存在性。     

Solution of approximate gluon evolution equation for small-x at fixed ‘ρ’

A simple method has been applied to solve the approximate gluon evolution equation for small-x at fixedρ(≡√ln(x ₀/x)/ln[ln(Q ²/Λ²)/ln(Q ₀ ² /Λ²)]. Numerical comparison is made with the predictions from ‘double asymptotic scaling’ and fit. Better agreement is found between our solution and fit nearρ=1. The solution gives approximate double scaling in this region having ‘hard’ pomeron with small contamination.