A new modified highly accurate Laplace-Fourier method for linear neutral delay differential equations

Gilbert Kerr,Gilberto Gonzalez-Parra
2024-03-23
Abstract:In this article, a new modified Laplace-Fourier method is developed in order to obtain the solutions of linear neutral delay differential equations. The proposed method provides a more accurate solution than the one provided by the pure Laplace method and the original Laplace-Fourier method. We develop and show the crucial modifications of the Laplace-Fourier method. As with the original Laplace-Fourier method, the new modified method combines the Laplace transform method with Fourier series theory. All of the beneficial features from the original Laplace-Fourier method are retained. The modified solution still includes a component that accounts for the terms in the tail of the infinite series, allowing one to obtain more accurate solutions. The Laplace-Fourier method requires us to approximate the formula for the residues with an asymptotic expansion. This is essential to enable us to use the Fourier series results that enable us to account for the tail. The improvement is achieved by deriving a new asymptotic expansion which minimizes the error between the actual residues and those which are obtained from this asymptotic expansion. With both the pure Laplace and improved Laplace-Fourier methods increasing the number of terms in the truncated series obviously increases the accuracy. However, with the pure Laplace, this improvement is small. As we shall show, with the improved Laplace-Fourier method the improvement is significantly larger. We show that the convergence rate of the new modified Laplace-Fourier solution has a remarkable order of convergence $O(N^{-3})$. The validity of the modified technique is corroborated by means of illustrative examples. Comparisons of the solutions of the new modified method with those generated by the pure Laplace method and the original/unmodified Laplace-Fourier approach are presented.
Numerical Analysis
What problem does this paper attempt to address?
This paper attempts to solve the accuracy problem in solving Linear Neutral Delay Differential Equations (NDDEs). Specifically, the author has developed a New Modified Laplace - Fourier Method, aiming to improve the accuracy of solutions of the existing pure Laplace method and the original Laplace - Fourier method. ### Main objectives of the paper: 1. **Improve the accuracy of solutions**: By introducing a new asymptotic expansion to minimize the error between the actual residual and the residual obtained from the asymptotic expansion, the accuracy of the solution is significantly improved. 2. **Improve the convergence rate**: The solution of the new method has a high convergence rate of \(O(N^{-3})\), while the convergence rate of the pure Laplace method is \(O(N^{-1})\). 3. **Retain the advantages of the original method**: The new method retains all the advantages of the original Laplace - Fourier method and further improves the accuracy by improving the method of calculating the residual. ### Key contributions of the paper: - **New asymptotic expansion**: A new asymptotic expansion is derived to approximate the residual more accurately, thereby improving the accuracy of the solution. - **Improved convergence rate**: Through theoretical analysis and numerical experiments, it is verified that the convergence rate of the new method is \(O(N^{-3})\). - **Examples of practical applications**: Through several specific examples of solving linear neutral delay differential equations, the superior performance of the new method in practical applications is demonstrated. ### Method overview: 1. **Laplace transform**: Perform a Laplace transform on the linear neutral delay differential equation to obtain an equation in the complex frequency domain. 2. **Fourier series**: Use Fourier series theory to handle the tail terms of the infinite series, thereby improving the accuracy of the solution. 3. **Asymptotic expansion**: Use a new asymptotic expansion to approximate the residual more accurately, thereby reducing the error. 4. **Numerical experiments**: Verify the effectiveness and superiority of the new method through several specific examples. ### Mathematical formulas: - **Characteristic equation**: \[ D(s)=s - a-(bs + c)e^{-s\tau}=0 \] - **Pole position**: \[ s_k=\frac{\ln(b)+2k\pi i}{\tau}, \quad k\in\mathbb{N} \] - **Residual formula**: \[ c_k=\frac{N(r_k)}{1+(br_k\tau - b + c\tau)e^{-r_k\tau}} \] - **New asymptotic expansion**: \[ c_a^k\approx-\frac{a_2}{\alpha_k^2}+i\frac{a_3}{\alpha_k^3}+O\left(\frac{1}{\alpha_k^4}\right) \] where, \[ \alpha_k=\frac{2k\pi}{\tau} \] \[ a_2=\frac{1}{\tau}\left[cH(-\tau)+aH(0)+bH'(-\tau)-H'(0)\right] \] \[ a_3=\frac{1}{\tau}\left[(a - 2\ln(b)/\tau)(cH(-\tau)+aH(0)+bH'(-\tau)-H'(0))+cH'(-\tau)+aH'(0)+bH''(-\tau)-H''(0)\right] \] Through these improvements, the new method shows higher accuracy and faster convergence speed when solving linear neutral delay differential equations.