Numerical Solutions of Fractional Variable Order Differential Equations via Using Shifted Legendre Polynomials

In this manuscript,an algorithm for the computation of numerical solutions to some variable order fractional differential equations(FDEs)subject to the boundary and initial conditions is developed.We use shifted Legendre polynomials for the required numerical algorithm to develop some operational matrices.Further,operational matrices are constructed using variable order differentiation and integration.We are finding the operationalmatrices of variable order differentiation and integration by omitting the discretization of data.With the help of aforesaid matrices,considered FDEs are converted to algebraic equations of Sylvester type.Finally,the algebraic equations we get are solved with the help of mathematical software like Matlab or Mathematica to compute numerical solutions.Some examples are given to check the proposed method’s accuracy and graphical representations.Exact and numerical solutions are also compared in the paper for some examples.The efficiency of the method can be enhanced further by increasing the scale level.

New Form of Legendre Polynomials Obtained by Virtue of Excited Squeezed State and IWOP Technique in Quantum Optics

Using the technique of integration within an ordered product of operators and the intermediate coordinatemomentum representation in quantum optics, as well as the excited squeezed state we derive a new form of Legendre polynomials.

Characterization of Optical Aberrations Induced by Thermal Gradients and Vibrations via Zernike and Legendre Polynomials

For every astronomical instrument, the operating conditions are undoubtedly different from those defined in a setup experiment. Besides environmental conditions, the drives, the electronic cabinets containing heaters and fans introduce disturbances that must be taken into account already in the preliminary design phase. Such disturbances can be identified as being mostly of two types: heat sources/sinks or cooling systems responsible for heat transfer via conduction, radiation, free and forced convection on one side and random and periodic vibrations on the other. For this reason, a key role already from the very beginning of the design process is played by integrated model merging the outcomes based on a Finite Element Model from thermo-structural and modal analysis into the optical model to estimate the aberrations. The current paper presents the status of such model, capable of analyzing the deformed surfaces deriving from both thermo-structural and vibrational analyses and measuring their effect in terms of optical aberrations by fitting them by Zernike and Legendre polynomial fitting respectively for circular and rectangular apertures. The independent contribution of each aberration is satisfied by the orthogonality of the polynomials and mesh uniformity.

Legendre Polynomial Kernel: Application in SVM

In machines learning problems, Support Vector Machine is a method of classification. For non-linearly separable data, kernel functions are a basic ingredient in the SVM technic. In this paper, we briefly recall some useful results on decomposition of RKHS. Based on orthogonal polynomial theory and Mercer theorem, we construct the high power Legendre polynomial kernel on the cube [-1,1]d. Following presentation of the theoretical background of SVM, we evaluate the performance of this kernel on some illustrative examples in comparison with Rbf, linear and polynomial kernels.

Legendre-Weighted Residual Methods for System of Fractional Order Differential Equations

The numerical approach for finding the solution of fractional order systems of boundary value problems (BPVs) is derived in this paper. The implementation of the weighted residuals such as Galerkin, Least Square, and Collocation methods are included for solving fractional order differential equations, which is broadened to acquire the approximate solutions of fractional order systems with differentiable polynomials, namely Legendre polynomials, as basis functions. The algorithm of the residual formulations of matrix form can be coded efficiently. The interpretation of Caputo fractional derivatives is employed here. We have demonstrated these methods numerically through a few examples of linear and nonlinear BVPs. The results in absolute errors show that the present method efficiently finds the numerical solutions of fractional order systems of differential equations.

Synchronization of two different chaotic systems using Legendre polynomials with applications in secure communications

In this study, a new controller for chaos synchronization is proposed. It consists of a state feedback controller and a robust control term using Legendre polynomials to compensate for uncertainties. The truncation error is also considered. Due to the orthogonal functions theorem, Legendre polynomials can approximate nonlinear functions with arbitrarily small approximation errors. As a result, they can replace fuzzy systems and neural networks to estimate and compensate for uncertainties in control systems. Legendre polynomials have fewer tuning parameters than fuzzy systems and neural networks. Thus, their tuning process is simpler. Similar to the parameters of fuzzy systems, Legendre coefficients are estimated online using the adaptation rule obtained from the stability analysis. It is assumed that the master and slave systems are the Lorenz and Chen chaotic systems, respectively. In secure communication systems, observer-based synchronization is required since only one state variable of the master system is sent through the channel. The use of observer-based synchronization to obtain other state variables is discussed. Simulation results reveal the effectiveness of the proposed approach. A comparison with a fuzzy sliding mode controller shows that the proposed controller provides a superior transient response. The problem of secure communications is explained and the controller performance in secure communications is examined.

