Exponential sum approximations for $t^{-β}$

William McLean

DOI: https://doi.org/10.1007/978-3-319-72456-0_40

2017-06-07

Abstract:Given $\beta>0$ and $\delta>0$, the function $t^{-\beta}$ may be approximated for $t$ in a compact interval $[\delta,T]$ by a sum of terms of the form $we^{-at}$, with parameters $w>0$ and $a>0$. One such an approximation, studied by Beylkin and Monzón, is obtained by applying the trapezoidal rule to an integral representation of $t^{-\beta}$, after which Prony's method is applied to reduce the number of terms in the sum with essentially no loss of accuracy. We review this method, and then describe a similar approach based on an alternative integral representation. The main difference is that the new approach achieves much better results before the application of Prony's method; after applying Prony's method the performance of both is much the same.

Numerical Analysis

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to provide an effective exponential sum approximation method for the function \( t^{-\beta} \) within a compact interval \([\delta, T]\). Specifically, the author studied two exponential sum approximation methods based on different integral representations and explored how to use the Prony method to reduce the number of terms in the approximation without losing precision.

### Main problems
1. **Efficient approximation**: For given \( \beta > 0 \) and \( \delta > 0 \), how to efficiently approximate the function \( t^{-\beta} \) within the compact interval \([\delta, T]\).
2. **Reduce computational complexity**: By representing the convolution kernel \( k(t)=t^{-\beta} \) in the form of an exponential sum, reduce the number of operations in numerical calculations, thereby avoiding the quadratically increasing computational complexity.
3. **Improve approximation accuracy**: Ensure that while applying the Prony method to reduce the number of terms, the high - precision of the approximation is maintained.

### Solutions
The author proposed two main methods to achieve these goals:

1. **Alternative method based on \( p = e^x \)**:
- Use the integral representation \( \frac{1}{t^\beta}=\frac{1}{\Gamma(\beta)}\int_{-\infty}^{\infty}\exp(-te^x+\beta x)\,dx \).
- Apply the infinite trapezoidal rule to obtain the approximation formula \( \frac{1}{t^\beta}\approx\frac{1}{\Gamma(\beta)}\sum_{n = -\infty}^{\infty}w_n e^{-a_n t} \), where \( a_n = e^{hn} \) and \( w_n = h e^{\beta nh} \).

2. **Alternative method based on \( p=\exp(x - e^{-x}) \)**:
- Use the integral representation \( \frac{1}{t^\beta}=\frac{1}{\Gamma(\beta)}\int_{-\infty}^{\infty}\exp(-\phi(x, t))(1 + e^{-x})\,dx \), where \( \phi(x, t)=t\exp(x - e^{-x})-\beta(x - e^{-x}) \).
- Similarly, apply the infinite trapezoidal rule to obtain a similar approximation formula.

### Key results
- **Application of the Prony method**: Through the Prony method, the number of terms in the approximation can be significantly reduced while keeping the relative error within an acceptable range.
- **Error analysis**: Provide detailed error estimates to ensure the effectiveness of the approximation under different parameter selections.
- **Performance comparison**: Compared with previous methods, the new method can achieve the same precision with fewer terms in some cases.

### Conclusion
Through these two methods, especially after applying the Prony method, the amount of calculation can be significantly reduced while ensuring precision, thereby improving the efficiency of numerical calculations. This is particularly important when dealing with large - scale vectors and long - time steps.