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.

Functional Integrals and Quantum Fluctuations on Two-Dimensional Noncommutative Space-Time

The generalized Thirring model with impurity coupling is defined on two-dimensional noncommutativespace-time,a modified propagator and free energy are derived by means of functional integrals method.Moreover,quantum fluctuations and excitation energies are calculated on two-dimensional black hole and soliton background.

An Eight Component Integrable Hamiltonian Hierarchy from a Reduced Seventh-Order Matrix Spectral Problem

We present an eight component integrable Hamiltonian hierarchy, based on a reduced seventh order matrix spectral problem, with the aim of aiding the study and classification of multicomponent integrable models and their underlying mathematical structures. The zero-curvature formulation is the tool to construct a recursion operator from the spatial matrix problem. The second and third set of integrable equations present integrable nonlinear Schrödinger and modified Korteweg-de Vries type equations, respectively. The trace identity is used to construct Hamiltonian structures, and the first three Hamiltonian functionals so generated are computed.

Chebyshev spectral variational integrator and applications

The Chebyshev spectral variational integrator(CSVI) is presented in this paper. Spectral methods have aroused great interest in approximating numerically a smooth problem for their attractive geometric convergence rates. The geometric numerical methods are praised for their excellent long-time geometric structure-preserving properties.According to the generalized Galerkin framework, we combine two methods together to construct a variational integrator, which captures the merits of both methods. Since the interpolating points of the variational integrator are chosen as the Chebyshev points,the integration of Lagrangian can be approximated by the Clenshaw-Curtis quadrature rule, and the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of configuration variables and the corresponding derivatives. The numerical float errors of the first-order spectral differentiation matrix can be alleviated by using a trigonometric identity especially when the number of Chebyshev points is large. Furthermore, the spectral variational integrator(SVI) constructed by the Gauss-Legendre quadrature rule and the multi-interval spectral method are carried out to compare with the CSVI, and the interesting kink phenomena for the Clenshaw-Curtis quadrature rule are discovered. The numerical results reveal that the CSVI has an advantage on the computing time over the whole progress and a higher accuracy than the SVI before the kink position. The effectiveness of the proposed method is demonstrated and verified perfectly through the numerical simulations for several classical mechanics examples and the orbital propagation for the planet systems and the Solar system.

Data acquisition,analysis and applications of multi-sensor integration

This paper presents some key techniques for multi-sensor integration system, which is applied to the intelligent transportation system industry and surveying and mapping industry, e.g. road surface condition detection, digital map making. The techniques are synchronization control of multi-sensor, space-time benchmark for sensor data, and multi-sensor data fusion and mining. Firstly, synchronization control of multi-sensor is achieved through a synchronization control system which is composed of a time synchronization controller and some synchronization sub-controllers. The time synchronization controller can receive GPS time information from GPS satellites, relative distance information from distance measuring instrument and send space-time information to the synchronization sub-controller. The latter can work at three types of synchronization mode, i.e. active synchronization, passive synchronization and time service synchronization. Secondly, space-time benchmark can be established based on GPS time and global reference coordinate system, and can be obtained through position and azimuth determining system and synchronization control system. Thirdly, there are many types of data fusion and mining, e.g. GPS/Gyro/DMI data fusion, data fusion between stereophotogrammetry and PADS, data fusion between laser scanner and PADS, and data fusion between CCD camera and laser scanner. Finally, all these solutions presented in paper have been applied to two areas, i.e. land-bone intelligent road detection and measurement system and 3D measurement system based on unmanned helicopter. The former has equipped some highway engineering Co. , Ltd. and has been successfully put into use. The latter is an ongoing resealch.

A new high-order spectral problem of the mKdV and its associated integrable decomposition

A new approach to formulizing a new high-order matrix spectral problem from a normal 2 × 2 matrix modified Korteweg-de Vries （mKdV） spectral problem is presented. It is found that the isospectral evolution equation hierarchy of this new higher-order matrix spectral problem turns out to be the well-known mKdV equation hierarchy. By using the binary nonlinearization method, a new integrable decomposition of the mKdV equation is obtained in the sense of Liouville. The proof of the integrability shows that r-matrix structure is very interesting,

A Spectral Integral Equation Solution of the Gross-Pitaevskii Equation

The Gross-Pitaevskii equation (GPE), that describes the wave function of a number of coherent Bose particles contained in a trap, contains the cube of the normalized wave function, times a factor proportional to the number of coherent atoms. The square of the wave function, times the above mentioned factor, is defined as the Hartree potential. A method implemented here for the numerical solution of the GPE consists in obtaining the Hartree potential iteratively, starting with the Thomas Fermi approximation to this potential. The energy eigenvalues and the corresponding wave functions for each successive potential are obtained by a spectral method described previously. After approximately 35 iterations a stability of eight significant figures for the energy eigenvalues is obtained. This method has the advantage of being physically intuitive, and could be extended to the calculation of a shell-model potential in nuclear physics, once the Pauli exclusion principle is allowed for.