Evaluate More General Integrals Involving Universal Associated Legendre Polynomials via Taylor's Theorem

Abstract A few important integrals involving the product of two universal associated Legendre polynomials Pl＇m＇, （x）,Pk＇n＇（x）and x2a（1-x2）-p-1,xb（1± x）-p-1and xc（1-x2）-p-1（1 ± x）axe evaluated using the operator form of Taylor＇s theorem and an integral over a single universal associated Legendre polynomial. These integrals are more general since the quantum numbers are unequal, i.e.l＇ ≠ k＇ and m＇≠ n＇ .Their selection rules are a/so given. We also verify the correctness of those integral formulas numerically.

Galerkin approximation with Legendre polynomials for a continuous-time nonlinear optimal control problem

We investigate the use of an approximation method for obtaining near-optimal solutions to a kind of nonlinear continuous-time(CT) system. The approach derived from the Galerkin approximation is used to solve the generalized Hamilton-Jacobi-Bellman(GHJB) equations. The Galerkin approximation with Legendre polynomials(GALP) for GHJB equations has not been applied to nonlinear CT systems. The proposed GALP method solves the GHJB equations in CT systems on some well-defined region of attraction. The integrals that need to be computed are much fewer due to the orthogonal properties of Legendre polynomials, which is a significant advantage of this approach. The stabilization and convergence properties with regard to the iterative variable have been proved.Numerical examples show that the update control laws converge to the optimal control for nonlinear CT systems.

Duality between Bessel Functions and Chebyshev Polynomials in Expansions of Functions

In expansions of arbitrary functions in Bessel functions or Spherical Bessel functions, a dual partner set of polynomials play a role. For the Bessel functions, these are the Chebyshev polynomials of first kind and for the Spherical Bessel functions the Legendre polynomials. These two sets of functions appear in many formulas of the expansion and in the completeness and (bi)-orthogonality relations. The analogy to expansions of functions in Taylor series and in moment series and to expansions in Hermite functions is elaborated. Besides other special expansion, we find the expansion of Bessel functions in Spherical Bessel functions and their inversion and of Chebyshev polynomials of first kind in Legendre polynomials and their inversion. For the operators which generate the Spherical Bessel functions from a basic Spherical Bessel function, the normally ordered (or disentangled) form is found.

SOME EXTREMAL PROPERTIES OF THE INTEGRAL OF LEGENDRE POLYNOMIALS

Some extremal properties of the integral of Legendre polynomials are given, which are of independent interest. Meanwhile they show that a conjecture of P. Erdos[1] is plausible and maybe provides some means to prove this conjecture.

A Chebyshev/Legendre polynomial interpolation approach for fingerprint orientation estimation smoothing and prediction

We introduce a novel coarse ridge orientation smoothing algorithm based on orthogonal polynomials, which can be used to estimate the orientation field (OF) for fingerprint areas of no ridge information. This method does not need any base information of singular points (SPs). The algorithm uses a consecutive application of filtering-and model-based orientation smoothing methods. A Gaussian filter has been employed for the former. The latter conditionally employs one of the orthogonal polynomials such as Legendre and Chebyshev type I or II, based on the results obtained at the filtering-based stage. To evaluate our proposed method, a variety of exclusive fingerprint classification and minutiae-based matching experiments have been conducted on the fingerprint images of FVC2000 DB2, FVC2004 DB3 and DB4 databases. Results showed that our proposed method has achieved higher SP detection, classification, and verification performance as compared to competing methods.

Two-Stream Approximation to the Radiative Transfer Equation:A New Improvement and Comparative Accuracy with Existing Methods

