Cyclic reduction and FACR methods for piecewise hermite bicubic orthogonal spline collocation

Cyclic reduction and Fourier analysis-cyclic reduction (FACR) methods are presented for the solution of the linear systems which arise when orthogonal spline collocation with piecewise Hermite bicubics is applied to boundary value problems for certain separable partial differential equations on a rectangle. On anN×N uniform partition, the cyclic reduction and Fourier analysis-cyclic reduction methods requireO(N ²log₂ N) andO(N ²log₂log₂ N) arithmetic operations, respectively.     

用频域法实现圆环形目标物的计数

研究了用频域法实现圆环形目标物计数的方法。对空域中含有圆环形目标物的二值图像，经傅里叶变换后在频域内进行滤波，再经傅里叶反变换后取一阈值进行处理，则对应于圆环形目标物即变为一个实心体，然后对实心体进行计数即可实现对圆环形目标物的计数。对该方法进行了数学描述，并分别用模拟图像和实际图像进行了实验。结果表明，该方法消除噪声的效果好，对部分封闭圆的检测同样有效，且可用于目标物有粘连的情况。该方法对圆环形目标物的尺寸相差不大的情况更为有效。     

The application of spectral distribution of product of two random matrices in the factor analysis

In the factor analysis model with large cross-section and time-series dimensions,we pro- pose a new method to estimate the number of factors.Specially if the idiosyncratic terms satisfy a linear time series model,the estimators of the parameters can be obtained in the time series model. The theoretical properties of the estimators are also explored.A simulation study and an empirical analysis are conducted.     

关于多元幂级数

引入了多元函数项级数的概念,给出了其收敛域及和函数的定义;通过详实的例子讨论了多元幂级数的收敛域、和函数及多元函数展开为多元幂级数的计算方法.     

Generators of some Ramanujan formulas

In this paper we prove some Ramanujan type formulas for 1/π but without using the theory of modular forms. Instead we use the WZ—method created by H. Wilf and D. Zeilberger and find some hypergeometric functions in two variables which are second components of WZ—pairs than can be certified using Zeilberger's EKHAD package. These certificates have an additional property which allows us to get generalized Ramanujan's type series which are routinely proven by computer. We call these second hypergeometric components of the WZ—pairs generators. Finding generators seems a hard task but using a kind of experimental research (explained below), we have succeeded in finding some of them. Unfortunately we have not found yet generators for the most impressive Ramanujan's formulas. We also prove some interesting binomial sums for the constant 1/π². Finally we rewrite many of the obtained series using pochhammer symbols and study the rate of convergence. 2000 Mathematics Subject Classification Primary—33C20     

Sparse generalized Fourier transforms

Block-diagonalization of sparse equivariant discretization matrices is studied. Such matrices typically arise when partial differential equations that evolve in symmetric geometries are discretized via the finite element method or via finite differences. By considering sparse equivariant matrices as equivariant graphs, we identify a condition for when block-diagonalization via a sparse variant of a generalized Fourier transform (GFT) becomes particularly simple and fast. Characterizations for finite element triangulations of a symmetric domain are given, and formulas for assembling the block-diagonalized matrix directly are presented. It is emphasized that the GFT preserves symmetric (Hermitian) properties of an equivariant matrix. By simulating the heat equation at the surface of a sphere discretized by an icosahedral grid, it is demonstrated that the block-diagonalization is beneficial. The gain is significant for a direct method, and modest for an iterative method. A comparison with a block-diagonalization approach based upon the continuous formulation is made. It is found that the sparse GFT method is an appropriate way to discretize the resulting continuous subsystems, since the spectrum and the symmetry are preserved. AMS subject classification (2000)  43A30, 65T99, 20B25     

Quadratic transformations and Guillera’s formulas for 1/π2

We prove two new series of Ramanujan type for 1/π².     

Optimal reconstruction of the solution of the wave equation from inaccurate initial data

In the present paper, we consider the problem of the optimal reconstruction of the solution of the wave equation from the approximate values of the Fourier coefficients of the function specifying the initial form of the string. For an operator defined on the weight space of vectors from l ₂, we present the solution of the more general problem of reconstruction from the approximate values of the coordinates of these vectors.     

On inequalities of Lieb-Thirring type

Applying the method proposed by Kashin for proving inequalities of Lieb-Thirring type for orthonormal systems, we prove a similar inequality in the multidimensional case.     

On multiple integrals of special form

Multiple integrals generalizing the iterated kernels of integral operators are expressed as single integrals in the case of a special representation of the kernel (this is our theorem). Besides integral equations, Markov processes involve these integrals as well. As a consequence of the theorem, we obtain transition probability densities of certain Markov processes. As an illustration, we consider nine examples.