Local polynomial fitting under association

We consider the estimation of multivariate regression functions r(x₁,…,x_d) and their partial derivatives up to a total order p1 using high-order local polynomial fitting. The processes {Y_i,X_i} are assumed to be (jointly) associated. Joint asymptotic normality is established for the estimates of the regression function r and all its partial derivatives up to the total order p. Expressions for the bias and variance/covariance matrix (of the asymptotic distribution) are given.     

Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

Linear least squares problems with data over incomplete grids

The paper addresses bivariate surface fitting problems, where data points lie on the vertices of a rectangular grid. Efficient and stable algorithms can be found in the literature to solve such problems. If data values are missing at some grid points, there exists a computational method for finding a least squares spline by fixing appropriate values for the missing data. We extended this technique to arbitrary least squares problems as well as to linear least squares problems with linear equality constraints. Numerical examples are given to show the effectiveness of the technique presented. AMS subject classification (2000)  65D05, 65D07, 65D10, 65F05, 65F20     

几种基于散乱数据拟合的局部插值方法

本文首先针对散乱数据拟合的Shepard方法,结合截断多项式、B样条基函数和指数函数来构造其权函数,使新的权函数具有更高的光滑度和更好的衰减性,并且其光滑性和衰减性可以根据实际需要自由调节,从而提高了曲面的拟合质量．同时还给出一种类似的局部插值方法．另外,本文还基于多重二次插值,结合多元样条的思想,给出了两个局部插值算法．该算法较好地继承了多重二次插值曲面的性质,从而保证了拟合曲面具有好地光顺性和拟合精度．曲面整体也具有较高的光滑性．     

Numerical analysis for the periodic response of nonsmooth dynamic systems of high-dimension

This paper deals with the numerical solution of high dimensional dynamic systems with nonsmooth characteristics, such as structural elastoplasticity, mechanical clearance and dry friction. As a stable periodic response of the system attracts its transient response, it may be extracted from the transient response somehow. The paper suggests a novel iteration scheme to extrapolate the periodic response from a short time history of its transient response by curve fitting. Compared with the current schemes such as shooting and incremental harmonic balancing, the present scheme makes full use of the information within the transient response and the system characters, then its computation efficiency has increased by an order and the numerical convergence depending on the initial iteration has greatly improved.The subject supported by the Chinese Foundation of Aeronautical Science     

灰色Verhulst模型的样条插值函数的残差修正

本用样条函数对灰色Verhulst模型的残差序列进行插值拟合，然后作用于二阶线性微分方程，并以此修正原模型，得到一种新的预测模型的数值解，提高了拟合的精确度。     

A statistical study of the analysis of congested spectra by total spectrum fitting

Heavily overlapped, or congested spectra often display much structure but few individual “lines.” Methods have been devised for analyzing such spectra through nonlinear least-squares fitting of the intensity as a function of wavelength or wavenumber. Such total spectrum fitting (TSF) methods are examined statistically for a simple diatomic model and compared with the standard “measure-assign-fit” (MAF) approach in use since the dawn of spectroscopy. Monte Carlo computations on typically 1000 synthetic spectra at a time verify that the predictions of the variance-covariance matrix are reliable under many circumstances. However in regions where the P and R branches double up, the predicted standard errors in the key spectroscopic constants rise sharply and greatly exceed estimates from the Monte Carlo ensemble statistics. In the same regions, the MAF method actually gives better precision. However, for imperfectly overlapped R and P branches, the MAF standard errors are typically three times larger than the TSF values; moreover, the MAF statistical errors are dwarfed by bias. The TSF approach, while clearly superior in these tests, has a practical drawback: it, too, can give significant bias if the spectra are analyzed with an incorrect model, as illustrated here through calculations employing the wrong function to describe the spectral lineshape.     

Least-Squares Fitting of Algebraic Spline Surfaces

We present an algorithm for fitting implicitly defined algebraic spline surfaces to given scattered data. By simultaneously approximating points and associated normal vectors, we obtain a method which is computationally simple, as the result is obtained by solving a system of linear equations. In addition, the result is geometrically invariant, as no artificial normalization is introduced. The potential applications of the algorithm include the reconstruction of free-form surfaces in reverse engineering. The paper also addresses the generation of exact error bounds, directly from the coefficients of the implicit representation.     

Likelihood estimation of soft-core interaction potentials for Gibbsian point patterns

The likelihood method is developed for the analysis of socalled regular point patterns. Approximating the normalizing factor of Gibbs canonical distribution, we simultaneously estimate two parameters, one for the scale and the other which measures the softness (or hardness), of repulsive interactions between points. The approximations are useful up to a considerably high density. Some real data are analyzed to illustrate the utility of the parameters for characterizing the regular point pattern.     

人工神经网络及其在分析化学中的应用

人工神经网络是一种新兴的计算方法，有着广阔的发展前途，目前在分析化学领域已经有了多方面的应用。本文简要介绍了人工神经网络的原理及其在分析化学中的应用。