第一百日（1）作屎的老貓（第6/8頁）

We have adopted an FFT because of its overwhelming speed， since the amount of numerical data to be decomposed into frequency components is terribly huge (several tens of Gbytes).

A typical example of the time–frequency map created by the above procedures is shown in a grey-scale diagram as Fig. 5， which shows the variation of periodicity in the eccentricity and inclination of Earth in N+2 integration. In Fig. 5， the dark area shows that at the time indicated by the value on the abscissa， the periodicity indicated by the ordinate is stronger than in the lighter area around it. We can recognize from this map that the periodicity of the eccentricity and inclination of Earth only changes slightly over the entire period covered by the N+2 integration. This nearly regular trend is qualitatively the same in other integrations and for other planets， although typical frequencies differ planet by planet and element by element.

4.2 Long-term exchange of orbital energy and angular momentum

We calculate very long-periodic variation and exchange of planetary orbital energy and angular momentum using filtered Delaunay elements L， G， H. G and H are equivalent to the planetary orbital angular momentum and its vertical component per unit mass. L is related to the planetary orbital energy E per unit mass as E=?μ2/2L2. If the system is completely linear， the orbital energy and the angular momentum in each frequency bin must be constant. Non-linearity in the planetary system can cause an exchange of energy and angular momentum in the frequency domain. The amplitude of the lowest-frequency oscillation should increase if the system is unstable and breaks down gradually. However， such a symptom of instability is not prominent in our long-term integrations.

In Fig. 7， the total orbital energy and angular momentum of the four inner planets and all nine planets are shown for integration N+2. The upper three panels show the long-periodic variation of total energy (denoted asE- E0)， total angular momentum ( G- G0)， and the vertical component ( H- H0) of the inner four planets calculated from the low-pass filtered Delaunay elements.E0， G0， H0 denote the initial values of each quantity. The absolute difference from the initial values is plotted in the panels. The lower three panels in each figure showE-E0，G-G0 andH-H0 of the total of nine planets. The fluctuation shown in the lower panels is virtually entirely a result of the massive jovian planets.

Comparing the variations of energy and angular momentum of the inner four planets and all nine planets， it is apparent that the amplitudes of those of the inner planets are much smaller than those of all nine planets: the amplitudes of the outer five planets are much larger than those of the inner planets. This does not mean that the inner terrestrial planetary subsystem is more stable than the outer one: this is simply a result of the relative smallness of the masses of the four terrestrial planets compared with those of the outer jovian planets. Another thing we notice is that the inner planetary subsystem may become unstable more rapidly than the outer one because of its shorter orbital time-scales. This can be seen in the panels denoted asinner 4 in Fig. 7 where the longer-periodic and irregular oscillations are more apparent than in the panels denoted astotal 9. Actually， the fluctuations in theinner 4 panels are to a large extent as a result of the orbital variation of the Mercury. However， we cannot neglect the contribution from other terrestrial planets， as we will see in subsequent sections.

4.4 Long-term coupling of several neighbouring planet pairs

Let us see some individual variations of planetary orbital energy and angular momentum expressed by the low-pass filtered Delaunay elements. Figs 10 and 11 show long-term evolution of the orbital energy of each planet and the angular momentum in N+1 and N?2 integrations. We notice that some planets form apparent pairs in terms of orbital energy and angular momentum exchange. In particular， Venus and Earth make a typical pair. In the figures， they show negative correlations in exchange of energy and positive correlations in exchange of angular momentum. The negative correlation in exchange of orbital energy means that the two planets form a closed dynamical system in terms of the orbital energy. The positive correlation in exchange of angular momentum means that the two planets are simultaneously under certain long-term perturbations. Candidates for perturbers are Jupiter and Saturn. Also in Fig. 11， we can see that Mars shows a positive correlation in the angular momentum variation to the Venus–Earth system. Mercury exhibits certain negative correlations in the angular momentum versus the Venus–Earth system， which seems to be a reaction caused by the conservation of angular momentum in the terrestrial planetary subsystem.

It is not clear at the moment why the Venus–Earth pair exhibits a negative correlation in energy exchange and a positive correlation in angular momentum exchange. We may possibly explain this through observing the general fact that there are no secular terms in planetary semimajor axes up to second-order perturbation theories (cf. Brouwer & Clemence 1961; Boccaletti & Pucacco 1998). This means that the planetary orbital energy (which is directly related to the semimajor axis a) might be much less affected by perturbing planets than is the angular momentum exchange (which relates to e). Hence， the eccentricities of Venus and Earth can be disturbed easily by Jupiter and Saturn， which results in a positive correlation in the angular momentum exchange. On the other hand， the semimajor axes of Venus and Earth are less likely to be disturbed by the jovian planets. Thus the energy exchange may be limited only within the Venus–Earth pair， which results in a negative correlation in the exchange of orbital energy in the pair.