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.

Associated Hermite Polynomials Related to Parabolic Cylinder Functions

In analogy to the role of Lommel polynomials ?in relation to Bessel functions Jv(z) the theory of Associated Hermite polynomials in the scaled form ?with parmeter v to Parabolic Cylinder functions Dv(z) is developed. The group-theoretical background with the 3-parameter group of motions M(2) in the plane for Bessel functions and of the Heisenberg-Weyl group W(2) for Parabolic Cylinder functions is discussed and compared with formulae, in particular, for the lowering and raising operators and the eigenvalue equations. Recurrence relations for the Associated Hermite polynomials and for their derivative and the differential equation for them are derived in detail. Explicit expressions for the Associated Hermite polynomials with involved Jacobi polynomials at argument zero are given and by means of them the Parabolic Cylinder functions are represented by two such basic functions.

Chebyshev Polynomials with Applications to Two-Dimensional Operators

A new application of Chebyshev polynomials of second kind Un(x) to functions of two-dimensional operators is derived and discussed. It is related to the Hamilton-Cayley identity for operators or matrices which allows to reduce powers and smooth functions of them to superpositions of the first N-1 powers of the considered operator in N-dimensional case. The method leads in two-dimensional case first to the recurrence relations for Chebyshev polynomials and due to initial conditions to the application of Chebyshev polynomials of second kind Un(x). Furthermore, a new general class of Generating functions for Chebyshev polynomials of first and second kind Un(x) comprising the known Generating function as special cases is constructed by means of a derived identity for operator functions f(A) of a general two-dimensional operator A. The basic results are Formulas (9.5) and (9.6) which are then specialized for different examples of functions f(x). The generalization of the theory for three-dimensional operators is started to attack and a partial problem connected with the eigenvalue problem and the Hamilton-Cayley identity is solved in an Appendix. A physical application of Chebyshev polynomials to a problem of relativistic kinematics of a uniformly accelerated system is solved. All operator calculations are made in coordinate-invariant form.

Operator Methods and SU(1,1) Symmetry in the Theory of Jacobi and of Ultraspherical Polynomials

Starting from general Jacobi polynomials we derive for the Ul-traspherical polynomials as their special case a set of related polynomials which can be extended to an orthogonal set of functions with interesting properties. It leads to an alternative definition of the Ultraspherical polynomials by a fixed integral operator in application to powers of the variable u in an analogous way as it is possible for Hermite polynomials. From this follows a generating function which is apparently known only for the Legendre and Chebyshev polynomials as their special case. Furthermore, we show that the Ultraspherical polynomials form a realization of the SU(1,1) Lie algebra with lowering and raising operators which we explicitly determine. By reordering of multiplication and differentiation operators we derive new operator identities for the whole set of Jacobi polynomials which may be applied to arbitrary functions and provide then function identities. In this way we derive a new “convolution identity” for Jacobi polynomials and compare it with a known convolution identity of different structure for Gegenbauer polynomials. In short form we establish the connection of Jacobi polynomials and their related orthonormalized functions to the eigensolution of the Schrödinger equation to Pöschl-Teller potentials.

Improving Video Watermarking through Galois Field GF(2^(4)) Multiplication Tables with Diverse Irreducible Polynomials and Adaptive Techniques

Video watermarking plays a crucial role in protecting intellectual property rights and ensuring content authenticity.This study delves into the integration of Galois Field(GF)multiplication tables,especially GF(2^(4)),and their interaction with distinct irreducible polynomials.The primary aim is to enhance watermarking techniques for achieving imperceptibility,robustness,and efficient execution time.The research employs scene selection and adaptive thresholding techniques to streamline the watermarking process.Scene selection is used strategically to embed watermarks in the most vital frames of the video,while adaptive thresholding methods ensure that the watermarking process adheres to imperceptibility criteria,maintaining the video's visual quality.Concurrently,careful consideration is given to execution time,crucial in real-world scenarios,to balance efficiency and efficacy.The Peak Signal-to-Noise Ratio(PSNR)serves as a pivotal metric to gauge the watermark's imperceptibility and video quality.The study explores various irreducible polynomials,navigating the trade-offs between computational efficiency and watermark imperceptibility.In parallel,the study pays careful attention to the execution time,a paramount consideration in real-world scenarios,to strike a balance between efficiency and efficacy.This comprehensive analysis provides valuable insights into the interplay of GF multiplication tables,diverse irreducible polynomials,scene selection,adaptive thresholding,imperceptibility,and execution time.The evaluation of the proposed algorithm's robustness was conducted using PSNR and NC metrics,and it was subjected to assessment under the impact of five distinct attack scenarios.These findings contribute to the development of watermarking strategies that balance imperceptibility,robustness,and processing efficiency,enhancing the field's practicality and effectiveness.

Probabilistic analysis of tunnel face seismic stability in layered rock masses using Polynomial Chaos Kriging metamodel

Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.

Three-dimensional pseudo-dynamic reliability analysis of seismic shield tunnel faces combined with sparse polynomial chaos expansion