Solving Large Scale Unconstrained Minimization Problems by a New ODE Numerical Integration Method

In reference [1], for large scale nonlinear equations , a new ODE solving method was given. This paper is a continuous work. Here has gradient structure i.e. , is a scalar function. The eigenvalues of the Jacobian of;or the Hessian of , are all real number. So the new method is very suitable for this structure. For quadratic function the convergence was proved and the spectral radius of iteration matrix was given and compared with traditional method. Examples show for large scale problems (dimension ) the new method is very efficient.

SPACE-TIME FINITE ELEMENT METHOD FOR SCHRDINGER EQUATION AND ITS CONSERVATION

Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.

A Finite-Dimensional Integrable System Related to the Complex 3 ×3 Spectral Problem and the Coupled Nonlinear Schrödinger Equation

The relation between the 3 × 3 complex spectral problem and the associated completely integrable system is generated. From the spectral problem, we derived the Lax pairs and the evolution equation hierarchy in which the coupled nonlinear Schr?dinger equation is included. Then, with the constraints between the potential function and the eigenvalue function, using the nonlineared Lax pairs, a finite-dimensional complex Hamiltonian system is obtained. Furthermore, the representation of the solution to the evolution equations is generated by the commutable flows of the finite-dimensional completely integrable system.

A New m × m Matrix Kaup-Newell Spectral Problem and Its Associated Integrable Decomposition

A new m × m matrix Kaup-Newell spectral problem is constructed from a normal 2 × 2 matrix Kaup-Newell spectral problem, a new integrable decomposition of the Kaup-Newell equation is presented. Through this process, we find the structure of the r-matrix is interesting.

Performance of Cross-Layer Design with Power Allocation in Space-Time Coded MIMO Systems

A cross-layer design(CLD)scheme with combination of power allocation,adaptive modulation(AM)and automatic repeat request(ARQ)is presented for space-time coded MIMO system under imperfect feedback,and the corresponding system performance is investigated in a Rayleigh fading channel.Based on imperfect feedback information,a suboptimal power allocation(PA)scheme is derived to maximize the average spectral efficiency(SE)of the system.The scheme is based on a so-called compressed SNR criterion,and has a closed-form expression for positive power allocation,thus being computationally efficient.Moreover,it can improve SE of the presented CLD.Besides,due to better approximation,it obtains the performance close to the existing optimal approach which requires numerical search.Simulation results show that the proposed CLD with PA can achieve higher SE than the conventional CLD with equal power allocation scheme,and has almost the same performance as CLD with optimal PA.However,it has lower calculation complexity.

A Pseudo-Spectral Scheme for Systems of Two-Point Boundary Value Problems with Left and Right Sided Fractional Derivatives and Related Integral Equations

We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the problem into an equivalent system of weakly singular integral equations.Then,a Legendre-based spectral collocation method is developed for solving the transformed system.Therefore,we can make good use of the advantages of the Gauss quadrature rule.We present the construction and analysis of the collocation method.These results can be indirectly applied to solve fractional optimal control problems by considering the corresponding Euler–Lagrange equations.Two numerical examples are given to confirm the convergence analysis and robustness of the scheme.

The Blue Stragglers and the Integrated Spectral Energy Distribution of Single Stellar Populations

Blue stragglers are a common observational fact for the Galactic clusters. Single Stellar Populations (SSPs) are basic to the studies of galaxy structre and evolution. SSPs are mainly based either on the observation of the integrated properties of star clusters, or on the theoretical understandings of single star evolution. Both of the two ways of making SSPs suffer from either observational uncertainties concerning field contaminations or lack of good models for close binary systems. Based on the photometry of the classical open cluster M67 and the thorough membership survey, we made a color-magnitude diagram (CMD) of high membership stars for the cluster. We will show that by including the contributions of the bright blue stragglers that is common to open clusters, the integrated properties of the clusters are quite different from tranditional SSP models. We further conclude that these blue light contributors are very important to SSP models, and may cast new lights on its applications in the studies of galaxies.

An accurate and efficient space-time Galerkin spectral method for the subdiffusion equation

