Improvement of Orbit Prediction Algorithm for Spacecraft Through Simplified Precession-Nutation Model Using Cubic Spline Interpolation Method

For the on-orbit flight missions,the model of orbit prediction is critical for the tasks with high accuracy requirement and limited computing resources of spacecraft.The precession-nutation model,as the main part of extended orbit prediction,affects the efficiency and accuracy of on-board operation.In this paper,the previous research about the conversion between the Geocentric Celestial Reference System and International Terrestrial Reference System is briefly summarized,and a practical concise precession-nutation model is proposed for coordinate transformation computation based on Celestial Intermediate Pole(CIP).The idea that simplifying the CIP-based model with interpolation method is driven by characteristics of precession-nutation parameters changing with time.A cubic spline interpolation algorithm is applied to obtain the required CIP coordinates and Celestial Intermediate Origin locator.The complete precession nutation model containing more than 4000 parameters is simplified to the calculation of a cubic polynomial,which greatly reduces the computational load.In addition,for evaluating the actual performance,an orbit propagator is built with the proposed simplified precession-nutationmodel.Compared with the orbit prediction results obtained by the truncated series of IAU2000/2006 precession-nutation model,the simplified precession-nutation model with cubic spline interpolation can significantly improve the accuracy of orbit prediction,which implicates great practical application value in further on-orbit missions of spacecraft.

Several Ways to Calculate the Universal Gravitational Constant <i>G</i>Theoretically and Cubic Splines to Verify Its Measured Value

In 1686, Newton discovered the laws of gravitation [1] and predicted the universal gravitational constant . In 1798, with a torsion balance, Cavendish [2] measured . Due to the low intensity of gravitation, it is difficult to obtain reliable results because they are disturbed by surrounding masses and environmental phenomena. Modern physics is unable to link G with other constants. However, in a 2019 article [3], with a new cosmological model, we showed that G seams related to other constants, and we obtained a theoretical value of . Here, we want to show that our theoretical value of G is the right one by interpreting measurements of G with the help of a new technique using cubic splines. We make the hypothesis that most G measurements are affected by an unknown systematic error which creates two main groups of data. We obtain a measured value of . Knowing that our theoretical value of G is in agreement with the measured value, we want to establish a direct link between G and as many other constants as possible to show, with 33 equations, that G is probably linked with most constants in the universe. These equations may be useful for astrophysicists who work in this domain. Since we have been able to link G with Hubble parameter H0 (which evolve since its reverse gives the apparent age of the universe), we deduce that G is likely not truly constant. It’s value probably slowly varies in time and space. However, at our location in the universe and for a relatively short period, this parameter may seem constant.

Dynamic GM（1,1） Model Based on Cubic Spline for Electricity Consumption Prediction in Smart Grid

Electricity demand forecasting plays an important role in smart grid expansion planning.In this paper,we present a dynamic GM(1,1) model based on grey system theory and cubic spline function interpolation principle.Using piecewise polynomial interpolation thought,this model can dynamically predict the general trend of time series data.Combined with low-order polynomial,the cubic spline interpolation has smaller error,avoids the Runge phenomenon of high-order polynomial,and has better approximation effect.Meanwhile,prediction is implemented with the newest information according to the rolling and feedback mechanism and fluctuating error is controlled well to improve prediction accuracy in time-varying environment.Case study using the living electricity consumption data of Jiangsu province in 2008 is presented to demonstrate the effectiveness of the proposed model.

A comparison of piecewise cubic Hermite interpolating polynomials,cubic splines and piecewise linear functions for the approximation of projectile aerodynamics

Modelling and simulation of projectile flight is at the core of ballistic computer software and is essential to the study of performance of rifles and projectiles in various engagement conditions.An effective and representative numerical model of projectile flight requires a relatively good approximation of the aerodynamics.The aerodynamic coefficients of the projectile model should be described as a series of piecewise polynomial functions of the Mach number that ideally meet the following conditions:they are continuous,differentiable at least once,and have a relatively low degree.The paper provides the steps needed to generate such piecewise polynomial functions using readily available tools,and then compares Piecewise Cubic Hermite Interpolating Polynomial(PCHIP),cubic splines,and piecewise linear functions,and their variant,as potential curve fitting methods to approximate the aerodynamics of a generic small arms projectile.A key contribution of the paper is the application of PCHIP to the approximation of projectile aerodynamics,and its evaluation against a set of criteria.Finally,the paper provides a baseline assessment of the impact of the polynomial functions on flight trajectory predictions obtained with 6-degree-of-freedom simulations of a generic projectile.

