The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the ...The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a dis- crete nonlinear programming （NLP） problem using the Chebyshev-Gauss pseudospectral method （CGPM）. Solutions of the NLP problem can be obtained using the sequential quadratic programming （SQP） algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss （CG） points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform （FFT）. To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and con- trol variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.展开更多
This paper proposes an optimal midcourse guidance method for dual pulse air-to-air missiles,which is based on the framework of the linear Gauss pseudospectral model predictive control method.Firstly,a multistage optim...This paper proposes an optimal midcourse guidance method for dual pulse air-to-air missiles,which is based on the framework of the linear Gauss pseudospectral model predictive control method.Firstly,a multistage optimal control problem with unspecified terminal time is formulated.Secondly,the control and terminal time update formulas are derived analytically.In contrast to previous work,the derivation process fully considers the Hamiltonian function corresponding to the unspecified terminal time,which is coupled with control,state,and costate.On the assumption of small perturbation,a special algebraic equation is provided to represent the equivalent optimal condition for the terminal time.Also,using Gauss pseudospectral collocation,error propagation dynamical equations involving the first-order correction term of the terminal time are transformed into a set of algebraic equations.Furthermore,analytical modification formulas can be derived by associating those equations and optimal conditions to eliminate terminal error and approach nonlinear optimal control.Even with their mathematical complexity,these formulas produce more accurate control and terminal time corrections and remove reliance on task-related parameters.Finally,several numerical simulations,comparisons with typical methods,and Monte Carlo simulations have been done to verify its optimality,high convergence rate,great stability and robustness.展开更多
Pseudospectral method is an efficient and high accuracy numerical method for simulating seismic wave propaga- tion in heterogeneous earth medium. Since its derivative operator is global, this method is commonly consid...Pseudospectral method is an efficient and high accuracy numerical method for simulating seismic wave propaga- tion in heterogeneous earth medium. Since its derivative operator is global, this method is commonly considered not suitable for parallel computation. In this paper, we introduce the parallel overlap domain decomposition scheme and give a parallel pseudospectral method implemented on distributed memory PC cluster system for modeling seismic wave propagation in heterogeneous medium. In this parallel method, the medium is decomposed into several subdomains and the wave equations are solved in each subdomain simultaneously. The solutions in each subdomain are connected through the transferring at the overlapped region. Using 2D models, we compared the parallel and traditional pseudospectral method, analyzed the accuracy of the parallel method. The results show that the parallel method can efficiently reduce computation time for the same accuracy as the traditional method. This method could be applied to large scale modeling of seismic wave propagation in 3D heterogeneous medium.展开更多
In view of generating optimal trajectories of Bolza problems, standard Chebyshev pseudospectral （PS） method makes the points＇ accumulation near the extremities and rarefaction of nodes close to the center of interv...In view of generating optimal trajectories of Bolza problems, standard Chebyshev pseudospectral （PS） method makes the points＇ accumulation near the extremities and rarefaction of nodes close to the center of interval, which causes an ill-condition of differentiation matrix and an oscillation of the optimal solution. For improvement upon the difficulties, a mapped Chebyshev pseudospectral method is proposed. A conformal map is applied to Chebyshev points to move the points closer to equidistant nodes. Condition number and spectral radius of differentiation matrices from both methods are presented to show the improvement. Furthermore, the modification keeps the Chebyshev pseudospectral method＇s advantage, the spectral convergence rate. Based on three numerical examples, a comparison of the execution time, convergence and accuracy is presented among the standard Chebyshev pseudospectral method, other collocation methods and the proposed one. In one example, the error of results from mapped Chebyshev pseudospectral method is reduced to 5% of that from standard Chebyshev pseudospectral method.展开更多
Aiming at increasing the calculation efficiency of the pseudospectral methods, a multiple- interval Radau pseudospectral method （RPM） is presented to generate a reusable launch vehicle （RLV） ＇s optimal re-entry t...Aiming at increasing the calculation efficiency of the pseudospectral methods, a multiple- interval Radau pseudospectral method （RPM） is presented to generate a reusable launch vehicle （RLV） ＇s optimal re-entry trajectory. After dividing the optimal control problem into many intervals, the state and control variables are approximated using many fixed- and low-degree Lagrange polyno- mials in each interval. Convergence of the numerical discretization is then achieved by increasing the number of intervals. With the application of the proposed method, the normal nonlinear program- ming （NLP） problem transcribed from the optimal control problem can avoid being dense because of the low-degree approximation polynomials in each interval. Thus, the NLP solver can easily compute a solution. Finally, simulation results show that the optimized re-entry trajectories satisfy the path constraints and the boundary constraints successfully. Compared with the single interval RPM, the multiple-interval RPM is significantly faster and has higher calculation efficiency. The results indicate that the multiple-interval RPM can be applied for real-time trajectory generation due to its high effi- ciency and high precision.展开更多
In this study,the problem of time-optimal reconnaissance trajectory design for the aeroassisted vehicle is considered.Different from most works reported previously,we explore the feasibility of applying a high-order a...In this study,the problem of time-optimal reconnaissance trajectory design for the aeroassisted vehicle is considered.Different from most works reported previously,we explore the feasibility of applying a high-order aeroassisted vehicle dynamic model to plan the optimal flight trajectory such that the gap between the simulated model and the real system can be narrowed.A highly-constrained optimal control model containing six-degree-of-freedom vehicle dynamics is established.To solve the formulated high-order trajectory planning model,a pipelined optimization strategy is illustrated.This approach is based on the variable order Radau pseudospectral method,indicating that the mesh grid used for discretizing the continuous system experiences several adaption iterations.Utilization of such a strategy can potentially smooth the flight trajectory and improve the algorithm convergence ability.Numerical simulations are reported to demonstrate some key features of the optimized flight trajectory.A number of comparative studies are also provided to verify the effectiveness of the applied method as well as the high-order trajectory planning model.展开更多
To improve the survivability during an emergency situation, an algorithm for aircraft forced landing trajectory planning is proposed. The method integrates damaged aircraft modelling and trajectory planning into an op...To improve the survivability during an emergency situation, an algorithm for aircraft forced landing trajectory planning is proposed. The method integrates damaged aircraft modelling and trajectory planning into an optimal control framework, in order to deal with the complex aircraft flight dynamics, a solving strategy based on Gauss pseudospetral method （GPM） is presented. A 3-DOF nonlinear mass-point model taking into account the wind is developed to approximate the aircraft flight dynamics after loss of thrust. The solution minimizes the forced landing duration, with respect to the constraints that translate the changed dynamics, flight envelope limitation and operational safety requirements. The GPM is used to convert the trajectory planning problem to a nonlinear programming problem （NLP）, which is solved by sequential quadratic programming algorithm. Simulation results show that the proposed algorithm can generate the minimum-time forced landing trajectory in event of engine-out with high efficiency and precision.展开更多
A pseudospectral method with symplectic algorithm for the solution of time-dependent Schrodinger equations （TDSE） is introduced. The spatial part of the wavefunction is discretized into sparse grid by pseudospectral...A pseudospectral method with symplectic algorithm for the solution of time-dependent Schrodinger equations （TDSE） is introduced. The spatial part of the wavefunction is discretized into sparse grid by pseudospectral method and the time evolution is given in symplectic scheme. This method allows us to obtain a highly accurate and stable solution of TDSE. The effectiveness and efficiency of this method is demonstrated by the high-order harmonic spectra of one-dimensional atom in strong laser field as compared with previously published work. The influence of the additional static electric field is also investigated.展开更多
To be close to the practical flight process and increase the precision of optimal trajectory, a six-degree-offreedom（6-DOF） trajectory is optimized for the reusable launch vehicle（RLV） using the Gauss pseudospectr...To be close to the practical flight process and increase the precision of optimal trajectory, a six-degree-offreedom（6-DOF） trajectory is optimized for the reusable launch vehicle（RLV） using the Gauss pseudospectral method（GPM）. Different from the traditional trajectory optimization problem which generally considers the RLV as a point mass, the coupling between translational dynamics and rotational dynamics is taken into account. An optimization problem is formulated to minimize a performance index subject to 6-DOF equations of motion, including translational and rotational dynamics. A two-step optimal strategy is then introduced to reduce the large calculations caused by multiple variables and convergence confinement in 6-DOF trajectory optimization. The simulation results demonstrate that the 6-DOF trajectory optimal strategy for RLV is feasible.展开更多
When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energ...When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.展开更多
In this paper, we propose to replace the Chebyshev series used in pseudospectral methods with the equivalent Chebyshev economized power series that can be evaluated more rapidly. We keep the rest of the implementation...In this paper, we propose to replace the Chebyshev series used in pseudospectral methods with the equivalent Chebyshev economized power series that can be evaluated more rapidly. We keep the rest of the implementation the same as the spectral method so that there is no new mathematical principle involved. We show by numerical examples that the new approach works well and there is indeed no significant loss of solution accuracy. The advantages of using power series also include simplicity in its formulation and implementation such that it could be used for complex systems. We investigate the important issue of collocation point selection. Our numerical results indicate that there is a clear accuracy advantage of using collocation points corresponding to roots of the Chebyshev polynomial.展开更多
Autonomous-rail rapid transit(ART)is a new medium-capacity rapid transportation system with punctuality,comfort and convenience,but low-cost construction.Combined velocity planning is a critical approach to meet the r...Autonomous-rail rapid transit(ART)is a new medium-capacity rapid transportation system with punctuality,comfort and convenience,but low-cost construction.Combined velocity planning is a critical approach to meet the requirements of energy-saving and punctuality.An ART velocity pre-planning and re-planning strategy based on the combination of punctuality dynamic programming(PDP)and pseudospectral(PS)method is proposed in this paper.Firstly,the longitudinal dynamics model of ART is established by a multi-particle model.Secondly,the PDP algorithm with global optimal characteristics is adopted as the pre-planning strategy.A model for determining the number of collocation points of the real-time PS method is proposed to improve the energy-saving effect while ensuring computation efficiency.Then the enhanced PS method is utilized to design the velocity re-planning strategy.Finally,simulations are conducted in the typical scenario with sloping roads,traffic lights,and intrusion of the pedestrian.The simulation results indicate that the ART with the proposed velocity trajectory optimization strategy can meet the punctuality requirement,and obtain better economy efficiency compared with the punctuality green light optimal speed advisory(PGLOSA).展开更多
This paper ix devoted to establishment of the Chebyshev pseudospectral domain de-composition scheme for solving two-dimensional elliptic equation. By the generalized equivalent variatiunal form, we can get the stabili...This paper ix devoted to establishment of the Chebyshev pseudospectral domain de-composition scheme for solving two-dimensional elliptic equation. By the generalized equivalent variatiunal form, we can get the stability and convergence of this new scheme.展开更多
A generalized Gauss-type quadrature formula is introduced, which assists in selection of collocation points in pseudospectral method for differential equations with two-point derivative boundary conditions. Some resul...A generalized Gauss-type quadrature formula is introduced, which assists in selection of collocation points in pseudospectral method for differential equations with two-point derivative boundary conditions. Some results on the related Jacobi interpolation are established. A pseudospectral scheme is proposed for the Kuramoto-Sivashisky equation. A skew symmetric decomposition is used for dealing with the nonlinear convection term. The stability and convergence of the proposed scheme are proved. The error estimates are obtained. Numerical results show the efficiency of this approach.展开更多
This paper is devoted to the Chebyshev pseudospectral domain decomposition method of one-dimensional elliptic problems,it is easily applied to complex geometry.The approximate accuracy can be increased by increasing t...This paper is devoted to the Chebyshev pseudospectral domain decomposition method of one-dimensional elliptic problems,it is easily applied to complex geometry.The approximate accuracy can be increased by increasing the order of approximation in fixed number of subdomains,rather than by resorting to a further partitioning.The stability and the convergence of this method are proved.展开更多
The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are est...The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are established, which play important roles in the related numerical methods for unbounded domains. As an example, the modified Laguerre spectral and pseudospectral methods are proposed for two-dimensional Logistic equation. The stability and convergence of the suggested schemes are proved. Numerical results demonstrate the high accuracy of these approaches.展开更多
The Chebyshev pseudospectral approximation of the homogeneous initial boundary value problem for a class of multi-dimensional generalized symmetric regularized long wave (SRLW) equations is considered. The fully discr...The Chebyshev pseudospectral approximation of the homogeneous initial boundary value problem for a class of multi-dimensional generalized symmetric regularized long wave (SRLW) equations is considered. The fully discrete Chebyshev pseudospectral scheme is constructed. The convergence of the approximation solution and the optimum error of approximation solution are obtained.展开更多
Results on the composite generalized Laguerre-Legendre interpolation in unbounded domains are established. As an application,a composite Laguerre-Legendre pseudospectral scheme is presented for nonlinear Fokker-Planck...Results on the composite generalized Laguerre-Legendre interpolation in unbounded domains are established. As an application,a composite Laguerre-Legendre pseudospectral scheme is presented for nonlinear Fokker-Planck equations on the whole line. The convergence and the stability of the proposed scheme are proved. Numerical results show the efficiency of the scheme and conform well to theoretical analysis.展开更多
We derive a new method for a coupled nonlinear Schr/Sdinger system by using the square of first-order Fourier spectral differentiation matrix D1 instead of traditional second-order Fourier spectral differentiation mat...We derive a new method for a coupled nonlinear Schr/Sdinger system by using the square of first-order Fourier spectral differentiation matrix D1 instead of traditional second-order Fourier spectral differentiation matrix D2 to approximate the second derivative. We prove the proposed method preserves the charge and energy conservation laws exactly. In numerical tests, we display the accuracy of numerical solution and the role of the nonlinear coupling parameter in cases of soliton collisions. Numerical experiments also exhibit the excellent performance of the method in preserving the charge and energy conservation laws. These numerical results verify that the proposed method is both a charge-preserving and an energy-preserving algorithm.展开更多
基金supported by the National Natural Science Foundation of China(No.11472058)
文摘The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a dis- crete nonlinear programming （NLP） problem using the Chebyshev-Gauss pseudospectral method （CGPM）. Solutions of the NLP problem can be obtained using the sequential quadratic programming （SQP） algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss （CG） points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform （FFT）. To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and con- trol variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.
基金Supported by Natural Science Basic Research Plan in Shaanxi Province of China(2014JQ8366)Fundamental Research Foundation of Northwestern Polytechnical University(JC20120210,JC20110238)Aeronautical Science Foundation of China(20120853007)
基金supported by the National Natural Science Foundation of China(No.62003019)the Young Talents Support Program of Beihang University,China(No.YWF-21-BJ-J-1180).
文摘This paper proposes an optimal midcourse guidance method for dual pulse air-to-air missiles,which is based on the framework of the linear Gauss pseudospectral model predictive control method.Firstly,a multistage optimal control problem with unspecified terminal time is formulated.Secondly,the control and terminal time update formulas are derived analytically.In contrast to previous work,the derivation process fully considers the Hamiltonian function corresponding to the unspecified terminal time,which is coupled with control,state,and costate.On the assumption of small perturbation,a special algebraic equation is provided to represent the equivalent optimal condition for the terminal time.Also,using Gauss pseudospectral collocation,error propagation dynamical equations involving the first-order correction term of the terminal time are transformed into a set of algebraic equations.Furthermore,analytical modification formulas can be derived by associating those equations and optimal conditions to eliminate terminal error and approach nonlinear optimal control.Even with their mathematical complexity,these formulas produce more accurate control and terminal time corrections and remove reliance on task-related parameters.Finally,several numerical simulations,comparisons with typical methods,and Monte Carlo simulations have been done to verify its optimality,high convergence rate,great stability and robustness.
基金National Natural Science Foundation of China (40474012 and 40521002)
文摘Pseudospectral method is an efficient and high accuracy numerical method for simulating seismic wave propaga- tion in heterogeneous earth medium. Since its derivative operator is global, this method is commonly considered not suitable for parallel computation. In this paper, we introduce the parallel overlap domain decomposition scheme and give a parallel pseudospectral method implemented on distributed memory PC cluster system for modeling seismic wave propagation in heterogeneous medium. In this parallel method, the medium is decomposed into several subdomains and the wave equations are solved in each subdomain simultaneously. The solutions in each subdomain are connected through the transferring at the overlapped region. Using 2D models, we compared the parallel and traditional pseudospectral method, analyzed the accuracy of the parallel method. The results show that the parallel method can efficiently reduce computation time for the same accuracy as the traditional method. This method could be applied to large scale modeling of seismic wave propagation in 3D heterogeneous medium.
基金supported by the National Natural Science Foundation of China (No.61203022)the Aeronautical Science Foundation of China (2012CZ51029)
文摘In view of generating optimal trajectories of Bolza problems, standard Chebyshev pseudospectral （PS） method makes the points＇ accumulation near the extremities and rarefaction of nodes close to the center of interval, which causes an ill-condition of differentiation matrix and an oscillation of the optimal solution. For improvement upon the difficulties, a mapped Chebyshev pseudospectral method is proposed. A conformal map is applied to Chebyshev points to move the points closer to equidistant nodes. Condition number and spectral radius of differentiation matrices from both methods are presented to show the improvement. Furthermore, the modification keeps the Chebyshev pseudospectral method＇s advantage, the spectral convergence rate. Based on three numerical examples, a comparison of the execution time, convergence and accuracy is presented among the standard Chebyshev pseudospectral method, other collocation methods and the proposed one. In one example, the error of results from mapped Chebyshev pseudospectral method is reduced to 5% of that from standard Chebyshev pseudospectral method.
文摘Aiming at increasing the calculation efficiency of the pseudospectral methods, a multiple- interval Radau pseudospectral method （RPM） is presented to generate a reusable launch vehicle （RLV） ＇s optimal re-entry trajectory. After dividing the optimal control problem into many intervals, the state and control variables are approximated using many fixed- and low-degree Lagrange polyno- mials in each interval. Convergence of the numerical discretization is then achieved by increasing the number of intervals. With the application of the proposed method, the normal nonlinear program- ming （NLP） problem transcribed from the optimal control problem can avoid being dense because of the low-degree approximation polynomials in each interval. Thus, the NLP solver can easily compute a solution. Finally, simulation results show that the optimized re-entry trajectories satisfy the path constraints and the boundary constraints successfully. Compared with the single interval RPM, the multiple-interval RPM is significantly faster and has higher calculation efficiency. The results indicate that the multiple-interval RPM can be applied for real-time trajectory generation due to its high effi- ciency and high precision.
文摘In this study,the problem of time-optimal reconnaissance trajectory design for the aeroassisted vehicle is considered.Different from most works reported previously,we explore the feasibility of applying a high-order aeroassisted vehicle dynamic model to plan the optimal flight trajectory such that the gap between the simulated model and the real system can be narrowed.A highly-constrained optimal control model containing six-degree-of-freedom vehicle dynamics is established.To solve the formulated high-order trajectory planning model,a pipelined optimization strategy is illustrated.This approach is based on the variable order Radau pseudospectral method,indicating that the mesh grid used for discretizing the continuous system experiences several adaption iterations.Utilization of such a strategy can potentially smooth the flight trajectory and improve the algorithm convergence ability.Numerical simulations are reported to demonstrate some key features of the optimized flight trajectory.A number of comparative studies are also provided to verify the effectiveness of the applied method as well as the high-order trajectory planning model.
基金supported by the National Key Basic Research Program of China(973 Program)(No.2011CB707002)
文摘To improve the survivability during an emergency situation, an algorithm for aircraft forced landing trajectory planning is proposed. The method integrates damaged aircraft modelling and trajectory planning into an optimal control framework, in order to deal with the complex aircraft flight dynamics, a solving strategy based on Gauss pseudospetral method （GPM） is presented. A 3-DOF nonlinear mass-point model taking into account the wind is developed to approximate the aircraft flight dynamics after loss of thrust. The solution minimizes the forced landing duration, with respect to the constraints that translate the changed dynamics, flight envelope limitation and operational safety requirements. The GPM is used to convert the trajectory planning problem to a nonlinear programming problem （NLP）, which is solved by sequential quadratic programming algorithm. Simulation results show that the proposed algorithm can generate the minimum-time forced landing trajectory in event of engine-out with high efficiency and precision.
基金Project supported by the National Natural Science Foundation of China （Grant Nos 10374119 and 10674154）, and The 0ne- Hundred-Talents Project of Chinese Academy of Science.Acknowledgments We gratefully acknowledge Professor Ding P Z and Professor Liu X S for their hospitality and help in symplectic algorithm.
文摘A pseudospectral method with symplectic algorithm for the solution of time-dependent Schrodinger equations （TDSE） is introduced. The spatial part of the wavefunction is discretized into sparse grid by pseudospectral method and the time evolution is given in symplectic scheme. This method allows us to obtain a highly accurate and stable solution of TDSE. The effectiveness and efficiency of this method is demonstrated by the high-order harmonic spectra of one-dimensional atom in strong laser field as compared with previously published work. The influence of the additional static electric field is also investigated.
基金supported by the National Basic Research Program of China(973 Program)(2012CB720003)the National Natural Science Foundation of China(10772011)
文摘To be close to the practical flight process and increase the precision of optimal trajectory, a six-degree-offreedom（6-DOF） trajectory is optimized for the reusable launch vehicle（RLV） using the Gauss pseudospectral method（GPM）. Different from the traditional trajectory optimization problem which generally considers the RLV as a point mass, the coupling between translational dynamics and rotational dynamics is taken into account. An optimization problem is formulated to minimize a performance index subject to 6-DOF equations of motion, including translational and rotational dynamics. A two-step optimal strategy is then introduced to reduce the large calculations caused by multiple variables and convergence confinement in 6-DOF trajectory optimization. The simulation results demonstrate that the 6-DOF trajectory optimal strategy for RLV is feasible.
文摘When there exists anisotropy in underground media elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.
文摘In this paper, we propose to replace the Chebyshev series used in pseudospectral methods with the equivalent Chebyshev economized power series that can be evaluated more rapidly. We keep the rest of the implementation the same as the spectral method so that there is no new mathematical principle involved. We show by numerical examples that the new approach works well and there is indeed no significant loss of solution accuracy. The advantages of using power series also include simplicity in its formulation and implementation such that it could be used for complex systems. We investigate the important issue of collocation point selection. Our numerical results indicate that there is a clear accuracy advantage of using collocation points corresponding to roots of the Chebyshev polynomial.
基金Supported by National Natural Science Foundation of China(Grant Nos.52072073 and 52025121)National Key R&D Program of China(Grant No.2018YFB1201602).
文摘Autonomous-rail rapid transit(ART)is a new medium-capacity rapid transportation system with punctuality,comfort and convenience,but low-cost construction.Combined velocity planning is a critical approach to meet the requirements of energy-saving and punctuality.An ART velocity pre-planning and re-planning strategy based on the combination of punctuality dynamic programming(PDP)and pseudospectral(PS)method is proposed in this paper.Firstly,the longitudinal dynamics model of ART is established by a multi-particle model.Secondly,the PDP algorithm with global optimal characteristics is adopted as the pre-planning strategy.A model for determining the number of collocation points of the real-time PS method is proposed to improve the energy-saving effect while ensuring computation efficiency.Then the enhanced PS method is utilized to design the velocity re-planning strategy.Finally,simulations are conducted in the typical scenario with sloping roads,traffic lights,and intrusion of the pedestrian.The simulation results indicate that the ART with the proposed velocity trajectory optimization strategy can meet the punctuality requirement,and obtain better economy efficiency compared with the punctuality green light optimal speed advisory(PGLOSA).
文摘This paper ix devoted to establishment of the Chebyshev pseudospectral domain de-composition scheme for solving two-dimensional elliptic equation. By the generalized equivalent variatiunal form, we can get the stability and convergence of this new scheme.
文摘A generalized Gauss-type quadrature formula is introduced, which assists in selection of collocation points in pseudospectral method for differential equations with two-point derivative boundary conditions. Some results on the related Jacobi interpolation are established. A pseudospectral scheme is proposed for the Kuramoto-Sivashisky equation. A skew symmetric decomposition is used for dealing with the nonlinear convection term. The stability and convergence of the proposed scheme are proved. The error estimates are obtained. Numerical results show the efficiency of this approach.
文摘This paper is devoted to the Chebyshev pseudospectral domain decomposition method of one-dimensional elliptic problems,it is easily applied to complex geometry.The approximate accuracy can be increased by increasing the order of approximation in fixed number of subdomains,rather than by resorting to a further partitioning.The stability and the convergence of this method are proved.
基金the Science Foundation of the Science and Technology Commission of Shanghai Municipality(No.075105118)the Shanghai Leading Academic Discipline Project(No.T0401)the Fund for E-institute of Shanghai Universities(No.E03004)
文摘The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are established, which play important roles in the related numerical methods for unbounded domains. As an example, the modified Laguerre spectral and pseudospectral methods are proposed for two-dimensional Logistic equation. The stability and convergence of the suggested schemes are proved. Numerical results demonstrate the high accuracy of these approaches.
文摘The Chebyshev pseudospectral approximation of the homogeneous initial boundary value problem for a class of multi-dimensional generalized symmetric regularized long wave (SRLW) equations is considered. The fully discrete Chebyshev pseudospectral scheme is constructed. The convergence of the approximation solution and the optimum error of approximation solution are obtained.
文摘Results on the composite generalized Laguerre-Legendre interpolation in unbounded domains are established. As an application,a composite Laguerre-Legendre pseudospectral scheme is presented for nonlinear Fokker-Planck equations on the whole line. The convergence and the stability of the proposed scheme are proved. Numerical results show the efficiency of the scheme and conform well to theoretical analysis.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11201169 and 11271195)the National Basic Research Program of China (Grant No. 2010AA012304)+1 种基金the Natural Science Foundation of Jiangsu Education Bureau,China (Grant Nos. 10KJB110001 and 12KJB110002)the Qing Lan Project of Jiangsu Province of China
文摘We derive a new method for a coupled nonlinear Schr/Sdinger system by using the square of first-order Fourier spectral differentiation matrix D1 instead of traditional second-order Fourier spectral differentiation matrix D2 to approximate the second derivative. We prove the proposed method preserves the charge and energy conservation laws exactly. In numerical tests, we display the accuracy of numerical solution and the role of the nonlinear coupling parameter in cases of soliton collisions. Numerical experiments also exhibit the excellent performance of the method in preserving the charge and energy conservation laws. These numerical results verify that the proposed method is both a charge-preserving and an energy-preserving algorithm.
