Abstract

Phase fitted and amplification fitted Numerov-type methods for periodic initial value problems with two frequencies are investigated. A one-frequency method and a two-frequency method are constructed. The two new methods both have algebraic order five and are dispersive of order six and dissipative of order five. The two-dimensional absolute stability region for the one-frequency method and the three-dimensional absolute stability region for the two-frequency method are plotted. Numerical experiments are reported to show the efficiency and competence of the two new methods.

1. Introduction

In this paper we are interested in the numerical integration of the initial value problem (IVP) of ordinary differential equations in the form whose solutions exhibit a pronounced oscillatory character. Such problems are frequently encountered in a variety of scientific fields and engineering applications. Since the analytical solutions of these equations are usually not available, the numerical solutions become very important. To our knowledge, most existing methods of fitted type (trigonometrically/exponentially fitted, phase fitted and amplification fitted, etc.) are adapted to one frequency (see [13]). We call them one-frequency methods. For example, by phase fitting, Van de Vyver [4] constructed an adapted explicit Numerov-type method for oscillatory problems with one frequency. Other related work can be found in [520]. However, for oscillatory problems with many frequencies in the solution, when we apply a one-frequency method, it is difficult to choose the most suitable fitting frequency. Contrast to the one-frequency methods, papers on integrators that are adapted to more than one frequency are few. In [21], Wang obtained a trigonometrically fitted Numerov-type method for periodic IVPs with two frequencies and provided some numerical examples to show that his method associated with two frequencies is more accurate and more efficient than the trigonometrically fitted Numerov methods with one frequency. More recently, Fang et al. [22] gave a trigonometrically fitted explicit Numerov-type method for the numerical integration of oscillatory problems with two frequencies. Inspired by Van de Vyver’s idea in [4], we investigate, in this paper, explicit Numerov-type methods for the numerical integration of periodic IVPs with two frequencies. The paper is organized as follows. In Section 2, the preliminaries of stability and phase properties of explicit Numerov-type methods are recalled. In Section 3, we construct two new phase fitted and amplification fitted explicit Numerov-type methods for periodic initial value problems with two frequencies. In Section 4, the stability and phase properties of the newly derived methods are analyzed. In Section 5, we present some numerical results. Section 6 is devoted to the conclusions.

2. Numerov-Type Methods

In this section, we recall some elementary notions about Numerov-type methods. The explicit Numerov-type methods considered here have the form with , , and can be expressed in the Butcher tableau as or equivalently by the triplet . The order conditions for this kind of Numerov-type methods can be found in [23]. A favorable property of these methods is that they need only function evaluations at each step. According to the idea of Lambert and Watson [24], the stability of the Numerov-type methods can be analyzed by the test equation Applying the explicit Numerov-type method (2.1) to (2.3), we obtain the following relation with Since the coefficient matrix satisfies , we have where is the -dimensional vector of units, and the vectors , and matrix are defined by the Butcher tableau of the scheme (2.1). The polynomial is called the stability polynomial of the difference equation (2.4) associated with the method (2.1).

Definition 2.1. The method (2.1) is said to have a periodicity interval if the roots of its stability polynomial, are conjugate complex and satisy , for all . The method is called P-stable if the periodicity interval .

The existence of a nonempty periodicity interval is equivalent to the fact that there exists a nonzero real number such that the coefficients of the stability polynomial satisfy the conditions Therefore, the -stable conditions are given by An interval is called the interval of absolute stability of the method (2.1) if (see [4])

Definition 2.2. The quantities (see [25]) are called the phase-lag and the dissipation, respectively. Accordingly, the method is said to be dispersive of order and dissipative of order if , and , respectively. If , the method is said to be phase fitted or zero-dispersive, and if , the method is said to be amplification fitted or zero-dissipative.

3. Construction of the New Methods

In this section, we construct explicit Numerov-type methods which are both phase fitted and amplification fitted. The construction is based on the following Butcher tableau of coefficients If we choose , the classical Numerov-type method is recovered in [26]. In the following, we derive the frequency-dependent weights of a one-frequency method and a two-frequency method by phase fitting and amplification fitting.

3.1. A One-Frequency Numerov-Type Method

Applying the method (3.1) to the test equation we obtain the following expressions for and Nullifying the dispersion and the dissipation in Definition 2.2 for the test equation (3.2) yields Choosing , we solve the linear systems (3.4) for the -values as follows where .

It is easy to verify the above coefficients satisfy the algebraic fifth-order conditions for the Numerov-type methods with one frequency (see [23]): This one-frequency phase fitted and amplification fitted fifth-order method is denoted as PFAFI.

3.2. A Two-Frequency Numerov-Type Method

For the method (3.1), we only need the weights . In order to solve the oscillatory problems with two frequencies, we consider another test equation With the same idea stated above, we find the following relation Solving (3.4) and (3.8) gives In which .

It is also easy to verify the algebraic fifth order conditions for the Numerov-type methods with two frequencies (see [23] for reference) The two-frequency fifth-order Numerov-type method with the weights in (3.9) is denoted as PFAFII.

4. Stability and Phase Analysis

In this section, we analyze the stability and phase properties of the two new methods derived in Section 3. The theory of Lambert and Watson was reconsidered by Coleman and Ixaru [27].