Displacement Response Reconstruction of Slender Flexible Structures Based on Cubic Spline Fitting Method

How to reconstruct a dynamic displacement of slender flexible structures is the key technology to develop smart structures and structural health monitoring(SHM), which are beneficial for controlling the structural vibration and protecting the structural safety. In this paper, the displacement reconstruction method based on cubic spline fitting is put forward to reconstruct the dynamic displacement of slender flexible structures without the knowledge of modeshapes and applied loading. The obtained strains and displacements are compared with the results calculated by ABAQUS to check the method's validity. It can be found that the proposed method can accurately identify the strains and displacement of slender flexible structures undergoing linear vibrations, nonlinear vibrations, and parametric vibrations. Under the concentrated force, the strains of slender flexible structures will change suddenly along the axial direction. With locally densified measurement points, the present reconstruction method still works well for the strain concentration problem.

A Finite Element Cable Model and Its Applications Based on the Cubic Spline Curve

For accurate prediction of the deformation of cable in the towed system,a new finite element model is presented that provides a representation of both the bending and torsional effects.In this paper,the cubic spline interpolation function is applied as the trial solution.By using a weighted residual approach,the discretized motion equations for the new finite element model are developed.The model is calculated with the computation program complier by Matlab.Several numerical examples are presented to illustrate the numerical schemes.The results of numerical simulation are stable and valid,and consistent with the mechanical properties of the cable.The model can be applied to kinematics analysis and the design of ocean cable,such as mooring lines,towing,and ROV umbilical cables.

AN INTERPOLATION METHOD BY BIVARIATE CUBIC SPLINES WITH C^2 -JOIN ON TYPE-II TRIANGULARS AT A RECTANGULAR DOMAIN

In this paper, an interpolating method for bivariate cubic splines with C2-join on type-II triangular at a rectangular domain is given, and the approximation degree, inter- polating existence and uniqueness of the cubic splines are studied.

Numerical solution of Poisson equation with wavelet bases of Hermite cubic splines on the interval

A new wavelet-based finite element method is proposed for solving the Poisson equation. The wavelet bases of Hermite cubic splines on the interval are employed as the multi-scale interpolation basis in the finite element analysis. The lifting scheme of the wavelet-based finite element method is discussed in detail. For the orthogonal characteristics of the wavelet bases with respect to the given inner product, the corresponding multi-scale finite element equation can be decoupled across scales, totally or partially, and suited for nesting approximation. Numerical examples indicate that the proposed method has the higher effciency and precision in solving the Poisson equation.

Generalized fairing algorithm of parametric cubic splines

Kjellander has reported an algorithm for fairing uniform parametric cubic splines. Poliakoff extended Kjellander’s algorithm to non-uniform case. However, they merely changed the bad point’s position, and neglected the smoothing of tangent at bad point. In this paper, we present a fairing algorithm that both changed point’s position and its corresponding tangent vector. The new algorithm possesses the minimum property of energy. We also proved Poliakoff’s fairing algorithm is a deduction of our fairing algorithm. Several fairing examples are given in this paper.

DEFICIENT CUBIC SPLINES WITH AVERAGE SLOPE MATCHING

We obtain a deficient cubic spline function which matches the functions with certain area matching over a greater mesh intervals, and also provides a greater flexibility in replacing area matching as interpolation. We also study their convergence properties to the interpolating functions.

Non-overshooting and Non-undershooting Cubic Spline Interpolation for Empirical Mode Decomposition

To suppress the overshoots and undershoots in the envelope fitting for empirical mode decomposition (EMD), an alternative cubic spline interpolation method without overshooting and undershooting is proposed. On the basis of the derived slope constraints of knots of a non-overshooting and non-undershooting cubic interpolant, together with "not-a-knot" conditions the cubic spline interpolants are constructed by replacing the requirement for equal second order derivatives at every knot with Brodlie's derivative formula. Analysis and simulation experiments show that this approach can effectively avoid generating new extrema, shifting or exaggerating the existing ones in a signal, and thus significantly improve the decomposition performance of EMD.

A polynomial smooth epsilon-support vector regression based on cubic spline interpolation

Regression analysis is often formulated as an optimization problem with squared loss functions.Facing the challenge of the selection of the proper function class with polynomial smooth techniques applied to support vector regression models,this study takes cubic spline interpolation to generate a new polynomial smooth function |x|_ε~2 in ε-insensitive support vector regression.Theoretical analysis shows that S_ε~2-function is better than p_ε~2-function in properties,and the approximation accuracy of the proposed smoothing function is two order higher than that of classical p_ε~2-function.The experimental data shows the efficiency of the new approach.

