共查询到20条相似文献，搜索用时 31 毫秒
1.
考虑周期解的数值延拓问题并提出基于Broyden拟牛顿法来延拓周期解的一种有效算法,先后以布鲁塞尔振子、平面圆型限制性三体问题(Planar Circular Restricted Three-Body Problem, PCRTBP)的周期解为例进行了验证.这里的Broyden方法包含线性搜索、正交三角分解求线性方程组的步骤.对一般的周期解,周期性条件方程组中含有周期作为待延拓参数,可用周期来决定积分时长,将解代入周期性条件得到积分型的非线性方程组,利用Broyden方法迭代延拓直至初值收敛.根据两次垂直通过一个超平面的轨道是对称周期轨道的性质,可采用插值的方法求得再次抵达超平面的解分量,得到周期性条件方程组,再用Broyden方法求解.结合哈密顿系统的对称性和PCRTBP周期轨道的一些分类,对2/1、3/1的内共振周期解族进行了数值研究.最后,对算法和计算结果做了总结和讨论. 相似文献
2.
In this paper, a new family of explicit and implicit multistep methods is presented both for the error-controlled and uncontrolled modes. The main concept is to replace the Newton interpolation with the Hermite interpolation, where the Hermite polynomial is fitted to the function values and its derivatives. This idea is very useful in the numerical solution of problems (e.g., orbit propagation problem) where higher-order derivatives can easily be computed. In addition to the theoretical concept, the stability regions of the proposed methods are determined. The new methods are more stable than the well-known multistep numerical integrators (i.e., Adams–Bashforth and Adams–Bashforth–Moulton) in the explicit, implicit, and predictor–corrector forms. Using the second-order derivatives gives smaller error constants in the proposed method. The new integrators are numerically tested for a few examples, and the solutions are compared with those of the well-known multistep methods. Moreover, the CPU time and absolute integration error are compared in the satellite orbit propagation problem using various integration methods. The CHAMP mission, i.e., a German small-satellite mission for geoscientific and atmospheric research and applications, is considered as a case study for comparing the achievable accuracy of the proposed method with the existing method for solving the two-body problem. 相似文献
3.
徐继鸿 《中国科学院上海天文台年刊》2002,(1)
概述了制约线性多步积分公式轨道积分状态的多种因素 ;提出了综合评价线性多步积分公式积分性能的两项新指标。建议在对数值计算有较高精度要求的科研项目中 ,应将构造并选择适合研究项目的线性多步积分公式以及高效的积分方式列为课题前期工作的重要部分 相似文献
4.
This paper describes a second-order upwind scheme for multidimensional magnetohydrodynamics, which uses a linear approximation for all Riemann problems except those involving strong rarefactions. This enables it to cope with initial data for which previously published schemes might fail. The condition ▽⊙ B = 0 is not enforced in multidimensions, but the numerical problems associated with this are dealt with by adding source terms to the equations, as suggested by Powell. We also show that there are advantages to adding second-order artificial dissipation at shocks. 相似文献
5.
Chris Koen 《Monthly notices of the Royal Astronomical Society》2005,361(3):887-896
The theory of low-order linear stochastic differential equations is reviewed. Solutions to these equations give the continuous time analogues of discrete time autoregressive time-series. Explicit forms for the power spectra and covariance functions of first- and second-order forms are given. A conceptually simple method is described for fitting continuous time autoregressive models to data. Formulae giving the standard errors of the parameter estimates are derived. Simulated data are used to verify the performance of the methods. Irregularly spaced observations of the two hydrogen-deficient stars FQ Aqr and NO Ser are analysed. In the case of FQ Aqr the best-fitting model is of second order, and describes a quasi-periodicity of about 20 d with an e-folding time of 3.7 d. The NO Ser data are best fitted by a first-order model with an e-folding time of 7.2 d. 相似文献
6.
A new second-order solution to the two-point boundary value problem for relative motion about orbital rendezvous in one orbit
period is proposed. First, nonlinear differential equations to describe the relative motion between a chaser and a target
are presented considering the second-order terms in the gravity. Then, by regarding the second-order terms as external accelerations,
we establish second-order state transition equations. Moreover, the J2 perturbations effects can also be considered in the
state transition equations. Last, the initial relative velocity to fulfill a rendezvous is determined by solving the state
transition equations. Numerical simulations show that the new second-order state transition equations are accurate. The second-order
solution to the two-point boundary value problem on eccentric orbits is valid even if the relative range is farther than 500 km. 相似文献
7.
The Integral Variation (IV) method is a technique to generate an approximate solution to initial value problems involving systems of first-order ordinary differential equations. The technique makes use of generalized Fourier expansions in terms of shifted orthogonal polynomials. The IV method is briefly described and then applied to the problem of near Earth satellite orbit prediction. In particular, we will solve the Lagrange planetary equations including the first three zonal harmonics and drag. This is a highly nonlinear system of six coupled first-order differential equations. Comparison with direct numerical integration shows that the IV method indeed provides accurate analytical approximations to the orbit prediction problem.Advanced Systems Studies; Bldg. 254EElectro-Optical Systems Laboratory; Bldg. 201. 相似文献
8.
This paper studies the relative motion of satellite formation flying in arbitrary elliptical orbits with no perturbation.
The trajectories of the leader and follower satellites are projected onto the celestial sphere. These two projections and
celestial equator intersect each other to form a spherical triangle, in which the vertex angles and arc-distances are used
to describe the relative motion equations. This method is entitled the reference orbital element approach. Here the dimensionless
distance is defined as the ratio of the maximal distance between the leader and follower satellites to the semi-major axis
of the leader satellite. In close formations, this dimensionless distance, as well as some vertex angles and arc-distances
of this spherical triangle, and the orbital element differences are small quantities. A series of order-of-magnitude analyses
about these quantities are conducted. Consequently, the relative motion equations are approximated by expansions truncated
to the second order, i.e. square of the dimensionless distance. In order to study the problem of periodicity of relative motion,
the semi-major axis of the follower is expanded as Taylor series around that of the leader, by regarding relative position
and velocity as small quantities. Using this expansion, it is proved that the periodicity condition derived from Lawden’s
equations is equivalent to the condition that the Taylor series of order one is zero. The first-order relative motion equations,
simplified from the second-order ones, possess the same forms as the periodic solutions of Lawden’s equations. It is presented
that the latter are further first-order approximations to the former; and moreover, compared with the latter more suitable
to research spacecraft rendezvous and docking, the former are more suitable to research relative orbit configurations. The
first-order relative motion equations are expanded as trigonometric series with eccentric anomaly as the angle variable. Except
the terms of order one, the trigonometric series’ amplitudes are geometric series, and corresponding phases are constant both
in the radial and in-track directions. When the trajectory of the in-plane relative motion is similar to an ellipse, a method
to seek this ellipse is presented. The advantage of this method is shown by an example. 相似文献
9.
A. R. Bestman 《Astrophysics and Space Science》1991,186(1):1-6
The two-dimensional ideal flow of a neutrino gas past a vertical hot flat plate subjected to temperature fluctuation, is analysed within the framework of the continuum theory. If the difference between the wall temperature and the free-stream temperature is small, the first-order perturbed set of equations is reduced to a set of well-posed coupled linear boundary value problems for the temperature and the transverse velocity. The axial velocity and number density may then be determined. 相似文献
10.
几类辛方法的数值稳定性研究 总被引：1，自引：0，他引：1
主要对一阶隐式Euler辛方法M1、二阶隐式Euler中点辛方法M2、一阶显辛Euler方法M3和二阶leapfrog显辛积分器M4共4种辛方法及一些组合算法进行了通常意义下的线性稳定性分析．针对线性哈密顿系统,理论上找到每个数值方法的稳定区,然后用数值方法检验其正确性．对于哈密顿函数为实对称二次型的情况,为了理论推导便利,特推荐采用相似变换将二次型的矩阵对角化来研究辛方法的线性稳定性．当哈密顿分解为一个主要部分和一个小摄动次要部分且二者皆可积时,无论是线性系统还是非线性系统,这种主次分解与哈密顿具有动势能分解相比,明显扩大了辛方法的稳定步长范围． 相似文献
11.
A. G. Nikoghossian 《Astrophysics》2004,47(1):104-116
This series of papers is devoted to multiple scattering of light in plane parallel, inhomogeneous atmospheres. The approach proposed here is based on Ambartsumyan's method of adding layers. The main purpose is to show that one can avoid difficulties with solving various boundary value problems in the theory of radiative transfer, including some standard problems, by reducing them to initial value problems. In this paper the simplest one dimensional problem of diffuse reflection and transmission of radiation in inhomogeneous atmospheres with finite optical thicknesses is considered as an example. This approach essentially involves first determining the reflection and transmission coefficients of the atmosphere, which, as is known, are a solution of the Cauchy problem for a system of nonlinear differential equations. In particular, it is shown that this system can be replaced with a system of linear equations by introducing auxiliary functions P and S. After the reflectivity and transmissivity of the atmosphere are determined, the radiation field in it is found directly without solving any new equations. We note that this approach can be used to obtain the required intensities simultaneously for a family of atmospheres with different optical thicknesses. Two special cases of the functional dependence of the scattering coefficient on the optical thickness, for which the solutions of the corresponding equations can be expressed in terms of elementary functions, are examined in detail. Some numerical calculations are presented and interpreted physically to illustrate specific features of radiative transport in inhomogeneous atmospheres. 相似文献
12.
A first-order minimum propellant guidance law is developed for multi-impulse trajectories in an inverse-square gravitational field. A second-order variational analysis is used to formulate the guidance problem as an accessory minimum problem, i.e. minimize a quadratic form (second-variation in propellant consumption) subject to linear constraints (variational equations of motion and deterministic boundary conditions). Solution of the accessory minimum problem provides the optimal guidance law in feedback form. It is emphasized that this analysis takes into account the nominal impulse programme when calculating the optimal guidance corrections. It is shown that for multi-impulse transfers it is in general, non-optimal to add impulses. All corrections to the trajectory should be made by a combination of small changes in timing, magnitude and direction of the nominal impulses. 相似文献
13.
《Chinese Astronomy and Astrophysics》2007,31(2):172-186
In this paper, we analyze the linear stabilities of several symplectic integrators, such as the first-order implicit Euler scheme, the second-order implicit mid-point Euler difference scheme, the first-order explicit Euler scheme, the second-order explicit leapfrog scheme and some of their combinations. For a linear Hamiltonian system, we find the stable regions of each scheme by theoretical analysis and check them by numerical tests. When the Hamiltonian is real symmetric quadratic, a diagonalizing by a similar transformation is suggested so that the theoretical analysis of the linear stability of the numerical method would be simplified. A Hamiltonian may be separated into a main part and a perturbation, or it may be spontaneously separated into kinetic and potential energy parts, but the former separation generally is much more charming because it has a much larger maximum step size for the symplectic being stable, no matter this Hamiltonian is linear or nonlinear. 相似文献
14.
Properties of dust-acoustic solitary waves in a warm dusty plasma are analyzed by using the hydrodynamic model for massive dust grains, electrons, ions, and streaming ion beam. For this purpose, Korteweg-de Vries (KdV) equation for the first-order perturbed potential and linear inhomogeneous KdV-type equation for the second-order perturbed potential have been derived and their analytical solutions are presented. In order to show the characteristics of the dust-acoustic solitary waves are influenced by the plasma parameters, the relevant numerical analysis of the KdV and linear inhomogeneous KdV-type equations are obtained. The dust-acoustic solitary waves, as predicted here, may be associated with the nonlinear structures caused by the interaction of polar jets with the interstellar medium, which is known as Herbig-Haro objects. 相似文献
15.
Unrestricted second-order tensor virial equations for linear oscillations of magnetic configurations
M. Goossens 《Astrophysics and Space Science》1979,60(2):401-421
This paper deals with the second-order tensor virial equations for the linear oscillations of a gaseous mass in the presence of a magnetic field. It is shown that the commonly used linearized versions of the tensor virial equations are restricted integral equations that incorporate the linearized equation of motion but not the boundary condition. These restricted equations only allow trial functions that fulfil the boundary condition and are of limited practical value.The unrestricted variational principle for the linear oscillations of a magnetic configuration is used to derive a more general formulation of the second-order tensor virial equations so that the linear trial function i =X
ij
x
j can be used to study the oscillations of a configuration with a magnetic field that extends in the exterior vacuum. The unrestricted virial equations have been applied to Ferraro's model and approximate results for the eigenfrequencies and eigenfunctions have been obtained for nine oscillation modes. 相似文献
16.
17.
《New Astronomy》2013
In this paper, we present a new second kind Chebyshev (S2KC) operational matrix of derivatives. With the aid of S2KC, an algorithm is described to obtain numerical solutions of a class of linear and nonlinear Lane–Emden type singular initial value problems (IVPs). The idea of obtaining such solutions is essentially based on reducing the differential equation with its initial conditions to a system of algebraic equations. Two illustrative examples concern relevant physical problems (the Lane–Emden equations of the first and second kind) are discussed to demonstrate the validity and applicability of the suggested algorithm. Numerical results obtained are comparing favorably with the analytical known solutions. 相似文献
18.
A new algorithm is presented for the numerical integration of second-order ordinary differential equations with perturbations
that depend on the first derivative of the dependent variables with respect to the independent variable; it is especially
designed for few-body problems with velocity-dependent perturbations. The algorithm can be used within extrapolation methods
for enhanced accuracy, and it is fully explicit, which makes it a competitive alternative to standard discretization methods. 相似文献
19.
20.
Numerical integration methods for orbital motion 总被引：1，自引：0，他引：1
O. Montenbruck 《Celestial Mechanics and Dynamical Astronomy》1992,53(1):59-69
The present report compares Runge-Kutta, multistep and extrapolation methods for the numerical integration of ordinary differential equations and assesses their usefulness for orbit computations of solar system bodies or artificial satellites. The scope of earlier studies is extended by including various methods that have been developed only recently. Several performance tests reveal that modern single- and multistep methods can be similarly efficient over a wide range of eccentricities. Multistep methods are still preferable, however, for ephemeris predictions with a large number of dense output points. 相似文献