对火星轨道变化问题的最后解释 (2 / 5)
        In Section 2 we briefly explain our dynamical model， numerical method and initial conditions used in our integrations. Section 3 is devoted to a description of the quick results of the numerical integrations. Very long-term stability of Solar system planetary motion is apparent both in planetary positions and orbital elements. A rough estimation of numerical errors is also given. Section 4 goes on to a discussion of the longest-term variation of planetary orbits using a low-pass filter and includes a discussion of angular momentum deficit. In Section 5， we present a set of numerical integrations for the outer five planets that spans ± 5 × 1010 yr. In Section 6 we also discuss the long-term stability of the planetary motion and its possible cause.

        2 Description of the numerical integrations

        （本部分涉及比较复杂的积分计算，作者君就不贴上来了，贴上来了起点也不一定能成功显示。）

        2.3 Numerical method

        We utilize a second-order Wisdom–Holman symplectic map as our main integration method (Wisdom & Holman 1991; Kino**a， Yoshida & Nakai 1991) with a special start-up procedure to reduce the truncation error of angle variables，‘warm start’(Saha & Tremaine 1992， 1994).

        The stepsize for the numerical integrations is 8 d throughout all integrations of the nine planets (N±1，2，3)， which is about 1/11 of the orbital period of the innermost planet (Mercury). As for the determination of stepsize， we partly follow the previous numerical integration of all nine planets in Sussman & Wisdom (1988， 7.2 d) and Saha & Tremaine (1994， 225/32 d). We rounded the decimal part of the their stepsizes to 8 to make the stepsize a multiple of 2 in order to reduce the accumulation of round-off error in the putation processes. In relation to this， Wisdom & Holman (1991) performed numerical integrations of the outer five planetary orbits using the symplectic map with a stepsize of 400 d， 1/10.83 of the orbital period of Jupiter. Their result seems to be accurate enough， which partly justifies our method of determining the stepsize. However， since the eccentricity of Jupiter (～0.05) is much smaller than that of Mercury (～0.2)， we need some care when we pare these integrations simply in terms of stepsizes.

        In the integration of the outer five planets (F±)， we fixed the stepsize at 400 d.

        We adopt Gauss' f and g functions in the symplectic map together with the third-order Halley method (Danby 1992) as a solver for Kepler equations. The number of maximum iterations we set in Halley's method is 15， but they never reached the maximum in any of our integrations.

        The interval of the data output is 200 000 d (～547 yr) for the calculations of all nine planets (N±1，2，3)， and about 8000 000 d (～21 903 yr) for the integration of the outer five planets (F±).

        Although no output filtering was done when the numerical integrations were in process， we applied a low-pass filter to the raw orbital data after we had pleted all the calculations. See Section 4.1 for more detail.

        2.4 Error estimation

        2.4.1 Relative errors in total energy and angular momentum

        According to one of the basic properties of symplectic integrators， which conserve the physically conservative quantities well (total orbital energy and angular momentum)， our long-term numerical integrations seem to have been performed with very small errors. The averaged relative errors of total energy (～10?9) and of total angular momentum (～10?11) have remained nearly constant throughout the integration period (Fig. 1). The special startup procedure， warm start， would have reduced the averaged relative error in total energy by about one order of magnitude or more.

        Relative numerical error of the total angular momentum δA/A0 and the total energy δE/E0 in our numerical integrationsN± 1，2，3， where δE and δA are the absolute change of the total energy and total angular momentum， respectively， andE0andA0are their initial values. The horizontal unit is Gyr.

        Note that different operating systems， different mathematical libraries， and different hardware architectures result in different numerical errors， through the variations in round-off error handling and numerical algorithms. In the upper panel of Fig. 1， we can recognize this situation in the secular numerical error in the total angular momentum， which should be rigorously preserved up to machine-ε precision.

  The content is not finished, continue reading on the next page

更多完整内容阅读登陆

《墨缘文学网，https://wap.mywenxue.org》