Mathematical modeling of the interaction between solar radiation and the Earth's atmosphere is formalized by the radiative transfer equation(RTE), whose resolution calls for two-stream approximations among other methods. This paper proposes a new two-stream approximation of the RTE with the development of the phase function and the intensity into a third-order series of Legendre polynomials. This new approach, which adds one more term in the expression of the intensity and the phase function, allows in the conditions of a plane parallel atmosphere a new mathematical formulation of γparameters. It is then compared to the Eddington, Hemispheric Constant, Quadrature, Combined Delta Function and Modified Eddington, and second-order approximation methods with reference to the Discrete Ordinate(Disort) method(δ –128 streams), considered as the most precise. This work also determines the conversion function of the proposed New Method using the fundamental definition of two-stream approximation(F-TSA) developed in a previous work. Notably,New Method has generally better precision compared to the second-order approximation and Hemispheric Constant methods. Compared to the Quadrature and Eddington methods, New Method shows very good precision for wide domains of the zenith angle μ 0, but tends to deviate from the Disort method with the zenith angle, especially for high values of optical thickness. In spite of this divergence in reflectance for high values of optical thickness, very strong correlation with the Disort method(R ≈ 1) was obtained for most cases of optical thickness in this study. An analysis of the Legendre polynomial series for simple functions shows that the high precision is due to the fact that the approximated functions ameliorate the accuracy when the order of approximation increases, although it has been proven that there is a limit order depending on the function from which the precision is lost. This observation indicates that increasing the order of approximation of the phase function of the RTE leads to a better precision in flux calculations. However, this approach may be limited to a certain order that has not been studied in this paper.

Massive Dirac particles based on gapped graphene with Rosen-Morse potential in a uniform magnetic field

We explore the gapped graphene structure in the two-dimensional plane in the presence of the Rosen-Morse potential and an external uniform magnetic field.In order to describe the corresponding structure,we consider the propagation of electrons in graphene as relativistic fermion quasi-particles,and analyze it by the wave functions of two-component spinors with pseudo-spin symmetry using the Dirac equation.Next,to solve and analyze the Dirac equation,we obtain the eigenvalues and eigenvectors using the Legendre differential equation.After that,we obtain the bounded states of energy depending on the coefficients of Rosen-Morse and magnetic potentials in terms of quantum numbers of principal n and spin-orbit k.Then,the values of the energy spectrum for the ground state and the first excited state are calculated,and the wave functions and the corresponding probabilities are plotted in terms of coordinates r.In what follows,we explore the band structure of gapped graphene by the modified dispersion relation and write it in terms of the two-dimensional wave vectors K_(x) and K_(y).Finally,the energy bands are plotted in terms of the wave vectors K_(x) and K_(y) with and without the magnetic term.

Geometric Covariance Modeling for Surface Variation of Compliant Parts Based on Hybrid Polynomial Approximation and Spectrum Analysis

Part variation characterization is essential to analyze the variation propagation in flexible assemblies. Aiming at two governing types of surface variation,warping and waviness,a comprehensive approach of geometric covariance modeling based on hybrid polynomial approximation and spectrum analysis is proposed,which can formulate the level and the correlation of surface variations accurately. Firstly,the form error data of compliant part is acquired by CMM. Thereafter,a Fourier-Legendre polynomial decomposition is conducted and the error data are approximated by a Legendre polynomial series. The weighting coefficient of each component is decided by least square method for extracting the warping from the surface variation. Consequently,a geometrical covariance expression for warping deformation is established. Secondly,a Fourier-sinusoidal decomposition is utilized to approximate the waviness from the residual error data. The spectrum is analyzed is to identify the frequency and the amplitude of error data. Thus,a geometrical covariance expression for the waviness is deduced. Thirdly,a comprehensive geometric covariance model for surface variation is developed by the combination the Legendre polynomials with the sinusoidal polynomials. Finally,a group of L-shape sheet metals is measured along a specific contour,and the covariance of the profile errors is modeled by the proposed method. Thereafter,the result is compared with the covariance from two other methods and the real data. The result shows that the proposed covariance model can match the real surface error effectively and represents a tighter approximation error compared with the referred methods.

Legendre expansion method for Helmholz equations

