Multilevel a posteriori error analysis for reaction–convection–diffusion problems

We present a new approach to the a posteriori error analysis of stable Galerkin approximations of reaction–convection–diffusion problems. It relies upon a non-standard variational formulation of the exact problem, based on the anisotropic wavelet decomposition of the equation residual into convection-dominated scales and diffusion-dominated scales. The associated norm, which is stronger than the standard energy norm, provides a robust (i.e., uniform in the convection limit) control over the streamline derivative of the solution. We propose an upper estimator and a lower estimator of the error, in this norm, between the exact solution and any finite dimensional approximation of it. We investigate the behaviour of such estimators, both theoretically and through numerical experiments. As an output of our analysis, we find that the lower estimator is quantitatively accurate and robust.     

Analysis of finite dimensional approximations to a class of partial differential equations

In this paper, a unified approach for analysing finite dimensional approximations to a class of partial differential equations boundary value problems (second‐kind Fredholm differential equations) is introduced. The approach is shown to be general despite of its extremely simple form. In particular, it is expected to be useful in the convergence analysis of finite element methods for solving PDE problems. Three specific examples are presented to illustrate the broad applicability of the approach. Copyright © 2004 John Wiley & Sons, Ltd.     

The evolution of a truncated gaussian probability density through time—modelling animal liveweights after selection

The form of the probability density derived from the evolution in time of a previously truncated frequency distribution of animal Liveweights is of interest in animal husbandry. Truncated frequency distributions arise when the heavier animals are sold for slaughter and the lighter animals retained. The demands of modern quality assurance schemes require that, given information on animal growth, the farmer is able to estimate the number of animals that would meet the specifications at some time in the future after truncation. Assuming that animal growth can be described by a linear stochastic differential equation, we derive an explicit expression for the probability density of animal Liveweights at any time after the truncation of an initial Gaussian density. It is shown that this probability density converges rapidly to a Gaussian density, so that after about 20 days of typical growth rates for lambs, the resulting density is practically indistinguishable from Gaussian.     

美式期权定价中非局部问题的有限元方法

在本文中 ,我们关心的是美式期权的有限元方法 .首先 ,根据 [4 ]我们对所讨论的问题引进一个新奇的实用的方法 ,它涉及到对原问题重新形成准确的数学公式 ,使得数值解的计算可以在非常小的区域上进行 ,从而该算法计算速度快精度高 .进而 ,我们利用超逼近分析技术得到了有限元解关于 L2 -模的最优估计 .     

DISTRIBUTIONS OF LIMIT CYCLES FOR POLYNOMIAL VECTOR FIELDS (P_2, Q_3)

The distributions of limit cycles of cubic vector fields (P2, Q3) are considered in this paper, where P2 and Q3 are polynomials of x and y of order two and three, respectively. It is possibly seven different distributions of limit cycles given in [1]. We now prove that in which three kinds of distributions are impossible and other four kinds all can be realized by concrete vector fields of (P2,Q3). Some related results are also given.     

Krylov subspaces and the analytic grade

Typical behaviour of the solution of a linear system of equations obtained iteratively by Krylov methods can be characterized by three stages. Initially the residual diminishes steadily; this is followed by stagnation and finally rapid convergence near the algebraic grade. This study examines this behaviour in terms of the concepts of approximately invariant subspace and what we have called the analytic grade of a Krylov sequence. It is shown how the small Ritz values play a vital role in the convergence and how this knowledge helps in the construction of an effective preconditioner. Copyright © 2004 John Wiley & Sons, Ltd.     

Space time-index plots for probing dynamical nonstationarity

We propose a simple method to efficiently probe dynamical nonstationarity in observed time series. In a space time-index plot, the density distributions as a function of normalized time-index are V-shaped due to nonstationarity. We show that this method is workable for short data sets and typical examples are illustrated.     

Compact Matrix Expressions for Generalized Wald Tests of Equality of Moment Vectors

Asymptotic chi-squared test statistics for testing the equality of moment vectors are developed. The test statistics proposed are generalized Wald test statistics that specialize for different settings by inserting an appropriate asymptotic variance matrix of sample moments. Scaled test statistics are also considered for dealing with nonstandard conditions. The specialization will be carried out for testing the equality of multinomial populations, and the equality of variance and correlation matrices for both normal and nonnormal data. When testing the equality of correlation matrices, a scaled version of the normal theory chi-squared statistic is proven to be an asymptotically exact chi-squared statistic in the case of elliptical data.     

The Nonexistence of Procedures with Bounded Performance Characteristics in Certain Parametric Inference Problems

This paper gives a condition which implies the nonexistence of parametric statistical procedures with bounded risk or error performance characteristics. Many examples for which such a condition is satisfied are considered.