On the convergence of the pseudospectral approximation of reproduction numbers for age-structured models

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: De Reggi, Simone, Scarabel, Francesca, Vermiglio, Rossana
Natura: Preprint
Pubblicazione: 2024
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866914934254731264
author De Reggi, Simone
Scarabel, Francesca
Vermiglio, Rossana
author_facet De Reggi, Simone
Scarabel, Francesca
Vermiglio, Rossana
contents We rigorously investigate the convergence of a new numerical method, recently proposed by the authors, to approximate the reproduction numbers of a large class of age-structured population models with finite age span. The method consists in reformulating the problem on a space of absolutely continuous functions via an integral mapping. For any chosen splitting into birth and transition processes, we first define an operator that maps a generation to the next one (corresponding to the Next Generation Operator in the case of R0). Then, we approximate the infinite-dimensional operator with a matrix using pseudospectral discretization. In this paper, we prove that the spectral radius of the resulting matrix converges to the true reproduction number, and the (interpolation of the) corresponding eigenvector converges to the associated eigenfunction, with convergence order that depends on the regularity of the model coefficients. Results are confirmed experimentally and applications to epidemiology are discussed.
format Preprint
id arxiv_https___arxiv_org_abs_2409_01520
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle On the convergence of the pseudospectral approximation of reproduction numbers for age-structured models
De Reggi, Simone
Scarabel, Francesca
Vermiglio, Rossana
Numerical Analysis
Dynamical Systems
34L16, 47B07, 65L15, 65L60, 65N12, 92D25
We rigorously investigate the convergence of a new numerical method, recently proposed by the authors, to approximate the reproduction numbers of a large class of age-structured population models with finite age span. The method consists in reformulating the problem on a space of absolutely continuous functions via an integral mapping. For any chosen splitting into birth and transition processes, we first define an operator that maps a generation to the next one (corresponding to the Next Generation Operator in the case of R0). Then, we approximate the infinite-dimensional operator with a matrix using pseudospectral discretization. In this paper, we prove that the spectral radius of the resulting matrix converges to the true reproduction number, and the (interpolation of the) corresponding eigenvector converges to the associated eigenfunction, with convergence order that depends on the regularity of the model coefficients. Results are confirmed experimentally and applications to epidemiology are discussed.
title On the convergence of the pseudospectral approximation of reproduction numbers for age-structured models
topic Numerical Analysis
Dynamical Systems
34L16, 47B07, 65L15, 65L60, 65N12, 92D25
url https://arxiv.org/abs/2409.01520