Abstract:We extend the piecewise orthogonal collocation method to computing periodic solutions of coupled renewal and delay differential equations. Through a rigorous error analysis, we prove convergence of the relevant finite-element method and provide a theoretical estimate of the error. We conclude with some numerical experiments to further support the theoretical results.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is how to calculate the periodic solutions of the system of coupled delay differential equations (DDEs) and renewal equations (REs). Specifically, the authors extended the piecewise orthogonal collocation method to deal with the periodic solutions of such equations and proved the convergence of this method through strict error analysis.
### Main problem decomposition
1. **Calculation of periodic solutions**:
- The paper focuses on how to effectively calculate the periodic solutions of the coupled DDEs and REs system. These equations are widely used in fields such as biology and epidemiology, for example, to describe the time - delay effects in population dynamics.
2. **Method extension**:
- The authors extended the piecewise orthogonal collocation method from individual DDEs or REs to the coupled system. This requires dealing with the coupling between two different types of equations and the different regularity of their respective defined spaces.
3. **Convergence proof**:
- To ensure the effectiveness of the method, the authors carried out strict error analysis and proved the convergence of this method. They provided a theoretical error estimate, indicating that the error decays as the grid division is refined.
4. **Numerical experiment verification**:
- Finally, the authors further supported their theoretical results through some numerical experiments, demonstrating the effectiveness of this method in practical applications.
### Overview of formulas and methods
- **Coupled system**:
\[
\begin{cases}
x(t)=F(x_t,y_t), \\
y'(t)=G(x_t,y_t),
\end{cases}
\]
where \(x_t\) and \(y_t\) represent the historical states at time \(t\) respectively.
- **Boundary value problem (BVP)**:
The periodic solution can be obtained by solving a boundary value problem in the form of:
\[
\begin{cases}
x(t)=F(x_t,y_t), & t\in[0,\omega], \\
y'(t)=G(x_t,y_t), & t\in[0,\omega], \\
x(0)=x(\omega), \\
y(0)=y(\omega), \\
p(x,y) = 0,
\end{cases}
\]
where \(p(x,y)=0\) is the phase condition used to remove translational invariance.
- **Piecewise orthogonal collocation method**:
This method numerically solves the above BVP by using piecewise polynomial to approximate the solution on each sub - interval and satisfying the equations at the selected collocation points.
### Conclusion
The main contribution of the paper is to provide the first numerical tool capable of calculating the periodic solutions of the coupled DDEs and REs system and prove its convergence through strict theoretical analysis. This provides strong support for studying complex systems such as physiologically structured populations.