A spectral method based on the Legendre polynomials for solving Helmholz equations was proposed. With an explicit formula for the Legendre polynomials in terms of arbitrary order of their derivatives, the successive integration of the Legendre polynomials was represented by the Legendre polynomials. Then the method was formulized for secondorder differential equations in one dimension and two dimensions. Numerical results indicate that the suggested method is significantly accurate and in satisfactory agreement with the exact solution.

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.

Kernel polynomial representation for imaginary-time Green's functions in continuous-time quantum Monte Carlo impurity solver

Inspired by the recently proposed Legendre orthogonal polynomial representation for imaginary-time Green s functions G（τ）,we develop an alternate and superior representation for G（τ） and implement it in the hybridization expansion continuous-time quantum Monte Carlo impurity solver.This representation is based on the kernel polynomial method,which introduces some integral kernel functions to filter the numerical fluctuations caused by the explicit truncations of polynomial expansion series and can improve the computational precision significantly.As an illustration of the new representation,we re-examine the imaginary-time Green＇s functions of the single-band Hubbard model in the framework of dynamical mean-field theory.The calculated results suggest that with carefully chosen integral kernel functions,whether the system is metallic or insulating,the Gibbs oscillations found in the previous Legendre orthogonal polynomial representation have been vastly suppressed and remarkable corrections to the measured Green＇s functions have been obtained.

Data fitting and modeling of regional geomagnetic field

The selection of the truncation level （TL） and the control of boundary effect （BE） are critical in regional geomagnetic field models that are based on data fitting. We combine Taylor and Legendre polynomials to model geomagnetic data over China's Mainland for years 1960, 1970, 1990, and 2000. To tackle the TL and BE problems, we first determine the range of TL by calculating the root-mean-square error （RMSE） of the models. Next, we determine the optimum TL using the Akaike information criterion （AIC） and the normalized root- mean-square error （NRMSE）. We use the regional anomaly addition （RAA） and the uniform addition （UA） method to add supplementary point outside the national boundary, and find that the intensities of extreme points gradually decrease and stabilize. The UA method better controls BEs over China, whereas the RAA method does a better job at smaller scales. In summary, we rely on a three-step method to determine the optimum TL and propose criteria to determine the optimum number of supplementary points.

An analysis of the characteristics of crustal magnetic anomaly in China based on CHAMP satellite data

Based on the observation data of CHAMP satellite from 2006 to 2009, a 2D crustal magnetic anomaly model in China is established to study the distribution characteristics of crustal magnetic anomaly. In this paper, the 2D anomaly model is derived from the Legendre polynomial expansion of harmonic term N =6-50. The result shows that many elaborate structures reflected in magnetic anomaly map well correspond to the geologic structures in China and its adjacent area. The magnetic anomaly at low satellite height behaves complexly, which is mainly caused by the magnetic disturbance of shallow rocks.In contrast, the magnetic field isolines at high satellite height are relatively sparse and only magnetic anomalies of deep crust are reflected. This fact implies that the 2D model of crustal magnetic anomaly provides an important method of the space prolongation of geomagnetic field, and is of theoretical and practice importance in geologic structure analysis and geophysical prospecting.

Critical velocity of sandwich cylindrical shell under moving internal pressure

Critical velocity of an infinite long sandwich shell under moving internal pressure is studied using the sandwich shell theory and elastodynamics theory. Propagation of axisymmetric free harmonic waves in the sandwich shell is studied using the sandwich shell theory by considering compressibility and transverse shear deformation of the core, and transverse shear deformation of face sheets. Based on the elastodynamics theory, displacement components expanded by Legendre polynomials, and position-dependent elastic constants and densities are introduced into the equations of motion. Critical velocity is the minimum phase velocity on the desperation relation curve obtained by using the two methods. Numerical examples and the finite element （FE） simulations are presented. The results show that the two critical velocities agree well with each other, and two desperation relation curves agree well with each other when the wave number k is relatively small. However, two limit phase velocities approach to the shear wave velocities of the face sheet and the core respectively when k limits to infinite. The two methods are efficient in the investigation of wave propagation in a sandwich cylindrical shell when k is relatively small. The critical velocity predicted in the FE simulations agrees with theoretical prediction.