To address the seismic face stability challenges encountered in urban and subsea tunnel construction,an efficient probabilistic analysis framework for shield tunnel faces under seismic conditions is proposed.Based on the upper-bound theory of limit analysis,an improved three-dimensional discrete deterministic mechanism,accounting for the heterogeneous nature of soil media,is formulated to evaluate seismic face stability.The metamodel of failure probabilistic assessments for seismic tunnel faces is constructed by integrating the sparse polynomial chaos expansion method(SPCE)with the modified pseudo-dynamic approach(MPD).The improved deterministic model is validated by comparing with published literature and numerical simulations results,and the SPCE-MPD metamodel is examined with the traditional MCS method.Based on the SPCE-MPD metamodels,the seismic effects on face failure probability and reliability index are presented and the global sensitivity analysis(GSA)is involved to reflect the influence order of seismic action parameters.Finally,the proposed approach is tested to be effective by a engineering case of the Chengdu outer ring tunnel.The results show that higher uncertainty of seismic response on face stability should be noticed in areas with intense earthquakes and variation of seismic wave velocity has the most profound influence on tunnel face stability.

A coupled Legendre-Laguerre polynomial method with analytical integration for the Rayleigh waves in a quasicrystal layered half-space with an imperfect interface

The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.

Generalized polynomial chaos expansion by reanalysis using static condensation based on substructuring

This paper presents a new computational method for forward uncertainty quantification(UQ)analyses on large-scale structural systems in the presence of arbitrary and dependent random inputs.The method consists of a generalized polynomial chaos expansion(GPCE)for statistical moment and reliability analyses associated with the stochastic output and a static reanalysis method to generate the input-output data set.In the reanalysis,we employ substructuring for a structure to isolate its local regions that vary due to random inputs.This allows for avoiding repeated computations of invariant substructures while generating the input-output data set.Combining substructuring with static condensation further improves the computational efficiency of the reanalysis without losing accuracy.Consequently,the GPCE with the static reanalysis method can achieve significant computational saving,thus mitigating the curse of dimensionality to some degree for UQ under high-dimensional inputs.The numerical results obtained from a simple structure indicate that the proposed method for UQ produces accurate solutions more efficiently than the GPCE using full finite element analyses(FEAs).We also demonstrate the efficiency and scalability of the proposed method by executing UQ for a large-scale wing-box structure under ten-dimensional(all-dependent)random inputs.

Jacobi-Sobolev Orthogonal Polynomialsand Spectral Methods for Elliptic Boundary Value Problems

Generalized Jacobi polynomials with indexes α,β∈ R are introduced and some basic properties are established. As examples of applications,the second- and fourth-order elliptic boundary value problems with Dirichlet or Robin boundary conditions are considered,and the generalized Jacobi spectral schemes are proposed. For the diagonalization of discrete systems,the Jacobi-Sobolev orthogonal basis functions are constructed,which allow the exact solutions and the approximate solutions to be represented in the forms of infinite and truncated Jacobi series. Error estimates are obtained and numerical results are provided to illustrate the effectiveness and the spectral accuracy.

Factorization of Cyclotomic Polynomials with Quadratic Radicals in the Coefficients

In this article we continue the consideration of geometrical constructions of regular n-gons for odd n by rhombic bicompasses and ruler used in [1] for the construction of the regular heptagon (n=7). We discuss the possible factorization of the cyclotomic polynomial in polynomial factors which contain not higher than quadratic radicals in the coefficients whereas usually the factorization of the cyclotomic polynomials is considered in products of irreducible factors with integer coefficients. In considering the regular heptagon we find a modified variant of its construction by rhombic bicompasses and ruler. In detail, supported by figures, we investigate the case of the regular tridecagon (n=13) which in addition to n=7 is the only candidate with low n (the next to this is n=769 ) for which such a construction by rhombic bicompasses and ruler seems to be possible. Besides the coordinate origin we find here two points to fix for the possible application of two bicompasses (or even four with the addition of the complex conjugate points to be fixed). With only one bicompass one has in addition the problem of the trisection of an angle which can be solved by a neusis construction that, however, is not in the spirit of constructions by compass and ruler and is difficult to realize during the action of bicompasses. As discussed it seems that to finish the construction by bicompasses the correlated action of two rhombic bicompasses must be applied in this case which avoids the disadvantages of the neusis construction. Single rhombic bicompasses allow to draw at once two circles around two fixed points in such correlated way that the position of one of the rotating points on one circle determines the positions of all the other points on the second circle in unique way. The known case n=17 embedded in our method is discussed in detail.

Flow Dynamics in Restricted Geometries: A Mathematical Concept Based on Bloch NMR Flow Equation and Boubaker Polynomial Expansion Scheme

