A second-order divided difference filter (SDDF) is derived for integrating line of sight measurement from vision sensor with acceleration and angular rate measurements of the follower to estimate the precise relative ...A second-order divided difference filter (SDDF) is derived for integrating line of sight measurement from vision sensor with acceleration and angular rate measurements of the follower to estimate the precise relative position,velocity and attitude of two unmanned aerial vehicles (UAVs).The second-order divided difference filter which makes use of multidimensional interpolation formulations to approximate the nonlinear transformations could achieve more accurate estimation and faster convergence from inaccurate initial conditions than standard extended Kalman filter.The filter formulation is based on relative motion equations.The global attitude parameterization is given by quarternion,while a generalized three-dimensional attitude representation is used to define the local attitude error.Simulation results are shown to compare the performance of the second-order divided difference filter with a standard extended Kalman filter approach.展开更多
A square-root version of the divided difference Rauch-Tung-Striebel (RTS) smoother is proposed in this paper. The square-root variant essentially propagates the square roots of the covariance matrices and can consiste...A square-root version of the divided difference Rauch-Tung-Striebel (RTS) smoother is proposed in this paper. The square-root variant essentially propagates the square roots of the covariance matrices and can consistently improve the numerical stability because all the resulting covariance matrices are guaranteed to stay positive semi-definite. Furthermore, the square-root form ensures reliable implementation in an embedded system with fixed or limited precision although it is algebraically equivalent to the standard form. The new smoothing algorithm is tested in a challenging two-dimensional maneuvering target tracking problem with unknown and time-varying turn rate, and its performance is compared with that of other de-facto standard filters and smoothers. The simulation results indicate that the proposed RTS smoother markedly outperforms the associated filters and gives slightly smaller error than an unscented-based RTS smoother.展开更多
A hierarchical robot control is proposed for robot soccer system. The Newton’s divided difference is utilized in robot path planning. This paper describes the problems encoutered, software design considerations, visi...A hierarchical robot control is proposed for robot soccer system. The Newton’s divided difference is utilized in robot path planning. This paper describes the problems encoutered, software design considerations, vision algorithm and controls of individual robots. The solutions to the problems implemented are simple and direct. It is observed that many of the ideas and solutions can be evolved based on simple theories and concepts. This paper focuses on software structure of multi agent controls, vision algorithm and simple path planning method.展开更多
The efficient and accurate approximate nonlinear filters have been widely used in the estimation of states and parameters of dynamical systems. In this paper, an adaptive divided difference filter is designed for prec...The efficient and accurate approximate nonlinear filters have been widely used in the estimation of states and parameters of dynamical systems. In this paper, an adaptive divided difference filter is designed for precise estimation of states and parameters of micromechanical gyro navigation system. Based on the investigation of nonlinear divided difference filter the adaptive divided difference filter(ADDF) was designed, which takes account of the incorrect time-varying noise statistics of dynamical systems and compensation of the nonlinearity effects neglected by linearization. And its performance is superior to that of DDF and extended Kalman filter(EKF). Simulation results indicate that the advantages of the proposed nonlinear filters make them attractive alternatives to the extended Kalman filter.展开更多
The aim of this work is to construct a new quadrature formula for Fourier Chebyshev coef ficients based on the divided differences of the integrand at points 1, 1 and the zeros of the n th Chebyshev polynomial o...The aim of this work is to construct a new quadrature formula for Fourier Chebyshev coef ficients based on the divided differences of the integrand at points 1, 1 and the zeros of the n th Chebyshev polynomial of the second kind. The interesting thing is that this quadrature rule is closely related to the well known Gauss Turn quadrature formula and similar to a recent result of Micchelli and Sharma, extending a particular case due to Micchelli and Rivlin.展开更多
New sigma point filtering algorithms, including the unscented Kalman filter (UKF) and the divided difference filter (DDF), are designed to solve the nonlinear filtering problem under the condition of correlated no...New sigma point filtering algorithms, including the unscented Kalman filter (UKF) and the divided difference filter (DDF), are designed to solve the nonlinear filtering problem under the condition of correlated noises. Based on the minimum mean square error estimation theory, the nonlinear optimal predictive and correction recursive formulas under the hypothesis that the input noise is correlated with the measurement noise are derived and can be described in a unified framework. Then, UKF and DDF with correlated noises are proposed on the basis of approximation of the posterior mean and covariance in the unified framework by using unscented transformation and second order Stirling's interpolation. The proposed UKF and DDF with correlated noises break through the limitation that input noise and measurement noise must be assumed to be uneorrelated in standard UKF and DDF. Two simulation examples show the effectiveness and feasibility of new algorithms for dealing with nonlinear filtering issue with correlated noises.展开更多
A new convergence theorem for the Secant method in Banach spaces based on new recurrence relations is established for approximating a solution of a nonlinear operator equation. It is assumed that the divided differenc...A new convergence theorem for the Secant method in Banach spaces based on new recurrence relations is established for approximating a solution of a nonlinear operator equation. It is assumed that the divided difference of order one of the nonlinear operator is Lipschitz continuous. The convergence conditions differ from some existing ones and are easily satisfied. The results of the paper are justified by numerical examples that cannot be handled by earlier works.展开更多
Based on the multiquadric trigonometric B-spline quasi-interpolant, this paper proposes a meshless scheme for some partial differential equations whose solutions are periodic with respect to the spatial variable. This...Based on the multiquadric trigonometric B-spline quasi-interpolant, this paper proposes a meshless scheme for some partial differential equations whose solutions are periodic with respect to the spatial variable. This scheme takes into ac- count the periodicity of the analytic solution by using derivatives of a periodic quasi-interpolant (multiquadric trigonometric B-spline quasi-interpolant) to approximate the spatial derivatives of the equations. Thus, it overcomes the difficulties of the previous schemes based on quasi-interpolation (requiring some additional boundary conditions and yielding unwanted high-order discontinuous points at the boundaries in the spatial domain). Moreover, the scheme also overcomes the dif- ficulty of the meshless collocation methods (i.e., yielding a notorious ill-conditioned linear system of equations for large collocation points). The numerical examples that are presented at the end of the paper show that the scheme provides excellent approximations to the analytic solutions.展开更多
Interpolatory subdivision algorithms for the generation of curves and surfaces play a veryimportant rule in shape design and modelling in CAD/CAM systems. In this paper, by using the dif-ference and divided difference...Interpolatory subdivision algorithms for the generation of curves and surfaces play a veryimportant rule in shape design and modelling in CAD/CAM systems. In this paper, by using the dif-ference and divided difference analysis, a systematic method to construct Cn (n≥ 0) interpolatorycurves by subdivision from given data is described and the mask (filter) of the algorithm is presentedexplicitly. This algorithm generates a Cn smooth curve which interpolates the initial control points.Control parameters are also provided so that the shape of the final curve can be adjusted according torequirements. An immediate generalisation of the method is the construction of smooth interpolatorysubdivision algorithms over uniform triangular networks (tensor product type data) in Rm. The mainresults of this algorithm for smooth interpolatory surface subdivision algorrthm are also included.AMS(MOS) : 65D05 , 65D15 , 65D17.展开更多
The n-divided difference of the composite function h := f o g of functions f, g at a group of nodes t0,t1,…,tn is shown by the combinations of divided differences of f at the group of nodes g(t0),g(t1),…,g(tm...The n-divided difference of the composite function h := f o g of functions f, g at a group of nodes t0,t1,…,tn is shown by the combinations of divided differences of f at the group of nodes g(t0),g(t1),…,g(tm) and divided differences of g at several partial group of nodes t0,t1,…,tn, where m = 1,2,…,n. Especially, when the given group of nodes are equal to each other completely, it will lead to Faà di Bruno's formula of higher derivatives of function h.展开更多
The remainder estimates of numerical divided difference formula are given for the functions of lower and higher smoothness, respectively. Then several divided difference formulas with super-convergence are derived wit...The remainder estimates of numerical divided difference formula are given for the functions of lower and higher smoothness, respectively. Then several divided difference formulas with super-convergence are derived with their remainder expressions.展开更多
The growth of meromorphic solutions of linear difference equations containing Askey-Wilson divided difference operators is estimated.Theφ-order is used as a general growth indicator,which covers the growth spectrum b...The growth of meromorphic solutions of linear difference equations containing Askey-Wilson divided difference operators is estimated.Theφ-order is used as a general growth indicator,which covers the growth spectrum between the logarithmic orderρlog(f)and the classical orderρ(f)of a meromorphic function f.展开更多
In this paper, we consider the higher divided difference of a composite function f(g(t)) in which g(t) is an s-dimensional vector. By exploiting some properties from mixed partial divided differences and multiva...In this paper, we consider the higher divided difference of a composite function f(g(t)) in which g(t) is an s-dimensional vector. By exploiting some properties from mixed partial divided differences and multivariate Newton interpolation, we generalize the divided difference form of Faà di Bruno's formula with a scalar argument. Moreover, a generalized Faà di Bruno's formula with a vector argument is derived.展开更多
In the paper, necessary and sufficient conditions are provided for a function involving the divided difference of two psi functions to be completely monotonic. Consequently, a class of inequalities for sums are presen...In the paper, necessary and sufficient conditions are provided for a function involving the divided difference of two psi functions to be completely monotonic. Consequently, a class of inequalities for sums are presented, the logarithmically complete monotonicity of a function involving the ratio of two gamma functions are derived, and two double inequalities for bounding the ratio of two gamma functions are discovered.展开更多
Newton's polynomial interpolation may be the favourite linear interpolation in the sense that it is built up by means of the divided differences which can be calculated recursively and produce useful intermediate res...Newton's polynomial interpolation may be the favourite linear interpolation in the sense that it is built up by means of the divided differences which can be calculated recursively and produce useful intermediate results. However Newton interpolation is in fact point based interpolation since a new interpolating polynomial with one more degree is obtained by adding a new support point into the current set of support points once at a time. In this paper we extend the point based interpolation to the block based interpolation. Inspired by the idea of the modern architectural design, we first divide the original set of support points into some subsets (blocks), then construct each block by using whatever interpolation means, linear or rational and finally assemble these blocks by Newton's method to shape the whole interpolation scheme. Clearly our method offers many flexible interpolation schemes for choices which include the classical Newton's polynomial interpolation as its special case. A bivariate analogy is also discussed and numerical examples are given to show the effectiveness of our method.展开更多
Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolatio...Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolation error of anisotropic elements based on Newton's formula of polynomial interpolation as well as its applications.展开更多
We present some matrix and determinant identities for the divided differences of the composite functions, which generalize the divided difference form of Faa di Bruno's formula. Some recent published identities can b...We present some matrix and determinant identities for the divided differences of the composite functions, which generalize the divided difference form of Faa di Bruno's formula. Some recent published identities can be regarded as special cases of our results.展开更多
By applying a Grobner-Shirshov basis of the symmetric group Sn, we give two formulas for Schubert polynomials, either of which involves only nonnegative monomials. We also prove some combinatorial properties of Schube...By applying a Grobner-Shirshov basis of the symmetric group Sn, we give two formulas for Schubert polynomials, either of which involves only nonnegative monomials. We also prove some combinatorial properties of Schubert polynomials. As applications, we give two algorithms to calculate the structure constants for Schubert polynomials, one of which depends on Monk's formula.展开更多
基金Sponsored by the Aerospace Technology Innovation Funding(Grant No. CASC0209)
文摘A second-order divided difference filter (SDDF) is derived for integrating line of sight measurement from vision sensor with acceleration and angular rate measurements of the follower to estimate the precise relative position,velocity and attitude of two unmanned aerial vehicles (UAVs).The second-order divided difference filter which makes use of multidimensional interpolation formulations to approximate the nonlinear transformations could achieve more accurate estimation and faster convergence from inaccurate initial conditions than standard extended Kalman filter.The filter formulation is based on relative motion equations.The global attitude parameterization is given by quarternion,while a generalized three-dimensional attitude representation is used to define the local attitude error.Simulation results are shown to compare the performance of the second-order divided difference filter with a standard extended Kalman filter approach.
基金the Fundamental Research Fund of Northwestern Polytechnical University( Grant No. JC20120210,JC20110238)
文摘A square-root version of the divided difference Rauch-Tung-Striebel (RTS) smoother is proposed in this paper. The square-root variant essentially propagates the square roots of the covariance matrices and can consistently improve the numerical stability because all the resulting covariance matrices are guaranteed to stay positive semi-definite. Furthermore, the square-root form ensures reliable implementation in an embedded system with fixed or limited precision although it is algebraically equivalent to the standard form. The new smoothing algorithm is tested in a challenging two-dimensional maneuvering target tracking problem with unknown and time-varying turn rate, and its performance is compared with that of other de-facto standard filters and smoothers. The simulation results indicate that the proposed RTS smoother markedly outperforms the associated filters and gives slightly smaller error than an unscented-based RTS smoother.
文摘A hierarchical robot control is proposed for robot soccer system. The Newton’s divided difference is utilized in robot path planning. This paper describes the problems encoutered, software design considerations, vision algorithm and controls of individual robots. The solutions to the problems implemented are simple and direct. It is observed that many of the ideas and solutions can be evolved based on simple theories and concepts. This paper focuses on software structure of multi agent controls, vision algorithm and simple path planning method.
文摘The efficient and accurate approximate nonlinear filters have been widely used in the estimation of states and parameters of dynamical systems. In this paper, an adaptive divided difference filter is designed for precise estimation of states and parameters of micromechanical gyro navigation system. Based on the investigation of nonlinear divided difference filter the adaptive divided difference filter(ADDF) was designed, which takes account of the incorrect time-varying noise statistics of dynamical systems and compensation of the nonlinearity effects neglected by linearization. And its performance is superior to that of DDF and extended Kalman filter(EKF). Simulation results indicate that the advantages of the proposed nonlinear filters make them attractive alternatives to the extended Kalman filter.
文摘The aim of this work is to construct a new quadrature formula for Fourier Chebyshev coef ficients based on the divided differences of the integrand at points 1, 1 and the zeros of the n th Chebyshev polynomial of the second kind. The interesting thing is that this quadrature rule is closely related to the well known Gauss Turn quadrature formula and similar to a recent result of Micchelli and Sharma, extending a particular case due to Micchelli and Rivlin.
基金Projects(61135001, 61075029, 61074155) supported by the National Natural Science Foundation of ChinaProject(20110491690) supported by the Postdocteral Science Foundation of China
文摘New sigma point filtering algorithms, including the unscented Kalman filter (UKF) and the divided difference filter (DDF), are designed to solve the nonlinear filtering problem under the condition of correlated noises. Based on the minimum mean square error estimation theory, the nonlinear optimal predictive and correction recursive formulas under the hypothesis that the input noise is correlated with the measurement noise are derived and can be described in a unified framework. Then, UKF and DDF with correlated noises are proposed on the basis of approximation of the posterior mean and covariance in the unified framework by using unscented transformation and second order Stirling's interpolation. The proposed UKF and DDF with correlated noises break through the limitation that input noise and measurement noise must be assumed to be uneorrelated in standard UKF and DDF. Two simulation examples show the effectiveness and feasibility of new algorithms for dealing with nonlinear filtering issue with correlated noises.
基金Supported by the National Natural Science Foundation of China (10871178)the Natural Science Foundation of Zhejiang Province of China (Y606154)Foundation of the Education Department of Zhejiang Province of China (20071362)
文摘A new convergence theorem for the Secant method in Banach spaces based on new recurrence relations is established for approximating a solution of a nonlinear operator equation. It is assumed that the divided difference of order one of the nonlinear operator is Lipschitz continuous. The convergence conditions differ from some existing ones and are easily satisfied. The results of the paper are justified by numerical examples that cannot be handled by earlier works.
基金supported by the Shanghai Guidance of Science and Technology,China(Grant No.12DZ2272800)the Natural Science Foundation of Education Department of Anhui Province,China(Grant No.KJ2013B203)the Foundation of Introducing Leaders of Science and Technology of Anhui University,China(Grant No.J10117700057)
文摘Based on the multiquadric trigonometric B-spline quasi-interpolant, this paper proposes a meshless scheme for some partial differential equations whose solutions are periodic with respect to the spatial variable. This scheme takes into ac- count the periodicity of the analytic solution by using derivatives of a periodic quasi-interpolant (multiquadric trigonometric B-spline quasi-interpolant) to approximate the spatial derivatives of the equations. Thus, it overcomes the difficulties of the previous schemes based on quasi-interpolation (requiring some additional boundary conditions and yielding unwanted high-order discontinuous points at the boundaries in the spatial domain). Moreover, the scheme also overcomes the dif- ficulty of the meshless collocation methods (i.e., yielding a notorious ill-conditioned linear system of equations for large collocation points). The numerical examples that are presented at the end of the paper show that the scheme provides excellent approximations to the analytic solutions.
文摘Interpolatory subdivision algorithms for the generation of curves and surfaces play a veryimportant rule in shape design and modelling in CAD/CAM systems. In this paper, by using the dif-ference and divided difference analysis, a systematic method to construct Cn (n≥ 0) interpolatorycurves by subdivision from given data is described and the mask (filter) of the algorithm is presentedexplicitly. This algorithm generates a Cn smooth curve which interpolates the initial control points.Control parameters are also provided so that the shape of the final curve can be adjusted according torequirements. An immediate generalisation of the method is the construction of smooth interpolatorysubdivision algorithms over uniform triangular networks (tensor product type data) in Rm. The mainresults of this algorithm for smooth interpolatory surface subdivision algorrthm are also included.AMS(MOS) : 65D05 , 65D15 , 65D17.
基金This work was supported by the National Science Foundation of China (Grant No.10471128).
文摘The n-divided difference of the composite function h := f o g of functions f, g at a group of nodes t0,t1,…,tn is shown by the combinations of divided differences of f at the group of nodes g(t0),g(t1),…,g(tm) and divided differences of g at several partial group of nodes t0,t1,…,tn, where m = 1,2,…,n. Especially, when the given group of nodes are equal to each other completely, it will lead to Faà di Bruno's formula of higher derivatives of function h.
基金supported by the National Natural Science Foundation of China(Grant No.10471128).
文摘The remainder estimates of numerical divided difference formula are given for the functions of lower and higher smoothness, respectively. Then several divided difference formulas with super-convergence are derived with their remainder expressions.
基金the support of the China Scholarship Council(Grant No.201806330120)supported by National Natural Science Foundation of China(Grant No.11771090)+1 种基金supported by the National Natural Science Foundation of China(Grants Nos.11971288 and 11771090)Shantou University SRFT(Grant No.NTF18029)。
文摘The growth of meromorphic solutions of linear difference equations containing Askey-Wilson divided difference operators is estimated.Theφ-order is used as a general growth indicator,which covers the growth spectrum between the logarithmic orderρlog(f)and the classical orderρ(f)of a meromorphic function f.
基金Acknowledgments. This work was supported by the National Science Foundation of China (Grant Nos. 10471128, 10731060).
文摘In this paper, we consider the higher divided difference of a composite function f(g(t)) in which g(t) is an s-dimensional vector. By exploiting some properties from mixed partial divided differences and multivariate Newton interpolation, we generalize the divided difference form of Faà di Bruno's formula with a scalar argument. Moreover, a generalized Faà di Bruno's formula with a vector argument is derived.
基金supported partially by the China Scholarship Council and the Science Foundation of Tianjin Polytechnic Universitysupported in part by the Natural Science Foundation Project of Chongqing,China(Grant No.CSTC2011JJA00024)+1 种基金the Research Project of Science and Technology of Chongqing Education Commission,China(Grant No.KJ120625)the Fund of Chongqing Normal University,China(Grant Nos.10XLR017 and 2011XLZ07)
文摘In the paper, necessary and sufficient conditions are provided for a function involving the divided difference of two psi functions to be completely monotonic. Consequently, a class of inequalities for sums are presented, the logarithmically complete monotonicity of a function involving the ratio of two gamma functions are derived, and two double inequalities for bounding the ratio of two gamma functions are discovered.
基金Project supported by the National Natural Science Foundation of China under Grant No.10171026 and No.60473114, and the Anhui Provincial Natural Science Foundation, China under Grant No.03046102.
文摘Newton's polynomial interpolation may be the favourite linear interpolation in the sense that it is built up by means of the divided differences which can be calculated recursively and produce useful intermediate results. However Newton interpolation is in fact point based interpolation since a new interpolating polynomial with one more degree is obtained by adding a new support point into the current set of support points once at a time. In this paper we extend the point based interpolation to the block based interpolation. Inspired by the idea of the modern architectural design, we first divide the original set of support points into some subsets (blocks), then construct each block by using whatever interpolation means, linear or rational and finally assemble these blocks by Newton's method to shape the whole interpolation scheme. Clearly our method offers many flexible interpolation schemes for choices which include the classical Newton's polynomial interpolation as its special case. A bivariate analogy is also discussed and numerical examples are given to show the effectiveness of our method.
基金the National Nutural Science Foundation of China(Grant Nos.10771198,10590353)
文摘Interpolation theory is the foundation of finite element methods.In this paper,after reviewing some existed interpolation theorems of anisotropic finite element methods,we present a new way to analyse the interpolation error of anisotropic elements based on Newton's formula of polynomial interpolation as well as its applications.
文摘We present some matrix and determinant identities for the divided differences of the composite functions, which generalize the divided difference form of Faa di Bruno's formula. Some recent published identities can be regarded as special cases of our results.
文摘By applying a Grobner-Shirshov basis of the symmetric group Sn, we give two formulas for Schubert polynomials, either of which involves only nonnegative monomials. We also prove some combinatorial properties of Schubert polynomials. As applications, we give two algorithms to calculate the structure constants for Schubert polynomials, one of which depends on Monk's formula.