带有临界Sobolev指数的奇异双调和问题

LetΩR~N be a smooth bounded domain such that 0∈Ω,N≥5,2~*:=(2N)/(N-4) is the critical Sobolev exponent,and f(x) is a given function.By using the variational methods, the paper proves the existence of solutions for the singular critical in the homogeneous problemΔ~u-μu/(|x|~4)=|u|~(2~*-2)u f(x) with Dirichlet boundary condition on Ωunder some assumptions on f(x) andμ.     

Tomographic reconstruction using heuristic Monte Carlo methods

A tomographic reconstruction method based on Monte Carlo random searching guided by the information contained in the projections of radiographed objects is presented. In order to solve the optimization problem, a multiscale algorithm is proposed to reduce computation. The reconstruction is performed in a coarse-to-fine multigrid scale that initializes each resolution level with the reconstruction of the previous coarser level, which substantially improves the performance. The method was applied to a real case reconstructing the internal structure of a small metallic object with internal components, showing excellent results.     

求解变延迟微分方程的一类线性多步方法的收缩性

This paper presents a sufficient condition on the contractivity of theoretical solution for a class of nonlinear systems of delay differential equations with many variable delays(MDDEs), which is weak,compared with the sufficient condition of previous articles.In addition,it discusses the numerical stability properties of a class of special linear nmltistep methods for this class nonlinear problems.And it is pointed out that not only the backwm‘d Euler method but also this class of linear multistep methods are GRNm-stable if linear interpolation is used.     

Predictor-Corrector Smoothing Methods for Monotone LCP

In this paper, we analyze the global and local convergence properties of two predictor-corrector smoothing methods, which are based on the framework of the method in [1], for monotone linear complementarity problems (LCPs). The difference between the algorithm in [1] and our algorithms is that the neighborhood of smoothing central path in our paper is different to that in [1]. In addition, the difference between Algorithm 2.1 and the algorithm in [1] exists in the calculation of the predictor step. Comparing with the results in [1],the global and local convergence of the two methods can be obtained under very mild conditions. The global convergence of the two methods do not need the boundness of the inverse of the Jacobian. The superlinear convergence of Algorithm 2.1‘ is obtained under the assumption of nonsingularity of generalized Jacobian of Φ(x,y) at the limit point and Algorithm 2.1 obtains superlinear convergence under the assumption of strict complementarity at the solution. The efficiency of the two methods is tested by numerical experiments.     

Strong tractability of multivariate integration using quasi-Monte Carlo algorithms

We study quasi-Monte Carlo algorithms based on low discrepancy sequences for multivariate integration. We consider the problem of how the minimal number of function evaluations needed to reduce the worst-case error from its initial error by a factor of depends on and the dimension . Strong tractability means that it does not depend on and is bounded by a polynomial in . The least possible value of the power of is called the -exponent of strong tractability. Sloan and Wozniakowski established a necessary and sufficient condition of strong tractability in weighted Sobolev spaces, and showed that the -exponent of strong tractability is between 1 and 2. However, their proof is not constructive.

In this paper we prove in a constructive way that multivariate integration in some weighted Sobolev spaces is strongly tractable with -exponent equal to 1, which is the best possible value under a stronger assumption than Sloan and Wozniakowski's assumption. We show that quasi-Monte Carlo algorithms using Niederreiter's -sequences and Sobol sequences achieve the optimal convergence order for any 0$"> independent of the dimension with a worst case deterministic guarantee (where is the number of function evaluations). This implies that strong tractability with the best -exponent can be achieved in appropriate weighted Sobolev spaces by using Niederreiter's -sequences and Sobol sequences.

     

Enhanced accuracy by post-processing for finite element methods for hyperbolic equations

We consider the enhancement of accuracy, by means of a simple post-processing technique, for finite element approximations to transient hyperbolic equations. The post-processing is a convolution with a kernel whose support has measure of order one in the case of arbitrary unstructured meshes; if the mesh is locally translation invariant, the support of the kernel is a cube whose edges are of size of the order of only. For example, when polynomials of degree are used in the discontinuous Galerkin (DG) method, and the exact solution is globally smooth, the DG method is of order in the -norm, whereas the post-processed approximation is of order ; if the exact solution is in only, in which case no order of convergence is available for the DG method, the post-processed approximation converges with order in , where is a subdomain over which the exact solution is smooth. Numerical results displaying the sharpness of the estimates are presented.

     

基于粗糙集理论的胸部X线数字图象增强方法的研究

采用一种基于粗糙集理论实现胸部X线数字图象增强的新方法,该方法根据粗糙集的条件属性,将一幅采用双能量算法获得的肺部软组织X线数字图象划分为不同的子图象,然后对各个子图象分别作衬比增强从结果中可看出,采用种方法处理后得到的图象效果较为理想,能给临床上的医学工作者提供清晰的图象,同时也为后续的计算机辅助诊断提供有效的依据.     

Transfer reactions and sub-barrier fusion in 40Ca + 90, 96Zr

The two systems ⁴⁰Ca + ^90,96Zr have been studied by measuring nucleon transfer reactions at two energies near the Coulomb barrier, thus complementing the available sub-barrier fusion cross-sections. Angular distributions for various transfer channels have been determined. Significantly larger neutron transfer cross-sections are found for the target ⁹⁶Zr that exhibits the larger enhancement in the sub-barrier fusion cross-sections. All data have been analyzed with a new model for heavy-ion collisions that calculates simultaneously transfer cross-sections, fusion excitation functions and barrier distributions. The model gives a good account of both transfer and fusion data. Received: 2 May 2002 / Accepted: 4 June 2002 / Published online: 26 November 2002 RID="a" ID="a"e-mail: montagnoli@pd.infn.it, Fax +39049 8277102, Tel. +39049 8277117. RID="b" ID="b"On leave from the China Institute for Atomic Energy, 102413 Beijing, China. Communicated by C. Signorini     

Skorokhod decomposition of reflected diffusions on bounded Lipschitz domains with singular non-reflection part

Let be a compact set with interior G. Let L ¹ (G,dx), >0 dx-a.e. on G, and m:=dx. Let A=(a _ij ) be symmetric, and globally uniformly strictly elliptic on G. Let be such that ; f, , is closable in L ² (G,m) with closure ( ^r ,D( ^r )). The latter is fulfilled if satisfies the Hamza type condition, or _i L ¹ _loc (G,dx), 1id. Conservative, non-symmetric diffusion processes X _t related to the extension of a generalized Dirichlet form where satisfies are constructed and analyzed. If G is a bounded Lipschitz domain, H ^1,1 (G), and a _ij D( ^r ), a Skorokhod decomposition for X _t is given. This happens through a local time that is uniquely associated to the smooth measure 1_{{ Tr ()>0}} d, where Tr denotes the trace and the surface measure on G.This research has been financially supported by TMR grant HPMF-CT-2000-00942 of the European Union. Mathematics Subject Classification (2000): 60J60, 60J55, 31C15, 31C25, 35J25     

多组态Dirac-Fock方法与准相对论组态相互作用方法的比较研究

利用多组态Dirac-Fock(MCDF)方法与准相对论组态相互作用方法,分别详细计算了低Z、中Z和高Z原子(离子)各个壳层上电子的束缚能、平均轨道半径、总束缚能、激发能、精细结构能级以及类Ne等电子系列离子的2p53s 1,3P1-2p6 1S0跃迁能,并且对这两种方法的结果进行了数值比较.