Stability Analysis of Spatial Cubic Spline Geometric Nonlinear Beam Element Considering the Second-Order Effect

To analyze the stability problem of spatial beam structure more accurately,a spatial cubic spline geometric nonlinear beam element was proposed considering the second-order effect.The deformation field was built with cubic spline function,and its curvature degree of freedom(DOF) was eliminated by static condensation method.Then we got the geometric nonlinear stiffness matrix of the new spatial two-node Euler-Bernoulli beam element.Several examples proved calculation accuracy of the critical load by meshing a bar to one element using the method of this paper was equivalent to mesh a bar to 3 or 4 traditional nonlinear beam elements.

AN INTEGRATION METHOD WITH FITTING CUBIC SPLINE FUNCTIONS TO A NUMERICAL MODEL OF 2ND-ORDER SPACE-TIME DIFFERENTIAL REMAINDER——FOR AN IDEAL GLOBAL SIMULATION CASE WITH PRIMITIVE ATMOSPHERIC EQUATIONS

In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.

Cubic Spline Interpolation on a Class of Triangulations

In this paper,we consider spaces of cubic C^1-spline on a class of trian- gulations.By using the inductive algorithm,the posed Lagrange interpolation sets are constructed for cubic spline space.It is shown that the class of triangulations considered in this paper are nonsingular for S(1/3) spaces.Moreover,the dimensions of those spaces exactly equal to L.L.Schumaker's low bounds of the dimensions.At the end of this paper,we present an approach to construct triangulations from any scattered planar points,which ensures that the obtained triangulations for S(1/3) space are nonsingular.

A simple piecewise cubic spline method for approximation of highly nonlinear data

Approximation methods are used in the analysis and prediction of data, especially laboratory data, in engineering projects. These methods are usually linear and are obtained by least-square-error approaches. There are many problems in which linear models cannot be applied. Because of that there are logarithmic, exponential and polynomial curve-fitting models. These nonlinear models have a limited application in engineering problems. The variation of most data is such that the nonlinearity cannot be approximated by the above approaches. These methods are also not applicable when there is a large amount of data. For these reasons, a method of piecewise cubic spline approximation has been developed. The model presented here is capable of following the local nonuniformity of data in order to obtain a good fit of a curve to the data. There is C1 continuity at the limits of the piecewise elements. The model is tested and examined with four problems here. The results show that the model can approximate highly nonlinear data efficiently.

An O(k²+kh²+h²) Accurate Two-level Implicit Cubic Spline Method for One Space Dimensional Quasi-linear Parabolic Equations

In this piece of work, using three spatial grid points, we discuss a new two-level implicit cubic spline method of O(k2 + kh2 + h4) for the solution of quasi-linear parabolic equation , 0 0 subject to appropriate initial and Dirichlet boundary conditions, where h > 0, k > 0 are grid sizes in space and time-directions, respectively. The cubic spline approximation produces at each time level a spline function which may be used to obtain the solution at any point in the range of the space variable. The proposed cubic spline method is applicable to parabolic equations having singularity. The stability analysis for diffusion- convection equation shows the unconditionally stable character of the cubic spline method. The numerical tests are performed and comparative results are provided to illustrate the usefulness of the proposed method.

Curve length estimation based on cubic spline interpolation in gray-scale images

This paper deals with a novel local arc length estimator for curves in gray-scale images.The method first estimates a cubic spline curve fit for the boundary points using the gray-level information of the nearby pixels,and then computes the sum of the spline segments’lengths.In this model,the second derivatives and y coordinates at the knots are required in the computation;the spline polynomial coefficients need not be computed explicitly.We provide the algorithm pseudo code for estimation and preprocessing,both taking linear time.Implementation shows that the proposed model gains a smaller relative error than other state-of-the-art methods.

Explicit Gaussian Quadrature Rules for C1 Cubic Splines with Non-uniform Knot Sequences

This paper provides the explicit and optimal quadrature rules for the cubic C1 spline space,which is the extension of the results in Ait-Haddou et al.(J Comput Appl Math 290:543–552,2015)for less restricted non-uniform knot values.The rules are optimal in the sense that there exist no other quadrature rules with fewer quadrature points to exactly integrate the functions in the given spline space.The explicit means that the quadrature nodes and weights are derived via an explicit recursive formula.Numerical experiments and the error estimations of the quadrature rules are also presented in the end.