求解垂直互补问题的参数牛顿法

给出了求解垂直互补问题的一种参数牛顿法,在较为温和的条件下证明了该方法的局部超线性收敛结果,并且给出了具体数值计算.     

Error estimates in S P for cubature formulas exact for haar polynomials in the two-dimensional case

On the spaces S _p , an upper estimate is found for the norm of the error functional δ _N (f) of cubature formulas possessing the Haar d-property in the two-dimensional case. An asymptotic relation is proved for $ \left\| {\delta _N (f)} \right\|_{S_p^* } On the spaces S _p , an upper estimate is found for the norm of the error functional δ _N (f) of cubature formulas possessing the Haar d-property in the two-dimensional case. An asymptotic relation is proved for with the number of nodes N ∼ 2 ^d , where d → ∞. For N ∼ 2 ^d with d → ∞, it is shown that the norm of δ _N for the formulas under study has the best convergence rate, which is equal to N ^−1/p . Original Russian Text ? K.A. Kirillov, M.V. Noskov, 2009, published in Zhurnal Vychislitel’noi Matematiki i Matematicheskoi Fiziki, 2009, Vol. 49, No. 1, pp. 3–13.     

High-Order Leap-Frog Based Discontinuous Galerkin Method for the Time-Domain Maxwell Equations on Non-Conforming Simplicial Meshes

A high-order leap-flog based non-dissipative discontinuous Galerkin timedomain method for solving Maxwell's equations is introduced and analyzed. The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements, with a Nth-order leap-frog time scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with highorder elements show the potential of the method.     

On the Caginalp system with dynamic boundary conditions and singular potentials

This article is devoted to the study of the Caginalp phase field system with dynamic boundary conditions and singular potentials. We first show that, for initial data in H ², the solutions are strictly separated from the singularities of the potential. This turns out to be our main argument in the proof of the existence and uniqueness of solutions. We then prove the existence of global attractors. In the last part of the article, we adapt well-known results concerning the Lojasiewicz inequality in order to prove the convergence of solutions to steady states.      

MATLAB GUI设计在几何作图中的应用

本文考虑MATLAB软件在几何作图中的应用.以在教学中遇到的两个几何作图为例,介绍了相应GUI设计界面的制作和演示过程.通过演示,得到了两个常微分方程的解的几何形状及其形成过程,从而对所得结论有更加直观和形象的认识,也加深了学生对问题本身的理解,进一步提高了教学效果.     

Riemann—Lebesgue定理及其应用

本文利用实变函数方法详细证明了Riemann-Lebesgue定理在开集上成立,并且给出了其应用.     

Asymptotic Convergence of Conjugate Gradient Methods for the Partial Symmetric Eigenproblem

Recently an efficient method (DACG) for the partial solution of the symmetric generalized eigenproblem A x = δB x has been developed, based on the conjugate gradient (CG) minimization of the Rayleigh quotient over successive deflated subspaces of decreasing size. The present paper provides a numerical analysis of the asymptotic convergence rate ρ_j of DACG in the calculation of the eigenpair λ_j, u _j, when the scheme is preconditioned with A⁻¹. It is shown that, when the search direction are A-conjugate, ρ_j is well approximated by 4/ξ_j, where ξ_j is the Hessian condition number of a Rayleigh quotient defined in appropriate oblique complements of the space spanned by the leftmost eigenvectors u ₁, u ₂,…, u _j−1 already calculated. It is also shown that 1/ξ_j is equal to the relative separation between the eigenvalue λ_j currently sought and the next higher one λ_j+1 and the next higher one λ_{j + 1}. A modification of DACG (MDACG) is studied, which involves a new set of CG search directions which are made M-conjugate, with M-conjugate, with M-conjugate, with M a matrix approximating the Hessian. By distinction, MDACG has an asymptotic rate of convergence which appears to be inversely proportional to the square root of ξ_j, in complete agreement with the theoretical results known for the CG solution to linear systems. © 1997 by John Wiley & Sons, Ltd.     

A nearly optimal preconditioning based on recursive red–black orderings

Considering matrices obtained by the application of a five-point stencil on a 2D rectangular grid, we analyse a preconditioning method introduced by Axelsson and Eijkhout, and by Brand and Heinemann. In this method, one performs a (modified) incomplete factorization with respect to a so-called ‘repeated’ or ‘recursive’ red–black ordering of the unknowns while fill-in is accepted provided that the red unknowns in a same level remain uncoupled. Considering discrete second order elliptic PDEs with isotropic coefficients, we show that the condition number is bounded by 𝒪(n ${}^{\text{½ log}{}_{\text{2}}\text{(√(5) −1)}}$ ) where n is the total number of unknowns (½ log₂(√(5) − 1) = 0.153), and thus, that the total arithmetic work for the solution is bounded by 𝒪(n^1.077). Our condition number estimate, which turns out to be better than standard 𝒪(log² n) estimates for any realistic problem size, is purely algebraic and holds in the presence of Neumann boundary conditions and/or discontinuities in the PDE coefficients. Numerical tests are reported, displaying the efficiency of the method and the relevance of our analysis. © 1997 John Wiley & Sons, Ltd.     

Three limit theorems for vacancy in multivariate coverage problems

Three limit theorems describing asymptotic distribution of vacancy in general multivariate coverage problems are proved, in which nk-dimensional spheres are distributed within a k-dimensional unit cube according to a density f. The first result (a central limit theorem) describes the case where the proportion of vacancy converges to a fixed constant lying between 0 and 1. The last two results treat the case where the proportion of vacancy tends to 1 as n → ∞. Results of this nature have hitherto been available only for restricted k and/or for f equal to the uniform density.