Computational techniques are invaluable to the continued success and development of Magnetic Resonance Imaging (MRI) and to its widespread applications. New processing methods are essential for addressing issues at each stage of MRI techniques. In this study, we present new sets of non-exponential generating functions representing the NMR transverse magnetizations and signals which are mathematically designed based on the theory and dynamics of the Bloch NMR flow equations. These signals are functions of many spinning nuclei of materials and can be used to obtain information observed in all flow systems. The Bloch NMR flow equations are solved using the Boubaker polynomial expansion scheme (BPES) and analytically connect most of the experimentally valuable NMR parameters in a simplified way for general analyses of magnetic resonance imaging with adiabatic condition.

Linear Functional Equations and Twisted Polynomials

A certain variety of non-switched polynomials provides a uni-figure representation for a wide range of linear functional equations. This is properly adapted for the calculations. We reinterpret from this point of view a number of algorithms.

Solving Some Problems and Elimination in Systems of Polynomial Equations

In a factorial ring, we can define the p.g.c.d. of two elements (defined to the nearest unit) and the notion of prime elements between them. More generally, Bezout’s identity characterizes two prime elements in a main ring. A ring that satisfies the property of the theorem is called a Bezout ring. We have given some geometry theorems that can be proved algebraically, although the methods of geometry and, in particular, of projective geometry are by far the most beautiful. Most geometric problems actually involve polynomial equations and can be translated into the language of polynomial ideals. We have given a few examples of a different nature without pretending to make a general theory.

Numerical Study for a Class of Variable Order Fractional Integral-differential Equation in Terms of Bernstein Polynomials

The aim of this paper is to seek the numerical solution of a class of variable order fractional integral-differential equation in terms of Bernstein polynomials.The fractional derivative is described in the Caputo sense.Four kinds of operational matrixes of Bernstein polynomials are introduced and are utilized to reduce the initial equation to the solution of algebraic equations after dispersing the variable.By solving the algebraic equations,the numerical solutions are acquired.The method in general is easy to implement and yields good results.Numerical examples are provided to demonstrate the validity and applicability of the method.

Approximation of Finite Population Totals Using Lagrange Polynomial

Approximation of finite population totals in the presence of auxiliary information is considered. A polynomial based on Lagrange polynomial is proposed. Like the local polynomial regression, Horvitz Thompson and ratio estimators, this approximation technique is based on annual population total in order to fit in the best approximating polynomial within a given period of time (years) in this study. This proposed technique has shown to be unbiased under a linear polynomial. The use of real data indicated that the polynomial is efficient and can approximate properly even when the data is unevenly spaced.

A “Hard to Die” Series Expansion and Lucas Polynomials of the Second Kind

We show how to use the Lucas polynomials of the second kind in the solution of a homogeneous linear differential system with constant coefficients, avoiding the Jordan canonical form for the relevant matrix.

Diophantine equations and Fermat's last theorem for multivariate(skew-)polynomials

Fermat’s Last Theorem is a famous theorem in number theory which is difficult to prove.However,it is known that the version of polynomials with one variable of Fermat’s Last Theorem over C can be proved very concisely.The aim of this paper is to study the similar problems about Fermat’s Last Theorem for multivariate(skew)-polynomials with any characteristic.

Sensitivity Analysis of Electromagnetic Scattering from Dielectric Targets with Polynomial Chaos Expansion and Method of Moments

In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is applied to accurately solve the electric field integral equation(EFIE)of electromagnetic scattering from homogeneous dielectric targets.Within the bistatic radar cross section(RCS)as the research object,the adaptive PCE algorithm is devoted to selecting the appropriate order to construct the multivariate surrogate model.The corresponding sensitivity results are given by the further derivative operation,which is compared with those of the finite difference method(FDM).Several examples are provided to demonstrate the effectiveness of the proposed algorithm for sensitivity analysis of electromagnetic scattering from homogeneous dielectric targets.

A Collocation Technique via Pell-Lucas Polynomials to Solve Fractional Differential EquationModel for HIV/AIDS with Treatment Compartment

In this study,a numerical method based on the Pell-Lucas polynomials(PLPs)is developed to solve the fractional order HIV/AIDS epidemic model with a treatment compartment.The HIV/AIDS mathematical model with a treatment compartment is divided into five classes,namely,susceptible patients(S),HIV-positive individuals(I),individuals with full-blown AIDS but not receiving ARV treatment(A),individuals being treated(T),and individuals who have changed their sexual habits sufficiently(R).According to the method,by utilizing the PLPs and the collocation points,we convert the fractional order HIV/AIDS epidemic model with a treatment compartment into a nonlinear system of the algebraic equations.Also,the error analysis is presented for the Pell-Lucas approximation method.The aim of this study is to observe the behavior of five populations after 200 days when drug treatment is applied to HIV-infectious and full-blown AIDS people.To demonstrate the usefulness of this method,the applications are made on the numerical example with the help of MATLAB.In addition,four cases of the fractional order derivative(p=1,p=0.95,p=0.9,p=0.85)are examined in the range[0,200].Owing to applications,we figured out that the outcomes have quite decent errors.Also,we understand that the errors decrease when the value of N increases.The figures in this study are created in MATLAB.The outcomes indicate that the presented method is reasonably sufficient and correct.