On the Evaluation of Oscillatory Kind with Cauchy Principal Value ()
1. Introduction
In this paper, we concerned about the Cauchy principal value in oscillatory kind function which
, specifically have the form like
where
is smooth and
, generally
is much large. Li [1] uses Levin method to solve the integral with Cauchy principal value. In that case, the principal value point is
and
, but the function
. And the Levin method uses Chebyshev-Gauss-Labatt nodes to make a matrix equation in order to get the approximate solution. Then Wang [2] has mentioned a technique of put the integral to complex plane so that we can use a loop integral. In this paper, the author chooses several paths to make it available tot Gauss-Laguerre quadrature. It is a great idea for us to solve Cauchy principal value problem and we are using this idea to solve problem in this paper. Furthermore, Daan [3] puts forward the form of integral form. But in his paper, the integral without Cauchy principal value. Importantly, we take his idea of choose path, which make it easy to have a good path to integral. Alfrddo [4] then put the question to a more general situation. With the function
not smooth, maybe has stationary points and so on, they use Gauss-Fraud quadrature to solve the path near the stationary points.
In this paper, we firstly transformed the equation into the form which we known before, then we take the path to get the solution with Gauss-Laguerre quadrature. Finally two numerical examples are given to prove the method. One is for real part, and the other is for imaginary part.
2. Integral with Gauss-Laguerre Quadrature
As we can see, the variable x has a different order. First we need to transform x to let them have the same order. The specific process is
In last equation, we use the variable substation of
to transform it. The equation can be calculated in some methods like Li [1] has put forward.
In such complex plane, the integral can be divided into several pieces which can be seen in above figure. Therefore, we only need to calculate the integral of path
so that we can have the integral of path
.
The path
and
have similar situation. have mentioned a method about select a path to infinity as following
where
is the path start in x and end in infinity,
is the function
, and i is the imaginary unit, which is
, p is the variable. Now we prove the start point and the end point. Because
is smooth, so
When
,
, when
,
. In this case, because of the known form of integral, the path can be written as following [3]
We need to notice that, p is the variable of
, and x is the start point of the path.
Observe that
where
in the last equation.
The same
We use polar coordinates to solve
. Let
, where r is the radio and θ is the variable.
When
, the equation above can be transformed into
As we discussed above, the integral () now can be written as following
Fortunately, the infinite integral can solve with Gauss-Laguerre quadrature with n-points. The
represents Gauss-Laguerre points and the
represents Gauss-Laguerre weights.
Remark 1
is the quadrature rule that we mentioned before with n points and weights. The approximation error behaves like
3. Numerical Example
In order to prove this method is effective, we will use two examples to show the approximate solution and exact solution. All calculation is computed in Matlab2019a, and exact solution is given by Mathematica 12.0 with 64-digit.
Example1. The integral
is analytic expression. Using above method, we can use
to approach it which have the form like following
Table 1 demonstrates the errors in real part. We can draw a conclusion that, with ω and n increasing, the accuracy of method is increasing. We take
as example, when
the results approach the max accuracy. When
, this process moves up when
.
Table 1. The real part error of
.
Table 2. The imaginary part error of
.
Example 2. The integral of
only have imaginary part, so in this example, we can see the performance in imaginary part. This time we take
, and we still take
to compare the efficiency about the method.
As we can see from Table 2, with n-points Gauss-Laguerre quadrature, the convergence is fast and only small n can obtain satisfying solution.
4. Conclusion
In this paper, we discuss the oscillatory kind function with Cauchy principal value which has a different form. The integral can be transformed into some known situation so that we can handle it easier. The method that puts the integral in complex plane in this paper comes from Wang [3] and the numerical solution is shown in Table 1 and Table 2. Within the two examples, we compare the real part and imaginary part between numerical solution and exact solution. And the absolute error represents the method for this integral is available.
Acknowledgements
The author would like to acknowledge to my tutor for teaching and useful discussion, and anonymous referees for suggestions.