4.1. Stability and Phase Analysis of PFAFI

We apply a four-stage one-frequency Numerov-type method to the test equation (2.3) and obtain the functions and in (2.4) as follows The stability and the phase properties depend on the bivariate functions and . In this situation, the interval of absolute stability is replaced by a two-dimensional region.

Definition 4.1. For the Numerov-type method with and given by (4.1), the quantities are called the phase-lag and the dissipation, respectively. If then the method is said to be dissipative of order and dissipative of order , respectively.

For convenience of analyzing the phase-lag and the dissipation, we denote the ratio . Then the phase-lag and the dissipation of the new method PFAFI are given by Thus the method PFAFI is dispersive of order six and dissipative of order five.

Definition 4.2. In the - plane, the region is said to be the region of absolute stability of the Numerov-type method PFAFI. Any closed curve defined by is a stability boundary of the method.

In Figure 1, we plot the region of the absolute stability of the new method PFAFI derived in Section 3.1 (a similar discussion can be found in [27, 28]).

It should be noted that when the exact value of the frequency of the problem is known (i.e., the fitted frequency equals the test frequency ), then and we have that and . Consequently, when the main frequency is known the interval of periodicity is except for . So the new method PFAFI is almost -stable.

4.2. Stability and Phase Analysis of Phase Fitted Numerov-Type Method with Two Frequencies

We apply a four-stage two-frequency Numerov-type method to the test equation (2.3) and obtain the functions and in (2.4) as follows

Definition 4.3. For the Numerov-type method with and given by (4.7), the quantities are called the phase-lag and the dissipation, respectively. If then this method is said to be dispersive of order and dissipative of order , respectively.

Denote and . Then the phase-lag and the dissipation of the new method PFAFII are given by Thus the new method PFAFII is dispersive of order six and dissipative of order five.

Following the main idea of Coleman and Ixaru in [27], the absolute stability of the two-frequency method PFAFII can be described by a three-dimensional region.

Definition 4.4. A three-dimensional region in the -- space is called the region of absolute stability of a Numerov-type method with the and given by (4.7), and any closed curve defined by is a stability boundary of the method.

In Figure 2, we plot the region of absolute stability for the method PFAFII, and Figure 3 shows the projection of the three-dimensional region onto the - plane. (A similar discussion can be found in [22, 29].)

It should be noted that, for the new method PFAFII and for a problem with one frequency, when the exact value of the frequency of the problem is known, that is, , then and we have that and . Consequently, the interval of periodicity is with the points deleted. That is, the new method PFAFII is also almost -stable.

5. Numerical Experiments

In this section, we compare the numerical performance of the new methods with some well-known methods from the literature. We select three periodic IVPs with two frequencies. The methods we choose for comparison are listed as follows. (i) EFRKN4F: the trigonometrically fitted RKN method of order 4 with the FSAL property given in [2]. (ii) VPHA: the phase fitted and amplification fitted Numerov-type method with one frequency given in [4]. (iii) TFNTONE: the trigonometrically fitted Numerov-type method obtained by Fang and Wu in [3]. (iv) PFAFI: the phase fitted and amplification fitted Numerov-type method with one frequency obtained in this paper. (v) PFAFII: the phase fitted and amplification fitted Numerov-type method with two frequencies obtained in this paper.

Problem 1. We consider the inhomogeneous equation whose exact solution is given by In our test we choose for VPHA, TFNTONE, EFRKN4F and PFAFI, and   for PFAFII. The problem is integrated on the interval [] and the numerical results are presented in Table 1.

Problem 2. Consider the famous Kramarz system (see [30]) whose exact solution is The eigenvalues of the matrix of coefficients of the equations for and are −1 and −2500. For this stiff equation, we choose for VPHA, TFNTONE, EFRKN4F and PFAFI, and for PFAFII. The problem is integrated on the interval [] and the numerical results are presented in Table 2.

Problem 3. Consider the linear stiff system in [31] whose analytic solution is given by In our test we choose for VPHA, TFNTONE, EFRKN4F and PFAFI, and for PFAFII. The problem is integrated on the interval [] and the numerical results are presented in Table 3 with , and in Table 4 with .

Problem 4. Consider the nonlinear system in [32] where . The analytic solution of the problem is

In our test we choose for VPHA, TFNTONE, EFRKN4F and PFAFI, and for PFAFII. The problem is integrated on the integral interval [] and the numerical results are stated in Table 5.

6. Conclusions

In this paper, we investigate Numerov-type methods adapted to oscillatory problems with two frequencies obtained by phase fitting and amplification fitting. The numerical stability and the phase properties of the new method are analyzed. Unlike the case of one-frequency methods, the absolute stability region of a two-frequency method is a region in the three-dimensional space. The results of numerical experiments on four oscillatory test problems with two frequencies, including non-stiff and stiff problems, show that the new methods are superior to some well-known high quality codes proposed in the recent scientific literature.

Acknowledgments

This research was partially supported by NSFC (no. 11101357, no. 11171155), the foundation of Shandong Outstanding Young Scientists Award Project (no. BS2010SF031), the foundation of Scientific Research Project of Shandong Universities (no. J11LG69), NSF of Shandong Province, China (no. ZR2011AL006) and the Fundamental Research Fund for the Central Universities (no. Y0201100265, no. KYZ201125).