In this paper, we design and analyze a space-time spectral method for the subdiffusion equation.Here, we are facing two difficulties. The first is that the solutions of this equation are usually singular near the initial time. Consequently, traditional high-order numerical methods in time are inefficient. The second obstacle is that the resulting system of the space-time spectral approach is usually large and time-consuming to solve. We aim at overcoming the first difficulty by proposing a novel approach in time, which is based on variable transformation techniques. Suitable ψ-fractional Sobolev spaces and a new variational framework are introduced to establish the well-posedness of the associated variational problem. This allows us to construct our space-time spectral method using a combination of temporal generalized Jacobi polynomials(GJPs) and spatial Legendre polynomials. For the second difficulty, we propose a fast algorithm to effectively solve the resulting linear system. The fast algorithm makes use of a matrix diagonalization in space and QZ decomposition in time. Our analysis and numerical experiments show that the proposed method is exponentially convergent with respect to the polynomial degrees in both space and time directions, even though the exact solution has very limited regularity.

Scattering Properties of Atmospheric Aerosols over Lanzhou City and Applications Using an Integrating Nephelometer

The data, measured by a three-wavelength Integrating Nephelometer over Lanzhou City during the winters of 2001/2002 and 2002/2003 respectively, have been analyzed for investigating the scattering properties of atmospheric aerosols and exploring their relationship and the status of air pollution. The aerosol particle volume distribution is inverted with the measured spectral scattering coe?cients. The results show that the daily variation of the aerosol scattering coe?cients is in a tri-peak shape. The average ratio of backscattering coe?cient to total scattering coe?cient at 550 nm is 0.158; there exists an excellent correlation between the scattering coe?cients and the concentration of PM10. The average ratio of the concentration of PM10 to the scattering coe?cients is 0.37 g m?2, which is contingent on the optical parameters of aerosol particles such as the size distribution, etc.; an algorithm is developed for inverting the volume distribution of aerosol particles by using the histogram and Monte-Carlo techniques, and the test results show that the inversion is reasonable.

An Integral Collocation Approach Based on Legendre Polynomials for Solving Riccati, Logistic and Delay Differential Equations

In this paper, we propose and analyze some schemes of the integral collocation formulation based on Legendre polynomials. We implement these formulae to solve numerically Riccati, Logistic and delay differential equations with variable coefficients. The properties of the Legendre polynomials are used to reduce the proposed problems to the solution of non-linear system of algebraic equations using Newton iteration method. We give numerical results to satisfy the accuracy and the applicability of the proposed schemes.

New Integrable Couplings of Generalized Kaup-Newell Hierarchy and Its Hamiltonian Structures

A new isospectral problem is firstly presented, then we derive integrable system of soliton hierarchy. Also we obtain new integrable couplings of the generalized Kaup-Newell soliton equations hierarchy and its Hamiltonian structures by using Tu scheme and the quadratic-form identity. The method can be generalized to other soliton hierarchy.

Sanxingdui Cultural Relics Recognition Algorithm Based on Hyperspectral Multi-Network Fusion

Sanxingdui cultural relics are the precious cultural heritage of humanity with high values of history,science,culture,art and research.However,mainstream analytical methods are contacting and detrimental,which is unfavorable to the protection of cultural relics.This paper improves the accuracy of the extraction,location,and analysis of artifacts using hyperspectral methods.To improve the accuracy of cultural relic mining,positioning,and analysis,the segmentation algorithm of Sanxingdui cultural relics based on the spatial spectrum integrated network is proposed with the support of hyperspectral techniques.Firstly,region stitching algorithm based on the relative position of hyper spectrally collected data is proposed to improve stitching efficiency.Secondly,given the prominence of traditional HRNet(High-Resolution Net)models in high-resolution data processing,the spatial attention mechanism is put forward to obtain spatial dimension information.Thirdly,in view of the prominence of 3D networks in spectral information acquisition,the pyramid 3D residual network model is proposed to obtain internal spectral dimensional information.Fourthly,four kinds of fusion methods at the level of data and decision are presented to achieve cultural relic labeling.As shown by the experiment results,the proposed network adopts an integrated method of data-level and decision-level,which achieves the optimal average accuracy of identification 0.84,realizes shallow coverage of cultural relics labeling,and effectively supports the mining and protection of cultural relics.

Numerical Solution of Mean-Square Approximation Problem of Real Nonnegative Function by the Modulus of Double Fourier Integral

A nonlinear problem of mean-square approximation of a real nonnegative continuous function with respect to two variables by the modulus of double Fourier integral dependent on two real parameters with use of the smoothing functional is studied. Finding the optimal solutions of this problem is reduced to solution of the Hammerstein type two-dimensional nonlinear integral equation. The numerical algorithms to find the branching lines and branching-off solutions of this equation are constructed and justified. Numerical examples are presented.